1111.zip_地震

ElCentro;

Accelerate= ElCentro(:,1)*9.8067;

N=length(Accelerate);

time=0:0.005:(N-1)*0.005;

Displace=zeros(1,N);

Velocity=zeros(1,N);

AbsAcce=zeros(1,N);

Damp=0.02;

TA=0.0:0.05:6;

NatualFrequency=2*pi/T;

DampFrequency=NatualFrequency*sqrt(1-Damp*Damp);

e_t=exp(-Damp*NatualFrequency*Dt);

s=sin(DampFrequency*Dt);

c=cos(DampFrequency*Dt);

A=zeros(2,2);

A(1,1)=e_t*(s*Damp/sqrt(1-Damp*Damp)+c);

A(1,2)=e_t*s/DampFrequency;

A(2,1)=-NatualFrequency*e_t*s/sqrt(1-Damp*Damp);

A(2,2)=e_t*(-s*Damp/sqrt(1-Damp*Damp)+c);

d_f=(2*Damp^2-1)/(NatualFrequency^2*Dt);

d_3t=Damp/(NatualFrequency^3*Dt);

B=zeros(2,2);

B(1,1)=e_t*((d_f+Damp/NatualFrequency)*s/DampFrequency+(2*d_3t+1/NatualFrequency^2)*c)-2*d_3t;

B(1,2)=-e_t*(d_f*s/DampFrequency+2*d_3t*c)-1/NatualFrequency^2+2*d_3t;

B(2,1)=e_t*((d_f+Damp/NatualFrequency)*(c-Damp/sqrt(1-Damp^2)*s)-(2*d_3t+1/NatualFrequency^2)*(DampFrequency*s+Damp*NatualFrequency*c))+1/(NatualFrequency^2*Dt);

B(2,2)=e_t*(1/(NatualFrequency^2*Dt)*c+s*Damp/(NatualFrequency*DampFrequency*Dt))-1/(NatualFrequency^2*Dt);

for i=1:(N-1)

if T==0.0

MaxA(1,t)=max(abs(Accelerate));

else

MaxA(1,t)=max(abs(AbsAcce));

end

Displace=zeros(1,N);

AbsAcce=zeros(1,N);

t=t+1;

End

close all