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

《实验三 离散信号与系统的连续频域分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验三离散信号与系统的连续频域分析一、实验目的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版本以上）四

2、、设计流程查找离散时间系统信号的幅度和相位函数→通过MATLAB帮助阅读函数的使用→编写程序→在MATLAB上调试→书写实验报告。五、程序清单参考程序附后。六、要求1．通过参考程序进行仿真，并理解程序；2．对重要语句进行解释，附在程序行后面；3．理解函数的含义及参数所表示的意义。本实验中如freqz函数、abs函数和angle函数。参考程序：1n=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);

3、subplot(2,1,1);plot(w/pi,magX);xlabel('w/pi');ylabel('幅度

4、X

5、');subplot(2,1,2);plot(w/pi,angX);xlabel('w/pi');ylabel('相位(rad/π)');2a=[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=angl

6、e(H);H1=freqz(b,a,w);magH1=abs(H1);angH1=angle(H1);subplot(2,2,2);plot(w/pi,angH/pi);grid;xlabel('w(frequencyinpiunits)');ylabel('相位(rad/π)');subplot(2,2,1);plot(w/pi,magH);grid;xlabel('w(frequencyinpiunits)');ylabel('幅度

7、H

8、');subplot(2,2,3);plot(w/pi,magH1);grid;xlabel('

9、w(frequencyinpiunits)');ylabel('幅度

10、H1

11、');subplot(2,2,4);plot(w/pi,angH1/pi);grid;xlabel('w(frequencyinpiunits)');ylabel('相位(rad/π)');axis([0,1,-0.8,0]);figure(2);zplane(b,a);