插值法VB源程序

1、第四章插值方法VB一、線性插值應(yīng)用實(shí)例VB1已知：169的平方根為13,196的平方根為14，編程求175的平方根。Private Sub Comma nd1_click()Dim x(1 To 2),y(1 To 2) As Si nglex(1)=169:y=13x(2)=196:y(2)=14x0=175yO=F(x,y,xO)a=Format$(yO,#.#)Text1.Text=yO= & Str$(a)End SubPrivate Fun cti on F(x,y,xO)F=y(1)+(y (2)-y(1)/(x (2)-x(1)*(x0-x(1)End Fun cti on線性插

2、值應(yīng)用實(shí)例VB2已知水的溫度與密度和關(guān)系如下:溫度5.010.015.020.025.030.0密度0.99990.9970.990.990.990.99試編程計(jì)算溫度為 6、13、16、19、23、26 oC時的密度。Private Sub Comma nd1_click()N = 6X = Array(0, 5, 10, 15, 20, 25, 30)Y = Array(0, 0.999965, 0.9997, 0.999099, 0.998203, 0.997044, 0.995646)xc = Array(0, 6, 13, 16, 19, 23, 26)List1.Addltem x

3、0 & y0For I = 1 To 6X0 = xc(I)Y0 = f(N, X, Y , X0)List1.AddItem Str$(X0) & Str$(Y0)Next IEnd SubPrivate Function f(N, X, Y , X0)For I = 1 To NIf X0 = X(I) The nW = IExit ForEnd IfNext IW = W - 1If W = N ThenW = N - 1: End Iff = Y(W) + (Y(W + 1) - Y(W) / (X(W + 1) - X(W) * (XO - X(W)End Fun cti on二、拉

4、格朗日三點(diǎn)插值實(shí)例例：實(shí)驗(yàn)測得25OC乙醇溶液的平均摩爾體積v與乙醇的摩爾分?jǐn)?shù)xi的關(guān)系數(shù)據(jù)如下表。試編程求摩爾數(shù)xi為0.3時的乙醇溶液的平均體積v。VB清單Dim X, v計(jì) Wv21.2222.1623.1824.3225.5726.95xi0.0890.1150.1430.1730.2060.242v28.4730.1532.0134.0736.37xi0.2810.3230.3690.4200.477找牯巾|日二點(diǎn)Jit.怛Private Sub Command1_click()N = 11X = Array(0, 0.0891, 0.1153, 0.1435, 0.1739, 0.

5、2068, 0.2424, 0.2811, 0.3234, 0.3697, 0.4207, 0.4771)v = Array(0, 21.22, 22.16, 23.18, 24.32, 25.57, 26.95, 28.47, 30.15, 32.01, 34.07, 36.37)M = 3X0 = Val(Text1.Text)If X0 X(N) ThenPrint errorElseEnd IfCall pt(N, M, XO, A, B)vO = F(XO, A, B)Listl. Additem a= & AListl. AddItem b= & BList1. AddItem X

6、0= & X0List1. AddItem v0= & Format$(v0, #.#)End SubPrivate Sub pt(N, M, XO, A, B)For I = 1 To NIf XO = X(I) ThenIf Abs(X0 - X(I - 1) = Abs(X0 - X(I) ThenA = I - Int(M + 1) / 2)Exit ForEnd IfElseA = I - Int(M / 2)End IfNext IIf A = N ThenB = N: End IfA = B - M + 1End SubPrivate Function F(X0, A, B)Fo

7、r I = A To BW = 1For J = A To BIf J I ThenW = W * (XO - X(J) / (X(I) - X(J): End IfNext Js = s + W * v(I)Next IF = sEnd Function三、拉格朗日全點(diǎn)插值VB試對上例編程對比利用三個、五個、七個、九個和全部節(jié)點(diǎn)時，摩爾分?jǐn)?shù) xi為0.30時乙醇溶液的平均摩爾體積v。Dim X, vPrivate Sub Comma nd1_click()N = 11X = Array 0.0891, 0.1153, 0.1435, 0.1739, 0.2068, 0.2424, 0.281

8、1, 0.3234, 0.3697, 0.4207, 0.4771)v = Array(0, 21.22, 22.16, 23.18, 24.32, 25.57, 26.95, 28.47, 30.15, 32.01, 34.07, 36.37)X0 = Val(Text1.Text)If X0 X(N) The nPrin t errorEndElseEnd IfList1.Addltem x0= & X0For M = 3 To N Step 2Call pt(N, M, X0, A, B)v0 = F(X0, A, B)List1.AddItem M= & M & ”A= & A & ”B= & B & ”v0 & Format$(v0, #.#)Next MEnd SubPrivate Sub pt(N, M, X0, A, B)For I = 1 To NIf X0 = X(l) The nIf Abs(X0 - X(l - 1) = Abs(X0 - X(l) The nA = I - I nt(M + 1) / 2)Exit ForEnd IfElseA = I - In t(M / 2)End IfNext IIf A = N The nB = N: End IfA = B - M + 1End SubPrivate Fu nction