实验三 离散信号与系统的连续频域分析.docVIP

实验三 离散信号与系统的连续频域分析 一、实验目的 1．离散时间信号的DTFT的MATLAB实现； 2．进行离散时间系统的DTFT分析； 3．理解系统函数和频率相应之间的关系。 二、实验内容 1．自定义一个长度为8点的信号，信号幅度值也由自己任意指定，对该信号作DTFT，分别画出幅度谱和相位谱； 2．已知离散时间系统差分方程为y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求出并画出其频率响应； 3．求该系统系统函数，并画极零点图，并通过freqz函数求频率响应。 三、实验平台 MATLAB集成系统（MATLAB6.5版本以上） 四、设计流程 查找离散时间系统信号的幅度和相位函数→通过MATLAB帮助阅读函数的使用→编写程序→在MATLAB上调试→书写实验报告。 五、程序清单 参考程序附后。 六、要求 1通过参考程序进行仿真，并理解程序； 2对重要语句进行解释，附在程序行后面； 3理解函数的含义及参数所表示的意义。本实验中如freqz函数、abs函数和 angle函数。 参考程序： 1 n=0:7; x=(0.9*exp(j*pi/3)).^n; w=-pi:pi/200:pi; X=x*(exp(-j*pi/4)).^(n*w); magX=abs(X);angX=angle(X); subplot(2,1,1);plot(w/pi,magX);xlabel(w/pi);ylabel(幅度|X|); subplot(2,1,2);plot(w/pi,angX);xlabel(w/pi);ylabel(相位(rad/π)); 2 a=[1,-0.5,0.06];b=[1,1,0]; m=0:length(b)-1;l=0:length(a)-1; w=0:pi/500:pi; num=b*exp(-j*m*w); den=a*exp(-j*l*w); H=num./den; magH=abs(H);angH=angle(H); H1=freqz(b,a,w); magH1=abs(H1);angH1=angle(H1); subplot(2,2,2);plot(w/pi,angH/pi);grid; xlabel(w (frequency in pi units));ylabel(相位(rad/π)); subplot(2,2,1);plot(w/pi,magH);grid;xlabel(w (frequency in pi units));ylabel(幅度|H|); subplot(2,2,3);plot(w/pi,magH1);grid;xlabel(w (frequency in pi units));ylabel(幅度|H1|); subplot(2,2,4);plot(w/pi,angH1/pi);grid; xlabel(w (frequency in pi units));ylabel(相位(rad/π));axis([0,1,-0.8,0]); figure(2);zplane(b,a);