功率谱估计MATLAB实现

BT methodРRxWin = Rx(1:M).*window_v; % smooth window multiply auto-correlation function РPw = abs(fft(RxWin, N));РРk = 0:floor(N/2 - 1);Рfigure;Рplot(k*Fs/N, 10*log10(Pw(k+1)/max(Pw)),'r');Рtitle('BT Method PSD Estimate');Рxlabel('Frequency(Hz)');ylabel('Normalized PSD(dB)'),grid on;Р%% Estimate PSD using Welch methodРwindow = 32; % the length of each segmentРnoverlap = 8; % overlap number for each segmentРnfft = pow2(nextpow2(length(x))); % nfft-points FFT for each segmentР[Pxx,f] = pwelch(x,window,noverlap,nfft,Fs); % call pwelch functionРfigure;Рplot(f,10*log10(Pxx/max(Pxx)),'g');Рtitle('Pwelch Method PSD Estimate');Рxlabel('Frequency(Hz)');ylabel('Normalized PSD(dB)'),grid on;РРР8