原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数值分析上机实验报告
姓名：
班级：
学号：
PAGE \* MERGEFORMAT
PAGE \* MERGEFORMAT 23
第一次上机：
1.题目：
分别用不动点法和牛顿法求解方程
x-exp(x)+4 = 0
Fixed point converge to negative root
The principal：
a number p is fixed point for a given function g if g(p) = p .And we can change The equation to fixed point form .as for fixed point if pk make g(pk) = pk .
We can change the form to :
g1 = exp(g0) -4
use the initial number as 1
The code:
g0 = 1;
g1 = exp(g0) -4;
while abs(g1 - g0)0.00001
g0 = g1;
g1 = exp(g0) -4;
end
g1
the results :
w1_fixed_point
g1 =
-3.9813
结果分析
通过控制迭代式的形式可以决定收敛到正跟或负根。
Fixed point converge to positive root
The principal：
We can change a new fixed point form to produce a positive root
We can change the form to :
g1 = log(4+g0);
use initial number as -3.9
The code:
g0 = -3.9;
g1 = log(4+g0);
while abs(g1 - g0)0.00001
g0 = g1;
g1 = log(g0+4);
end
g1
the results :
w1_fixed_point_positive
g1 =
1.7490
Newton’s method to find negative roots
The principal：
We derive Newton’s method ,based on taylor polynomials .newton’s method is a kind of fixed point method which has two order convergence.
We can change the form to :
g1 = g0-(g0-exp(g0)+4)/(1-exp(g0) )
use initial nuber as -3
The code:
g0 = -3;
g1 = g0-(g0-exp(g0)+4)/(1-exp(g0));
while abs(g1 - g0)0.00001
g0 = g1;
g1 = g0-(g0-exp(g0)+4)/(1-exp(g0))
end
the results :
w1_newton
g1 =
-3.981342639636226
g1 =
-3.981339370911418
结果分析：
收敛速度明显快于不动点（结果太长没打出）。
Newton’s method to find positive roots
The principal：
Use initial number as 1
The code:
g0 = 1;
g1 = g0-(g0-exp(g0)+4)/(1-exp(g0));
while abs(g1 - g0)0.00001
g0 = g1;
g1 = g0-(g0-exp(g0)+4)/(1-exp(g0));
end
g1
the results :
w1_newton
g1 =
1.749031386012702
Multiple roots
结果分析：
通过确定初始值得正负来控制收敛到哪个根。
2.p100#9题
题目：
Use each of the following methods to find a solution in[0.1,1] accurate to within 10^-4 for
600*x^4 -550*x^3+200*x^2-20*x-1=0
a.bisection b.newton’s method
c.secant method d.meth
文档评论（0）