说明：  function C=capacity(SNR,Nt,Nr) %SNR=0; %定义信噪比为0HZ A=10^(SNR/10); %信噪比的单位转换关系式 Im=eye(min(Nr,Nt)); %这里的"1"为发射天线和接收天线两者数目少的根数 sum=0; %初始化sum为0 for i=1:10000 %程序循环10000次 H=randn(Nr,Nt); %产生高斯信号源 if Nr<Nt Q=H*H'; %求威沙特矩阵 else Q=H'*H; end sum=sum+log2(det(Im+A*Q/Nt)); %10000次循环得到的信道容量总和 end C=sum/10000 %10000次循环的平均值
(function C=capacity(SNR,Nt,Nr) %SNR = 0;% defines SNR as 0Hz A = 10 ^ (SNR / 10);% s / N unit conversion relation Im = eye (min (NR, NT));% here "1" is the number of the less number of transmitting antenna and receiving antenna Sum = 0;% initialize sum to 0 For I = 1:10000% program cycles 10000 times H = randn (NR, NT);% produces Gaussian signal source if Nr<Nt Q = h * H ';% to find the weissard matrix else Q=H'*H; End Sum = sum + log2 (DET (IM + A * q / NT));% sum of channel capacity from 10000 cycles End C = sum / 10000% average of 10000 cycles)