質問編集履歴
2
手で解いて実装しました
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1 +1,25 @@
|
|
|
1
1
|
与えられたデータについて2次曲線ax^2 +bでフィッティングしたいのですが、(一次の部分はなし)何か方法はないでしょうか。polyfitではできないようです。
|
|
2
|
+
|
|
3
|
+
|
|
4
|
+
|
|
5
|
+
皆様ありがとうございました。最小二乗法を手で解いてきました。以下のコードでです。大変お世話になりました。
|
|
6
|
+
|
|
7
|
+
|
|
8
|
+
|
|
9
|
+
|
|
10
|
+
|
|
11
|
+
```
|
|
12
|
+
|
|
13
|
+
# x, yはarray. A, Bはy = A*x**2 + Bの値
|
|
14
|
+
|
|
15
|
+
n = len(x)
|
|
16
|
+
|
|
17
|
+
A = ( -n*np.dot(x**2,y)+ np.sum(y)*np.sum(x**2) )/( (np.sum(x**2))**2 - n*np.sum(x**4) )
|
|
18
|
+
|
|
19
|
+
B = ( -np.sum(x**2)*np.dot(x**2,y) + np.sum(x**4)*np.sum(y) )/( n*np.sum(x**4) - (np.sum(x**2))**2.0 )
|
|
20
|
+
|
|
21
|
+
|
|
22
|
+
|
|
23
|
+
|
|
24
|
+
|
|
25
|
+
```
|
1
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
与えられたデータについて2次曲線ax^2 +bでフィッティングしたいのですが、(一次の部分はなし)何か方法はないでしょうか。polyfitではできないようです
|
|
1
|
+
与えられたデータについて2次曲線ax^2 +bでフィッティングしたいのですが、(一次の部分はなし)何か方法はないでしょうか。polyfitではできないようです。
|