- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PAGE \* MERGEFORMAT16
本科生实验报告
实验课程 数值计算方法
学院名称 信息科学与技术学院
专业名称 计算机科学与技术
学生姓名
学生学号
指导教师
实验地点
实验成绩
二〇 一六 年 五 月 二〇一六 年 五 月
实验一 非线性方程求根
1.1问题描述
实验目的:掌握非线性方程求根的基本步骤及方法,。
实验内容:试分别用二分法、简单迭代法、Newton迭代法、弦截法(割线法、双点弦法),求 x5-3x3+x-1= 0 在区间 [-8,8]上的全部实根,误差限为10-6。
要求:讨论求解的全过程,对所用算法的局部收敛性,优缺点等作分析及比较,
第2章 算法思想
2.1二分法
思想:在函数的单调有根区间内,将有根区间不断的二分,寻找方程的解。
步骤: 1.取中点mid=(x0+x1)/2
2.若f(mid)=0,则mid为方程的根,否则比较与两端的符号,若与f(x0)
异号,则根在[x0,mid]之间,否则在[mid,x1]之间。
3并重复上述步骤,直达达到精度要求,则mid为方程的近似解。
2.2 简单迭代法
思想:迭代法是一种逐次逼近的方法,它是固定公式反复校正跟的近似值,使之逐步精确,最后得到精度要求的结果。
步骤:1.构造迭代公式f(x),迭代公式必须是收敛的。
2.计算x1,x1=f(x0).
3.判断|x1-x0|是否满足精度要求,如不满足则重复上述步骤。
4.输出x1,即为方程的近似解。
2.3 Newton迭代法
思想:设r是的根,选取作为r的初始近似值,过点做曲线的切线L,L的方程为,求出L与x轴交点的横坐标,称x1为r的一次近似值。过点做曲线
的切线,并求该切线与x轴交点的横坐标,称为r的二次近似值。重复以上过程,得r的近似值序列,其中,称为r的
次近似值
步骤:1.计算原函数的导数f’(x);构造牛顿迭代公式
2.计算 ,若f’(x0)=0,退出计算,否则继续向下迭代。
3.若|x1-x0|满足精度要求,x1即为方程的近似解。
2.4弦截法
思想:为加速收敛,改用两个端点都在变动的弦,用差商替代牛顿迭代公式的导数f’(x)。
步骤: 1.构造双??弦法的公式
2.计算x2=x1-f(x1)(x1-x0)/f(x1)-f(x0);
3.判断f(x2)是否满足精度要求,若没有则按照上述步骤继续迭代,否则输出x2.x2即为方程的近似解。
第3章 测试结果及分析
测试结果
函数图像
函数 Y=x5-3x3+x-1
二分法(表1-1,1-2,1-3)
[-1.6,-1.3]
kxkkxkkxk0-1.455-1.5015610-1.504931-1.5256-1.5039111-1.5052-1.48757-1.5050812-1.505043-1.506258-1.5044913-1.505064-1.496889-1.5047914-1.50507 表1-1
区间[-1.2,-0.9]
kxkkxkkxk0-1.055-0.99843710-1.000051-0.9756-1.0007811-0.9999762-1.01257-0.99960912-1.000013-0.993758-1.000213-0.9999944-1.003129-0.99990214-1表1-2
区间[1.5,1.8]
kxkkxkkxk01.6571.69102141.6902911.72581.69043151.6902921.687591.69014161.6902931.70625101.69028171.6902841.69687111.69036181.6902851.69219121.6903261.68984131.6903表1-
文档评论(0)