数值分析-练习题.docx

《数值分析-练习题.docx》由会员分享，可在线阅读，更多相关《数值分析-练习题.docx（8页珍藏版）》请在第一文库网上搜索。

1、JTY1. 5, 2上的二次最小二乘拟合多项式.函数/(x) = Sin,给出该函数在点0, 0. 5, 1,10给出公式中的参数,使其代数精度到达最高，并指明代数精度，用它计算积分JIogK)Xeh1 Mdx =等(/(。) + /() +(/()- fb). a写出下面线性代数方程组的Jacobi迭代格式,并说明它是否收敛用Euler公式求解下面的常微分方程初值问题,dy1- = x+V-I1 dxy(0) = l, 0xl取步长 = 0.2:利用Cholesky分解计算方程组的解251用牛顿插值法计算满足下面条件的三次多项式Ml) = 5/6,以2) = 17/3,2= 31/2,以4)

2、 = 97/3.下面题目为编程题，请勿手算，请只用MatIab语言.七、(10分)切比雪夫函数4(%)定义如下:(x) = l, 7；(X) = %,Tn+1(x) = 2xTn(x)Tnxy编程画出线(x)到4(%)在区间1, 1上的图像有两个Matlab函数文件定义如下： 文件1function V = f(x) v = - cos (x);文件2function Z = g(x)z = x*lOg(X)-1;请编写一个二分法求根的程序，给定初始区间，给出求解这两个函数根的恰当的调用方式. 解：一、给出一个收敛的迭代方法计算非线性方程XInX = 2的根，要求得到的根的近似值至少有5位有效

3、 数字.二、利用LU分解求解下面线性代数方程组-22-22三、求函数/(X) = J方在区间T, 1上的二次最正确平方逼近多项式. 1 + x四、函数y = Iog2 有函数值表如下1/21241012利用三次插值方法近似计算log2 3.五、确定下面积分公式中的参数，使其代数精度到达最高，并指明代数精度：2 /氏=%/(1) + o() +。-2下面题目为编程题，请勿手算，请只用MatIab语言.六、(16分)微分方程初值问题如下：Ay=cos j + log x, Ax，y=2,lx3.(1) 利用Matlab中的函数求出该微分方程初值问题的数值解，并画出其函数图像；(2) 利用上面小题的

4、计算结果，求函数y = y(x)在区间1, 3上的复合梯形积分的积分值.矩阵的维数由命令行输七、编写一个Gauss-Seidel迭代方法求解下面线性代数方程组AX = Z;的解, 入：(-2.I，I = J1z y,与=In J , lij乙 +kL(Xz 乙)+ /(Xz+1, Vil)计算常微分方程初值问题曳= (l + 2)y,O%l dxy(0) = 1,的数值解，取步长九=02数据13579%22112试用最小二乘法求形如5(%)= Q +匕欠的拟合曲线.以下两题为Matlab编程题：七、(10分)用牛顿法计算方程在炉-3-l = 在XO =2附近的根，要求满足精度Ixk +1 fi

5、d V Ie - 7八、12分用乘幕法计算矩阵4O21512973413342621的按模最大特征值与对应的特征向量.-3 2 113设A= 2 2 1 , b= 2 ,将A进行ChOliSky分解A = ZZr,其中，：下三角矩阵，并由1 1 1J1此求解线性方程组Ax =从对线性方程组2Xi 一 + % = 1 x1-x3+ 5x4 = 6 x2 + 4x3 -x4 = 8设法导出使Jacobi迭代和Gauss-Seidel迭代收敛的格式，要求分别写出格式，并说明 收敛的理由。确定系数。/,c和d,使得以下函数S(X)= lx22x33(x-l) + 2(x-l)2-(x-l)3,a +

6、b(x 2) + C(X 2)2 + d(x 2)3,是一个三次样条，且J(x)公最小。4、15分)求拟合以下数据的线性最小二乘多项式解。确定以下带权求积公式中的系数AB使其代数精度最高，指明代数精度: f (X)Gdx 4f(0.289949) + (0.821162)并用上述公式近似计算积分 e* GdX o推导以下数值微分公式f 4(X + z) - 3/(%) - J(X + 2切并证明它的误差项为C) 以下面的函数头完成牛顿法求方程的根的程序：function x,it = newton(f, g,x, to),其 中：/为非线性方程所对应的函数，g为/的一阶导数，为迭代初始值，力/

7、为允许精度。 并编写一段调用程序来求解/(X) = /+2泥，+产在_2,2内的一个根，要求初始值取XO=L 迭代精度控制国-/一 I10%以下面的函数头完成4阶Runge-Kutta方法解常微分方程初值问题的程序：function x, y = oderk4(f,a,b,h),其中：/为常微分方程的右端项，为求解区间的左右端点，Zz为自变量的步长，x,y分别是计算完成时的自变量取值和对应点上的函数值,并编写一段调用程序来求解以下常微分方程初值问题。步长力= 0.01=xy2, 0xl dxJ(O) = 1此处的4阶Runge-Kutta方法的公式如下：%+1 =+ = (+ 2 左2 + 2

8、 左3 + 左4) OK =好(%, %)=V(+p+y)k4=hf(xn+h,yn+k3)4-3 3设 A= 3 2-61 -5 34b= 7 ,将A进行LU分解，其中心单位下三角矩阵，U ： 3上三角矩阵，并由此求解线性方程组Ax = O.对线性方程组3jV 1OX2 = 7 9x1 - 4x2 = 51)问用JaCObi迭代和GaussSeidel迭代求解此方程组是否收敛？并说出理由.2)设法导出使JaCObi迭代和GaUSSSeidel迭代收敛的格式，要求分别写出格式，并 说明收敛的理由.确定系数。/,。和d,使得以下函数3(x-l) + 2(x-l)2-(x-l)3,lx2CI +

9、b(x 2) + C(X 2)2 + d (X 2)3,2 x 3是一个三次样条，同时具有性质SQ) =b(3).求以下数据表的最小二乘三次逼近多项式.X10123于(X)-31339161确定以下求积公式中的系数AI使其代数精度最高，指明代数精度:f(x)dx+并完成以下两个问题：(1) 用上述公式近似计算积分j e2dx.11.5(2) 将1, L 5分成两个区间，用它的复合求积公式近似计算积分J e*dx.1计算以下表达式中的系数ABC，/(x)4f()+w1)+()使得它对以下三个多项式能精确成立：L % -期,(X -/, 其中 XO %2, 再一 %o=z, x2-xl= ah.以下面的函数头完成牛顿法求方程的根的程序：function x,it = newton(f, g,x, tol),其 中：/为非线性方程所对应的函数，g为/的一阶导数，升为迭代初始值，加/为允许精度,X是迭代完成时的解，z1(%) = 3+22+oX-20在1,2内的一个根，要求初始值取XO=L迭代精度控制匕-以下面的函数头完成反塞法求实矩阵的按模最小的特征值和其对应的特征向量的程序： function x,lam,it = anpow(A,x,tol),其中：A为输入的矩阵，/为初始迭代特征向量, tol 为精度控制参数，XjammA -xm2 106.