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网站征稿正在进行时,欢迎踊跃投稿啦 | |
【有奖活动】智能可穿戴设备AR/VR如何引领科技新潮流! | |
奖!发布技术笔记,技术评测贴换取您心仪的礼品 | |
【有奖活动】震撼来袭!这场直播将直击工程师的心灵! |
打赏帖 | |
---|---|
如何实现基于NXPiMX.RT1021的BH1730采集被打赏50分 | |
【换取手持数字示波器】AHT10温度检测分享被打赏40分 | |
【换取手持数字示波器】ACM32F070开发板点亮LCD屏和触控按键,串口,ADC被打赏40分 | |
【换取手持数字示波器】ACM32F070LCD屏和触控按键功驱动蜂鸣器分享被打赏40分 | |
【换取手持数字示波器】国民技术PWM功能知识分享被打赏40分 | |
【换取手持数字示波器】放大器运放知识分享被打赏40分 | |
【分享评测,赢取加热台】+开关电源AC输入知识分享被打赏30分 | |
【换取手持数字示波器】N32G430点亮数码管与串口打印、ADC被打赏40分 | |
老胖子聊电路--分析一个反馈电路被打赏50分 | |
【换取手持数字示波器】+点亮WS2812RGB灯被打赏50分 |