EEPW论坛

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');


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

赞~谢谢分享

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