bartlett.m

传输参考TH－PPM调制超宽带信号的功率谱密度计算

%figure;
fc = 50e9;
Stx = transmitter_TR_TH_2PPM;
N = length(Stx);
L = 100;%%%%%%%%L越大，PSD方差越小
M = N/L;
P = zeros(1,M);
%计算每段数据功率谱
for i=1:L
    si = Stx((i-1)*M+1:i*M);%第i段数据功率谱
    Si = fft(si,M);%%%FFT
    Pi = (abs(Si)).^2/M;
    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('功率谱 /dBW'),title('TR-PPM PSD');
%grid;