- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析实验2new.doc
《数值分析》课程实验报告
学院:数学学院
专业:信息与计算科学
班级:10信计
学号:
姓名
实验名称 用牛顿迭代法求解非线性方程 成绩
一、实验目的
1.了解求解非线性方程的解的常见方法。
2.编写牛顿迭代法程序求解非线性方程。
二、实验内容
分别用初值,和求,要求精度为。
三、算法描述
牛顿法实质是一种线性化方法,其基本思想是将非线性方程逐步归结为某种线性方程来求解。
对于牛顿法来说它的一般步骤为:
选定初始近似值x0,计算f0=f(x0),f0=f(x0)
按公式想x(k+1=)=x(k)-f(xk)/f(xk)进行迭代下去
使得所求的根在精度的范围内
如果迭代次数达到预先指定的次数或f1=0,则方法不可以进行下去,否则依次进行迭代下去,直到结果出现为止。
四、实验步骤
1.根据实验题目,给出题目的程序。
上机输入和调试自己所编的程序。
function [k,xk,yk,piancha,xdpiancha]=newtonqx(x0,tol,ftol,gxmax)
x(1)=x0;
for i=1: gxmax
x(i+1)=x(i)-fnq(x(i))/(dfnq(x(i))+eps); piancha=abs(x(i+1)-x(i));
xdpiancha= piancha/( abs(x(i+1))+eps); i=i+1;
xk=x(i);yk=fnq(x(i)); [(i-1) xk yk piancha xdpiancha]
if (abs(yk)ftol)((pianchatol)|(xdpiancha tol))
k=i-1; xk=x(i);[(i-1) xk yk piancha xdpiancha]
return;
end
end
if igxmax
disp(??×¢òa£oμü′ú′?êy3?1y???¨μ?×?′ó?μgxmax?£)
k=i-1; xk=x(i);[(i-1) xk yk piancha xdpiancha]
return;
end
[(i-1),xk,yk,piancha,xdpiancha];
M文件:
1.function y=fnq(x)
y=x^2-113;
2. function y=dfnq(x)
y=2*x;
运行结果:
1.初值为0.01:
k
xk
Yk
piancha
xdpiancha
1
0.0006
3.19
0.0006
0
2
0.0028
7.98
0.0028
0
3
0.0014
2
0.0014
0
4
0.0071
4.99
0.0071
0
5
0.0035
1.25
0.0035
0
6
0.0177
3.11
0.0176
0.0001
7
0.0887
7.76
0.0881
0.001
8
0.045
1.91
0.0437
0.001
9
23.751
451
21.2395
0.8943
10
14.254
90.2
9.4967
0.6662
11
11.091
10
3.1635
0.2852
12
10.64
0.2
0.4512
0.0424
13
10.63
0
0.0096
0.0009
14
10.63
0
0
0
2.初值为10时
k
xk
yk
piancha
xdpiancha
1
10.65
0.42
0.65
0.061
2
10.63
0
0.0198
0.0019
3
10.63
0
0
0
4
106301
0
0
0
3.初值为300时:
k
xk
yk
piancha
xdpiancha
1
0.015
2.2444
0.015
0.0001
2
0.0755
5.5828
0.0747
0.001
3
0.0385
1.368
0.037
0.001
4
20.7101
315.9065
17.7738
0.8582
5
13.0832
58.1694
7.6269
0.583
6
10.8601
4.942
2.2231
0.2047
7
10.6326
0.0518
0.2275
0.0214
8
10.6301
0
0.0024
0.0002
9
10.6301
0
0
0
2)开方方法:
function [k,xk,yk,piancha,xdpi
文档评论(0)