长江大学信号与系统matlab实验答案

0Magnitude (dB)Bode Diagram-50-100-150180Phase (deg)900-90-21010-110Frequency (rad/s)0101102

Eg3.29(a) 代码如下：

b=[0 1]; n=[1 3];m=[1 6]; c=conv(n,m); a=c; figure(1)

w=linspace(0,200,200); zpplot(b,a) freresp(b,a,w) figure(2) t=0:0.02:5; subplot(1,2,1); impulse(b,a,t); subplot(1,2,2); step(b,a,t); figure(3) bode(b,a);

结果：

零极点分布图0.0660.0540.042幅频特性图0-20-40-60相频特性图相位(度)0100?2000j)||F(?0.03-80-100-120-140-160-20.02-40.01-60-505-1800100?200

Impulse ResponseStep Response0.090.080.060.050.070.06Amplitude0.04Amplitude0.050.040.030.020.030.020.010.01001234Time (seconds)5001234Time (seconds)5

-20Magnitude (dB)Bode Diagram-40-60-80-100-1200Phase (deg)-45-90-135-180-11010010Frequency (rad/s)1102103

实验4离散信号波形绘制和迭代法计算

上机内容

? 用MATLAB计算习题4-4。

用MATLAB计算习题4-6。并显示前15项数值，画出其波形。 Eg4.4

代码如下：

function [y]=f(k) %定义函数f(k)

y=-((k-3)/3).*(u(k)-u(k-4))+delta(k+1)+delta(k+2); end

k=-5:6; y=f(k); y1=f(k+2);

y2=f(k+2).*u(-k-2); y3=f(-k+2);

y4=f(-k+2).*u(k-1);

subplot(5,1,1),ylabel('f(k)') mystem(k,y),title('信号f(k)的波形图'); subplot(5,1,2),ylabel('f(k+2)')

mystem(k,y1),title('信号f(k+2)的波形图'); subplot(5,1,3),ylabel('f(k+2)ε(-k-2)')