数值分析上机题目详解解读.docx

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第一章 一、题目 设,其精确值为。 1) 编制按从大到小的顺序,计算SN的通用程序。 2) 编制按从小到大的顺序,计算SN的通用程序。 3) 按两种顺序分别计算,并指出有效位数。(编制程序时用单精度) 4) 通过本次上机题,你明白了什么? 二、通用程序 N=input(Please Input an N (N1):); AccurateValue=single((0-1/(N+1)-1/N+3/2)/2); Sn1=single(0); for a=2:N; Sn1=Sn1+1/(a^2-1); end Sn2=single(0); for a=2:N; Sn2=Sn2+1/((N-a+2)^2-1); end fprintf(The value of Sn (N=%d)\n,N); fprintf(Accurate Calculation %f\n,AccurateValue); fprintf(Caculate from large to small %f\n,Sn1); fprintf(Caculate from small to large %f\n,Sn2); disp(____________________________________________________) 三、结果 从结果可以看出有效位数是6位。 感想:可以得出,算法对误差的传播有一定的影响,在计算时选一种好的算法可以使结果更为精确。从以上的结果可以看到从大到小的顺序导致大数吃小数的现象,容易产生较大的误差,求和运算从小数到大数所得到的结果才比较准确。 第二章 一、题目 (1)给定初值及容许误差,编制Newton法解方程f(x)=0的通用程序。 (2)给定方程,易知其有三个根 a) 由Newton方法的局部收敛性可知存在当时,Newton迭代序列收敛于根x2*。试确定尽可能大的。 b)试取若干初始值,观察当时Newton序列的收敛性以及收敛于哪一个根。 (3)通过本上机题,你明白了什么? 二、通用程序 1、定义函数和导函数 %% 定义函数f(x) function Fx=fx(x) Fx=x^3/3-x; ---------------------------------- %% 定义导函数df(x) function Fx=dfx(x) Fx=x^2-1; 2、寻找最大的 clear flag=1; k=1; x0=0; while flag==1 delta=k*10^-6; x0=delta; k=k+1; m=0; flag1=1; while flag1==1 m=10^3 x1=x0-fx(x0)/dfx(x0); if abs(x1-x0)10^-6 flag1=0; end m=m+1; x0=x1; end if flag1==1||abs(x0)=10^-6 flag=0; end end fprintf(The maximun delta is %f\n,delta); 3、Newton法求方程的根 clear ef=10^-6; %% 给定容许误差10^-6 k=0; x0=input(Please input initial value Xo:); disp(k Xk); fprintf(0 %f\n,x0); flag=1; while flag==1 k=10^3 x1=x0-fx(x0)/dfx(x0); if abs(x1-x0)ef flag=0; end k=k+1; x0=x1; fprintf(%d %f\n,k,x0); end 三、结果 寻找最大的δ值: 在题目给出的不同区间内进行牛顿迭代: 综上所述:首先得到最大的δ值为0.774597。在各区间任意取初值进行迭代得到(-∞,-1)区间收敛于-1.73205,(-1,-δ)区间局部收敛于0, (-δ,δ)区间收敛于0,(δ,1)区间类似于(-1,δ)区间,收敛于0,(1,∞)收敛于1.73205。 感想:通过本上机题,了解到用Newton法求多根方程的根时,迭代序列收敛于某一个根有一定的区间限制。如果不清楚这个限制随意取值的话,会出现在一个区间上局部收敛于不同的根的情况。如下面的迭代: 第三章 一、

文档评论(0)

shuwkb + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档