MIMO注水功率分配MATLAB程序

%收发为圆形阵列的MIMO信道容量分析，注水功分
%Transfer function\
clc;
clear all;
close all;
belta=1;
phi=0;% 取值小于2*pi/Nrx;
BW=1;%1hz
lambda=0.125;% f=2.4Ghz
rx=1*lambda;
tx=1*lambda;  %半径的影响比较大
Nrx=12;Ntx=12;
H_d=zeros（NrxNtx）;
D_rei=lambda;%2*（rx^2）/lambda;%瑞利距离
cnt_row=1;
res=5;%阵列距离为1个瑞利距离
%%
Sdb=0;Edb=30;cntPoint=300;  %满功率信噪比范围
Capt=zeros（1cntPoint）;
D=res*D_rei;% distance between the transmitter-reciever array
for p=1:Ntx
    for q=1:Nrx
        theta=2*pi*（q/Nrx-p/Ntx）+phi;
        d_theta=sqrt（D^2+rx^2+tx^2-2*rx*tx*（cos（theta））^2）;
        H_d（qp）=belta*lambda/（4*pi*d_theta）*exp（-1j*2*pi*d_theta/lambda）;
    end
end
S=svd（H_d）;%return the signular value vector
for nuo=10.^（linspace（SdbEdbcnt