bartlett.m

来自「传输参考TH-PPM调制超宽带信号功率谱密度」· M 代码 · 共 23 行

figure;
fc = 50e9;
Stx = TRANSMITTER_TR_TH_2PAM;
N = length(Stx);
L = 100;
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-PAM PSD');
%grid;