gusucode.com > m4k20功率分配图 mimo功率优化,OFDM典型仿真程序,贪婪算法仿真程序 > code/15883844m4k20power/m4k20功率分配图xin/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('信道情况较好时功率分配图');