?? ofdm_double_different_number.m
字號:
%雙層次相位序列優化算法（不同的分塊取值）
clear;
clc;
v=16;
nl=2;
for n=1:nl;
b1=ones(1,v);
x12=zeros(1,512);
x14=zeros(v,512);
x16=zeros(256,512);
a=rand(1,256)>0.5;
C=(a*2-1)/sqrt(2);
sreal=C(1:2:256);
simage=C(2:2:256);
D=sreal+j.*simage;
% oversample 4
z1=[D(1:128/2),zeros(1,3*128),D(128/2+1:128)];
z2=ifft(z1)*512;
l1=length(z1);
% divided into v sets
for k1=1:v;
z(k1,:)=[zeros(1,(k1-1)*128/v),D((k1-1)*128/v+1:k1*128/v),zeros(1,(l1-(k1-1)*128/v-128/v))];
iz(k1,:)=ifft(z(k1,:));
iz(k1,:)=iz(k1,:)*512;
x12(1,:)=x12(1,:)+b1(k1)*iz(k1,:);
end
x13=abs(x12).^2;
m10=mean(x13);
v10=max(x13);
papr5(n)=10*log10(v10/m10);
papr6(n)=papr5(n);
%次優化搜索
index=1;
while index<v+1
b1(index)=-1;
for m=1:v;
x14(index,:)=x14(index,:)+b1(m)*iz(m,:);
m+1;
end
x15=abs(x14(index,:)).^2;
m1=mean(x15);
v1=max(x15);
papr0000(n)=10*log10(v1/m1);
if papr0000(n)<papr5(n)
papr5(n)=papr0000(n);
else
b1(index)=1;
end
index=index+1;
end
%全搜索算法
jjj=1;
for i1=-1:2:1;
for i2=-1:2:1;
for i3=-1:2:1
for i4 = -1:2:1
for i5=-1:2:1
for i6=-1:2:1
for i7=-1:2:1
for i8=-1:2:1
for i9=-1:2:1;
for i10=-1:2:1;
for i11=-1:2:1
for i12= -1:2:1
for i13=-1:2:1
for i14=-1:2:1
for i15=-1:2:1
for i16=-1:2:1
A(jjj,1:16)=[i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,i11,i12,i13,i14,i15,i16];
jjj=jjj+1;
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
for i=2:256;
for k2=1:v;
x16(i,:)=x16(i,:)+A(i,k2)*iz(k2,:);
end
x17=abs(x16(i,:)).^2;
v5=max(x17);
m5=mean(x17);
papr000(n)=10*log10(v5/m5);
if papr000(n)<papr6(n);
papr6(n)=papr000(n);
end
end
% divided into v sets and v sets divided into b sets
t=1;
d=2; % 2,4,8
while d<9
z=zeros(d,512);
x2=zeros(1,128);
x4=zeros(v,128);
x6=zeros(1,l1);
x8=zeros(d,l1);
x10=zeros(1,l1);
w=zeros(1,128/d);
b=ones(1,v/d);
c=ones(1,d);
for k=1:d;
q=zeros(v/d,d);
p=zeros(v/d,d);
z(k,:)=[zeros(1,(k-1)*128/d),D((k-1)*128/d+1:k*128/d),zeros(1,(l1-(k-1)*128/d-128/d))];
w=z(k,(k-1)*128/d+1:k*128/d);%長度32
for m=1:v/d;
s(m,:)=[zeros(1,(m-1)*128/v),w((m-1)*128/v+1:m*128/v),zeros(1,(128-(m-1)*128/v-128/v))];
is(m,:)=ifft(s(m,:));
is(m,:)=is(m,:)*128;
x2(1,:)=x2(1,:)+b(m)*is(m,:);
end
x3=abs(x2).^2;
m0=mean(x3);
v0=max(x3);
papr0(n)=10*log10(v0/m0);
index=1;
while index<v/d+1
b(index)=-1;
for m=1:v/d;
x4(index,:)=x4(index,:)+b(m)*is(m,:);
m+1;
end
x5=abs(x4(index,:)).^2;
m1=mean(x5);
v1=max(x5);
papr1(n)=10*log10(v1/m1);
if papr1(n)<papr0(n)
papr0(n)=papr1(n);
q(index,k)=-1;
else
b(index)=1;
q(index,k)=1;
end
index=index+1;
end
iz(k,:)=ifft(z(k,:));
iz(k,:)=iz(k,:)*512;
x6(1,:)=x6(1,:)+c(k)*iz(k,:);
end
x7=abs(x6).^2;
m1=mean(x7);
v1=max(x7);
papr00(n)=10*log10(v1/m1);
%次優化搜索
index=1;
while index<d+1
c(index)=-1;
for m=1:d;
x8(index,:)=x8(index,:)+c(m)*iz(m,:);
m+1;
end
x9=abs(x8(index,:)).^2;
m1=mean(x9);
v1=max(x9);
papr2(n)=10*log10(v1/m1);
if papr2(n)<papr00(n)
papr00(n)=papr2(n);
c(1,index)=-1;
else
c(index)=1;
c(1,index)=-1;
end
index=index+1;
end
%矩陣相乘再計算峰均比
p(:,k)=q(:,k)*c(k);
P=reshape(p,1,v);
% oversample 4
for w=1:v;
z3(w,:)=[zeros(1,(w-1)*128/v),D((w-1)*128/v+1:w*128/v),zeros(1,(l1-(w-1)*128/v-128/v))];
iz(w,:)=ifft(z3(w,:));
iz(w,:)=iz(w,:)*512;
x10(1,:)=x10(1,:)+P(w)*iz(w,:);
end
x11=abs(x10).^2;
v3=max(x11);
m3=mean(x11);
papr3(n)=10*log10(v3/m3);
if t==1
papr7(n)=papr3(n);
elseif t==2
papr8(n)=papr3(n);
else
papr9(n)=papr3(n);
end
t=t+1;
d=2*d;
end
end
ccdf0=zeros(1,131);
ccdf2=zeros(1,131);
ccdf4=zeros(1,131);
ccdf6=zeros(1,131);
ccdf8=zeros(1,131);
NN=0:.1:13;
for n=1:nl
for l=1:131;
if papr5(n)>NN(l);
ccdf0(l)=ccdf0(l)+1;
end
if papr6(n)>NN(l);
ccdf2(l)=ccdf2(l)+1;
end
if papr7(n)>NN(l);
ccdf4(l)=ccdf4(l)+1;
end
if papr8(n)>NN(l);
ccdf6(l)=ccdf6(l)+1;
end
if papr9(n)>NN(l);
ccdf8(l)=ccdf8(l)+1;
end
end
end
ccdf1=ccdf0./nl;%迭代
ccdf3=ccdf2./nl;%全搜索
ccdf5=ccdf4./nl;%d=2
ccdf7=ccdf6./nl;%d=4
ccdf9=ccdf8./nl;%d=8
figure(1)
semilogy(NN,ccdf1,'gO-',NN,ccdf3,'bX-',NN,ccdf5,'cs-',NN,ccdf7,'r*-',NN,ccdf9,'kd-')
title('CCDF曲線比較')
xlabel('PAPR( in dB)'),ylabel('CCDF Pr(PAPR>PAPR0')
legend('PTS次優化算法','全搜索算法','雙層次相位序列優化v=16 d=2','v=16 d=4','v=16 d=8');
grid on?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -