![d2a61ad76080150d9def7eec428f3fb3.png](https://i-blog.csdnimg.cn/blog_migrate/a53134c8e8577cdf120ad4a0da53ab6a.jpeg)
两种一般的方法:
- 根据傅里叶变换的定义:
![fd32540c094865eab76d5e2c47b670de.png](https://i-blog.csdnimg.cn/blog_migrate/7faaafe8faa526b85e50fe3ab26662e6.jpeg)
典型程序:
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)')
参考资料:
![0dca332e3ddcbe7a9701ad9e2e91822b.png](https://i-blog.csdnimg.cn/blog_migrate/2f18109a81b34f0fb96bf479ee66d840.png)