質問編集履歴
1
誤字
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -38,8 +38,6 @@
|
|
|
38
38
|
|
|
39
39
|
i = 4
|
|
40
40
|
|
|
41
|
-
# うまくいかないポイント
|
|
42
|
-
|
|
43
41
|
for o, g in zip(name, initial):
|
|
44
42
|
|
|
45
43
|
|
|
@@ -58,13 +56,9 @@
|
|
|
58
56
|
|
|
59
57
|
|
|
60
58
|
|
|
61
|
-
# 水頭
|
|
62
|
-
|
|
63
59
|
a = data[o].iloc[i:i+400,1]
|
|
64
60
|
|
|
65
61
|
|
|
66
|
-
|
|
67
|
-
# 体積含水率
|
|
68
62
|
|
|
69
63
|
d = data[o].iloc[i:i+400,2]
|
|
70
64
|
|
|
@@ -78,15 +72,11 @@
|
|
|
78
72
|
|
|
79
73
|
|
|
80
74
|
|
|
81
|
-
# データにおける最大含水量
|
|
82
|
-
|
|
83
75
|
m = data[o].iloc[5:,2].max()
|
|
84
76
|
|
|
85
77
|
s = data[o].iloc[5:,2].min()
|
|
86
78
|
|
|
87
79
|
|
|
88
|
-
|
|
89
|
-
#フィッティング関数の指定
|
|
90
80
|
|
|
91
81
|
def fitfunc(h, x, z, n ):
|
|
92
82
|
|
|
@@ -96,11 +86,9 @@
|
|
|
96
86
|
|
|
97
87
|
|
|
98
88
|
|
|
99
|
-
para_ini =[g, 0.5, 1.5]
|
|
89
|
+
para_ini =[g, 0.5, 1.5]
|
|
100
90
|
|
|
101
91
|
|
|
102
|
-
|
|
103
|
-
# 拘束条件
|
|
104
92
|
|
|
105
93
|
para_bounds=([0,-np.inf,1],[s,1,np.inf])
|
|
106
94
|
|
|
@@ -110,8 +98,6 @@
|
|
|
110
98
|
|
|
111
99
|
|
|
112
100
|
|
|
113
|
-
#得られたフィッティングパラメータの表示
|
|
114
|
-
|
|
115
101
|
print ("Fitted x =", str(para_opt[0]))
|
|
116
102
|
|
|
117
103
|
print ("Fitted z =", str(para_opt[1]))
|
|
@@ -119,8 +105,6 @@
|
|
|
119
105
|
print ("Fitted n =", str(para_opt[2]))
|
|
120
106
|
|
|
121
107
|
|
|
122
|
-
|
|
123
|
-
# for plot
|
|
124
108
|
|
|
125
109
|
x=para_opt[0]
|
|
126
110
|
|