Please check following matlab codes for following equation is correct

clc;
clear all;
close all;
K=5;                  
N=16;
aa = 1000;
ba = 5000;
r = (ba-aa).*rand(5,1) + aa;
qq=10/10;
q=10^qq;
s = 0 + q.*randn(K,1);
v = s/10;
si=10.^v;
u=1;
a=3.8;  
%r=real(r);
b=(u.*si)/(r.^a);
D=diag(b(1:K));
X=sqrt(D);
H = rand(N,K) + i*rand(N,K);
G=H * X;