实验三数字信号处理

《实验三数字信号处理》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、07电信35号林烟泉1.假设系统用下面差分方程描述：y(n)=x(n)+ay(n-1)假设a=0.7,0.8,0.9，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。解：B=1;A=[1,-0.7];subplot(3,3,3);zplane(B,A);xlabel('实部Re→');ylabel('虚部Im→');title('y(n)=x(n)-0.7y(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');subplot(3,3,2);plot(w/pi,abs(H),'linewidth',2);gridon;

2、xlabel('omega/pi');ylabel('

3、H(e^j^omega)

4、');title('幅频响应特性');axis([0,2,0,6]);subplot(3,3,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('omega/pi');ylabel('phi(omega)');title('相频响应特性');B=1;A=[1,-0.8];subplot(3,3,6);zplane(B,A);xlabel('实部Re→');ylabel('虚部Im

5、→');title('y(n)=x(n)-0.8y(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');subplot(3,3,5);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('

6、H(e^j^omega)

7、');title('幅频响应特性');axis([0,2,0,6]);subplot(3,3,4);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xl

8、abel('omega/pi');ylabel('phi(omega)');title('相频响应特性');B=1;A=[1,-0.9];subplot(3,3,9);zplane(B,A);xlabel('实部Re→');ylabel('虚部Im→');title('y(n)=x(n)-0.9y(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');subplot(3,3,8);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('

9、H

10、(e^j^omega)

11、');title('幅频响应特性');axis([0,2,0,6]);subplot(3,3,7);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('omega/pi');ylabel('phi(omega)');title('相频响应特性');图形如下所示：07电信35号林烟泉07电信35号林烟泉2.假设系统用下面差分方程描述：y(n)=x(n)+ax(n-1)假设a=0.7,0.8,0.9，分别在三种情况下分析系统的频率特性，并打印幅度特性

12、曲线。解：B=[1,0.7];A=1;subplot(3,3,3);zplane(B,A);xlabel('实部Re→');ylabel('虚部Im→');title('y(n)=x(n)+0.7x(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');subplot(3,3,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('

13、H(e^j^omega)

14、');title('幅频响应特性');axis([0,2,0,2]);subplo

15、t(3,3,1);plot(w/pi,angle(H),'linewidth',2);gridon;axis([-0.1,2.1,-3,3]);xlabel('omega/pi');ylabel('phi(omega)');title('相频响应特性');B=[1,0.8];A=1;subplot(3,3,6);zplane(B,A);xlabel('实部Re→');ylabel('虚部Im→');title('y(n)=x(n)+0.8x(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');sub