11/5のGoogleトップページはAICの赤池さん、ということで多項式などのフィッティングができるVBAプログラムをUP!(Numerical Recipes移植版）

トップページがAICで有名な赤池さんになっていて驚いた。

そこでAICの計算をカシオの高精度計算サイトにでもUPしようか、と思ったけど多数のデータを読み込ませるのができない、、、

じゃあExcel VBAででもやるか、と思ったけどその前段階のフィッティングする関数とかのほうが便利かということでそちらを紹介。

昔、Numerical Recipes in Cのlfit関数をVBAに移植していたのだ。

こちら。

「lfitLib.bas」をダウンロード

これを標準モジュールとして読み込む。

使い方は、例えば、Excelではできない8次多項式のフィッティングをやってみる。

データはExcelのセルに100個入力されているとして、、、

Option Explicit
Option Base 1

Private Sub CommandButton1_Click()
    Dim x() As Double, y() As Double, sig() As Double
    Dim a() As Double, ia() As Integer
    Dim n As Integer, m As Integer, i As Integer
    Dim covar() As Double, chisq As Double
    
    n = 100
    m = 9
    
    ReDim x(n), y(n), sig(n)
    ReDim a(m), ia(m), covar(m, m)
    
    For i = 1 To n
        x(i) = CDbl(Worksheets("Sheet1").Cells(i + 1, 4))
        y(i) = CDbl(Worksheets("Sheet1").Cells(i + 1, 5))
        sig(i) = 1#
    Next i
    
    For i = 1 To m
        ia(i) = 1
    Next i
    
    
    Call lfit(x, y, sig, n, a, ia, m, covar, chisq)
    
    For i = 1 To m
        Worksheets("Sheet1").Cells(i + 3, 2) = a(i)
    Next i
    
    For i = 1 To n
        
        x(i) = CDbl(Worksheets("Sheet1").Cells(i + 1, 4))

    Next i
    
End Sub

のように使う。引数の詳しい説明は、Numerical Recipesを読んでください。

無償で読めます。

http://apps.nrbook.com/c/index.html

の15章、Modeling of Dataのところ。

残念ながらVBAでは関数を引数に取れないので、フィッティング関数はfuncsに決め打ちです、、、ここを変えたら多項式からより一般の

y(x)=Σa_k*X_k(x)

のフィッティングができます。

その他も移植してますが、またおいおい、、、