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)');
共6条
1/1 1 跳转至页


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');
一个求平均功率谱的程序,该程序从示波器采样的数据中读取数据。

共6条
1/1 1 跳转至页
回复
有奖活动 | |
---|---|
【EEPW电子工程师创研计划】技术变现通道已开启~ | |
发原创文章 【每月瓜分千元赏金 凭实力攒钱买好礼~】 | |
【EEPW在线】E起听工程师的声音! | |
“我踩过的那些坑”主题活动——第001期 | |
高校联络员开始招募啦!有惊喜!! | |
【工程师专属福利】每天30秒,积分轻松拿!EEPW宠粉打卡计划启动! | |
送您一块开发板,2025年“我要开发板活动”又开始了! | |
打赏了!打赏了!打赏了! |
打赏帖 | |
---|---|
【我踩过的那些坑】STM32的硬件通讯调试过程的“坑”被打赏50分 | |
【我踩过的那些坑】晶振使用的问题被打赏100分 | |
【我踩过的那些坑】电感选型错误导致的处理器连接不上被打赏50分 | |
【我踩过的那些坑】工作那些年踩过的记忆深刻的坑被打赏10分 | |
【我踩过的那些坑】DRC使用位置错误导致的问题被打赏100分 | |
我踩过的那些坑之混合OTL功放与落地音箱被打赏50分 | |
汽车电子中巡航控制系统的使用被打赏10分 | |
【我踩过的那些坑】工作那些年踩过的记忆深刻的坑被打赏100分 | |
分享汽车电子中巡航控制系统知识被打赏10分 | |
分享安全气囊系统的检修注意事项被打赏10分 |