dft对称性 matlab实验,数字信号处理实验13MATLAB

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

一、实验目的

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

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

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

二、实验原理

1、线性性质

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

y(n)=ax1(n)+b x2(n) (a,b均为常数)

则该y(n)的N点DFT为

Y(k)=DFT[y(n)]=aX1(k)+bX2(k) 0≤k≤N-1

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

例13-1 已知x1(n) =[0,1,2,4]、x2(n)=[1,0,1,0,1]，求：

(1)y(n)= 2x1(n)+3x2(n),再由y(n)的N点DFT获得Y(k);

(2) 由x1(n)和x2(n)求X1(k)和X2(k)，再求)Y(k)= 2X1(k)+3X2(k)。

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

解MATLAB程序如下：

xn1=[0,1,2,4]; %建立xn1序列

xn2=[1,0,1,0,1]; %建立xn2序列

N1=length(xn1);

N2=length(xn2);

N=max(N1,N2); %确定N

if N1>N2 xn2=[xn2,zeros(1,N1-N2)]; %对长度短的序列补0

elseif N2>N1 xn1=[xn1,zeros(1,N1-N2)];

end

yn=2*xn1+3*xn2; %计算yn

n=0:N-1;

k=0:N-1;

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

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

Yk2=xn2*(exp(-j*2*pi/N)).^(n’*k);

Yk2=2*Xk1+3*Xk2; %由Xk1，Xk2求Yk

以上程序作图部分省略。

用两种方法求得Y(k)结果一致，如下所所示：

Yk=

23.0000 -7.5902+1.5388i 3.5902-0.3633i

3.5902+0.3633i -7.5902-1.5388i

运行结果如图13-1所示。