q06,08,10修正

Author: kuu8902

kkoiso/chapter05/q06.py

@@ -23,7 +23,7 @@
 x2 = np.roll(s, 10) + epsilon
 x3 = np.roll(s, 20) + epsilon
 
-x = (x1 + x2 + x3) / 3
+x = x1 + x2 + x3
 plt.figure(figsize=(10, 6))
 plt.plot(t[: int(0.01 * fs)], x[: int(0.01 * fs)])
 plt.xlabel("Time [s]")

kkoiso/chapter05/q08.py

@@ -27,8 +27,6 @@ def plot_beam_pattern(w_f, p, fs):
 theta = 45
 F = 1000
 linear_coords = [((m - 1) * d, 0) for m in range(M)]
-
-
 fs = 16000
 
 

kkoiso/chapter05/q10.py

@@ -2,6 +2,7 @@
 import matplotlib.pyplot as plt
 from scipy.signal import stft
 from q04 import spatial_correlation_matrix
+from q01 import array_manifold_vector_linear
 
 
 def compute_spatial_spectrum(z, p, fs):
@@ -12,7 +13,7 @@ def compute_spatial_spectrum(z, p, fs):
     tau = np.array([0, 10 / fs, 20 / fs])
 
     for i, theta in enumerate(angles):
-        w = 1 / 3 * np.exp(-1j * 2 * np.pi * np.cos(np.deg2rad(theta)) * tau[:, None])
+        w = 1 / 3 * np.exp(-1j * 2 * np.pi * theta * tau)
         for f_bin in range(20, 31):
             P[f_bin, i] = np.abs(np.conj(w.T).dot(R_z[f_bin]).dot(w))