Ergodic Capacity Rayelig Fading Channel

[RajaKashif]

Ergodic Capacity Rayeligh Fading Channel

Hi,

Is the following correct to obtain the ergodic capacity in a Rayleigh Fading Channel.

******************************************************************************

SNRdB=2;
SNR = 10 ^ (SNRdB/10);

nsample_h=2000;
N0=1; % Noise power density

cap_inst=zeros(1,nsample_h);

for ChannelLoop = 1:nsample_h
h_11 = sqrt( N0 /2)* ( randn + sqrt(-1) * randn ) ;
cap_inst(ChannelLoop)=log2 (1 + SNR * abs(h_11)^2 );
clear h_11;
end

cap_ergodic=mean ( cap_inst );

 

[abhijithbg]

Check out the following code. It give Ergodic capacity for Rayleigh channel
clear all;close all;clc;
global SNRVal;
BW = 1;
SNRdB = -15:5:20;
x0 = 1; % Make a starting guess at the solution
GammaArr = zeros(1,length(SNRdB));
for ii = 1:length(SNRdB)
SNRVal = 10^(SNRdB(ii)/10);
[x,fval] = fsolve(@myfun,x0);
GammaArr(ii)= x;
end
Crayleigh = BW*(log(1+GammaArr));

%%--------------------Calling function:myfun.m------
function F = myfun(x)
global SNRVal;
% F = [2*x(1) - x(2) - exp(-x(1)); -x(1) + 2*x(2) - exp(-x(2))];
F = x-log(1+x)-SNRVal;
end