ifft_papr.m

课程作业

%***************************************************
% 用IFFT法计算OFDM符号PARP, 采用BPSK调制 
% N 载波数  符号数为2^N
% M 过采样因子  可选2^n,IFFT后每符号实际输出N*M个样点
% pt 时域信号的包络功率
% papr 每符号内峰均功率比
%---------------------------------------------------
clc; clear
N = 4; M = 4;
pt = zeros(1,(2^N)*N*M); papr = zeros(1,2^N);
for i = 0:(2^N-1) %对每个符号循环
    symbol = dec2bin(i,N);
    temp = zeros(1,N*M);
    for j = 1:N   %对符号的每一位循环
        if (symbol(N-j+1)=='0') %字符串尾为最低位 
            temp(j) = -1;
        else
            temp(j) = 1;
        end
    end
    temp = abs(ifft(temp)).^2
    pt(i*N*M+1:(i+1)*N*M) = temp;                %计算功率
    papr(i+1) = 10*log10(max(temp)/mean(temp));  %计算PAPR
end
figure(1);
subplot(2,1,1);plot(pt);xlim([0 2^N*N*M]);set(gca,'Xtick',[0:N*M:2^N*N*M]);
grid;xlabel('时域采样');ylabel('包络功率(W)');
subplot(2,1,2);plot([0:1:2^N-1],papr);set(gca,'Xtick',[0:1:2^N]);xlim([-0.5  2^N-0.5]);
grid;xlabel('码字');ylabel('峰平功率比(dB)');