質問編集履歴

2021/10/01 11:49

投稿

スコア9

title CHANGED Viewed

File without changes

body CHANGED Viewed

@@ -9,21 +9,17 @@
 ### リンクの周波数空間でのカットオフのコード
 ```python
-def LPF_MAM(x,times,tau = 0.01):
+def  LPF_CF(x,times,fmax):
-    k = np.round(tau /(times[1] - times[0])).astype(int)
+    freq_X = np.fft.fftfreq(times.shape[0],times[1] - times[0])
+    X_F = np.fft.fft(x)
+    X_F[freq_X>fmax] = 0
+    X_F[freq_X<-fmax] = 0
+#     虚数は削除
-    x_mean =  np.zeros(x.shape)
+    x_CF  = np.fft.ifft(X_F).real
-    N = x.shape[0]
-    for i in range(N):
-        if  i-k//2 <0 :
-            x_mean[i]  = x[: i - k//2 +k].mean()
-        elif i - k//2 +k>=N:
-            x_mean[i]  = x[i - k//2 :].mean()
-        else :
-            x_mean[i]  = x[i - k//2 : i - k//2 +k].mean()
-    return x_mean
+    return x_CF
-#tau = 0.035(sin wave), 0.051(step)
+#fmax = 5(sin wave), 13(step)
-x_MAM = LPF_MAM(x,times,tau)
+x_CF = LPF_CF(x,times,fmax)
 ```
 ###使用環境