MATLAB的答案

为什么周期图随数据长度?

3视图(30天)
维尼
维尼 2013年10月14日
回答: 维尼2013年10月14日
你好,
我想更好地理解使用用纯正弦信号周期图是如何工作的。
Fs = 1000;
t = 0:1 / Fs: 1 - 1 / f;
x =罪(2 *π* 100 * t);
[psdestx, Fxx] =周期图(x,[],长度(x), Fs);
情节(Fxx psdestx);网格在;
包含(“赫兹”);
标题(的周期图估计功率谱密度);
我把最上面的代码的MATLAB示例文档中心: //www.tianjin-qmedu.com/help/signal/ug/psd-estimate-using-fft.html
0.5响应,峰值在100赫兹,我看起来是正确的,因为理论正弦信号的平均功率(^ 2)/ 2。这也是同样的答案时,我得到类型:
意思是(x ^ 2)。
我怀疑出现增加数据长度,从代码的第二行:
t = 0:1 / Fs: 2 - 1 / f;
为什么这个改变我的PSD估计?因为我处理一个周期信号时,不应该平均功率保持不变,尽管数据长度?
谢谢你,维尼

在回答这个问题。

接受的答案

韦恩王
韦恩王 2013年10月14日
编辑:韦恩王 2013年10月14日
因为周期图严格PSD估计,除非你指定“权力”选项。增加一倍的时间间隔你积分区间的两倍。
如果指定“权力”的选项,你会得到你所期望的0.5
Fs = 1000;
t = 0:1 / Fs: 2 - 1 / f;
x =罪(2 *π* 100 * t);
[psdestx, Fxx] =周期图(x,[],长度(x), Fs,“权力”);
情节(Fxx psdestx);网格在;
包含(“赫兹”);
标题(“权力周期图估计”);

更多的答案(3)

维尼
维尼 2013年10月14日
谢谢你的回复韦恩。
“权力”选项不工作在我的版本的MATLAB (2012 a),但我可以取代它的“女士”。你知道如果它是一样的吗?我看不出它所描述的任何文档。
也,你能澄清价值观的意义我认为情节与缺省选项?例如,假设我看到0.5的峰值在100赫兹和2的峰值在200赫兹,否则“干净”的周期图。我认为这意味着,重建原始信号,我需要两个正弦曲线,与100和200赫兹的频率,每个平均功率为0.5和2,分别调整(+阶段)。那是不正确的吗?
最后,是默认选项的情节更有意义,或者更多的使用?,为什么?
韦恩王
韦恩王 2013年10月14日
在R2012a,您可以使用频谱。周期图然后msspectrum
Fs = 1000;
t = 0:1 / Fs: 2 - 1 / f;
x = 0.5 * cos(2 *π* 100 * t) + 2 *罪(2 *π* 200 * t);
惠普= spectrum.periodogram;
惠普。Windowname =“平顶”;
hms = msspectrum(惠普x,“Fs”Fs);
情节(hms.Frequencies hms.Data)
看情节(这是一个单边功率谱),intepretation是这样的,在100赫兹是0.125所以cos(忽略阶段)的振幅
sqrt (2 * 0.125)
功率估计在200赫兹是2有一个正弦信号的振幅
sqrt (2 * 2)
维尼
维尼 2013年10月14日
太好了,这正是我的理解!
我仍然不明白情节与默认选项是什么意思。为什么是默认?第二个选项似乎更符合逻辑,很容易解释。

类别

找到更多的在谱估计在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!