资源描述：
《数值分析课件典型例题与习题.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、《数值分析》典型例题I一、二章内容提要典型例题分析例题与练习题实验题介绍化大为小化繁为简化难为易核心的概念误差算法的构造与分析收敛性稳定性复杂度(时间与空间)等有效数字概念若近似值x的绝对误差限是某一位上的半个单位,该位到x的第一位非零数字一共有n位,则称近似值x有n位有效数字。从左向右看第一个非零数误差限不超过该位的半个单位n位有效数字如果x具有n位有效数字,则相对误差满足:其绝对误差满足:如果一个规格化浮点数则称近似数x具有n位有效数字。迭代法思想:收敛性收敛速度Iterate:Tosayordoagainoragainandagain例1.经过四舍五入得出x1=6.1025和x2
2、=80.100,试问它们分别具有几位有效数字?解：例2.已知近似数x有两位有效数字,试求其相对误差限。解：
3、er(x)
4、<5*10-2例3.如下近似值的绝对误差限均为0.005,问各近似值有几位有效数值x1=1.38,x2=-0.0312,x3=0.00086。例4.二次方程x2–16x+1=0,取求使具有4位有效数。解:直接计算x1≈8–7.937=0.063修改算法4位有效数计算出的x1具有两位有效数例5.采用迭代法计算,取x0=7(k=0，1，2，……)若xk具有n位有效数字,求证xk+1具有2n位有效数字。Ex1：对是否都有这一性质?例6.序列{yn}满足递推关系yn=10yn-1–1
5、(n=1,2,·····)若取y0=21/2≈1.41(三位有效数字)。递推计算y10时误差有多大？思考:这个计算过程稳定吗？例7.设y0=28,按递推公式yn=yn-1–(783)1/2/100(n=1,2,·····)计算到y100。若取(783)1/2≈27.982(5位有效数字),试问计算y100将有多大的误差？例8.设计算球体V允许其相对误差限为1%,问测量球半径R的相对误差限最大为多少?解：由球体计算公式分析误差传播规律故当球体V的相对误差限为1%时，测量球半径R的相对误差限最大为0.33%。相对误差传播规律例9.利用级数可计算出无理数的近似值。由于交错级数的部分和数列Sn在
6、其极限值上下摆动,试分析为了得到级数的三位有效数字近似值应取多少项求和。解:由部分和只需n>1000时,Sn有三位有效数。例10.在计算机上对调和级数逐项求和计算当n很大时，Sn将不随n的增加而增加。试分析原因。例11.证明方程1-x-sinx=0在区间[0,1]上有一根,使用二分法求误差不大于0.5*10-4的根需要二分多少次？提示:f(0)=1,f(1)=-sin1<0。且f′(x)=-1-cosx在区间(0,1]严格单调递减。例12.构造求ex+10x-2=0根的迭代法。提示:故迭代法算法一阶收敛。例13.应用牛顿迭代法于方程x3–a=0,导出求立方根的迭代公式,并讨论其收敛阶。解：令
7、f(x)=x3–a，则牛顿迭代公式故立方根迭代算法二阶收敛例14.设a为正实数,试建立求1/a的牛顿迭代公式,要求在迭代公式中不含有除法运算,并考虑迭代公式的收敛。xn+1=xn(2–axn)，(n=0，1，2……)所以,当
8、1–ax0
9、<1时,迭代公式收敛。解:建立方程利用牛顿迭代法，得1–axn+1=(1–axn)2整理，得例15.证明对于C>0,迭代格式例16.解:Ex2.若x*是f(x)=0的m重根,试证明修正的牛顿迭代法至少为二阶收敛。[f(x)]1/m或f(x)/f′(x)单根Ex3对于复变量z=x+iy的复值函数f(z)应用牛顿迭代公式时为避开复数运算,令zn=xn+iynf(z
10、n)=An+iBn，f′(zn)=Cn+iDn证明例17.提示:取初值x1=21/2,考虑序列单调有界,则该序列必有极限。例18.例19.已知方程x3-x2-1=0在x0=1.5附近有根,试判断下列迭代格式的收敛性。例20.证明由迭代格式xn+1=xn/2+1/xn产生的迭代序列{xn},对任意的x0>0,均收敛于21/2。牛顿迭代法的收敛域问题:用牛顿迭代法求解方程zd–1=0的复根。例如d=3时,方程在复平面上三个根分别是z1=1选择中心位于坐标原点，边长为2的正方形内的任意点作初始值，进行迭代，把收敛到三个根的初值分为三类，并分别标上不同颜色(例如红、绿和蓝)。对充分多的初始点进行实验，
11、绘出牛顿迭代法对该方程的收敛域彩色图。%%PerformNewtoniterationsfork=1:maxIter;Z=Z-(f(Z,d)./fprime(Z,d));endfunctiony=f(x,d);y=(x.^d)-1;endfunctiony=fprime(x,d);y=d*(x.^(d-1));end代码片段1:%%Finddrootsofunity,andthemaskforj=
