https://blog.csdn.net/weixin_44225182/article/details/102488649
dt=0.02;%采样间隔
N=512;%采样点数
t=0:dt:(N-1)*dt;%采样时刻
fs=1/dt;%采样频率,与才采样间隔互为倒数
n=0:1:N-1;
f=(fs/N).*n;%X轴每个点对应的频率
x=0.5*sin(2*pi*3*t)+cos(2*pi*10*t);%信号
figure(1)
plot(t,x)
y=fft(x);%傅里叶变换得到一个复数
Ay=abs(y);%取模
Ayy=Ay*2/N;%转换成实际的幅值
figure(2)
plot(f(1:N/2),Ayy(1:N/2))
f1=8;
f2=15;
yy=zeros(1,length(y));
for m=0:N-1
if(m*(fs/N)>f1&m*(fs/N)&&(fs-f2)&m*(fs/N)<(fs-f1));%将奈奎斯特之后的频率也滤除点掉
yy(m+1)=0;
else
yy(m+1)=y(m+1);
end
end %将频率为8Hz-12Hz的信号的幅值置0
yyi=abs(yy);
figure(3)
plot(f(1:N/2),yyi(1:N/2))
yi=ifft(yy);
figure(4)
plot(t,real(yi))
如何理解离散傅里叶变换https://blog.csdn.net/xfgryujk/article/details/78906502