Is the unit of Cxy from cpsd in decibel?

Question

Kay 2011-10-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/18915-is-the-unit-of-cxy-from-cpsd-in-decibel

Hi,

I want to use the Signal Processing Toolbox to compute the Cross spectrum gain and phase for my sequences X and Y. I wondered the Cxy that cpsd gives me is in the units of dB? I saw the example and it says the unit is dB/rad/freq, how can I convert it to dB?

Thank you so much!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Wayne King 2011-10-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/18915-is-the-unit-of-cxy-from-cpsd-in-decibel#answer_25287

在 MATLAB Online 中打开

Hi, no, the output is the cross spectrum, it is complex-valued. If you want the output in db. Use

10*log10(abs(Pxy));

For example:

   Fs = 1000;  
   t = 0:1/Fs:.296;
   x = cos(2*pi*t*200)+randn(size(t)); 
   y = cos(2*pi*t*100)+randn(size(t));  
   [Pxy,F] = cpsd(x,y,[],[],[],Fs);  
   plot(F,10*log10(abs(Pxy))); grid on;
   % compare to 
   figure;
   cpsd(x,y,[],[],[],Fs);

4 个评论
显示 2更早的评论隐藏 2更早的评论

Wayne King 2011-10-20

yes, angle(Pxy) gives you the phase. The difference between 10*log10() and 20*log10() is the difference between magnitude and magnitude squared, they are both in dB per unit frequency.

Kay 2011-10-20

is dB same as dB/Hz?

if not, how do I get rid it?

请先登录，再进行评论。

Answer 2

Kay 2011-12-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/18915-is-the-unit-of-cxy-from-cpsd-in-decibel#answer_30020

I'm getting some offset of the result after cpsd, do you know if the algorithm of the cpsd would give me a larger set of numbers as a result?

thanks, Kay

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Wayne King 2011-12-1

offset after what result?

请先登录，再进行评论。

Answer 3

Giandomenico 2013-1-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/18915-is-the-unit-of-cxy-from-cpsd-in-decibel#answer_73019

编辑：Giandomenico 2013-1-31

在 MATLAB Online 中打开

I have a problem with cpsd amplitude. I generate a 20s signal with two harmonic components at 5 and 10Hz, both with an amplitude of 10, then I perform the cpsd I have the two components at the right frequencies, but both with an amplitude of 1000??? it should be 100, right?

fs=512;         % sampling frequency
L=10240;        % samples
t=(0:L-1)'/fs;
x=10*sin(2*pi*10*t)+10*sin(2*pi*5*t);
[Gxx, f]=cpsd(x,x,rectwin(L),0,L,fs);
plot(f,abs(Gxx),'.r')

Where I am wrong?