４章前半を追加

Author: 153hashimoto

khashimoto/chapter04/q01.py

@@ -5,21 +5,27 @@
 
 def zeropad(L, S, x):
     zero1 = np.zeros(L - S)
-    x_pad = np.concatenate([zero1, x, zero1])  # 1.
+    # x_pad = np.concatenate([zero1, x, zero1])  # 1.
+    x_pad = np.concatenate([zero1, x])
     zero2 = x_pad.size % S
     if zero2 != 0:
         x_pad = np.concatenate([x_pad, np.zeros(S - zero2)])  # 2.
+    x_pad = np.concatenate([x_pad, np.zeros(L - S)])
     return x_pad
 
 
 """
 # 確認
 x = np.ones(5)
 print(x)
-x_pad = zeropad(4, 3, x)
+x_pad = zeropad(4, 2, x)
 print(x_pad)
 
 
 # [1. 1. 1. 1. 1.]
 # [0. 1. 1. 1. 1. 1. 0. 0. 0.]
+
+# [1. 1. 1. 1. 1.]
+# [0. 0. 1. 1. 1. 1. 1. 0. 0. 0.]
+# 不足分が出ないように零詰め
 """

khashimoto/chapter04/q02.py

@@ -18,7 +18,8 @@ def frame(L, S, x):
 # 確認
 x = np.ones(5)
 print(x)
-x_2 = frame(4, 3, x)
+print(zeropad(4, 2, x))
+x_2 = frame(4, 2, x)
 print(x_2)
 
 

khashimoto/chapter04/q04.py

@@ -17,7 +17,7 @@
 
 
 # プロット
-T = np.arange(len(x_stft)) * 0.1  # 時間軸
+T = np.arange(len(x_stft))  # 時間軸
 fr = np.arange(x_stft.shape[1]) / x_stft.shape[1] * fs / 2  # 周波数軸
 
 plt.subplot(1, 2, 1)

khashimoto/chapter04/q05.py

@@ -1,6 +1,7 @@
 # 合成窓
 
 import numpy as np
+import matplotlib.pyplot as plt
 
 
 def window(S, w):
@@ -18,8 +19,10 @@ def window(S, w):
 
 """
 # 確認
-w = np.ones(6)
-print(window(2, w))
+L = 1024
+w = 0.54 - 0.46 * np.cos(2 * np.pi * np.arange(L) / (L - 1))  # Hamming窓
+plt.plot(window(256, w))
+plt.show()
 
 
 # [0.33333333 0.33333333 0.33333333 0.33333333 0.33333333 0.33333333]