数字信号处理实验-matlab版-离散傅里叶变换的性质.pdfVIP

实验13离散傅里叶变换的性质

(完美格式版，本人自己完成，所有语句正确，不排极个别错误，特别适用于山大，勿用冰点等

工具下载，否则下载之后的word格式会让很多部分格式错误，谢谢)

XXXX学号姓名处XXXX

一、实验目的

1加深对离散傅里叶变换(DFT)基本性质的理解。

2了解有限长序列傅里叶变换(DFT)性质的研究方法c

3掌握用MATLAB语言进行离散傅里叶变换性质分析时程序编写的方法。

二、实验内容

1线性性质。

2循环移位性质。

3循环折叠性质。

4时域和频域循环卷积特性。

5循环对称性。

三、实验环境

MATLAB7.0

四、实验原理

1线性性质

如果两个有限长序列分别为xl(n)和x2(n),长度分别为N1和N2,且

y(n)=axl(n)+bx2(n)(a、b均为常数)

则该y(n)的N点DFT为

Y(k)-DFTCy(n)]-aXl(k)+bX2(k)OWkWN-1

其中：N=max[Nl,N2],Xl(k)和X2(k)分别为xl(n)和x2(n)的N点DFT。

例13・1已知x1(n)=[0,1,2,4],x2(n)=[1,0,I,0,1],求：

(l)y(n)=2xl(n)+3x2(n),再由y(n)的N点DFT获得Y(k)；

(2)由xl(n)、x2(n)求Xl(k)、X2(k),再求Y(k)=2Xl(k)+3X2(k)。

用图形分别表示以上结果，将两种方法求得的Y(k)进行比较，由此验证有限长序列傅里叶变

换(DFT)的线性性质。

解MATLAB程序如下：

»xnl=[0,l,2,4J；%建立xnl序列

»xn2=[l,0JAH；%建立xn2序列

»N1=lcngth(xnl);N2=lcngth(xn2);

»N=max(NI,N2);%确定N

»ifNlN2xn2=[xn2,zeros(l,Nl-N2)];%对长度短的序列补0

»elseifN2NIxnl=[xnl,zeros(l,N2-Nl)];

»end

»yn=2*xnl+3*xn2;%计算yn

»n=0:N-l;k=0:N-l;

A

»Yk1=yn*(exp(-j*2*pi/N)).(n**k);%求yn的N点DFT

A

»Xk1=xn1*(exp(-i*2*pi/N)).(n**k);%求xnl的N点DFT

A

»Xk2=xn2*(exp(-j*2*pi/N)).(n*k);%求xn2的N点DFT

»Yk2=2*Xkl+3*Xk2;%©Xkl、Xk2求Yk

»subplot(4,2,1),stem(n,