信号与系统实验四 离散时间LTI系统分析实验报告资料

结论：该因果系统的极点全部在单位圆内，故系统是稳定的。

（5）试用MATLAB命令画出该系统的频率响应曲线。 程序：

a=[1 0.2 -0.24]; b=[1 1 0];

[H,w]=freqz(b,a,400,'whole'); Hm=abs(H); Hp=angle(H);

subplot(211);plot(w,Hm);

xlabel('\\omega(rad/s)');ylabel('幅度'); title('离散系统幅频特性曲线');

subplot(212);plot(w,Hp);

xlabel('\\omega(rad/s)');ylabel('相位'); title('离散系统相频特性曲线');

运行结果：

5

2、已知系统函数H(z)?zz?k(k为常数)，求系统的频率响应，k?0,0.5,1三种情况下系统的幅度响应和相位响应。

程序：

k=0; %k=0,0.5,1; a=[1 k]; b=[1 0];

[H,w]=freqz(b,a,400,'whole'); Hm=abs(H); Hp=angle(H);

subplot(211);plot(w,Hm);

xlabel('\\omega(rad/s)');ylabel('幅度'); title('k=0时幅度响应');

subplot(212);plot(w,Hp);

xlabel('\\omega(rad/s)');ylabel('相位'); title('k=0时相位响应');

运行结果：

6

并画出

7

3、已知某离散LTI系统的差分方程为：

1y(n)?y(n?1)?x(n)

3??1?n?1?n?（1）若系统的零状态响应为y(n)?3???????u(n)，求出并画出激励信号

???2??3???x(n)；

程序：

y=3*((1/2).^n-(1/ 3).^n).*uDT(n);

y1=3*((1/2).^(n-1)-(1/ 3).^(n-1)).*uDT(n); x=y-1/3*y1; stem(n,x,'fill');

运行结果： ans =

(3*(1/2)^n)/2 - 4*(1/3)^n + 3*kroneckerDelta(n, 0)

8