实验九离散信号分析电子信息.pdfVIP

实验九离散信号的分析

一、实验目的

1.学习离散信号的分析方法；

2.加深对离散傅氏变换的理解和认识。

二、实验内容

1.自行编写函数构建成为为的离散信号的DFT(DiscreteFourierTransform)矩J

2

阵，并与MATLAB自带函数dftmtx产生的DFT矩阵对比。

EditorWindow

clc,clear,closeall;

%%通过调用自定义函数和系统自带函数产生DFT矩阵

N=10;

WNnk1=dftmtx(2^N);

WNnk2=dft1(N);

%%调用imagesc对比两个结果

figure

subplot(121)

imagesc(real(WNnk1));colorbar;title(系统自带函数的DFT矩阵实

部)

subplot(122)

imagesc(real(WNnk2));colorbar;title(自编写函数的DFT矩阵实部

)

figure

subplot(121)

imagesc(imag(WNnk1));colorbar;title(系统自带函数的DFT矩阵虚

部)

subplot(122)

imagesc(imag(WNnk2));colorbar;title(自编写函数的DFT矩阵虚部

)

自编写函数：

functionoutput=dft1(m)

N=2^m;

n=[0:1:N-1];%矩阵行

k=[0:1:N-1];%矩阵列

WN=exp(-j*2*pi/N);

WNnk=WN.^(n*k);

output=WNnk;





2.利用1中的DFT矩阵，计算ftsintut10ut10在t1时得



10



到的等间隔采样序列fn的变换Fk，并与MATLAB自带函数fft得到



的变换结果对比，画出各自的幅度谱和相位谱。

EditorWindow

clc,clear,closeall;

%%参数

dt=1;

t_min=-10;

t_max=10;

t=t_min:dt:t_max;

n=length(t);

%%得到序列并计算

f_n=sin(pi/10*t).*(u(t+10)-u(t-10));

mtx=dft2(n);

FK_fn=dft2(n)