按我的理解你只希望PSD情节60赫兹,而不是围绕它的情节你生成你看到一个在光谱中传播。你得到这个传播在光谱频率步(Fs / N)并不是一个合适的你的目标频率因子。在你的情况中一步频率(32000/8192)和利益是60赫兹的频率,如果这是你的问题然后试着改变Fs值32768,在上面的代码,你仍能看到一个更遍布60 Hz,但它会给你正确的大小值。我们看到一个传播因为我们试图计算FFT有限窗口的时间信号,希望这回答了你的查询
策划periodgram pwelch /其他PSD函数模拟数据
2视图(30天)
显示旧的评论
我开发一个1 v peake-to-peak 60赫兹频率的正弦波。我这样做为了groundtruth确定测量声音的方法。我遇到了麻烦在使用和策划开发一个PSD的各种方法。我有下面的代码。运行这个情节罚款正弦波和RMS值的数据点(好!)。我希望看到一个PSD的绘图区,带有1 v 60 hz,没有能源外,乐队。我没有得到,我不知道为什么。供参考,我是一个生物学家——这是所有希腊给我但我尽最大的努力。代码:
% %时间规范:
Fs = 32000;%样本每秒
dt = 1 / f;%秒每样
StopTime = 5;%秒
t = (0: dt: StopTime-dt) ';%秒
% %正弦波:
Fc = 60;%赫兹
x = cos(2 *π* Fc * t);
%画出信号与时间:
图;
情节(t, x);
包含(的时间(以秒为单位));
标题(“信号与时间”);
NFFT = 8192;%
NOVERLAP =圆(0.75 * NFFT);
w =汉宁(NFFT);
[sensor_spectrum,频率]= pwelch (x, w, NOVERLAP NFFT, fs);
图;
情节(频率、sensor_spectrum);
