MATLAB: Is the unit of Cxy from cpsd in decibel

cpsdcxydecibelgainphasesignal processingSignal Processing Toolbox

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!

Best Answer

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);

Related Solutions

MATLAB: I want to represent the PSD of three signals on a waterfall plot. How can this be done

First of all, do not name a variable the same as a function! Change the the first psd variable to

psd0 = psd(Hs,xn,'Fs',fs,'NFFT',1024,'SpectrumType','onesided');

then make the call to waterfall like this:

waterfall([psd0.data psd1.data psd2.data]')

MATLAB: Evaluate the system impose response function and the FRF with given input and output data

HELLO

had to find another way than using the ID Toolbox (as I don't have it)

so I basically tried to fit an IIR model to the FRF and then you can create the impulse response

data = importdata ('beam_experiment.mat');
x = transpose(data.x); %input
y = transpose(data.y); %output
fs = data.fs;   % sampling frequency
T1 = 20;
NFFT = 2048;
NOVERLAP = 0.75*NFFT;
[Txy,F] = tfestimate(x,y,hanning(NFFT),NOVERLAP,NFFT,fs);
% IIR model
W = linspace(0,pi,length(F));
ind = find(F>5 & F <80); % frequency weighting ; data reliable between 5 and 80 Hz
WT = zeros(size(W));
WT(ind) = 1;
NB = 6;
NA = NB+2;
[B,A] = invfreqz(Txy,W,NB,NA,WT,1e4);
% check bode plots
[H,WW] = freqz(B,A,length(F));
figure(1),
subplot(2,1,1),plot(F,20*log10(abs(Txy)),'b',F,20*log10(abs(H)),'r');grid
subplot(2,1,2),plot(F,180/pi*(angle(Txy)),'b',F,180/pi*(angle(H)),'r');grid
% Impulse response 
[IR,X] = dimpulse(B,A);
samples = length(IR);
time = 1/fs*(1:samples);
figure(2),
plot(time,IR,'r');grid
title('Impulse response')
xlabel('Time (s)');
ylabel('Amplitude')

Best Answer

Related Solutions

MATLAB: I want to represent the PSD of three signals on a waterfall plot. How can this be done

MATLAB: Evaluate the system impose response function and the FRF with given input and output data

Related Question