离散傅里叶变换

存个档，以后就用这版吧

fs = 1000;            %采样频率
L= 20000;           %数据长度
t = (0:L-1)/fs                  %时间
y = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);  %采样频率>2倍信号最大频率
Y= fft(y);%默认做L点傅里叶变换
Y0= abs(Y)/L; %双边谱，做L点傅里叶变换，频率幅值要除以L             
Y1= Y0(1:L/2); %单边谱，关于L/2对称
Y0(2:end) = 2*Y0(2:end);   %1是直流
Y1(2:end) = 2*Y1(2:end);
F0= fs*(0:(L-1))/L;%频谱分辨率，共有L个点，每两点之间的频率是fs/L
F1= fs*(0:(L/2)-1)/L;
subplot(2,1,1)
plot(F0,Y0)
title('双边频谱')
xlabel('f(Hz)')
ylabel('Y0')
subplot(2,1,2)
plot(F1,Y1)
title('单边频谱')
xlabel('f(Hz)')
ylabel('Y1')
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容