function p=lnorrandWN(M,N,mdB,sdB)
%function p=lnorrandWN(M,N,mdB,sdB)
C=10*log10(exp(1)); 
s=sdB/C;
m=10^(mdB/10);
p=m*exp(randn(M,N)*s);
%mean value exp(log(m)+s^2/2); %variance exp(2*log(m)+s^2)*(exp(s^2)-1);