original_papr.m

OFDM Source is very easy. so you can easily understand.

function [PAPR] = original_PAPR(M, N, factor);

D = N * 2;                       % 涝仿 单捞磐 肮荐
Fs=4;                       % Oversampling 
Fd=1;                       % Data frequency 

%=================================================%
%-data 惯积-%
[sig_in]= gen_qam_qpsk(M, D);
%=================================================%
%================================================================
%=========================== OFDM ==============================%
sp=vec2mat(conj(sig_in), N)';   %S/Pa
os=[sp(1: N /2,:); zeros(N * Fs / Fd - N,size(sp,2)); sp(N /2+1: N,:);];
ifft_out=ifft(os);                    %IFFT 函炼
ps=vec2mat(ifft_out',1)';       %P/S
y=ps*Fs/Fd * N;                    
%============================= CCDF ==============================%
PAPR=10*log10(max(abs(y).^2)/(norm(y)^2/length(y)));