dt=0.005;
N=2000;
global w;
load w1.dat -ascii;
w=w1;
%初始值
t0=0.00001;
y0=[0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001]';

[t,y]=LK2('LK1',t0,dt,N-1,y0);

%繪時域圖
figure(1);
subplot(2,2,1);
y1=y(1,:);
ts=0:N-1;
plot(dt*ts,y1);
xlabel('時間  t / s');
ylabel('路面位移輸入I1 / m');
subplot(2,2,2);
y2=y(2,:);
plot(dt*ts,y2);
xlabel('時間  t / s');
ylabel('路面位移輸入I2 / m');
subplot(2,2,3);
y3=y(3,:);
plot(dt*ts,y3);
xlabel('時間  t / s');
ylabel('路面位移輸入I3 / m');
subplot(2,2,4);
y4=y(4,:);
plot(dt*ts,y4);
xlabel('時間  t / s');
ylabel('路面位移輸入I4 / m');