计算方法实验3（Computational method experiment 3）.doc

计算方法实验3（Computational method experiment 3） 拉格朗日插值法 # include stdio.h Void main () { The double L (double x); Double a, X; Scanf (% lf, x); A = l (x); Printf (% F,); } The double L (double x) {double the [5] = {0.4,0.5,0.6,0.7,0.8}; Double B [5] = {- 0.916291, 0.693147, 0.510826, 0.356675, -0.223144}; Double n = 0, l, k = 1, s = 1; Int w I; For (I = 0; I = 4; I + +) { For (w = 0, w = 4; W + +) { If (w = = I) Continue; K = k * (x - [w]); S = s * (a [i] - [w]); } L = K / s * B [i]; S = 1; K = 1; N = n + L; } The return (n); } 牛顿法 # include stdio.h # include stdlib.h # include iostream.h Typedef struct date { Float X; Float y; } / / 变量x和函数值y的结构 date; Date d [20]; / / 最多二十组数据 Float f (int s, int t) / / 牛顿插值法, 用以返回插商 { If (t = = S + 1) The return (d [t].Y - D [S].Y) / (d [t].X - D [S].X); Else The return (f (s + 1, t) - f (s, t) - 1) / (d [t].X - D [S].X); } Newton (float x, float, int count) { Int n; While (1) { Cout 请输入n值 (即n次插值):; / / 获得插值次数 CIN n; If (n = count - 1) / / 插值次数不得大于count－1次 Break; Else System (CLS); } / / 初始化t, y, YT. Float t = 1.0; Float y = D [0].Y; Float YT = 0.0; / / 计算y值 For (int j = 1; J = n; j + +) { T = (x - D [J-1].X) * t; YT = f (0, J * t; / / cout f (0, J) endl; Y = y + YT; } Return y; } Lagrange (float x, float, int count) { Float y = 0.0; For (int k = 0; K count; K + +) / / 这儿默认为count－1次插值 { Float P = 1.0; / / 初始化p For (int j = 0, J count; j + +) {/ / 计算p的值 If (k = J) continue; / / 判断是否为同一个数 P = P * (x - D [J].X) / (d [k].X - D [J].X); } Y = y + p * d [k].Y; / / 求和 } Return y; / / 返回y的值 } Void main () { Float x, y; Int count; While (1) { Cout 请输入x [i], y [i] 的组数, 不得超过20组:; / / 要求用户输入数据组数 CIN count; If (count = 20) Break; / / 检查输入的是否合法 System (CLS); } / / 获得各组数据 For (int i = 0; I count; I + +) { Cout 请输入第 I + 1 组x的值:; CIN d [i].X; Cout 请输入第 I + 1 组y的值:; CIN d [i].Y; System (CLS); } Cout 请输入x的值:; / / 获得变量x的值 CIN x; While (1) { Int choice = 3; Cout 请您选择使用哪种插值法计算: endl; Cout (0): 退出 endl; Cou