I want to plot PSD for my input data(meters)
2 views (last 30 days)
Show older comments
I want to plot PSD for my input data (meters) if the input in volts the the PSD in dB/Hz what if the input in meters then what are the units of PSD please some body help me Ts=30; Fs =1/Ts; y=MP1; %input data(meters) h = spectrum.welch; Hpsd=psd(h,y,'Fs',Fs); what are units of y axis after the plot
0 Comments
Accepted Answer
Dr. Seis
on 29 May 2012
If you are measuring displacement in meters in the time domain, then your frequency domain amplitudes (after performing the FFT) should be in units of meters / Hz. Squaring those amplitudes should give you meters-squared / Hertz-squared. However... since Matlab generates their power spectra by normalizing their FFT result by the number of samples (N) instead of the sampling frequency (Fs), their output might have PSD units associated with meters-squared / samples-squared (which is completely meaningless to me).
2 Comments
Dr. Seis
on 29 May 2012
It will take all of them... or you can re-parameterize your "h = spectrum.welch" to operate on a sub-set (or window) of your data with pre-defined overlap amounts. If you want the PSD of the full data, then you might have to do Hpsd = psd(h,y,'NFFT',650,'Fs',Fs) since in my test (running without 'NFFT',650) it performs a 256-point FFT. Check out the documentation for more specific examples of setting up/paramertizing these functions.
More Answers (1)
Wayne King
on 29 May 2012
Is your sampling interval in meters, or are you really saying that your measured variable is in meters?
I'm assuming you just mean that your data is sampled spatially. What is actually being measured?
See Also
Categories
Find more on Parametric Spectral Estimation in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!