f14_6.m

matlab时频分析配套源码第14章的程序

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');