求高手帮助 用 matlab 做 OFDM的限幅法降低PAPR值 的源码 得出CCDF曲线

我也在做这个 很是头疼 我找的一个程序 你看能不能帮弯悉歼助你吧 希望可以交流

clear all

close all

c1=128

n1=10000

Fs=100

MM=1:.1:13 % CCDF的门限值埋冲

ccdf1=zeros(1,121)

ccdf2=zeros(1,121)

for i=1:n1

x(:,1)=randsrc(c1,1,[+1 -1 +3 -3])

x(:,2)=randsrc(c1,1,[+1 -1 +3 -3])

y1=squeeze(x)

a=amodce(y1,Fs,'qam') % QAM modulated information

y2=a.'

z1=[y2(1:c1/2),zeros(1,3*c1),y2(c1/2+1:c1)] % 4 oversample

w1=ifft(z1)

w1=w1*c1

x2=(abs(w1)).^2

m1=mean(x2)

v1=max(x2)

papr(i)=10*log10(v1/m1)

y3=abs(w1)

CR=1.4

A=CR*sqrt(2)*std(w1) % CR=Amax/sqrt(Pin) Pin:the input power of the ofdm signal before clipping

w2=w1

for h=1:4*c1

if y3(h)>A

w2(h)=A*w1(h)/y3(h)

end

end

x3=(abs(w2)).^2

m2=mean(x3)

v2=max(x3)

papr0(i)=10*log10(v2/m2)

for j=1:121

if papr(i)>MM(j)

ccdf1(j)=ccdf1(j)+1

end

if papr0(i)>MM(j)

ccdf2(j)=ccdf2(j)+1

end

end

end

[c,d]=butter(6,0.5)

w3=filter(c,d,w2)

MM=1:.1:13

ccdf0=ccdf1./n1

ccdf3=ccdf2./n1

semilogy(MM,ccdf0,'b',MM,ccdf3,'r'陆亩)

xlabel('PAPR门限值/dB')ylabel('CCDF')

legend('clipping前','clipping后',2)

title('clipping的CCDF曲线比较')

grid on

% clipping的CCDF曲线比较 clipping_ccdf.m16QAM调制 子载波数为128

% papr of the ofdm signals

clear all

close all

c1=128

n1=10000

Fs=100

MM=1:.1:10 % CCDF的门限值

ccdf1=zeros(1,91)

ccdf2=zeros(1,91)

for i=1:n1

x(:,1)=randsrc(c1,1,[+1 -1 +3 -3])

x(:,2)=randsrc(c1,1,[+1 -1 +3 -3])

y1=squeeze(x)

a=amodce(y1,Fs,'qam'枯伍) % QAM modulated information

y2=a.'

z1=[y2(1:c1/2),zeros(1,3*c1),y2(c1/2+1:c1)] % 4 oversample

w1=ifft(z1)

w1=w1*c1

x2=(abs(w1)).^2

m1=mean(x2)

v1=max(x2)

papr(i)=10*log10(v1/m1)

y3=abs(w1)

CR=1.4

A=CR*sqrt(2)*std(w1) % CR=Amax/sqrt(Pin) Pin:the input power of the ofdm signal before clipping

w2=w1

for h=1:4*c1

if y3(h)>A

w2(h)=A*w1(h)/y3(h)

end

end

x3=(abs(w2)).^2

m2=mean(x3)

v2=max(x3)

papr0(i)=10*log10(v2/m2)

for j=1:91

if papr(i)>MM(j)

ccdf1(j)=ccdf1(j)+1

end

if papr0(i)>MM(j)

ccdf2(j)=ccdf2(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

semilogy(MM,ccdf0,'b'运缺,MM,ccdf3,'r')

xlabel('PAPR门限值/dB')ylabel('CCDF')

legend('clipping前','clipping后',2)

title('clipping的旁败辩CCDF曲线比较')

grid on

---