将复杂的AWGN信号添加到bin的2D矩阵中

Question

我正在尝试将复数awgn噪声添加到复数的2D矩阵中，我想将该2D矩阵用作箱，以便稍后估计噪声的pmf

import numpy as np    
x = np.array([(1+ 0j),(-0.50+.87j), (-0.5-0.87j)]) 
Nx = len(x)
card_Y = 80
###############Calculating histogram bins
yr = np.linspace(-2, +2, card_Y)
yj = np.linspace(-2, 2, card_Y)
yimag, yreal = np.meshgrid(yr, yj, indexing='ij')
yimag = yimag*-1j
histogram_bins = yreal +yimag
#histogram_bins = np.reshape(histogram_bins, (1,card_Y*card_Y), order='C')
###############Calculate complex AWGN noise noise
samples = 1000000
SNR_dB = -10

SNR = 10**(SNR_dB / 10)
Esym = np.sqrt(2)
N0 = Esym / SNR
N0_db = 10 * np.log10(Esym / N0)

noiseSigma = np.sqrt(N0)
n = noiseSigma * (np.random.randn(1, samples) + 1j* np.random.randn(1, samples))

由于无法弄清楚如何对复杂数字使用直方图函数，我感到困惑。也许我的方法在估计pmf方面是错误的，如果有人可以提供帮助，我将不胜感激。