一．作业：②
用弦线法或牛顿法求方程lg 1x x ?=的实根。求解区间[]2,3x ∈。
二．程序:
1.弦线法：
#include
#include
using namespace std;
double f(double x){
return x*log10(x)-1;
}
int main(){
double x1(2),x2(3),x3,s;
cout<<"Please input the acceptable error: ";
cin>>s;
x3=x1-(x2-x1)*f(x1)/(f(x2)-f(x1));
while(abs(f(x3))>s){
if(f(x3)*f(x1)<0)
x2=x3;
else
x1=x3;
x3=x1-(x2-x1)*f(x1)/(f(x2)-f(x1));
}
cout<<"the root of the equation between 2 and 3 is : "< return 0; } 运行结果： 2.牛顿法 #include #include using namespace std; double f(double x){ return x*log10(x)-1; } double f1(double x){ return log10(x)+1/log(10); } int main(){ double x1(2),x2(3),x3,s; cout<<"Please input the acceptable error: "; cin>>s; x3=x1-f(x1)/f1(x1); while(abs(f(x3))>s&&x3<=x2){ x1=x3; x3=x1-f(x1)/f1(x1); } if(x3==x2&&abs(f(x2))>s){ cout<<"There is no root in [2,3]."< return 0; } cout<<"The root of the equation is: "< return 0; } 运行结果： 《MATLAB 程序设计实践》课程考核 1、 编程实现以下科学计算方法，并举一例应用之（参考书籍《精通MATLAB 科学计算》，王正林著，电子工业出版社，2009年）“弦截法非线性方程求解” 算法说明： （1） 过两点))(,(a f a ,))(,(b f b 作一直线，它与x 轴有一个交点，记为1x ； （2） 如果f(a)f(x1)<0,过两点))(,(a f a ，))(,(11x f x 作一直线，它与x 轴的交点记为2x ， 否则过两点))(,(b f b ，))(,(11x f x 作一直线，它与x 轴的交点记为2x ； （3） 如此下去，直到 ε<--1n n x x 就可以认为n x 为0)(=x f 在区间],[b a 上的一个根。 （4） k x 的递推公式为： ??? ???? >---=<---=------0)()(),()()(0)()(),()()(111111k k k k k k k k x f a f b f b f x f b x b x x f a f a f a f x f a x a x 且) ()(1a f b f a b a x --- =。 在MATLAB 中编程实现的弦截法的函数为：Secant. 功能：用弦截法求函数在某个区间的一个零点。 调用格式：root=Secant(f,a,b,eps). 其中，f 为函数名； a 为区间左端点； b 为区间左端点； eps 为根的精度； root 为求出的函数零点。 流程图： 源代码： function root=Secant(f,a,b,eps) %弦截法求函数在区间[a,b]上的一个零点 %函数名：f %区间左端点：a %区间右端点：b %根的精度：eps %求出函数的函数零点：root if (nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if (f1==0) 学科前沿讲座论文 班级：工程力学13-1班姓名：陆树飞 学号：02130827 牛顿法求非线性方程的根 一 实验目的 （1）用牛顿迭代法求解方程的根 （2）了解迭代法的原理，了解迭代速度跟什么有关 题目：用Newton 法计算下列方程 （1） 013=--x x ， 初值分别为10=x ，7.00=x ，5.00=x ； （2） 32943892940x x x +-+= 其三个根分别为1,3,98-。当选择初值02x =时 给出结果并分析现象,当6510ε-=?，迭代停止。 二 数学原理 对于方程f(x)=0，如果f(x)是线性函数，则它的求根是很容易的。牛顿迭代法实质上是一种线性化方法，其基本思想是将非线性方程f(x)=0逐步归结为某种线性方程来求解。 设已知方程f(x)=0有近似根x k (假定k f'(x )0≠) ，将函数f(x)在点x k 进行泰勒展开，有 k k k f(x)f(x )+f'(x )(x-x )+≈??? 于是方程f(x)=0可近似的表示为 k k k f(x )+f'(x )(x-x )=0 这是个线性方程，记其根为x k+1，则x k+1的计算公式为 k+1k ()x =x -'() k k f x f x ，k=0,1,2,… 这就是牛顿迭代法。 三 程序设计 （1）对于310x x --=，按照上述数学原理，编制的程序如下 program newton implicit none real :: x(0:50),fx(0:50),f1x(0:50)!分别为自变量x ，函数f(x)和一阶导数f1(x) integer :: k write(*,*) "x(0)=" read(*,*) x(0) !输入变量：初始值x(0) open(10,file='1.txt') do k=1,50,1 fx(k)=x(k-1)**3-x(k-1)-1 f1x(k)=3*x(k-1)**2-1 x(k)=x(k-1)-fx(k)/f1x(k) !牛顿法 write(*,'(I3,1x,f11.6)') k,x(k) !输出变量：迭代次数k 及x 的值 write(10,'(I3,1x,f11.6)') k,x(k) if(abs(x(k)-x(k-1))<1e-6) exit !终止迭代条件 end do stop end （2）对于32943892940x x x +-+=，按照上述数学原理，编制的程序如下 program newton implicit none 非线性方程数值解法及其应用 摘要:数值计算方法主要研究如何运用计算机去获得数学问题的数值解的理论和算法。 本文主要介绍非线性方程的数值解法以及它在各个领域的应用。是直接从方程出发，逐步缩小根的存在区间，或逐步将根的近似值精确化，直到满足问题对精度的要求。我将从二分法、Steffensen 加速收敛法、Newton 迭代法、弦截法来分析非线性方程的解法及应用。 关键字：非线性方程;二分法；Steffensen 加速收敛法；代数Newton 法；弦截法 一、前言 随着科技技术的飞速发展，科学计算越来越显示出其重要性。科学计算的应用之广已遍及各行各业，例如气象资料的分析图像，飞机、汽车及轮船的外形设计，高科技研究等都离不开科学计算。因此经常需要求非线性方程 f(x) = O 的根。方程f(x) = O 的根叫做函数f(x)的零点。由连续函数的特性知：若f(x)在闭区间[a ，b]上连续，且f(a)·f(b) 《MATLAB 程序设计实践》课程考核 ---第37-38页 题１ ： 编程实现以下科学计算算法，并举一例应用之。（参考书籍《精 通ＭＡT ＬＡＢ科学计算》，王正林等著，电子工业出版社，２００９ 年） “牛顿法非线性方程求解” 弦截法本质是一种割线法，它从两端向中间逐渐逼近方程的根；牛顿法本质上是一种切线法，它从一端向一个方向逼近方程的根，其递推公式为： - =+n n x x 1) ()(' n n x f x f 初始值可以取)('a f 和)('b f 的较大者，这样可以加快收敛速度。 和牛顿法有关的还有简化牛顿法和牛顿下山法。 在MATLAB 中编程实现的牛顿法的函数为：NewtonRoot 。 功能：用牛顿法求函数在某个区间上的一个零点。 调用格式：root=NewtonRoot )(```eps b a f 其中，f 为函数名； a 为区间左端点； b 为区间右端点 eps 为根的精度； root 为求出的函数零点。 ， 牛顿法的matlab程序代码如下： function root=NewtonRoot(f,a,b,eps) %牛顿法求函数f在区间[a,b]上的一个零点%函数名:f %区间左端点:a %区间右端点:b %根的精度:eps %求出的函数零点:root if(nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if (f1==0) root=a; end if (f2==0) root=b; end if (f1*f2>0) disp('两端点函数值乘积大于0 ！'); return; else tol=1; fun=diff(sym(f)); %求导数 fa=subs(sym(f),findsym(sym(f)),a); fb=subs(sym(f),findsym(sym(f)),b); dfa=subs(sym(fun),findsym(sym(fun)),a); dfb=subs(sym(fun),findsym(sym(fun)),b); if(dfa>dfb) %初始值取两端点导数较大者 root=a-fa/dfa; else root=b-fb/dfb; end while(tol>eps) r1=root; fx=subs(sym(f),findsym(sym(f)),r1); dfx=subs(sym(fun),findsym(sym(fun)),r1); %求该点的导数值 root=r1-fx/dfx; %迭代的核心公式 tol=abs(root-r1); end end 例：求方程3x^2－exp（x）＝0的一根 解：在MATLAB命令窗口输入： >> r=NewtonRoot('3*x^2-exp(x)',3,4) 输出结果： X=3.7331 中南大学MATLAB程序设计实践学长有爱奉献，下载填上信息即可上交，没有下载券的自行百度。所需m文件照本文档做即可，即新建（FILE）→脚本（NEW-Sscript）→复制本文档代码→运行（会跳出保存界面，文件名默认不要修改，保存）→结果。第一题需要把数据文本文档和m文件放在一起。全部测试无误，放心使用。本文档针对做牛顿法求非线性函数题目的同学，当然第一题都一样，所有人都可以用。←记得删掉这段话 班级： ？ 学号： 姓名： 一、《MATLAB程序设计实践》Matlab基础 表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散 空间函数值来表示取向分布函数，是三维取向分布函数的一个实例。 由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一 个matlab程序画出如下的几种图形来分析其取向分布特征： （1）用Slice函数给出其整体分布特征； " ~ （2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况(需要用到subplot函数)； (3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。 ( 备注：数据格式说明 解： （1）（ （2）将文件内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如 下： fid=fopen(''); for i=1:18 tline=fgetl(fid); end phi1=1;phi=1;phi2=1;line=0; f=zeros(19,19,19); [ while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1;数据说明部分，与 作图无关此方向表示f随着 φ1从0,5,10,15, 20 …到90的变化而 变化 此方向表示f随着φ 从0,5,10,15, 20 … 到90的变化而变化 表示以下数据为φ2＝0的数据，即f（φ1，φ，0） (1)二分法求解非线性方程： #include 第七章 非线性方程求根 教学目的与要求： 理解二分法求根的思想；掌握二分法求解过程；了解二分法的优点和缺点。了解迭代法的基本思想，迭代法的收敛条件以及局部收敛性的定义；理解基本迭代法的迭代思路，收敛条件的产生与求证过程；掌握基本迭代法的迭代格式，收敛条件的应用以及局部收敛定理。 重点和难点：迭代法的基本思想，迭代法的收敛性 ■ 教学内容： 基本概念： 的零点； 的m 重零点。 )(x f )(x f 非线性方程的求根通常分为两个步骤：一是对根的搜索，二是根的精确化，求得根的足够精确的近似值。 求方程的有根区间有如下方法： (1)描图法。画出的简图，从曲线与)(x f y =x 轴交点的位置确定有根区间。 (2)解析法。根据函数的连续性、介值定理以及单调性等寻找有根区间。 § 1 二分法 分析二分法的基本原理 例1 用二分法求方程的一个正根，要求误差不超过. 01)(6=??=x x x f 2105.0?ק 2 迭代法及其收敛性 一、迭代法的定义 二、基本迭代法 定义：将方程改写成以下等价形式() x x ?=取定初始值0x ，由迭代公式1() (0,1,2,)n n x x n ?+==L 产生迭代序列{}n x 。显然，若{}n x 收敛于*x ,()x ?在*x 处连续，就有** 1lim lim ()()n n n n x x x ??+→∞→∞ ===x 即*x 是方程() x x ?=的解，从而也是0)(=x f 的解。故当充分大时，可取作为方程根的近似值。用迭代格式求得方程近似根的方法称为基本迭代法，n n x )(x ?称为迭代函数。由于收敛点*x 满足*()* x x ?=，故称*x 为)(x ?的不动点 例 求方程的一个实根，要求精确到六位小数。 032)(3 =??=x x x f 注意：把此方程转换成三种等价形式 ,32)(31+==x x x ?)3(2 1)(32?= =x x x ?, 3)(33??==x x x x ?三、迭代法的收敛条件 牛顿迭代法求解非线性方程组 非线性方程组如下： 221122121210801080 x x x x x x x ?-++=??+-+=?? 给定初值()00.0T x =，要求求解精度达到0.00001 1.首先建立函数()F X ,方程编程如下，将F.m 保存到工作路径中： function f=F(x) f(1)=x(1)^2-10*x(1)+x(2)^2+8; f(2)=x(1)*x(2)^2+x(1)-10*x(2)+8; f=[f(1),f(2)] ; 2.建立函数()DF X ,用于求方程的jacobi 矩阵,将DF.m 保存到工作路径中： function df=DF(x) df=[2*x(1)-10,2*x(2);x(2)^2+1,2*x(1)*x(2)-10]; %jacobi 矩阵是一阶偏导数以一定方式排列成的矩阵。 3.编程牛顿迭代法解非线性方程组，将newton.m 保存在工作路径中： clear,clc; x=[0,0]'; f=F(x); df=DF(x); fprintf('%d %.7f %.7f\n',0,x(1),x(2)); N=4; for i=1:N y=df\f'; x=x-y; f=F(x); df=DF(x); fprintf('%d %.7f %.7f\n',i,x(1),x(2)); if norm(y)<0.0000001 break; else end end ezplot('x^2-10*x+y^2+8',[-6,6,-6,6]); hold on ezplot('x*y^2+x-10*y+8',[-6,6,-6,6]); 运行结果如下： 0 0.0000000 0.0000000 1 0.8000000 0.8800000 2 0.9917872 0.9917117 一.实验目的： 通过本节实验课的学习，要求我们理解并掌握二分法、不动点迭代、牛顿切线法及弦截法解非线性方程求根的原理，掌握相应的算法原理，通过计算机解决实验问题 二.实验内容： 1、用对分区间法方程 1-x-sinx=0 在区间[0,1]上的误差小于10^(-4)的一个根，并记录对方区间的次数。 2、用不动点迭代法求解方程下 x-log(x)=2(x>1) 要求相对误差容限e=10^(-8)。 3、用Newton法求方程 x^3-x-1=0 在区间[-3,3]上的误差不大于10^(-5)的根，分别取初值x0=1.5, x0=0, x0=-1进行计算，比较他们的迭代次数。 三. 实验方案（程序设计说明）[包括算法设计思路，必要的流程图，界面设计 说明、使用模块及变量的说明等。] 1、二分法是对区间收索法的一种改进，具体做法为：先求一区间的中点，并计算其 函数值，若恰好有函数值为0，就是方程的根，若不为0，在判断此点的函数值与两端的函数值乘积的情况，取小于0的那个端点在进行上述对分，直到满足要求为止。 2、迭代法分为两种，一种是从任何可取的初值出发都能保证收敛，称之为大范围收 敛的方法。另一类称之为局部收敛法，即为了保证收敛必须选取初值充分接近于所要求的解。迭代法的基本思想是一种逐渐逼近的方法，首先给定一个粗造的初值，然后用一个迭代公式，反复矫正这个初值，直到满足预先给出的精确要求为止。 3、双点弦接法与Newton法不同，两者有本质的区别，它分为两步，不属于不动点 迭代法。 四. 实验步骤或程序（经调试后正确的源程序） （填写主要步骤与程序代码等，不够可附页） 1、f=inline('x+sin(x)-1'); a=0; b=1; dlt=1.0e-4; k=1; while abs(b-a)>dlt c=(a+b)/2; if f(c)==0 break; elseif f(c)*f(b)<0 a=c; C++实现牛顿迭代解非线性方程组（二元二次为例） 求解{0=x*x-2*x-y+0.5; 0=x*x+4*y*y-4; }的方程 #include 2012-2013（1）专业课程实践论文弦截法求非线性方程的根 姓名：秦基超，0818180121，R数学08-1班 一、算法理论 弦截法是在牛顿法的基础上得出的求解非线性方程的一种十分重要的插值方法。用牛顿方法求解非线性方程的根时，每一步都要计算导数值，如果函数f(x)比较复杂时，计算f(x)的导数往往比较困难。而弦截法使用差商代替牛顿法中的导数值进行迭代，避免了计算函数的导数值，并且收敛速度很快。 设k x ,1-k x 是)(x f =0的近似根，利用)(k x f ,)(1-k x f 构造一次插值多项式)(1x p ,并用)(1x p =0的根作为)(x f =0的新的近似根1+k x 。 由于)(1x p =)(k x f +)()()(1 1k k k k k x x x x x f x f ----- 因此有1+k x =k x -)() ()()(11----k k k k k x x x f x f x f 上述公式可以看做牛顿公式1+k x =k x - ) ()(k k x f x f '中导数)(k x f '用差商11)()(----k k k k x x x f x f 取代的结果。 由上面1+k x 式中可以得到满足公式的点的序列，且点的序列收敛于根。 弦截法的算法步骤： （1）确定初始值以及精度要求。 （2）计算函数值。 （3）利用公式计算方程新的近似根。 （4）若满足精度时，便可停止迭代，此时的x 作为方程的近似根，计算结果。否则以),(11f x ，),(22f x 分别代替),(00f x ，),(11f x 然后继续迭代。 牛顿迭代法c++程序设计 求解{0=x*x-2*x-y+0.5; 0=x*x+4*y*y-4; }的方程 #include 基于Matlab 实现牛顿迭代法解非线性方程组 已知非线性方程组如下 2211221212 10801080x x x x x x x ?-++=??+-+=?? 给定初值0(0,0)T x =，要求求解精度达到0.00001 首先建立函数F(x)，方程组编程如下，将F.m 保存到工作路径中： function f=F(x) f(1)=x(1)^2-10*x(1)+x(2)^2+8; f(2)=x(1)*x(2)^2+x(1)-10*x(2)+8; f=[f(1) f(2)]; 建立函数DF(x),用于求方程组的Jacobi 矩阵，将DF.m 保存到工作路径中： function df=DF(x) df=[2*x(1)-10,2*x(2);x(2)^2+1,2*x(1)*x(2)-10]; 编程牛顿迭代法解非线性方程组，将newton.m 保存到工作路径中： clear; clc x=[0,0]'; f=F(x); df=DF(x); fprintf('%d %.7f %.7f\n',0,x(1),x(2)); N=4; for i=1:N y=df\f'; x=x-y; f=F(x); df=DF(x); fprintf('%d %.7f %.7f\n',i,x(1),x(2)); if norm(y)<0.0000001 break ; else end end 运行结果如下： 0 0.0000000 0.0000000 1 0.8000000 0.8800000 2 0.9917872 0.9917117 3 0.9999752 0.9999685 4 1.0000000 1.0000000 Newton迭代法求解非 线性方程 一、 Newton 迭代法概述 构造迭代函数的一条重要途径是用近似方程来代替原方程去求根。因此，如果能将非线性方程f (x )=0用线性方程去代替，那么，求近似根问题就很容易解决，而且十分方便。牛顿(Newton)法就是一种将非线性方程线化的一种方法。 设k x 是方程f (x )=0的一个近似根，把如果)(x f 在k x 处作一阶Taylor 展开，即： )x x )(x ('f )x (f )x (f k k k -+≈ (1-1) 于是我们得到如下近似方程： 0)x x )(x ('f )x (f k k k =-+ (1-2) 设0)('≠k x f ，则方程的解为： x ?=x k +f (x k ) f (x k )? (1-3) 取x ~作为原方程的新近似根1+k x ，即令： ) x ('f ) x (f x x k k k 1k -=+, k=0,1,2,… (1-4) 上式称为牛顿迭代格式。用牛顿迭代格式求方程的根的方法就称为牛顿迭代法，简称牛顿法。 牛顿法具有明显的几何意义。方程： )x x )(x ('f )x (f y k k k -+= (1-5) 是曲线)x (f y =上点))x (f ,x (k k 处的切线方程。迭代格式(1-4)就是用切线式(1-5)的零点来代替曲线的零点。正因为如此，牛顿法也称为切线法。 牛顿迭代法对单根至少是二阶局部收敛的，而对于重根是一阶局部收敛的。一般来说，牛顿法对初值0x 的要求较高，初值足够靠近*x 时才能保证收敛。若 要保证初值在较大范围内收敛，则需对)x (f 加一些条件。如果所加的条件不满足，而导致牛顿法不收敛时，则需对牛顿法作一些改时，即可以采用下面的迭代格式： ) x ('f ) x (f x x k k k 1k λ -=+, ?=,2,1,0k (1-6) 上式中，10<λ<，称为下山因子。因此，用这种方法求方程的根，也称为牛顿下山法。 牛顿法对单根收敛速度快，但每迭代一次，除需计算)x (f k 之外，还要计算 )x ('f k 的值。如果)x (f 比较复杂，计算)x ('f k 的工作量就可能比较大。为了避免计算导数值，我们可用差商来代替导数。通常用如下几种方法： 1. 割线法 如果用 1 k k 1k k x x ) x (f )x (f ----代替)x ('f k ，则得到割线法的迭代格式为： )x (f ) x (f )x (f x x x x k 1k k 1 k k k 1k --+---= (1-7) 2. 拟牛顿法 如果用 ) x (f )) x (f x (f )x (f k 1k k k ---代替)x ('f k ，则得到拟牛顿法的迭代格式为： )) x (f x (f )x (f ) x (f x x 1k k k k 2k 1k -+--- = (1-8) 3. Steffenson 法 如果用 ) x (f ) x (f ))x (f x (f k k k k -+代替)x ('f k ，则得到拟牛顿法的迭代格式为： ) x (f ))x (f x (f ) x (f x x k k k k 2k 1 k -+- =+ Newton法解非线性方程组 一．题目重述：编程实现非线性方程组的牛顿解法，并求解如下方程组。 3x1?cos x2x3?0.5=0 x12?81x2+0.12+sin x3+1.06=0 e?x1x2+20x3+10π?3 3 =0 二．算法： 非线性方程组的牛顿法为：给定初始解向量x(0),对于k≥1生成 x(k)=x(k?1)?J x k?1?1F(x(k?1)). 三．编程实现： 这里用MATLAB程序实现，建立三个文件如下： 1.函数F(X)文件 function F =F( X) F(1,1)=3*X(1)-cos(X(2)*X(3))-0.5; F(2,1)=X(1)^2-81*(X(2)+0.1)^2+sin(X(3))+1.06; F(3,1)=exp(-X(1)*X(2))+20*X(3)+(10*pi-3)/3; end 2.J(X) 函数（即Jacobian矩阵）文件 function F1= F1(X ) F1(1,:)=[3,sin(X(1)*X(2))*X(3),sin(X(1)*X(2))*X(2)]; F1(2,:)=[2*X(1),-162*(X(2)+0.1),cos(X(3))]; F1(3,:)=[exp(-X(1)*X(2))*(-X(2)),exp(-X(1)*X(2))*(-X(1)),20]; end 3.解题脚本文件 文件名zu %% 牛顿法解非线性方程组 clear; X0=[0.1;0.1;-0.1]; for i=1:200 X=X0-F1(X0)\F(X0); %这里采用MATLAB的左除方法，避免算逆矩阵X0=X; end X Chap3 非线性方程求解的相关考题 一、判断题 1.（2008）（2011）若()()0f a f b <，则方程()0f x =在区间(,)a b 内至少有一个根 ( ) 2.（2009）设*x 是方程0)(=x f 的根，则求*x 的Newton 迭代法至少是平方收敛的。 ( ) 3.（2013）设A 是n 阶非奇异方阵，则解方程组A =x b 的迭代法收敛的充要条件是A 的谱半径()1A ρ<。 ( ) 4.（2016）设*x 是函数()x ?的不动点，且()1x ?'<，则由迭代格式1() k k x x ?-=(1,2,)k =L 产生的序列{}k x 收敛到* x 。 ( ) 二、填空题 1.（2008）（2010A ）（2010B ）设*x 是方程0)(=x f 的3重实根，则改进的Newton 迭代公式是 . 2.（2012）设*3x =是方程327159x x x -+=的2重实根，则求*x 的改进的Newton 迭代公式为 . 三、解答题 1、（2008）(本题满分14分) 已知方程3210x x --=在区间]2,1[内有一个实根 *x . (1)建立一个收敛的迭代公式求*x （不要用Newton 迭代法），并证明其收敛性。(4分) (2)取初值0 1.4x =，用Newton 迭代法求*x （只迭代两次）。(5分) (3)取初值011.4, 1.5x x ==，用弦截法求*x （只迭代两次）。(5分) 2、（2009）（2010A ）（2010B ）（2012）(本题满分10分) 已知方程10x xe -=在 00.5x =附近有一个实根*x . 实验名称： 牛顿法解非线性方程组 1 引言 我们已经知道，线性方程组我们可以采取Jacobi 迭代法，G-S 迭代法以及SOR 迭代方法求解。而在科学技术领域里常常提出求解非线性方程组的问题，例如，用非线性函数拟合实验数据问题、非线性网络问题，用差分法求解非线性微分方程问题等。 我们在解非线性方程组时，也考虑用迭代法求解，其思路和解非线性方程式一样，首先要将F(x)=0转化为等价的方程组 12(,,,),(1,2, )i i n x g x x x i n == 或者简记为x =g （x ），其中:,:n n n i g R R g R R →→ 112 2()()(),()n n n g x g x g R g x ???? ????????==∈???? ???????????? x x x x x 迭代法：首先从某个初始向量(0)x 开始，按下述逐次代入方法构造一向量序列(){}k x ： (1)()() 1(,,),(1,2,,)k k k i i n x g x x i n +== 其中，()()() ()12 (,,,)k k k k T n x x x =x 。 或写成向量形式：(1)()(),(0,1,2,)k k g k +==x x 如果()*lim k k →∞ ≡x x （存在），称(){}k x 为收敛。且当()i g x 为连续函数时，可得 *()*(lim )()k k g g →∞ ==x x x 说明*x 为方程组的解。又称为x =g （x ）的不动点。 本实验中采用牛顿迭代法来求解非线性方程组。 2 实验目的和要求 运用matlab 编写一个.m 文件，要求用牛顿法非线性方程组： 12(0)(1)()3211 cos 02,(取(0,0),要求10)1sin 0 2 T k k x x x x x x x +-∞ ?-=??=- 中南大学 MATLAB程序设计实践学长有爱奉献，下载填上信息即可上交，没有下载券 的自行百度。所需m文件照本文档做即可，即新建（FILE）→脚本（NEW-Sscript）→复制本文档代码→运行（会跳出 保存界面，文件名默认不要修改，保存）→结果。第 一题需要把数据文本文档和m文件放在一起。全部测 试无误，放心使用。本文档针对做牛顿法求非线性函 数题目的同学，当然第一题都一样，所有人都可以用。 ←记得删掉这段话 班级： 学号： 姓名： 一、《MATLAB程序设计实践》Matlab基础 表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散空间函数值来表示取向分布函数，是三维取向分布函数的一个实例。由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一个matlab程序画出如下的几种图形来分析其取向分布特征：（1）用Slice函数给出其整体分布特征； （2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况(需要用到subplot函数)； (3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。 备注：数据格式说明 解： （1）将文件内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如下： fid=fopen(''); for i=1:18 tline=fgetl(fid); end phi1=1;phi=1;phi2=1;line=0; f=zeros(19,19,19); while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1; if mod(line,20)==1 phi2=(data/5)+1; phi=1; 数据说明部分，与作图无关 此方向表示f 随着φ1从0,5,10,15, 20 …到90的变化而变化 此方向表示f 随着φ从0,5,10,15, 20 …到90的变化而变化 表示以下数据为φ2＝0的数据，即f （φ1，φ，0） matlab实现牛顿迭代法求解非线性方程组 已知非线性方程组如下 3*x1-cos(x2*x3)-1/2=0 x1^2-81*(x2+0.1)^2+sin(x3)+1.06=0 exp(-x1*x2)+20*x3+(10*pi-3)/3=0 求解要求精度达到0.00001 ———————————————————————————————— 首先建立函数fun 储存方程组编程如下将fun.m保存到工作路径中: function f=fun(x); %定义非线性方程组如下 %变量x1 x2 x3 %函数f1 f2 f3 syms x1 x2 x3 f1=3*x1-cos(x2*x3)-1/2; f2=x1^2-81*(x2+0.1)^2+sin(x3)+1.06; f3=exp(-x1*x2)+20*x3+(10*pi-3)/3; f=[f1 f2 f3]; ———————————————————————————————— 建立函数dfun 用来求方程组的雅克比矩阵将dfun.m保存到工作路径中: function df=dfun(x); %用来求解方程组的雅克比矩阵储存在dfun中 f=fun(x); df=[diff(f,'x1');diff(f,'x2');diff(f,'x3')]; df=conj(df'); ———————————————————————————————— 编程牛顿法求解非线性方程组将newton.m保存到工作路径中: function x=newton(x0,eps,N); con=0; %其中x0为迭代初值eps为精度要求N为最大迭代步数con用来记录结果是否收敛for i=1:N; f=subs(fun(x0),{'x1' 'x2' 'x3'},{x0(1) x0(2) x0(3)}); df=subs(dfun(x0),{'x1' 'x2' 'x3'},{x0(1) x0(2) x0(3)}); x=x0-f/df; for j=1: length(x0); il(i,j)=x(j); end if norm(x-x0) 牛顿法求解非线性方程组matlab源程序Newton-Raphson 求解非线性方程组matlab源程序 matlab程序如下： function hom [P,iter,err]=newton('f','JF',[;; ],,,1000); disp(P); disp(iter); disp(err); function Y=f(x,y,z) Y=[x^2+y^2+z^2-1; 2*x^2+y^2-4*z; 3*x^2-4*y+z^2]; function y=JF(x,y,z) f1='x^2+y^2+z^2-1'; f2='2*x^2+y^2-4*z'; f3='3*x^2-4*y+z^2'; df1x=diff(sym(f1),'x'); df1y=diff(sym(f1),'y'); df1z=diff(sym(f1),'z'); df2x=diff(sym(f2),'x'); df2y=diff(sym(f2),'y'); df2z=diff(sym(f2),'z'); df3x=diff(sym(f3),'x'); df3y=diff(sym(f3),'y'); df3z=diff(sym(f3),'z'); j=[df1x,df1y,df1z;df2x,df2y,df2z;df3x,df3y,df3z]; y=(j); function [P,iter,err]=newton(F,JF,P,tolp,tolfp,max) %输入P为初始猜测值，输出P则为近似解 %JF为相应的Jacobian矩阵 %tolp为P的允许误差 %tolfp为f(P)的允许误差 %max：循环次数 Y=f(F,P(1),P(2),P(3)); for k=1:max弦截法非线性方程求解2
牛顿法求非线性方程的根
非线性方程数值解法及其应用
牛顿法非线性方程求解
matlab程序设计实践-牛顿法解非线性方程
二分法和牛顿法求解非线性方程(C语言)
非线性方程求根
牛顿迭代法求解非线性方程组的代码
求解非线性方程实验报告
C++实现 牛顿迭代 解非线性方程组
弦截法求非线性方程的根
c++求解非线性方程组的牛顿顿迭代法
基于Matlab的牛顿迭代法解非线性方程组
Newton迭代法求解非线性方程
Newton 法解非线性方程组
chap3 非线性方程的求解
牛顿法解非线性方程组实验报告
matlab程序设计实践-牛顿法解非线性方程
matlab实现牛顿迭代法求解非线性方程组
牛顿法求解非线性方程组matlab源程序
相关主题
文本预览