inhermitemainmodule.bas

几种常用的内插数值计算方法

Attribute VB_Name = "MainModule"
Sub Main()
    Dim i As Integer
    Dim t As Double, z As Double, x(10) As Double, y(10) As Double, dy(10) As Double
    Dim s As String
        
    ' 给定结点
    x(1) = 0.1
    x(2) = 0.15
    x(3) = 0.3
    x(4) = 0.45
    x(5) = 0.55
    x(6) = 0.6
    x(7) = 0.7
    x(8) = 0.85
    x(9) = 0.9
    x(10) = 1#
    
    y(1) = 0.904837
    y(2) = 0.860708
    y(3) = 0.740818
    y(4) = 0.637628
    y(5) = 0.57695
    y(6) = 0.548812
    y(7) = 0.496585
    y(8) = 0.427415
    y(9) = 0.40657
    y(10) = 0.367879
    
    ' 求示例函数的一阶导数值
    For i = 1 To 10
        dy(i) = -y(i)
    Next i
    
    ' 插值位置
    t = 0.356
    ' 插值
    z = INHermite(10, x, y, dy, t)
    s = s & "f(" & t & ") = " & z & Chr(13)
    
    MsgBox s
    
End Sub