Pwelch只返回单个值吗?

4视图(30天)
Shayma阿尔阿里
Shayma阿尔阿里 2022年2月25日
评论道: 马修•诺伊 2022年3月8日
我试图使用pwelch方法以计算两个信号之间的相干幅度平方的置信区间。然而,pwelch函数只返回单个值。我做错了什么?
我也附上我的代码中使用的变量。
我的代码:
%计算pwelch
fs = 1 / i_ph1.spacing_l;
[pwelch_1, f1, pxx1] = pwelch (i_ph1.l,汉宁(长度(i_ph1.l)),(长度(i_ph1.l)) / 2, fs,“ConfidenceLevel”,0.95);
[pwelch_2, f2, pxx2] = pwelch (i_ph2.l,汉宁(长度(i_ph2.l)),(长度(i_ph2.l)) / 2, fs,“ConfidenceLevel”,0.95);
[pwelch_3, f3, pxx3] = pwelch (i_ph3.l,汉宁(长度(i_ph3.l)),(长度(i_ph3.l)) / 2, fs,“ConfidenceLevel”,0.95);
[pwelch_12、f12 pxx12] = pwelch (xcorr (i_ph1.l i_ph2.l),汉宁(长度(i_ph1.l)),(长度(i_ph1.l)) / 2, fs,“ConfidenceLevel”,0.95);
[pwelch_23, f23 pxx23] = pwelch (xcorr (i_ph2.l i_ph3.l),汉宁(长度(i_ph2.l)),(长度(i_ph2.l)) / 2, fs,“ConfidenceLevel”,0.95);
[pwelch_31、f13 pxx13] = pwelch (xcorr (i_ph3.l i_ph1.l),汉宁(长度(i_ph3.l)),(长度(i_ph3.l)) / 2, fs,“ConfidenceLevel”,0.95);
%计算CMS
% Pxy ^ 2 / Pxx * Pyy组
CMS12 = (pwelch_12)。^ 2. / (pwelch_1。* pwelch_2);
CMS23 = (pwelch_23)。^ 2. / (pwelch_2。* pwelch_3);
CMS13 = (pwelch_31)。^ 2. / (pwelch_3。* pwelch_1);
%计算CMS的conf间隔
CMS_c13 = (pxx12)。^ 2. / (pxx1。* pxx2);
CMS_c23 = (pxx23)。^ 2. / (pxx2。* pxx3);
CMS_c13 = (pxx13)。^ 2. / (pxx1。* pxx3);
1评论
马修•诺伊
马修•诺伊 2022年3月8日
你好
你有试过 mscohere吗?
帮助mscohere
mscohere平方级相干估计。
Cxy = mscohere (X, Y)估计,平方级相干估计
系统的输出输入X和Y使用韦尔奇的平均,
改进的周期图法。连贯性是频率的函数
值在0和1之间,说明如何输入X对应
在每个频率输出Y。大小的平方连贯性,
Cxy由Cxy = (abs (Pxy) ^ 2)。/ (Pxx。* Pyy组)Pxx和肽yy
功率谱密度(PSD)估计的X和Y,分别;和
Pxy是Cross-PSD(运行CPSD)估计的X和y看到“帮助pwelch”
“帮助运行cpsd”的完整细节。

登录置评。

答案(0)

社区寻宝

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

开始狩猎!