matlab 离散傅里叶变换_Matlab中的傅里叶变换

两种一般的方法：

• 根据傅里叶变换的定义： ​

​典型程序：

t=0:1000;
xt=sin(t).*exp(-0.02*t); %随时间变化的信号 （时间空间）
subplot(2,1,1);
plot(t,xt); %时间图
N=length(xt); %傅里叶变换公式中的n，N和k
k=0:N-1;
n=0:N-1;
xk=xt*exp(-1j*2*pi/N).^(n'*k); %傅里叶变换公式
subplot(2,1,2);
plot(n,xk); %傅里叶变换后的信号（频率空间）
title('x(k)');

---

• FFT命令实现

典型程序：

clear all
f0 = 100;
fs = 500;
Ts = 1/fs;
n=1:1:1000;
N = length(n);
y = sin(2*pi*f0*n*Ts);
plot(n,y)
y_fft=fft(y);
P2_y_fft =abs(y_fft/N);
P1_y_fft = P2_y_fft(1:N/2+1);
P1_y_fft(2:end-1) = 2*P1_y_fft(2:end-1)
f = fs*(0:N/2)/N;
figure
plot(f,P1_y_fft)
xlabel('f (Hz)')

参考资料：

matlab关于fft的应用（一） - 你好，我是洪大帅~ - CSDN博客​blog.csdn.net