MATLAB功率谱函数psd与pwelch的使用

博客探讨了在MATLAB中使用psd和pwelch函数计算功率谱密度时出现的30dB差异问题。原因是psd函数在新版本中需要结合spectrum使用,而pwelch直接计算功率谱密度。通过调整psd函数的结果,可以使其与pwelch得到相同的结果。此外,文章还对比了两个函数的用法,并提供了等价的代码示例,展示两者在功率谱计算上的等效性。
摘要由CSDN通过智能技术生成

目录

 

一、遇到的问题

二、功率谱函数psd与pwelch的用法与区别


一、遇到的问题

在网上看到这样一个问题,然后搬过来记一下。

大家好,我用matlab的psd函数和pwelch函数求功率谱密度,但是得到的结果为什么下差大约30db,这个是怎么回事?

摘自:https://www.ilovematlab.cn/thread-339674-1-1.html

clear;
Fs=1000;
n=0:1/Fs:1;
y=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
Nfft=1000;
window=hanning(100);   
noverlap=20;   
dflag='none';  
[Pxx,f]=psd(y,Nfft,Fs,window,noverlap);   
figure
plot(f,10*log10(Pxx));  
xlabel('frequency/Hz');ylabel('PSD/dB');
title('PSD-Welch method')
hold on;
[Pxx1,f1]=pwelch(y,window,noverlap,N
MATLAB中的功率谱密度函数有多种实现方式。在较早版本中,可以使用PSD函数来计算功率谱密度。该函数具有以下特点:(1)功率谱密度函数P(ω)是实数;(2)功率谱函数是非负的,即P(ω)≥0;(3)功率谱密度函数的逆傅里叶变换是信号的自相关函数C(τ);(4)功率谱密度对频率的积分给出信号f(t)的方差。\[1\] 然而,在现在的新版本中,PSD函数的用法已经改变了,它结合了spectrum函数一起使用。在较早版本中,PSD函数得到的是功率谱,而不是功率谱密度,因此与pwelch函数计算的功率谱密度结果会有30dB的差别。\[2\] 另外,还有另一种常用的求功率谱密度的函数pwelchpwelch函数可以直接计算信号的功率谱密度,并且结果可以直接使用。它的等价关系可以表示为:\[px,fx\]=psd(x,nfft,fs,Windows,noverlap); px=px/fs*2; px(\[1,end\])=px(\[1,end\])/2;其中px为功率谱密度,fx为对应的频率。\[3\] 因此,如果你想在MATLAB中计算功率谱密度,可以使用pwelch函数或者结合spectrum函数使用PSD函数。 #### 引用[.reference_title] - *1* [MATLAB功率谱](https://blog.csdn.net/CV_Jason/article/details/84996400)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* *3* [MATLAB功率谱函数psdpwelch使用](https://blog.csdn.net/qq_41687938/article/details/117869445)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

子木呀

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值