原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PAGE
PAGE 2
《MATLAB语言及应用》
大作业
姓名：
学号：
学院：
班级：
题目编号：
- PAGE 4 -
割线法求解非线性方程
割线法的数学理论
割线法原理：牛顿迭代法的收敛速度快，但是每迭代一次，除需计算飞f(x)的值外，还要计算f’(x)的值外。如果f(x)比较复杂，计算f(x)’的工作量就可能很大。为了避免计算导数值，我们用差商来代替导数。?
设经过k次迭代后，与求xk-1。用f(x)在xk,两点的差商来代替牛顿迭代公式中的导数值f’(xk),于是我们得到如下迭代公式：
二、割线法的算法和流程图
f’(x)=0 迭代次数大于N
f’(x)=0 迭代次数大于N
算法失败
输入：N, x0 ，ε η
开始
输出：x
输出：x1
三、割线法的Matlab实现
编辑m文件：Secant.m
function[pl,err,k,y]=secant(f,p0,p1,delta,max)
%f是给定的非线性函数
%p0,p1为初始值
%delta为给定误差界
%max1是迭代次数的上限
%p1为所求得的方程的近似解
%err为p1-p0的绝对值
%k为所需要的迭代次数
%y=f(pl)
k=0,p0,pl,feval(‘f’,p0),feval(‘f’,pl)
for k=1:max1
p2=p1-feval(‘f’,pl)*(pl-p0)/feval(‘f’,pl)-feval(‘f’,p0);
err=abs(p2-p1);
p0=p1;
p1=p2;
k,pl,err,y=feval(‘f’,pi);
if(err,delta)｜(y==0),break,end
end
四、割线法的算例实现
例: 非线性方程 y=x3+2x?0.205061208 .
建立gexianfa.m源程序，源程序代码为:
function?x=gexian(f,x0,x1,e)
if?nargin4,e=1e-4;end
y=x0;x=x1;
while?abs(x-y)e
z=x-(feval(f,x)*(x-y))/(feval(f,x)-feval(f,y));?????????? y=x;
x=z;
end
在matlab软件中执行下列语句
?clear
?fun=inline(x^3+2*x-0.205061208);
?gexianfa(fun,0,1,1e-6),format?short
n?=???
5?
ans?=?
???0.1020
文档评论（0）