计算方法上机作业集合.docxVIP

精品word学习资料可编辑 名师归纳总结——欢迎下载 第一次 其次次上机作业 上机作业： 1.在 Matlab 上执行： 5.1-5-0.1 和 1.5-1-0.5 精品word学习资料可编辑 名师归纳总结——欢迎下载 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 给出执行结果，并简要分析一下产生现象的缘由；解：执行结果如下： 在 Matlab 中，小数值很难用二进制进展描述； 由于运算精度的影响，相近两数相减会显现误差； 2.〔课本 181 页第一题〕 解：〔1〕n=0 时，积分得 =ln6-ln5 ，编写如如下图代码 精品word学习资料可编辑 名师归纳总结——欢迎下载 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 从 以 上 代 码 显 示 的 结 果 可 以 看 出 ， 的 近 似 值 为 0.012712966517465 〔2〕 = dx,可得 dx dx dx,得 ,如此有 , 精品word学习资料可编辑 名师归纳总结——欢迎下载 取 ,以此逆序估算 ；代码段与结果如如下图所示 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 结 果 是 从 逆 序 输 出 至 ， 所 以 得 到 的 近 似 值 为 0.088392216030227 ； 〔3〕从 估量的过程更为牢靠；第一依据积分得表达式是可知，被积函数随着 n 的增大，其所围面积应当是逐步减小 的，即积分值应是随着 n 的递增二单调减小的， 〔1〕中输出的值不满意这一条件， 〔 2〕满意；设 表示 的近似值， - = ( 依据递推公式可以导出此式 ),可以看 精品word学习资料可编辑 名师归纳总结——欢迎下载 出，随着 n 的增大，误差也在增大，所以次序估量时，算法不稳固性逐步增大，逆序估量情形如此刚好相反，误差不断减小，算法逐步趋于稳固； 2.〔课本 181 页其次题〕 精品word学习资料可编辑 名师归纳总结——欢迎下载 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 〔1〕上机代码如以下图 求得近似根为 0.09058 〔2〕上机代码如以下图 精品word学习资料可编辑 名师归纳总结——欢迎下载 得近似根为 0.09064 得近似根为 0.09064； 〔3〕牛顿法上机代码如下 |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 运算所得近似解为 0.09091 第三次上机作业 精品word学习资料可编辑 名师归纳总结——欢迎下载 上机作业 上机作业 181 页第四题 线性方程组为 = 精品word学习资料可编辑 名师归纳总结——欢迎下载 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载.  ） 次序消元法 A=[1.1348,3.8326,1.1651,3.4017;0.5301,1.7875,2.5330,1.5435;3.4129,4. 9317,8.7643,1.3142;1.2371,4.9998,10.6721,0.0147]; b=[9.5342;6.3941;18.4231;16.9237]; 上机代码〔函数局部〕如下 function [b] = gaus( A,b ) %用b 返回方程组的解B=[A,b]; n=length(b); RA=rank(A); RB=rank(B); dif=RB-RA; if dif0 disp( 此方程组无解 ); return end if RA==RB 精品word学习资料可编辑 名师归纳总结——欢迎下载 if RA==n format long; disp( 此方程组有唯独解 ); for p=1:n-1 for k=p+1:n 精品word学习资料可编辑 名师归纳总结——欢迎下载 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. m=B(k,p)/B(p,p); B(k,p:n+1)=B(k,p:n+1)-m*B(p,p:n+1); end end %次序消元形成上三角矩阵b=B(1:n,n+1); A=B(1:n,1:n); b(n)=b(n)/A(n,n); for q=n-1:-1: