这些小活动你都参加了吗?快来围观一下吧!>>
电子产品世界 » 论坛首页 » 嵌入式开发 » MCU » 使用FFT进行谱分析的MATLAB源程序

共6条 1/1 1 跳转至

使用FFT进行谱分析的MATLAB源程序

工程师
2011-06-03 11:55:57     打赏
    t=0:1:65536;
    Fs=1e6;
    Data=cos(2*pi*200e3*t/Fs);
    FFT_L = 2^(floor(log2(length(Data))));
    FFT_Data = Data(1:FFT_L);
    FFT_Data_Wn = (hanning(FFT_L))';
    FFT_Out = fft(FFT_Data_Wn.*FFT_Data,FFT_L);
    figure;
    fig_x = linspace(-Fs/2,Fs/2,FFT_L);
    fig_y = 10*log10(abs([FFT_Out(FFT_L/2+1:end),FFT_Out(1:FFT_L/2)]/FFT_L).^2);
    plot(fig_x,fig_y,'b');
    axis([-Fs/2,Fs/2,-100,10]);
    title('PSD of input Data');
    xlabel('Frequency');
    ylabel('Power(dB)');



关键词: 使用     进行     分析     MATLAB     源程序    

专家
2011-06-03 17:24:11     打赏
2楼
不错,表扬一个!

高工
2011-11-17 19:19:52     打赏
3楼

filename = “示波器数据”
Fs = 1/csvread(filename,6,1,[6,1,6,1]);
data_all =  csvread(filename,16,1,[16,1,8192/2*18,1]);
average_num = 16;
N = 8192;
sp = zeros(N,1);
fre =[0:N-1]*Fs/N;

for i = 1:1:average_num                                                         
    sp = sp + abs(fft(data_all(N/2*(i - 1)+1:N/2*(i + 1)),N)).*2/N;
end
sp = sp/average_num;
plot(fre(1:floor(length(fre)/2)),sp(1:floor(length(sp)/2)),'k');


一个求平均功率谱的程序,该程序从示波器采样的数据中读取数据。


菜鸟
2012-02-13 20:53:54     打赏
4楼
赞一个!!!!!!

菜鸟
2012-04-19 04:59:10     打赏
5楼

赞~谢谢分享


助工
2012-05-29 21:08:01     打赏
6楼

不用matlab好多年了,不过还能看懂一些,楼主的程序大概就是生成一个汉宁窗函数(正弦吧?)然后做傅里叶出频谱,再转成分贝,画成曲线,有点不明白就是转分贝的时候为什么要平方一下啊?


共6条 1/1 1 跳转至

回复

匿名不能发帖!请先 [ 登陆 注册 ]