数值算法的稳定性教材.pptVIP

有效数字 如果|e| = |x* - x|? 1/2 ?10m-n 称近似数x准确到小数点后第n位,从这小数点后第n位数字直到最左边非零数字之间的所有数字都称为有效数字. 截断误差用我们原有的数学思维方式就比较容易理解的，而舍入误差则是本课程引入的一个新概念。要真正理解舍入误差，特别是它在计算中的传播及最终对计算结果的影响，是初步具备科学计算能力的重要标志。希望大家在完成实验后，认真仔细去体会截断误差和舍入误差的含义及对计算结果的影响。 图 0 10-7 10-5 10-5 10-4 10-3 10-2 10-1 10-3 10-1 误差 步长h 实验1.2（误差传播与算法稳定性）（P13 12题） 实验目的：体会算法的稳定性在选择算法中的地位 问题提出：考虑一个简单的积分序列 n=1,2,… 显然In0,n=1,2,…当n=1时，得： 当n≥2时，由分部积分可得： n=2,3,… 另外，还有： 实验内容：由递推关系In=1-nIn-1，可得计算积分序列{In}的两种算法： 算法一、直接使用递推公式得： In=1-nIn-1 n=2,3… 算法二、利用递推公式变形得： 实验要求：用上述两种算法分别在计算中采用5位、6位和7位有效数字，请判断哪种算法给出的结果更精确. 实验分析： 两种算法的优劣可能与你的第一感觉完全不同。 设算法一中I1的误差为e1，由I1递推计算In的误差为en 算法二中IN的误差为εN，由IN向前递推计算In (nN)的误差为εn 如果假定上述两种算法在后面的计算中都不再引入其他误差，则有： 由此可见，算法一中的e1可能很小，但在计算中它的影响急剧扩散，结果真实的数据很快被淹没了；而算法二中的εN可能相对比较大，但在计算中误差影响不扩散，即某一步产生误差后，该误差对后面的影响不断衰减。 误差扩散的算法是不稳定的，是我们所不期望的；误差衰减的算法是稳定的，是我们努力寻求的，也是贯穿本课程始终的目标。 熟悉几个必备的matlab命令： 1、help (帮助查找已知命令) 例如：help plot3 (显示命令plot3的功能) 2、lookfor(帮助查找未知命令) 例如：lookfor contour (显示含有contour的所有函数) 3、exist(检验函数) exist(‘A’) %A是变量 返回0，代表变量不存在 返回其它数值，变量存在 第一次实验内容 熟悉几个必备的matlab命令： 1、help (帮助查找已知命令) 例如：help plot3 (显示命令plot3的功能) 2、lookfor(帮助查找未知命令) 例如：lookfor contour (显示含有contour的所有函数) 3、exist(检验函数) exist(‘A’) %A是变量 返回0，代表变量不存在 返回其它数值，变量存在 作业：P13 9、12 实验内容 实验1.1、实验 1.2（P1312） 作业：P12—13 1、3、 4、6、8 老教材 实验1.1（数值微分精度与步长的关系） 实验目的：数值计算中误差是不可避免的，要求通过本实验初步认识数值分析中两个重要概念：截断误差和舍入误差，并认真体会误差对计算结果的影响。 问题提出：设一元函数f :R→R，则f在x0的导数定义为: 实验内容:根据不同的步长h可设计两种算法,计算f在x0处的导数 计算一阶导数的算法: 请给出两个计算高阶导数的近似算法,并完成如下工作: (1) (2) 1、对同样的h，比较(1)式和(2)式的计算结果； 2、针对计算高阶导数的算法，比较h取不同值时(1)式和(2)式的计算结果； 实验要求：选择有代表性的函数f (x)（最好多选择几个），利用Matlab提供的绘图工具画出该函数在某区间的导数曲线f (s)(x)，再将数值计算的结果用Matlab 画出来，认真思考实验的结果。 function y=dsh1(fu,x,h) y=(feval(fu,x+h)-feval(fu,x))/h; y 函数文件(扩展名为.m),可以被其命令调用 格式: function输出参数=函数名(输入参数) 可以是多个 函数名 一个 注意: 1、保存的文件名与函数名要一致，这样才能保证调用成功； 2、调用时的输入输出参数要与函数中的一样； 3、养成良好的注释习惯，便于自己或用户调用定义的函数。 注释语句(前面有%) 程序语句 function y=dsh2(fu,x,h) y=(feval(fu,x+h)-feval(fu,x-h))/(2*h); y 调用格式：dsh2(‘sin,1,0.01) 调用格式：dsh1(