资源描述:
《离散傅立叶变换的性质及应用》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、贵州大学实验报告姓名学号实验组实验时间2015.11指导教师成缋实验项0名称实验一:离散傅立叶变换的性质及应用实验i的1、了解DFT的性质及其应用。2、熟悉MATLAB编程特点。实验内容二、实验内容1、用三种不同的DFT程序实现一维数字信号的傅立叶变换。2、利用DFT实现两序列的卷积运©,并研究DFT点数与混叠的关系。3、研究高密度频谱与高分辨率频谱。4、序列的内插和抽取时所对应的傅立叶变换。实验原理KDFT5換正变換:AM反变换:2,序列卷积设序列4^(?^的长度为N、序列的长度为M。则分别对两个咧作£>点的on■得對⑷和则两序例
2、的线1注卷积;等于瓜(叫。即时域卷积频域为相莱关系。实验仪器计算机一台;Matlab软件学院:大数据与信息工程学院专业:通信工程班级:通信132六、实验步豫1、用三种不同的OFT程序计算x(«)=Jes(«)的傅立叶变换XOt),并比较三种程序计算机的运行时间。(1)««用以循环语句的M函数文件dftkm,用褊环变屋逐点计算X(幻;(2)»写用MATLAB矩阵运箅的M函数文件dft2.m,完成下列矩阵运貝:X(Q)X()x(0)x(l)(3)调用FFT库函数,吉接计算T(幻;(4)分别利用上述三种不同方式编写的DFT程序计算序列x
3、OO的傅立叶变换JTOt),并画出相应的幅频和相频特性,再比较各个程序的计算机运行时间。2、利用DFT实现两序列的卷积运S,并研究DFT点數与混©的关系。给定*(«)=«‘(>:),利用圆周卷积计見线性卷积(快速卷积,頻域);计莒不同DFT点數下的圆周卷积箱出:并用函敦5tem(n,y)画出相应图形。选择不同的DFT点数进行对比,观禁其混瞾效应。3、研究高密度频谱与高分辨率频谱设有连续信号y/Z)=cos(2^x6.5x103f)+cos(2^x7x103r)+cos(2>tx9x10"Z)以采样频率乂=32妍2对该信号采样,分析下
4、列三种恼况的幅频特性。(1)采策数据长度、V=16点,做.V=16点的DFT,并画出幅频特性。(2)采集数据长度,V=16点,补芩到256点的DFT,并画出幅频特性(3)采集数据长度,V=256点,做#=256点的DFT,并画出幅频特性观察三幅不同频率特性图,分析和比较它们的特点以及形成的原囚。4、实现序列的内插和抽取所对应的傅立叶变换。给定序列x⑻=
5、cos^^?/j+co<⑻,做128点的W立叶变换,并求0n荠4k//=4k为整数六、实验步豫1、用三种不同的OFT程序计算x(«)=Jes(«)的傅立叶变换XOt),并比较三种程序
6、计算机的运行时间。(1)««用以循环语句的M函数文件dftkm,用褊环变屋逐点计算X(幻;(2)»写用MATLAB矩阵运箅的M函数文件dft2.m,完成下列矩阵运貝:X(Q)X()x(0)x(l)(3)调用FFT库函数,吉接计算T(幻;(4)分别利用上述三种不同方式编写的DFT程序计算序列xOO的傅立叶变换JTOt),并画出相应的幅频和相频特性,再比较各个程序的计算机运行时间。2、利用DFT实现两序列的卷积运S,并研究DFT点數与混©的关系。给定*(«)=«‘(>:),利用圆周卷积计見线性卷积(快速卷积,頻域);计莒不同DFT点數
7、下的圆周卷积箱出:并用函敦5tem(n,y)画出相应图形。选择不同的DFT点数进行对比,观禁其混瞾效应。3、研究高密度频谱与高分辨率频谱设有连续信号y/Z)=cos(2^x6.5x103f)+cos(2^x7x103r)+cos(2>tx9x10"Z)以采样频率乂=32妍2对该信号采样,分析下列三种恼况的幅频特性。(1)采策数据长度、V=16点,做.V=16点的DFT,并画出幅频特性。(2)采集数据长度,V=16点,补芩到256点的DFT,并画出幅频特性(3)采集数据长度,V=256点,做#=256点的DFT,并画出幅频特性观察三幅
8、不同频率特性图,分析和比较它们的特点以及形成的原囚。4、实现序列的内插和抽取所对应的傅立叶变换。给定序列x⑻=
9、cos^^?/j+co<⑻,做128点的W立叶变换,并求0n荠4k//=4k为整数实骑代码Dftl:function[Am,pha]=dftl(x)N=length(x);w=exp(-j*2*3.1415926/N);fork=l:Nsum=0;forn=l:Nsum=sum+x(n)*wA((k一1)*(n-1));endAm(k)=abs(sum〉;pha(k)=angle(sum);endenddft2:funct
10、ion[Am,pha]=dft2(x)N=length(x);n=[0:N-1];k=[0:N-l];w=exp(-j*2*3.1415926/N);nk=n1*k;wnk=w•(nk);Xk=x*wnk;Am=abs(Xk);pha