原创力文档- 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）