welch.m

传输参考TH-PPM调制超宽带信号功率谱密度

%采用Welch法计算正弦信号功率谱,Hamming Window;

figure;
fc = 50e9;
Stx = TRANSMITTER_TR_TH_2PAM;
N = length(Stx);
Li = 100;%Bartlett法的分段数
M = N/Li;%每段数据长度
P = zeros(1,M);

%%%%%%%%%%%窗函数%%%%%%%%%%%%%%
n = 1:1:M;
w = zeros(1,M);
w(n) = 0.54-0.46*cos(2*pi*(n-1)/M);
square = w.^2;
totle = sum(square)/M;
%%%%%%%%%%%%%%%%%%%%%%%%%
L = (N-M/2)/(M/2);%真正分段数
P = zeros(1,M);
s1 = Stx(1:M);
s1_w = s1.*w;
S1_w = fft(s1_w,M);
P = (abs(S1_w)).^2/M/totle;
%计算每段数据功率谱
for i=1:L-1
    si = Stx(i*M/2+1:(i+2)*M/2);%第i段数据功率谱
    si_w = si.*w;
    Si_w = fft(si_w,M);%%%FFT
    Pi = (abs(Si_w)).^2/M/totle;
    P = P+Pi;
end
P = P/L;
%横坐标，赫兹；纵坐标，dB
t = 0:round(M/2-1);
k = t*fc/M;
P_dB = 10*log10(P(t+1));
plot(k,P_dB);
xlabel('频率 /Hz'),ylabel('功率谱密度 /dB'),title('TR TH-PAM PSD');