2014-03-29 11 views
-1

여기에 W가 내 데이터 인 microvolt^2/Hz의 내 PSD 프로그램이 있습니다.전력 스펙트럼은 마이크로 볼트^2/Hz 단위로 나타 납니까?

Fs=128 
x1 =W; 
L = length(x1); 
NFFT = 2^nextpow2(L); 
f = Fs/2*linspace(0,1,NFFT/2+1); 
n = 0:L-1; 
pxx1 = pwelch(x1,L); 
plot(pxx); 

마이크로 볼트^2/Hz에서 PSD를 제작하려면 내 프로그램이 맞는지 틀린 지 알아야합니다.

+1

테스트해야합니다. 알려진 데이터에 대해 올바른 값을 제공합니까? –

답변

0

실제로 pwelch을 사용하면 Hz 당^2 단위의 전력 스펙트럼 밀도를 얻는 것이 옳습니다. 그러나 pwelch에 몇 가지 인수를 더 제공 할 수 있습니다.

% input is a time series x 

fs=128     % sample rate (Hz) 
bw = 0.1;     % desired spectrum resolution (Hz) 
nfft = fs/bw;    % length of fft (samples) 
nfft = 2^nextpow2(nfft); % round to nearest power of 2 
bw = fs/nfft;   % actual resolution 

% Estimate the power spectral density 
[Pxx, f] = pwelch(x, hanning(nfft), nfft/2, nfft, fs); 

% Plot the results 
loglog(f, Pxx) 
xlabel('frequency [Hz]'); 
ylabel('units^2/Hz');