my_original_papr.asv

OFDM峰均比经典算法matlab程序

clc;
clear;
%参数设置
ifftn = 128;  %FFT点数
symbol_num = 1000; %OFDM符号数
data_length = ifftn*symbol_num;  %数据长度
j = 0;

%原始数据产生/BPSK调制
data_source = rand(1,data_length);
data_source(find(data_source<0.5)) = -1;
data_source(find(data_source >=0.5)) = 1;


MM=1:.1:10;
%OFDM 调制
data4ifft = reshape(data_source,ifftn,symbol_num);  %S/P
% data4ifft = data4ifft';
ccdf1=zeros(1,91);
for i=1:symbol_num;
%      z1=[data4ifft(i)(1:ifftn/2),zeros(1,3*ifftn),data4ifft(i)(ifftn/2+1:ifftn)]; %4 oversample
    w1(:,i)=ifft(data4ifft(:,i));            
    w1(:,i)=w1(:,i)*ifftn;
    x2=(abs(w1(:,i))).^2;%抽样的信号功率
    m1=mean(x2);
    v1=max(x2);
    papr(i)=10*log10(v1/m1);

    for j=1:91;
        if papr(i)>MM(j);
            ccdf1(j)=ccdf1(j)+1;
        end
    end
end
% [c,d]=butter(6,0.5);
% w3=filter(c,d,w2);
% MM=1:.1:10;
%  ccdf0=ccdf1./n1;
%  ccdf3=ccdf2./n1;
ccdf1 = ccdf1./symbol_num;
 semilogy(MM,ccdf1,'b-*');%,MM,ccdf3,'b-');
 title('clipping的CCDF曲线比较')
xlabel('papr(dB)'),ylabel('ccdf')
legend('clipping前','clipping后',2)
grid on;
hold on;