信号与线性系统实验报告(17)

x=[n>=0];%与x=stepseq(0,-N,k);作用相同 %x=stepseq(0,-N,k);

yzs(N+1:end)=filter(b,a,x(N+1:end));

yzi(N+1:end)=filter([0 0],a,x(N+1:end),zic); y(N+1:end)=filter(b,a,x(N+1:end),zic); figure

subplot(3,1,1) stem(n,x)

title('激励 x(n)'); subplot(3,1,2) stem(n,h)

title('冲激响应 h(n)'); subplot(3,1,3) hold on

stem(n,yzs,'g') stem(n,yzi,'r') stem(n,y) hold off

legend('yzs','yzs','yzi','yzi','y','y') title('输出响应 y(n)'); text={... ''

' 单位冲激响应 h=' ''

[' '*ones(1,12),num2str(h)] ''

' 零状态响应yzs=' ''

[' '*ones(1,12),num2str(yzs)] ''

' 零输入响应yzi=' ''

[' '*ones(1,12),num2str(yzi)] ''

' 全响应 y=' ''

[' '*ones(1,12),num2str(y)] ''};

textwin('差分方程数值解',text)