>> 0.5*log2(det(H*Rxx*H'+1))

ans =

    6.3219

>> H

H =

    80    50

>> Rxx

Rxx =

    0.9996         0
         0    0.0004

>> 0.5*log2(det(H*Rxx*H'+1))

ans =

    6.3219

>> Rbinv=H'*H+inv(Rxx)

Rbinv =

  1.0e+003 *

    6.4010    4.0000
    4.0000    5.0020

>> Gbar=chol(Rbinv)

Gbar =

   80.0063   49.9961
         0   50.0239

>> G=inv(diag(diag(Gbar)))*Gbar

G =

    1.0000    0.6249
         0    1.0000

>> S0=diag(diag(Gbar))*diag(diag(Gbar))

S0 =

  1.0e+003 *

    6.4010         0
         0    2.5024

>> b=0.5*log2(det(Rxx*S0))

b =

    6.3219

>> SNR=Rxx*S0

SNR =

  1.0e+003 *

    6.3984         0
         0    0.0010

>> Wunb=SNR*inv(SNR-eye(2))*inv(S0)*inv(G')*H'

Wunb =

    0.0125
    0.0200

>> Gunb=eye(2)+SNR*inv(SNR-eye(2))*(G-eye(2))

Gunb =

    1.0000    0.6250
         0    1.0000