结论:该因果系统的极点全部在单位圆内,故系统是稳定的。
(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