Galileo_ideal_sig_gen_导航信号_GALILEOE5_IDEAL_信号仿真_Galileo_

function [Pxx,fx,SigPow]=PWelch(SignalCalibration,Fs,Nfft,overlap,BandWid) w = 0.5*(1 - cos(2*pi*(1:Nfft)/(Nfft+1))); n = length(SignalCalibration); nwind = length(w); k = fix((n-overlap)/(nwind-overlap)); index = 1:nwind; KMU = k*norm(w)^2; % Normalizing scale factor ==> asymptotically unbiased % Calculate PSD Spec = zeros(1,Nfft); for i=1:k xw = w.*SignalCalibration(index); index = index + (nwind - overlap); Xx = abs(fft(xw,Nfft)).^2; Spec = Spec + Xx; end Pxx = (Spec / (KMU *Fs)); f = 0:(2*pi/Nfft):(2*pi-2*pi/Nfft); f = f*Fs/(2*pi); % freqvector wanted in Hz Pxx = fftshift(Pxx); MidIndex = Nfft/2+1; fx = (f-Fs/2)/1e6; IndOffset = round(BandWid/(Fs/Nfft)*0.5); SigPow = 10*log10(sum(Pxx((MidIndex-IndOffset:MidIndex+IndOffset)))*Fs/Nfft); Pxx = 10*log10((Pxx*Fs/Nfft*10));