Civil Engineering Reference
In-Depth Information
psi=sqrt(sigma-1.0);
else
u=(ti-bk)/aj;
rnum=sin(sigma*pi*u)-sin(pi*u);
den=u*pi*sqrt(sigma-1.0);
psi=rnum/den;
end
w=w+zz(i)*psi*tsm;
end
baj=sqrt(abs(aj));
wvcfacc(ibk,j+18)=w/baj;
end
end
save(‚wvcfacccmpsps.mat','wvcfacc'); %2047 rows, 22 columns
%
% Calculation of wavelet coefficients finished
%
% Ground motion characterization starts now (vide Section
3.1)
%
for i=1:2047
for j=1:22
wvcfaccsq(i,j)=wvcfacc(i,j)^2;
end
end
bb=0.2;
cpsi=2.0*log(sigma)/(2.0*(sigma-1.0)*pi);
sk=2.0/(4.0*pi*cpsi)*(sigma-1.0/sigma);
skk=sk/((sigma-1.0)*pi);
for jj=-17:4
j=jj+18;
tj=(4.0*(sigma^jj))/(1.0+sigma);
jjj=tj/tsm;
im=1;
sum=0.0;
is=0;
for i=im:2047
is=is+1;
sum=sum+wvcfaccsq(i,j);
if((is<jjj)&&(i==2047))
jjj=is;
wvs=sum/jjj;
ik=i-jjj+1;
for ie=ik:i
wv(ie,j)=wvs;
end
