?? zperiodogram.m
字號:
%Periodogram函數:功率譜估計
%用FFT直接法計算功率譜
Fs=2000;
NFFT=1024;
n=0:1/Fs:1;
x=sin(2*pi*100*n)+4*sin(2*pi*500*n)+randn(size(n));
figure(1)
window=boxcar(length(x));
periodogram(x,window,NFFT,Fs);
%對噪聲序列進行計算
X=fft(x,NFFT);
Pxx=abs(X).^2/length(n); %求解PSD
t=0:round(NFFT/2-1);
k=t*Fs/NFFT;
P=10*log10(Pxx(t+1)); %縱坐標的單位為dB
figure(2)
plot(k,P);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -