error_main.m

ofdm系统功率自适应功率分配算法研究的matlab程序

clc
clear  
K=20;
 global Mm;
 Mm=4;
 M=Mm;
for i=1:1:3
Pmax=20;
a=2*(sqrt(M)-1)/(sqrt(M)*log2(sqrt(M)));
b=3/(M-1);
x1=zeros(1,20);
x2=zeros(1,20);
x3=zeros(1,20);
y1=zeros(1,20);
y2=zeros(1,20);
y3=zeros(1,20);
for k=1:20
aa=2*k*linspace(0.1,1,20);
pp=optimal_pow_alloc(aa,K);
p=suboptimal_alloc(aa,K);

%plot(aa,pp,'r-x',aa,p,'g-*');%看功率分配情况时的图
%ylabel('功率分配');
%xlabel('信道情况');
%title('红x为最佳算法；绿*为次佳算法');
grid
%pause
r_opt=aa.*pp;%信噪比
r_sub=aa.*p;
x1(k)=x1(k)+sum(r_opt)/20;%平均信噪比
x2(k)=x2(k)+sum(r_sub)/20;
x3(k)=x3(k)+sum(aa)/20;
t1(k)=10*log10(x1(k));
t2(k)=10*log10(x2(k));
t3(k)=10*log10(x3(k));
e_opt=a*Q(sqrt(b.*aa.*pp));%误比特率矩阵
e_sub=a*Q(sqrt(b.*aa.*p));
e_equ=a*Q(sqrt(b.*aa));
y1(k)=y1(k)+sum(e_opt)/20;
y2(k)=y2(k)+sum(e_sub)/20;
y3(k)=y3(k)+sum(e_equ)/20;
end
semilogy(t1,y1,t2,y2,'r-*',t3,y3,'g-*')
title('曲线为最佳算法；红x为次最佳算法；绿*为等功率分配算法');
text(5,0.01,'\fontsize{14}\itM=4')
text(8,0.03,'\fontsize{14}\itM=16')
text(12,0.25,'\fontsize{14}\itM=64')
ylabel('功率EBR');
xlabel('平均SNR（db）');
grid on
hold on
M=M*4
Mm=M;
end
%end

%title('信道情况较好时功率分配图');