質問編集履歴
1
回答のように変更
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -13,6 +13,7 @@
|
|
|
13
13
|
#それを時系列とともにplotする。
|
|
14
14
|
###############################
|
|
15
15
|
import numpy as np
|
|
16
|
+
import scipy
|
|
16
17
|
from scipy import integrate
|
|
17
18
|
import matplotlib.pyplot as plt
|
|
18
19
|
from sympy import *
|
|
@@ -77,12 +78,13 @@
|
|
|
77
78
|
def integrate(vlf,lf,hf,F_abs_amp,fq):
|
|
78
79
|
|
|
79
80
|
|
|
80
|
-
LF= integrate.quad(F_abs_amp, (fq,vlf,lf))
|
|
81
|
+
LF= scipy.integrate.quad(F_abs_amp, (fq,vlf,lf))
|
|
81
|
-
HF= integrate.quad(F_abs_amp,(fq,lf,
|
|
82
|
+
HF= scipy.integrate.quad(F_abs_amp, (fq,vlf,lf))
|
|
82
83
|
LF_HF = float(LF) / HF
|
|
84
|
+
print(LF_HF)
|
|
83
85
|
|
|
86
|
+
return LF_HF
|
|
84
87
|
|
|
85
|
-
return LF_HF
|
|
86
88
|
|
|
87
89
|
|
|
88
90
|
def FFT(data_input, dt, window_F):
|
|
@@ -156,7 +158,7 @@
|
|
|
156
158
|
vlf = 0.04
|
|
157
159
|
lf = 0.15
|
|
158
160
|
hf = 0.4
|
|
159
|
-
|
|
161
|
+
|
|
160
162
|
FFT_main(t, x, dt, split_t_r, overlap, window_F,output_FN, y_label, y_unit)
|
|
161
163
|
```
|
|
162
164
|
|