Complex normal Gaussian noise
You should use the the standard formula:

s = randn(m, n) + 1i*randn(m, n);


And as pointed out by MBaz, the output should be scaled accordingly by $\frac{1}{\sqrt{2}}$

s = s/sqrt(2);


More on that topic you can find for example in a book Digital Media Processing by Hazarathaiah Malepati:



It should generate exactly what you want.