阶跃响应.m

来自「系统辨识与参数估计源代码。可以对一个未知模型做出估计。」· M 代码 · 共 18 行

clear
L=3000;
for i=1:L
   u(i)=1;
  
end
figure(1);
subplot(3,1,1);
stem(u),grid on;
v=randn(1,L);
subplot(3,1,2);
plot(v),grid on;
y(1)=0;y(2)=0;y(3)=0;
for k=4:L
    y(k)=1.5*y(k-1)-0.7*y(k-2)+u(k-2)+0.5*u(k-3)+v(k)+0.9*v(k-1)+0.7*v(k-2);
    subplot(3,1,3);
    plot(y),grid on;
end