MATLAB的答案

复制“pspectrum”使用的周期图

13个视图(30天)
埃德温
埃德温 2023年1月12日
编辑: Nayan 2023年4月13日
你好,
我试图理解pspectrum是如何工作的。阅读文档之后,我想我知道它使用kaiser窗和发展以下代码来繁殖的结果。但不幸的是我还似乎少了什么。
同时,好奇为什么权力的大小也是不同的吗?像一个~ 500倍。任何见解也会有帮助。
非常感谢看到的准确代码如何复制和其他预处理需要匹配两种方法。提前谢谢你!
欢呼,
埃德温
%方法1,使用pspectrum的功率谱。
[p f] = pspectrum (x, fs,“漏”,0.5);
%方法2,使用周期图的功率谱
路= 0.5;
M =长度(x);
g = kaiser (M, 40 * (1-lk));
[p f] =周期图(x, g,长度(x), fs);
1评论
埃德温
埃德温 2023年1月13日
调查后的代码“pspectrum”我相信我现在看到它是如何通过Matlab计算能力 信号分析仪 。
文档信号分析仪说:
  1. 把信号分成等长段。段必须足够短,信号的频率含量在一段没有明显变化。段可能会或可能不会重叠。
  2. 窗口每一部分和计算得到的光谱短时傅里叶变换。
  3. 显示segment-by-segment每个光谱分贝的力量。描绘了并排大小作为一个形象magnitude-dependent colormap。
现在这个有意义等“pspectrum”看到其他选项“NumPowerBins”和“OverlapPercent”。所以,因此,它似乎不能使用周期图复制pspectrum但需要使用STFT的方法。可以找到其他有用的信息 在这里 。
干杯

登录置评。

在回答这个问题。

答案(1)

Nayan
Nayan 2023年4月13日
编辑:Nayan 2023年4月13日
跑:
就像我 理解它 你想 理解 ”的工作 pspectrum (x) ”。我建议你 要注意 为更好地了解以下几点。
计算信号光谱, pspectrum 发现光谱分辨率之间的妥协实现信号的整个长度和性能的局限性而导致的计算大型fft算法:
  • 如果可能的话,这个函数计算单个修改整个信号的周期图使用Kaiser窗。
  • 如果不可能计算一个改进的周期图在合理的时间内,该函数计算韦尔奇周期图:它将信号划分为重叠部分,windows使用Kaiser窗每一部分,和平均周期图的部分。
因此,有两种不同的技术 pspectrum 根据评估情况吗 的手。我 建议你 去 通过下面的MATLAB文档链接 一个 详细的了解 //www.tianjin-qmedu.com/help/releases/R2023a/signal/ref/pspectrum.html?s_tid=doc_ta d124e149267
我建议你 分享 功率谱的确切数据你正在评估使用“ pspectrum(x) ”和“ 周期图(x窗口,f, f) ” 比较 (因为该方法用于评价的权力 范围 很大程度上取决于数据)
你可以通过以下链接了解” 周期图(x窗口,f, f) ”, 它的 工作。
//www.tianjin-qmedu.com/help/releases/R2023a/signal/ref/periodogram.html?s_tid=doc_ta btuc2ow-1
还发现的代码片段,产生相同的功率谱 方法 。
%生成一个测试信号采样频率100赫兹和长度为1000
Fs = 100;%采样频率(赫兹)
t = (0:999) / Fs;%的时间向量(sec)
x = cos(2 *π* 10 * t) + 0.5 *罪(2 *π* 20 * t) + randn (1000);%测试信号
%计算PSD使用周期图方法
[pxx f] =周期图(x [] [], Fs);
%可视化PSD
情节(f, 10 * log10 (pxx))
包含(的频率(赫兹))
ylabel (“功率/频率(dB / Hz)”)
标题(“周期图PSD估计”)
%放大显示只有50赫兹频率
xlim (50 [0])
%生成一个测试信号采样频率100赫兹和长度为1000
Fs = 100;%采样频率(赫兹)
t = (0:999) / Fs;%的时间向量(sec)
x = cos(2 *π* 10 * t) + 0.5 *罪(2 *π* 20 * t) + randn (1000);%测试信号
%计算PSD使用“pspectrum”功能和默认设置
[pxx f] = pspectrum (x, Fs);
%可视化PSD
情节(f, 10 * log10 (pxx))
包含(的频率(赫兹))
ylabel (“功率/频率(dB / Hz)”)
标题(“使用pspectrum PSD估计”)
%放大显示只有50赫兹频率
xlim (50 [0])
希望这可以帮助!

类别

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

标签

社区寻宝

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

开始狩猎!