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 跳转至页
回复
| 打赏帖 | |
|---|---|
| 【S32K146】S32DS watchdog 配置使用被打赏20分 | |
| 【Zephyr】使用 IAR 调试 Zephyr 镜像被打赏20分 | |
| 【Zephyr】MCXN947 Zephyr 开发入门适配shell被打赏20分 | |
| 【我要开发板】6.联合MATLAB记录数据被打赏50分 | |
| 【瑞萨RA2E1开发板】：使用ADC功能实现位移传感器采集方案被打赏20分 | |
| 【nRF7002DK】基于sht30的温湿度计被打赏20分 | |
| 【nRF7002DK】日志打印被打赏20分 | |
| 【换取手持示波器】RGB屏幕移植ARM-2D库被打赏35分 | |
| 【分享开发笔记，赚取电动螺丝刀】分享一下如何解决瑞萨RA2E1使用printf编译报错问题被打赏27分 | |
| rtthread硬件加密-5hash加密分析被打赏10分 | |
STM32
MCU
通讯及无线技术
物联网技术
电子DIY
板卡试用
基础知识
软件与操作系统
我爱生活
小e食堂