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网站征稿正在进行时，欢迎踊跃投稿啦 | |
| 奖！发布技术笔记，技术评测贴换取您心仪的礼品 | |
| 打赏了！打赏了！打赏了！ | |
| 打赏帖 | |
|---|---|
| 与电子爱好者谈读图二被打赏50分 | |
| 【FRDM-MCXN947评测】Core1适配运行FreeRtos被打赏50分 | |
| 【FRDM-MCXN947评测】双核调试被打赏50分 | |
| 【CPKCORRA8D1B评测】---移植CoreMark被打赏50分 | |
| 【CPKCORRA8D1B评测】---打开硬件定时器被打赏50分 | |
| 【FRDM-MCXA156评测】4、CAN loopback模式测试被打赏50分 | |
| 【CPKcorRA8D1评测】--搭建初始环境被打赏50分 | |
| 【FRDM-MCXA156评测】3、使用FlexIO模拟UART被打赏50分 | |
| 【FRDM-MCXA156评测】2、rt-thread MCXA156 BSP制作被打赏50分 | |
| 【FRDM-MCXN947评测】核间通信MUTEX被打赏50分 | |
STM32
MCU
通讯及无线技术
物联网技术
电子DIY
板卡试用
基础知识
软件与操作系统
我爱生活
小e食堂