cpsd output vector differs from plotted result

3 次查看（过去 30 天）
DLR
DLR 2024-8-1
评论： Voss 2024-9-5
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')
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!

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

采纳的回答

Voss
Voss 2024-8-1
编辑：Voss 2024-8-1
Ran in:
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);
Warning: Rounding order to nearest integer.
% 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);
Warning: Rounding order to nearest integer.
legend('CPSD output','CPSD direct plot')
% zoom in to see that the lines match
xlim([0 0.01])

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息

标签

产品


版本

R2024a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by