gusucode.com > m4k20功率分配图 mimo功率优化,OFDM典型仿真程序，贪婪算法仿真程序 > code/m4k20power/m4k20功率分配图/error_main.m

    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
    nn=2*k*linspace(0.1,1,20);

aa=2*k*linspace(0.1,1,20);
tt=2*k*linspace(0.1,1,20);
pp=optimal_pow_alloc(aa,K);
p=suboptimal_alloc(aa,K);
c=equ(tt,K)
n=dengxin(nn,K);

plot(nn,n,'r^',tt,c,'rd',aa,pp,'rp',aa,p,'rh');%看功率分配情况时的图
ylabel('功率分配');
xlabel('信道情况');
title('五角最佳；六角次佳；三角等信躁；菱形等功率');
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('信道情况较好时功率分配图');