资源描述:
《实验2离散时间傅里叶变换》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、文档电子科技大学实验报告学生姓名:项阳学号:2010231060011指导教师:邓建一、实验项目名称:离散时间傅里叶变换二、实验目的:熟悉序列的傅立叶变换、傅立叶变换的性质、连续信号经理想采样后进行重建,加深对时域采样定理的理解。三、实验内容:1.求下列序列的离散时间傅里叶变换(a)(b)2.设画出并观察其周期性。3.设画出并观察其共轭对称性。4.验证离散时间傅里叶变换的线性、时移、频移、反转(翻褶)性质。5.已知连续时间信号为,求:(a)的傅里叶变换;(b)采样频率为5000Hz,绘出,用理想内插函数重建,并对结果进行讨论;
2、(c)采样频率为1000Hz,绘出,用理想内插函数重建,并对结果进行讨论。四、实验原理:1.离散时间傅里叶变换(DTFT)的定义:文档2.周期性:是周期为的函数3.对称性:对于实值序列,是共轭对称函数。4.线性:对于任何,有5.时移6.频移7.反转(翻褶)五、实验器材(设备、元器件):PC机、WindowsXP、MatLab7.1六、实验步骤:本实验要求学生运用MATLAB编程产生一些基本的离散时间信号,并通过MATLAB的几种绘图指令画出这些图形,以加深对相关教学内容的理解,同时也通过这些简单的函数练习了MATLAB的使用。
3、七、实验源代码:文档1.(a)w=[0:1:500]*pi/500;x=exp(j*w)./(exp(j*w)-0.5*ones(1,501));magx=abs(x);angx=angle(x);realx=real(x);imagx=imag(x);subplot(2,2,1);plot(w/pi,magx);gridxlabel('frequencyinpiunits');title('MagnitudePart');ylabel('Magnitude')subplot(2,2,3);plot(w/pi,angx);gr
4、idxlabel('frequencyinpiunits');title('AnglePart');ylabel('Radians')subplot(2,2,2);plot(w/pi,realx);gridxlabel('frequencyinpiunits');title('RealPart');ylabel('Real')subplot(2,2,4);plot(w/pi,imagx);gridxlabel('frequencyinpiunits');title('ImaginaryPart');ylabel('Imagin
5、ary')1.(b)n=-1:3;x=1:5;k=0:500;w=(pi/500)*k;X=x*(exp(-j*pi/500)).^(n'*k);magX=abs(X);angX=angle(X);realX=real(X);imagX=imag(X);subplot(2,2,1);plot(k/500,magX);gridxlabel('frequencyinpiunits');title('magnitudePart')subplot(2,2,3);plot(k/500,angX);gridxlabel('frequenc
6、yinpiunits');title('AnglePart')subplot(2,2,2);plot(k/500,realX);gridxlabel('frequencyinpiunits');title('RealPart')subplot(2,2,4);plot(k/500,imagX);gridxlabel('frequencyinpiunits');title('ImaginaryPart')2.n=0:10;x=(0.9*exp(j*pi/3)).^n;k=-200:200;w=(pi/100)*k;文档X=x*(e
7、xp(-j*pi/100)).^(n'*k);magX=abs(X);angX=angle(X);subplot(2,1,1);plot(w/pi,magX);gridxlabel('frequencyinunitsofpi');ylabel('
8、x
9、')title('MagnitudePart')subplot(2,1,2);plot(w/pi,angX/pi);gridxlabel('frequencyinunitsofpi');ylabel('radians/pi')title('AnglePart')3.subplot
10、(1,1,1)n=-5:5;x=(-0.9).^n;k=-200:200;w=(pi/100)*k;X=x*(exp(-j*pi/100)).^(n'*k);magX=abs(X);angX=angle(X);subplot(2,1,1);plot(w/pi,magX);gr