实验三数字信号处理.doc
《实验三数字信号处理.doc》由会员分享,可在线阅读,更多相关《实验三数字信号处理.doc(5页珍藏版)》请在装配图网上搜索。
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)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,3,2);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(0,2,0,6);subplot(3,3,1);plot(w/pi,angle(H),linewidth,2);grid on;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);title(y(n)=x(n)-0.8y(n-1)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,3,5);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(0,2,0,6);subplot(3,3,4);plot(w/pi,angle(H),linewidth,2);grid on;axis(-0.1,2.1,-3,3);xlabel(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)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,3,8);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(0,2,0,6);subplot(3,3,7);plot(w/pi,angle(H),linewidth,2);grid on;axis(-0.1,2.1,-3,3);xlabel(omega/pi);ylabel(phi(omega);title(相频响应特性); 图形如下所示:2. 假设系统用下面差分方程描述: y(n) = x(n) +ax(n-1)假设a=0.7, 0.8, 0.9 , 分别在三种情况下分析系统的频率特性,并打印幅度特性曲线。解: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)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,3,2);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(0,2,0,2);subplot(3,3,1);plot(w/pi,angle(H),linewidth,2);grid on;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)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,3,5);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(0,2,0,2);subplot(3,3,4);plot(w/pi,angle(H),linewidth,2);grid on;axis(-0.1,2.1,-3,3);xlabel(omega/pi);ylabel(phi(omega);title(相频响应特性);B=1,0.9;A=1;subplot(3,3,9);zplane(B,A);xlabel(实部Re);ylabel(虚部Im);title(y(n)=x(n)+0.9x(n-1)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,3,8);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(0,2,0,2);subplot(3,3,7);plot(w/pi,angle(H),linewidth,2);grid on;axis(-0.1,2.1,-3,3);xlabel(omega/pi);ylabel(phi(omega);title(相频响应特性); 图形如下所示:3. 假设系统函数用下式描述: y(n)=1.273y(n-1)-0.81y(n-2)+x(n)+x(n-1)试分析它的频率特性,要求打印其幅度特性曲线,并求出峰值频率和谷值频率。解:B=1,1;A=1,-1.273,0.8;subplot(3,1,3);zplane(B,A);xlabel(实部Re);ylabel(虚部Im);title(y(n)=1.273y(n-1)-0.81y(n-2)+x(n)+x(n-1)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,1,1);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(-0.1,2.1,0,15);subplot(3,1,2);plot(w/pi,angle(H),linewidth,2);grid on;axis(-0.1,2.1,-3,3);xlabel(omega/pi);ylabel(phi(omega);title(相频响应特性); 图形输出如下:4. 假设系统函数用下式描述:(另加的一个) y(n)=0.66y(n-7)-0.77y(n-8)+0.88x(n)-0.99x(n-5)试分析它的频率特性,要求打印其幅度特性曲线。解:B=0.88,0,0,0,0,-0.99;A=1,0,0,0,0,0,0,-0.66,0.77;subplot(3,1,3);zplane(B,A);xlabel(实部Re);ylabel(虚部Im);title( y(n)=0.66y(n-7)-0.77y(n-8)+0.88x(n)-0.99x(n-5)传输函数零、极点分布);grid onH,w=freqz(B,A,whole);subplot(3,1,1);plot(w/pi,abs(H),linewidth,2);grid on;xlabel(omega/pi);ylabel(|H(ejomega)|);title(幅频响应特性);axis(-0.1,2.1,0,15);subplot(3,1,2);plot(w/pi,angle(H),linewidth,2);grid on;axis(-0.1,2.1,-3,3);xlabel(omega/pi);ylabel(phi(omega);title(相频响应特性);- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 数字信号 处理
装配图网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文