对单频信号做DFT处理matlab,数字信号处理DFTMATLAB程序

实验三 频域信号处理

1. 实验目的

(1) 学习信号DFT 变换的matlab 实现；

(2) 学习fft 的matlab 实现；

(3) 验证DFT 的相关性质。

2. 思考题

(1) 若()()()sin 8sin 4x n n n ππ=+是一个128点的有限长序列，求其128点DFT

结果；

程序如下：

求DFT 变换矩阵A ：

clc;

clear;

N=128;

A=dftmtx(N)

Ai=conj(dftmtx(N));

n=0:(N-1);

k=0:(N-1);

nk=n'*k;

Wn=(sin(pi/8)+sin(pi/4)).^nk

Wk=conj(Wn)/N;

求128点的DFT (分别用FFT 函数和dftmtx 函数)

clc;

clear;

N=128;

n=0:N-1;

x=sin(pi/8*n)+sin(pi/4*n);

subplot(3,1,1)

plot(n,x);

grid on

title('原图')

y1=fft(x,N);

A=dftmtx(N);

y2=x(1:N)*A;

subplot(3,1,2)

plot(n,y1)

grid on

title('FFT')