From 42040f6c0d168a3ebe0f3b6a2037d4c7bfd08fa0 Mon Sep 17 00:00:00 2001 From: Marc Lichtman Date: Mon, 22 Jan 2024 01:54:24 -0500 Subject: [PATCH] fix error in MVDR section --- content/doa.rst | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/content/doa.rst b/content/doa.rst index 1a5bf1c..e372269 100644 --- a/content/doa.rst +++ b/content/doa.rst @@ -416,8 +416,7 @@ We can implement the equations above in Python fairly easily: Rinv = np.linalg.pinv(R) # 3x3. pseudo-inverse tends to work better than a true inverse w = 1/(a.conj().T @ Rinv @ a) # MVDR equation! denominator is 1x3 * 3x3 * 3x1 - metric = metric[0,0] # convert the 1x1 matrix to a Python scalar, it's still complex though - metric = np.abs(metric) # take magnitude + metric = np.abs(w[0,0]) # take magnitude metric = 10*np.log10(metric) # convert to dB so its easier to see small and large lobes at the same time results.append(metric)