SNR=5;                             %定義信噪比為0HZ
A=10^(SNR/10);                     %信噪比的單位轉(zhuǎn)換關(guān)系式
Im=eye(4);                         %這里的"1"為發(fā)射天線和接收天線兩者數(shù)目少的根數(shù)
sum=0;                             %初始化sum為0
for i=1:10000                      %程序循環(huán)10000次
    H=randn(4,30);                  %產(chǎn)生高斯信號(hào)源
      Q=H*H';                      %求威沙特矩陣
 sum=sum+log2(det(Im+A*Q/30));      %10000次循環(huán)得到的信道容量總和
end
C=sum/10000                         %10000次循環(huán)的平均值