How to do PSD measuremnet for the given signal
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Hi, I am doing PSD (ms2) analysis for the frequency spectrum attached in the figure. I am able to find PSD but the results are not in PSD(ms2).
Here is the code:
figure
[pxx,f] = pwelch(rwave_decimate_1,[],[],[],fs/N);
%the mean
mean = 0;
for i=1:length(pxx)
mean = mean + pxx(i);
end
M = mean/length(pxx);
%the Varaince
var = 0;
for j=1:length(pxx)
var = var + (pxx(j)-M)^2;
end
V = var/length(pxx);
plot(f,pxx)
It plots the PSD as in the attached figure. however, I want PSD in (ms2 or s2) on yaxis as in the figure 2, to make proper analysis on LF/HF region. Also, how can I highlight the LF and HF region with some shadded areas as in the figure2 ?
:- ms2 = mean-variance square
it would be appreciated if someone give solution.
Thanks in the advacnce.
2 comentarios
dpb
el 5 de Jul. de 2021
M=mean(pxx);
V=var(pxx);
Don't use mean as a variable that aliases the builtin mean function in MATLAB.
Respuestas (1)
Arthi Sathyamurthi
el 19 de Jul. de 2021
Hello Sarfaraz,
The units of the PSD estimate from the pwelch function are in squared magnitude units of the time series data per unit frequency. Thus, if the input data, rwave_decimate is in ms and the sample rate given as fs/N is in Hz, the PSD estimate which in your case pxx, ms^2/Hz.
3 comentarios
Arthi Sathyamurthi
el 20 de Jul. de 2021
If the input data, rwave_decimate is in ms and the sample rate given as fs/N is in Hz, the PSD estimate pxx will be in ms^2/Hz.
Ver también
Categorías
Más información sobre Multirate Signal Processing en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!