《计算方法》课程论文.doc

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
《计算方法》课程论文《计算方法》课程论文《计算方法》课程论文

PAGE PAGE 4 《计算方法》课程论文 《计算方法》课程论文 牛顿迭代法的应用 班 级: 计算机班 学 号: 姓 名: 授课教师: 论文成绩: 20XX年X月 牛顿迭代法的应用 摘 要:针对小流域设计洪水计算中,采用推理公式求解洪峰流量Qm及汇流时间τ的方法繁琐,查算精度较低等问题,提出用牛顿迭代法求解,该方法较原图解试算法收敛快,迭代次数少,计算速度快,精度可靠。计算实例表明,用简单迭代法与牛顿迭代法计算设计洪峰流量,其结果与图解试算法完全一致。 关键词:推理公式;牛顿迭代法;简单迭代法;图解试算法 。 1 原理简介   迭代法是计算方法中的一种基本方法,是用来求解非线性方程(代数方程或超越方程)根的近似算法,迭代法的计算优劣取决于迭代序列的收敛性且依赖于迭代函数的构造.构造迭代函数的一条重要途径,是用近似方程来代替原方程. 1.1 简单迭代法   将方程f(x)=0化为一个等价(同解)的方程:x=φ(x),给定一个初值x0,代入右端可算得一个x1=φ(x0),再以x1代入右端,又可得x2=φ(x1)……如此继续下去,会得到一个序列{xk},其中: ??? xk+1=φ(xk)  k=0,1,2,…,n(1)     {xk}称为迭代序列,φ(x)称为迭代函数,式(1)称为迭代格式.如果迭代序列是收敛的,且收敛于X*,则当φ(x)连续时,必有:   这说明,只要迭代序列收敛,一般总收敛于原方程的解.实际计算当然不能做无穷多步,迭代到一定程度,就取xk+1作为原方程根的近似值.这种求根法称为简单迭代法,或称逐次逼近法. 1.2 牛顿迭代法 ??? 设xk是f(x)=0的一个近似根,把f(x)在xk处作泰勒展开:   显然是f(x)=0的同解方程,所以迭代函数为   故在α的邻域R内,对任意初值x0,由式(4)得到的迭代序列收敛于α.迭代式(4)所确定的方法称为牛顿迭代法.   牛顿迭代法有明显的几何意义.由式(4)知xk+1是点(xk,f(xk))处y=f(x)的切线:   与X轴交点的横坐标如图1所示.也就是说,新的近似值xk+1是用代替曲线y=f(x)的切线与x轴相交得到的.继续取点(xk+1,f(xk+1)),再作切线与X轴相交,又可得xk+2……由图1可见,只要初值取得充分靠近α,这个序列就会很快收敛于α.   由于牛顿迭代法的局部收敛性,又对初值要求较高,只有初值取得充分靠近α,才能保证序列收敛. 2 算 例   某河流流域面积F=142.1 km2,河长L=40.45km,河道坡度J=0.4%,汇流参数m=1.2,流域损失参数μ=2.5 mm/h,百年一遇最大24 h暴雨量436.7 mm,暴雨衰减指数n=0.75,要求用推理公式计算该河百年一遇的设计洪峰流量. 2.1 简单迭代法 将1428.6作为初值进行迭代,迭代计算结果如下:当Qm初值为1428.6,1330,1311,1307,1306时,Qm分别为1330,1311,1307,1306,1306. ??? 于是百年一遇设计洪峰流量为1306 m3/s,对应的τ值为 故经检验,产流模式符合原假定tc>τ. 2.2 牛顿迭代法   如前所述,牛顿迭代法只是在构造迭代函数方面与简单迭代法有所不同,在求解推理公式时,将式(7)进行变换可得下式:   假定初始值Qm(1)=1 428.6,迭代计算结果如下:当Qm初值为1428.6,1 307,1306时,Qm分别为1307,1306,1306.故Qm=1306 m3/s. 结束语 简单迭代法与牛顿迭代法的计算结果与图解试算法完全一样,且牛顿迭代法收敛较快,迭代次数少,计算速度快,精度可靠. 参考文献: [1]SDJ22-79,水利水电工程设计规范[S]. [2]徐德龙.推理公式计算方法的探讨[J].水文水资源,2000(2):12~13.

文档评论(0)

优秀文档 + 关注
实名认证
内容提供者

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

1亿VIP精品文档

相关文档