数值计算课堂提问含参考结果解析.docVIP

课堂提问及参考答案 第1章 误差与误差分析 1、 即 为避免相近数相减，从而丧失大量的有效数字，另一个根的计算可写成如下形式： 第2章 非线性方程的数值解法 1、证明1-x-sin(x)=0在区间[0,1]内有一个根，若使用二分法求误差不大于 0.5*10^（-4）的根要二分多少次?若取，能否用不动点迭代法求根？ 答：令f(x)= 1-x-sin(x)，显见f(x)为连续函数。f(x)=-1-cos(x),当x在区间[0,1]时，0cos(1)=cos(x)=1，则f(x)0，表明其单调递减。同时f(0)=10,f(1)=-sin(1)0，即f(0)*f(1)0，可得区间[0,1]内有且只有一个实根。 根据题意:a=0,b=1,ε=0.5*10^（-4） 由公式(2.4):可得 ，取k=14，即需要二分14次。 取迭代公式，则， 即当x在区间(0,1]时，， 根据定理2.4(p23)，该迭代公式在取迭代初值在(0,1]时，收敛，可求出根。 2、用不动点迭代法求根号131的数值解。 答：本题的重点是构造迭代公式ψ(x)： 令，为避免开方运算，变换为 ，即，但是注意到其导数远大于1，很可能发散。因此，可做如下变换： ， 即取，此时，因为根号131的值可以估计在区间(11,12)，故此其导数小于1，有可能收敛。 注意：此迭代公式对于根在（0,12）之间才可能有效，因此不具有一般性。实际上，使用加权加速或松弛加速均可以得到类似的迭代公式。（教材，公式2.19或2.23）以加权加速为例： ，由于根在在区间(11,12)，取 ，根据公式2.19，迭代公式为： 迭代的加速还有更好的解决办法， 比如：牛顿迭代法构造迭代公式 ，即取 %代码如下： A=131; fun1=@(x)x^2+x-A; %x=x^2+x-A fun2=@(x)-(x^2)/24+x+A/24; % 类似加权加速的思想 fun3=@(x)x/2+A/(2*x); %x=x/2+A/(2*x) 牛顿迭代 [x_star,k]=Iterate1(fun2,11) %取初值11，迭代公式fun2 [x_star,k]=Iterate1(fun3,11) %取初值11，迭代公式fun3 hold on,grid on fplot(fun2,[11,12]) fplot(fun3,[11,12],r) 3. 本章课后作业题1,利用迭代法证明?a)x=5/9*x+5*A/(9*x^2)-A*A/(9*x^5) b)x=A/(3*x^2)+2*x/3 哪种格式收敛较快？为什么？ 答：a) 令 已知，带入上式得 ，同样带入，得 ，带入，得 所以a)格式是3阶收敛的。 b) 令 已知，带入上式得 所以b)格式是2阶收敛的，b)格式其实是牛顿迭代法的迭代格式，牛顿迭代法对于单根是平方收敛的。 综上，a)格式收敛快。 下面以求为例，画出不同迭代格式对应的曲线图形，可以看到 A=125; fun1=@(x)x; %对角线 fun2=@(x)5/9*x+5*A/(9*x^2)-A*A/(9*x^5); %三阶收敛，红色线 fun3=@(x)A/(3*x^2)+2*x/3; %牛顿迭代，平方收敛，黑色 fun4=@(x)-x^3/48+x+A/48; %加权迭代，线性收敛，绿色 hold on,grid on fplot(fun1,[4.5,5.5]) fplot(fun2,[4.5,5.5],r) fplot(fun3,[4.5,5.5],k) fplot(fun4,[4.5,5.5],g) 5、以下不属于不动点迭代法的有(B、E) A、加权加速法； B、埃特金加速法 C、牛顿下山法； D、割线法 E、二分法 F、松弛法 G、牛顿平行弦法 H、斯蒂芬森加速法 I、单点割线法 6、设用牛顿法求方程在区间[a,b]内的单根，已知其收敛，则它具有（ B ）敛速。 A. 线性 B. 平方 C. 超线性 D. 三阶 7.设函数,写出解该方程的牛顿迭代格式，并求出此格式的收敛阶。 答：牛顿迭代格式： 所以，即 已知，则，所以对于重根，此时牛顿迭代法线性收敛。 8.若把上式中的导数改为，其几何含义是什么？ 答：导数用差商形式，即是单点割线法的迭代形式。 第3章 线性方程组的数值解法 1、线性方程组Ax=b，有解的条件是什么？ 答：R(A)=R(Ab)，即矩阵A和增广矩阵Ab秩相等时，该线性方程组有解；若R(A)=R(Ab)=n，即矩阵A和增广矩阵Ab满秩时，该线性方程组有唯一的解向量。 2、已知向