Topic starter
12/15/2022 15:06
PolynomialFitがうまくいかないことがある。下記が例。
- y = a0 + a1 * x + a2 * x2 + a3 * x3 + a4 * x4 + a5 * x5 + a6 * x6
- a0 = 2、a1 = 0.1、a2 = 0.02、a3 = 0.003、a4 = 0.0004、a5 = 0.00005、a6 = 0.000006
- というリストを作って、x=0,2,4,6,...20 で11個の点で多項式フィットする。
dat ={ {0,2.00000000000E+00}, {2,2.31238400000E+00}, {4,3.09017600000E+00}, {6,5.15513600000E+00}, {8,1.04656640000E+01}, {10,2.30000000000E+01}, {12,4.99159040000E+01}, {14,1.02986816000E+02}, {16,2.00314496000E+02}, {18,3.68318144000E+02}, {20,6.44000000000E+02}}; Print[PolynomialFit[dat,6]];
- メイントランク版(PFサーバ上の古いものを含む)のSADの結果(PF計算機サーバ、KEKBサーバ) → よく収束しない
Welcome to SAD Ver.1.0.10.10.2k64 built at 2016-09-15 00:08:51 +0900 {{ 1.9065249634, .4170309801, -.1378238690, .0330817610, -.0022752014, 1.619098253E-4, 4.218062455E-6},{Residual-> .0828379198}} Welcome to SAD Ver.1.0.10.5.14a5 built at 2017-08-08 17:49:32 +0900 {{ 1.9065249634, .4170309801, -.1378238690, .0330817610, -.0022752014, 1.619098253E-4, 4.218062455E-6},{(Residual-> .0828379198)}} Welcome to SAD Ver.1.0.10.10.1a built at 2021-10-06 13:54:01 +0900 Today: 15:04:05 Thursday 12/15/2022 {{ 1.9065249634, .4170309801, -.1378238690, .0330817610, -.0022752014, 1.619098253E-4, 4.218062455E-6},{(Residual-> .0828379198)}}
- 生出さんバージョンの結果(PF計算機サーバ) → 良く収束している。
Welcome to SAD Ver.1.1.9.0.5k64 built at 2020-09-09 16:04:26 +0900 {{ 2.0000000000, .1000000000, .0200000000, .0030000000, 4.000000000E-4, 5.000000000E-5, 6.000000000E-6},{Residual-> 1.24370785E-13}}