上理机械计算机辅助设计课程线性、抛物线、拉格朗日插值代码.docVIP

上理机械计算机辅助设计课程线性、抛物线、拉格朗日插值代码 Public Class Form1 Function INLagrn(ByVal n As Integer, ByVal x() As Double, ByVal y() As Double, ByVal t As Double) As Double Dim i As Integer, j As Integer, k As Integer, m As Integer Dim z As Double, s As Double z = 0.0# If (n 1) Then INLagrn = z Exit Function End If If (n = 1) Then z = y(1) INLagrn = z Exit Function End If If (n = 2) Then z = (y(1) * (t - x(2)) - y(2) * (t - x(1))) / (x(1) - x(2)) INLagrn = z Exit Function End If i = 0 While ((x(i) t) And (i = n)) i = i + 1 End While k = i - 4 If (k 0) Then k = 0 m = i + 3 If (m n - 1) Then m = n - 1 For i = k To m s = 1.0# For j = k To m If (j i) Then s = s * (t - x(j + 1)) / (x(i + 1) - x(j + 1)) 拉格朗日插值公式 Next j z = z + s * y(i + 1) Next i INLagrn = z End Function Dim t As Double, z As Double Dim a(0 To 9) As Double Dim k(0 To 9) As Double Private Sub Form1_Load(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles MyBase.Load a(0) = 6.58 a(1) = 5.99 a(2) = 5.3 a(3) = 4.92 a(4) = 4.57 a(5) = 4.0 a(6) = 3.54 a(7) = 3.23 a(8) = 2.85 a(9) = 2.1 k(0) = 4.23 k(1) = 4.0 k(2) = 3.56 k(3) = 3.21 k(4) = 3.02 k(5) = 2.56 k(6) = 2.35 k(7) = 1.95 k(8) = 1.55 k(9) = 1.04 End Sub Private Sub TextBox3_TextChanged(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles TextBox3.TextChanged End Sub Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button1.Click If TextBox1.Text = Then Text