【正文】 f4=zeros(1,length(Capacity))。 if q==1 f41(q,1)=0。 elseif q~=1 f41(q,1)=count41(1)/N。 end for k=2:(length(Capacity)) f4(k)=count41(k)/N。 f41(q,k)=f41(q,k1)+f4(k)。 end end %圖形繪制 f1=f41(1,:)。 f2=f41(2,:)。 f3=f41(3,:)。 f4=f41(4,:)。 f5=f41(5,:)。 %這里 q的最大值是 5，如果取更大自己添加 11 plot(Capacity,f1,39。39。,Capacity,f2,39?！? Capacity,f3,39。o39。,Capacity,f4,39。*39。,Capacity,f5,39。+39。) axis([0 30 0 1]) xlabel(39。C bit/s/Hz39。) ylabel(39。MIMO容量小于橫坐標概率分布 39。) title(39。發(fā)射信噪比 SNR=10dB39。) legend(39。1x139。,39。3x339。,39。5x539。,39。7x739。,39。9x939。) grid on 2 固定天線數(shù)目時，發(fā)射信噪比為變量， MIMO 容量累積分布 %MIMO_indB,CSI未知情況下平均分配功率 %慢瑞利衰落信道。 t=input(39。置信區(qū)間寬度 39。)。 SNR_indB_max=input(39。信噪比 maxdB39。)。 SNR_indB_min=input(39。信噪比 mindB39。)。 z=input(39。輸入信噪比步長 39。)。 SNR_indB=SNR_indB_min:z:SNR_indB_max。 Capacity=0:t:40。 N=input(39。迭帶次數(shù) 39。)。 Nt=input(39。輸入 發(fā)射天線數(shù)目 39。)。 Nr=input(39。輸入接收天線數(shù)目 39。)。 q1=waitbar(0,39。MIMO39。)。 for x=1:length(SNR_indB) waitbar(x/length(SNR_indB))。 SNR(x)=10^(SNR_indB(x)/10)。 for i=1:N H=raylrnd(,Nr,Nt)。 a=abs(H)。 for m=1:Nr for n=1:Nt b(m,n)=a(m,n)^2。 end end c=sum(b)。 u=0。 for m=1:Nt u=u+log2(1+SNR(x)*c(m)/Nt)。 end C4(x,i)=u。 end count41=panduan(C4(x,:),N,length(Capacity)1,t)。 f4=zeros(1,length(Capacity)1)。 f41(x,1)=count41(1)/N。 for k=2:(length(Capacity)1) f4(k)=count41(k)/N。 f41(x,k)=f41(x,k1)+f4(k)。 12 end end close(q1)。 for k=1:(length(Capacity)1) w1(k)=(k1)*t。 end f1=f41(1,:)。 f2=f41(2,:)。 f3=f41(3,:)。 f4=f41(4,:)。 f5=f41(5,:)。 f6=f41(6,:)。 %信噪比更大時自己添加 plot(w1,f1,39。39。,w1,f2,39。.39。,w1,f3,39。39。,w1,f4,39。*39。,w1,f5,39。o39。,w1,f6,39。+39。) axis([0 40 0 1]) xlabel(39。C bit/s/Hz39。) ylabel(39。MIMO總?cè)萘扛怕史植?CDF39。) title(39。天線數(shù)目為 :4x439。) legend(39。0dB39。,39。5dB39。,39。10dB39。,39。15dB39。,39。20dB39。,39。25dB39。) grid on 3 MIMO 容量與天線的關系 （ 1）、 %MIMO_天線 ,CSI未知情況下平均分配功率 %平坦瑞利衰落環(huán)境 %中斷概率為 %，即容量等級概率為 99% SNR_indB=0:5:25。 %發(fā)射信噪比 N=input(39。迭帶次數(shù) 39。)。 Na=2:2:20。 %發(fā)射天線等于接收天線 q1=waitbar(0,39。MIMO天線 39。)。 for x=1:length(SNR_indB) SNR=10^(SNR_indB(x)/10)。 for y=1:length(Na)。 for i=1:N H=raylrnd(,Na(y),Na(y))。 a=abs(H)。 b=a.*a。 c=sum(b)。 u=0。 snr=SNR*c/Na(y)。 for m=1:Na(y) u=u+log2(1+snr(m))。 end C4(y,i)=u。 end f41(x,y)=*max( C4(y,:))。 %取容量等級 end waitbar(x/length(SNR_indB))。 13 end close(q1)。 f1=f41(1,:)。 f2=f41(2,:)。 f3=f41(3,:)。 f4=f41(4,:)。 f5=f41(5,:)。 f6=f41(6,:)。 plot(Na,f1,39。39。,Na,f2,39。.39。,Na,f3,39。39。