error_main.m

ofdm系统中功率自适应分配算法研究的matlab源代码（运行出误比特图形）

clc
clear  
K=20;
M=4;
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,'r-x',t2,y2,'g-*',t3,y3)

title('红x为最佳算法；绿*为次佳算法;蓝线为相等功率分配')
ylabel('平决误比特率');
xlabel('平决信道增益（db）');
grid