clc;
clear?all;
for?mr=1:4
????if（mr==1）?
???????mt=1;
????else
???????mt=2;
????end
????ESNO=-5:5:30;
????a=length（ESNO）;
????cequ=zeros（1a）;?
????for?x=1:length（ESNO）
???????esno=ESNO（x）;
???????esno=10^（esno/10）;
???????sum=0;
???????for?m=1:1000
??????????hr=randn（mrmt）/sqrt（2）;
??????????hi=randn（mrmt）/sqrt（2）;
??????????h=hr+j*hi;
??????????imr=eye（mr）;
??????????c=log2（det（imr+esno*h*（h‘）/mt））;
??????????sum=sum+c;
???????end
???????cequ（x）=sum/1000;