输入额e(k)=1;输出为u(k);k=1:1:100;
u(k)=173.2e(k)-346.3e(k-1)+173e(k-2)+1.951u(k-1)-0.9512u(k-2)
初始时刻应该都是取0吧?e(0)=e(-1)=u(0)=u(-1)=0;e(k)表示单位阶跃输入;采样时间ts=0.0001;
time(k)=k*ts;
plot(time,e,'k',time,u,'k');
我不管你的采样时间,只讨论离散部分。
e = zeros(102, 1); % 把-1和0也算进去,共102项