Forum

PolynomialFit

1 Posts
1 Users
0 Reactions
302 Views
Member
Joined: 2 years ago
Posts: 4
Topic starter  

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}}

   
ReplyQuote

Leave a reply

Author Name

Author Email

Title *

Preview 0 Revisions Saved