数值分析验指导书.docVIP

PAGE PAGE 23 《数值分析》 实 验 指 导 书 潍坊学院数学与信息科学学院 2012年04月 目 录 TOC \o 1-1 \h \z \u HYPERLINK \l _Toc323327426 目录 PAGEREF _Toc323327426 \h I HYPERLINK \l _Toc323327427 实验一 插值与曲线拟合的最小二乘法 PAGEREF _Toc323327427 \h 1 HYPERLINK \l _Toc323327428 实验二 数值积分 PAGEREF _Toc323327428 \h 4 HYPERLINK \l _Toc323327429 实验三 解线性方程组的直接法 PAGEREF _Toc323327429 \h 9 HYPERLINK \l _Toc323327430 实验四 解线性方程组的迭代法 PAGEREF _Toc323327430 \h 11 HYPERLINK \l _Toc323327431 实验五 非线性方程的数值解法 PAGEREF _Toc323327431 \h 13 HYPERLINK \l _Toc323327432 实验六 常微分方程数值解法 PAGEREF _Toc323327432 \h 17 实验一 插值与曲线拟合的最小二乘法 一、实验目的： 1．了解拉格朗日插值法、牛顿插值法、曲线拟合最小二乘法的基本原理和方法； 2．掌握拉格朗日插值多项式牛顿插值多项式的用法； 3．掌握最小二乘原理，会求拟合函数及超定方程组的最小二乘解。 二、实验内容： 1．用拉格朗日插值公式和牛顿插值公式确定函数值； 2．对函数f (x)进行拉格朗日插值和牛顿插值； 3．利用Polyfit拟合幂函数，利用Polyfit拟合多项式。 三、实验过程： 1．给定函数四个点的数据如下: SKIPIF 1 0 ， 试用插值公式确定函数在 SKIPIF 1 0 处的函数值 SKIPIF 1 0 。 MATLAB程序如下： X=[1.1,2.3,3.9,5.1]; Y =[3.877,4.726,4.651 ,2.117]; p1=poly(X(1)); p2=poly(X(2)); p3=poly(X(3)); p4=poly(X(4)); l01= conv ( conv (p2, p3), p4)/(( X(1)- X(2))* ( X(1)- X(3)) * ( X(1)- X(4))), l11= conv ( conv (p1, p3), p4)/(( X(2)- X(1))* ( X(2)- X(3)) * ( X(2)- X(4))), l21= conv ( conv (p1, p2), p4)/(( X(3)- X(1))* ( X(3)- X(2)) * ( X(3)- X(4))), l31= conv ( conv (p1, p2), p3)/(( X(4)- X(1))* ( X(4)- X(2)) * ( X(4)- X(3))), l0=poly2sym (l01), l1=poly2sym (l11),l2=poly2sym (l21), l3=poly2sym (l31), P = l01* Y(1)+ l11* Y(2) + l21* Y(3) + l31* Y(4), 运行后输出的基函数l0，l1，l2和l3为 l0 =-1/24*x^3+1/8*x^2-1/12*x，l1 =1/4*x^3-1/4*x^2-x+1 l2 =-1/3*x^3+4/3*x，l3 =1/8*x^3+1/8*x^2-1/4*x 输入程序 L=poly2sym (P),x=2.101; Y = polyval(P,x) 运行后输出插值多项式和插值为 L=-629/5376*x^3+31433/53760*x^2-63029765850741/281474976710656*x+2010616283501353/562949953421312 Y =4.5969 输入程序 L=poly2sym (P),x=4.234; Y = polyval(P,x) 运行后输出插值多项式和插值为 L=-629/5376*x^3+31433/53760*x^2-63029765850741/281474976710656*x+2010616283501353/562949953421312 Y =4.2244 L=145616387951645/9007199254740992*x^3-2517512191700115/4503599627370496*x^2+14477/6000*x+2