%***************************************************
% 基于編碼方法和SLM方法的OFDM信號PAR抑制算法的仿真
% 課程：移動通信中的關鍵技術
% 作者：劉熹
% 日期：2008/05/30
%***************************************************
clc; clear
%***************************************************
% 用時域法計算OFDM符號PARP, BPSK調制, 符號數為2^N
% N 載波數
% T 符號周期
% M 每符號采樣次數
% pt 時域信號的包絡功率
% papr 每符號內峰均功率比
%---------------------------------------------------
N = 4; T = 1; M = 4;
pt = zeros(1,(2^N)*M); papr = zeros(1,2^N);
E = (exp(i*2*pi*([0:1:N-1]'/T)*[0:1/M:1-1/M]))/sqrt(N);
for i = 0:(2^N-1) %對每個符號循環
    symbol = dec2bin(i,N);
    for j = 1:N
        temp = zeros(1,2^N);
        if (symbol(j)=='0') %對符號的每一位循環求和
            temp = temp - E(j,:);
        else
            temp = temp + E(j,:);
        end
    end
    temp = abs(temp).^2   %計算功率
    pt(i*M+1:(i+1)*M) = temp;
    papr(i+1) = 10*log10(max(temp)/mean(temp))
end
figure(1);plot(pt);grid;xlabel('符號');ylabel('包絡功率(W)');
plotyy(X1,Y1,X2,Y2,'function1','function2')