数值分析上机实习题.docVIP

2019-2020 第1学期数值分析上机实习题 总目标：会算，要有优化意识。（以下程序要求以附件1例题代码格式给出） 1. 对给定的线性方程组进行迭代求解。 （1）给出Jacobi迭代的通用程序。 （2）给出Gauss-Seidel迭代的通用程序。 调用条件：系数矩阵，右端项，初值，精度要求。 输出结果：方程组的近似解。 给定线性方程组，和，取初值为0， 分别利用Jacobi迭代和G-S迭代进行求解，观察并解释其中的数学现象。 2. 利用紧凑格式（即直接分解法或逐框运算法）对给定的矩阵进行Doolittle分解，并用其求线性方程组的解。 调用条件：矩阵。 输出结果：单位下三角矩阵和上三角矩阵。 给定矩阵，利用以下算法： 1）将A作Doolittle分解， 2）令，并对作Doolittle分解， 3）重复2）的过程令，并对作Doolittle分解，， 观察，，的变化趋势，思考其中的数学现象。 3. 给定函数，取，用等距节点 对原函数进行多项式插值和五次多项式拟合，试画出插值和拟合曲线，并给出数学解释。 4. 给出迭代法求非线性方程的根的程序。 调用条件：迭代函数，初值 输出结果：根的近似值和迭代次数 给定方程，用迭代格式计算附近的根，要使计算结果具有四位有效数字，利用估计式，或估计式来判断需要的迭代次数，分别需要迭代多少次？两者是否有冲突？ 5. 利用数值求积算法计算。 （1）利用Romberg算法计算积分。 调用条件：被积函数，精度要求。 输出结果：定积分的近似值。 （2）将区间进行等分，每个小区间上利用点高斯型求积公式计算。 调用条件：被积函数，区间等分数，每个子区间的高斯点数。 输出结果：定积分的近似值。 给定，利用上述方法进行计算，并与准确值进行比较。 6. 给定常微分方程， （1）给出Runge-Kutta 4阶算法的通用程序。 （2）给出Adams显式4阶公式的通用程序（初值由Runge-Kutta 4阶算法提供）。 调用条件：，区间，初值，步长。 输出结果：节点处函数的近似值。 针对初值问题，分别取步长计算各点（）的值，并与准确值比较，解释其中的数学现象。 请各位同学在规定的时间内按老师要求提交给老师。 具体要求如下： 1．应提交一份完整的实习作业。具体要求如下： （1）排版应美观漂亮；要标明姓名、学号、专业和联系电话； （2）要有目录，指明题目、程序、计算结果，图表和分析等内容所在位置，做到信息简明而完全； （3）尽量清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分的功能和作用。 2．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。 4．拷贝或抄袭他人结果是不良行为，将视为不合格。 5．请按任课老师要求的时间和载体（电子或纸质）提交给任课老师。 6. 提交电子邮件请务必按下述格式命名邮件主题和文档名，以方便统计 学号姓名数值分析1920第1学期计算作业 如：2019201597张帅哥数值分析1920第1学期计算作业 附件1 例题1. 给出二分法求函数零点的程序，要求如下： 调用条件：输入函数表达式f(x)，隔根区间[a,b] 输出结果：零点的值x和精度e. 给出调用说明和例子。 提交的合格答案如下： 二分法求函数零点的程序（函数型M文件）代码如下： function [c,err,yc]=bisect(f,a,b,delta) %Input - f is the function input as a stringf % - a and b are the left and right end points % - delta is the tolerance %0utput- c is the zero point % - yc=f(c) % - err is the error estimate for c ya=feval(f,a); yb=feval(f,b); if ya*yb0,break,end max1=1+round((log(b-a)-log(delta))/log(2)); for k=1:max1 c=(a+b) /2; yc=feval(f,c); if yc==0,a=c;b=c; elseif yb*yc0,b=c;yb=yc; else a=c;ya=yc; end if b-a delta, break,end end c=(a+b)/2; er