cpsd output vector differs from plotted result
30 次查看(过去 30 天)
显示 更早的评论
When using cpsd to calculate the cross power spectral density, I've noticed that the plot produced by directly running the command with no output arguments is slightly different from what I get when I run it with an output argument and plot that output myself.
The following code (attached along with example signal file) produces the figure below.
signalx = randn(2e4,1);
signaly = randn(2e4,1);
window = length(signalx)/3;
noverlap = round(0.9*window);
fs = 2074;
[Pxy,F] = cpsd(signalx,signaly,window,noverlap,[],fs);
plot(F./1000,10.*log10(real(Pxy)))
hold on
cpsd(signalx,signaly,window,noverlap,[],fs);
legend('CPSD output','CPSD direct plot')
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1744356/image.png)
I'm trying to run cpsd column-wise on a pair of very large matrices. Obviously I can't just plot them one by one, but I don't know how to interpret this difference and I'm not sure if I trust the output argument. Does anyone know what's going on here? Thanks!
0 个评论
回答(1 个)
Voss
2024-8-1,4:13
编辑:Voss
2024-8-1,4:29
Looks like cpsd with no outputs plots the magnitude, abs(Pxy), rather than the real part, real(Pxy).
signalx = randn(2e4,1);
signaly = randn(2e4,1);
window = length(signalx)/3;
noverlap = round(0.9*window);
fs = 2074;
[Pxy,F] = cpsd(signalx,signaly,window,noverlap,[],fs);
% use a thick line so it can be seen underneath the thin line on top
plot(F./1000,10.*log10(abs(Pxy)),'--','LineWidth',2)
hold on
cpsd(signalx,signaly,window,noverlap,[],fs);
legend('CPSD output','CPSD direct plot')
% zoom in to see that the lines match
xlim([0 0.01])
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!