求助：反应时程曲线总是不对，程序如下

fenghuayu

高手帮忙看下我这个matlab小程序，输入地震波求解结构的位移和加速度反应，但是输出的位移和加速度反应时程曲线总是不对，不知道哪里出错了。专家帮忙看看呗？

clear all
close all
load taft.dat
elcen=taft;
dt=elcen(1);nn=length(elcen);
for i=1:nn-1
    e(2,i)=elcen(i+1);
end
w=4*e(2,:)/max(abs(e(2,:)));
%dt=e(1,2)-e(1,1);
wno=1;
t=[0:dt:(length(w)-1)*dt];
n=1;
k=9.3436e4;
c=1367;
m=3126.8e3;
%a=[0 1;-inv(m)*k -inv(m)*c];
%b=[0;-1];
%e=[1 0;-inv(m)*k -inv(m)*c];
%d=[0;-1];
a=[zeros(n,n)  eye(n);  -inv(m)*k  -inv(m)*c];
b=[zeros(n,n);-1];
d=[eye(n)  zeros(n,n);-inv(m)*k  -inv(m)*c];
e=[zeros(n,n);-1];
[y,x]=lsim(a,b,d,e,w,t);

dd1=max(abs(y(:,1)))*1000    %位移单位mm
aa1=max(abs(y(:,2)))*100     %加速度单位cm/s2,即gal
ttt=t';
y1=[ttt y(:,1)*1000];
y2=[ttt y(:,2)*100];
end_out=[dd1,aa1]

liuyang5299

能否发一下输出的结果，另外建议对一些自定义的变量作一些说明，这样也好帮你分析。而且时程分析有好几种方法，不知你依据什么公式计算的。