Matlab实现CoMP多用户注水算法在最最基础的注水算法的基础上

clc clear all; nt=2; nr=3; SNR_V_db = [10:1:20]; SNR_V_db_low=[0:1:10]; SNR_V = 10.^(SNR_V_db/10); SNR_V_low = 10.^(SNR_V_db_low/10); N0 = 1e-4; B = 1; for(i = 1 : length(SNR_V)) for(j = 1 : length(SNR_V_low)) Pt=N0; for(k = 1 : 1e3) H = random('rayleigh',1,nr,nt); [S V D] = svd(H); segama(:,k) = diag(V); H = segama(:,k).^2/(B*N0); landas(1,k)=SNR_V(i)*H(1); landas(2,k)=SNR_V_low(j)*H(1); landas(3,k)=SNR_V_low(j)*H(2); landas(4,k)=SNR_V(i)*H(2); [Capacity(k) PowerAllo] = WaterFilling_Self(Pt,landas(:,j),B,N0); end MeanCapacity(i,j)=mean(Capacity); end end [x,y]=meshgrid(SNR_V_db,SNR_V_db_low); mesh(x,y,MeanCapacity); set(h,'EdgeColor','r','FaceColor','r','MarkerEdgecolor','r','MarkerFacecolor','r');