matlab拟合反应谱,由地震波拟合反应谱的matlab程序调试

function pars(ethwv,tmitvl,dmpr)

ethwv1;

tmitvl1;

dmpr1;

% function pars(ethwv,tmitvl,dmpr)

% Purpose: To draw pseudo-acceleration response spectrum figure in 6 s

% variable names:

% ethwv acceleration of earthquake wave record in one column form

% tmitvl time interval of the earthquake wave record

% dmpr damping ratio

t = 0:tmitvl:(length(ethwv)-1)*tmitvl;

Prd = 0:0.01:6; Prd = Prd;

omgar = 2*pi./Prd;

for i = 1:length(omgar)

SYS = tf(1,[1,2*dmpr*omgar(i),omgar(i)*omgar(i)]);

disp(i) = max(abs(lsim(SYS,ethwv,t)));

spd(i) = 2*pi.*disp(i)/Prd(i);

acc(i) = 2*pi.*spd(i)/Prd(i);

end

plot(Prd,acc);