《非线性方程数值解法的探讨与MATLAB实现(DOC 16页)》由会员分享，可在线阅读，更多相关《非线性方程数值解法的探讨与MATLAB实现(DOC 16页)（16页珍藏版）》请在金锄头文库上搜索。
1、 非线性方程数值解法的探讨与MATLAB实现摘 要非线性方程在实际问题中经常出现，如科学与工程计算，因此研究和探讨非线性方程求解的方法具有重要的实际意义. 摘要字数：150-500，内容至少分两段：第一段是背景叙述、第二段论文研究的方法与结论等 本文首先介绍了非线性方程的概念及相关背景；其次描述了非线性方程常用的一些数值方法：二分法，牛顿迭代法，割线法；然后针对各种方法编写MATLAB程序，并对非线性方程的实例进行数值计算，比较各种算法的优劣；最后介绍了割线法在实际生活中的应用.关键词：非线性方程，二分法，牛顿迭代法，割线法关键词的个数至少3个，用逗号分隔、结尾无符号Discussion an
2、d MATLAB Realization of Nonlinear Numerical Solution 题目实词的首写字母大写，其余均小写 Author: Feng Wenghao Tutor: Zhang RongyanAbstractIn practical problems, nonlinear equations often appear, such as scientific and engineering computing. Therefore, to study and explore methods of nonlinear equations has important
3、practical significance.此部分与中文摘要一致；标点符号均为英文状态下的符号. 其它格式与中文摘要一致In this paper, we describe the concept and background of nonlinear equations first; second, it describes the commonly used numerical methods of nonlinear equations: Dichotomy, Newton iteration method, Secant method. Then we use MATLAB to w
4、rite computer program aimed at these methods，do numerical calculations in some practical problems of nonlinear equations and make comparison of various methods. Finally，the application of secant method in real life.Keywords: Nonlinear equation, Dichotomy, Newton iteration method, Secant method（空一行、小
5、五号、1倍行距）目 录目录中间空7个字符，居中、三号黑体、1.5倍行距（空一行、三号、1.5倍行距）1 绪论11.1 数值方法背景简介11.2 非线性方程简介21.2.1 非线性方程的背景21.2.2 非线性方程的研究内容21.2.3 根的存在性定理32 非线性方程的数值解法42.1 引言42.2 二分法42.2.1 二分法简介42.2.2 二分法的原理42.3 牛顿迭代法52.3.1 牛顿迭代法的简介52.3.2 牛顿迭代法的原理52.3.3 牛顿迭代法的几何意义62.4 割线法72.4.1 割线法简介72.4.2 割线法的原理83 非线性方程的MATLAB实现93.1 二分法93.1.1
6、二分法的MATLAB程序93.1.2 应用举例103.2 牛顿迭代法123.2.1 牛顿迭代法的MATLAB程序123.2.2 应用举例133.3 割线法143.3.1 割线法的MATLAB程序143.3.2 应用举例154 方法的分析与对比174.1 构造非线性方程迭代公式174.2 计算迭代公式175 实际应用215.1 引言215.2 问题提出215.3 模型建立225.4 模型求解23结论25致谢26参考文献271、1.5倍行距2、小四号、3、一级标题黑体字、其余为宋体4、注意各级标题和页码要与文章中的标题和页码一致，更改后及时更新1 绪论一级标题：段前、段后25磅，两倍行距，三号、黑
7、体，数字新罗马体1.1 数值方法背景简介二级、三级等：段前、段后0.5行，1.5倍行距，小四号、黑体，数字新罗马体数值分析又称数值计算方法或数值方法，是研究数值问题算法的一门学问，它研究如何借用计算工具求得数学问题的数值解答这里的数学问题是指数值问题，所谓数值问题，是给出一组数值型数据，通常是一些实数，称为初始数据，去求解另外一组数值型数据，问题的本身反映了这两组数据之间的一定关系函数计算、方程求根就是数值问题的典型例子除此之外，数学中还存在着大量的非数值问题，如定理的证明、几何的作图和组态的枚举等问题都不是研究的对象同时数值分析(Numerical Analysis)是研究数值求解各类数学问
8、题的方法和相应的数学理论研究的对象是数学问题，所用的方法是数学方法，因此也称数值数学(Numerical Mathematics)研究的内容可划分为以下几个主要方面：正文：1.5倍行距、小四号、汉字宋体、数字新罗马体首行缩进两个字符。标点符号：句号用英文状态全角“”，其余标点符号用中文状态下。(1) 数值代数：主要包括线性方程组和非线性方程组的数值解法、特征值与特征向量的数值计算等内容；(2) 数值逼近；(3) 常微分方程和动力系统的数值解法；(4) 偏微分方程的数值解法；(5) 最优化理论与方法；(6) 误差理论数值方法的特点，就是说求解一个数学问题的数值方法是要给出该问题的一个近似的数值结
9、果因此，首先数值结果要能算得出来其次结果应有一定的精度，满足实际问题的要求，一般要求误差满足指定的值第三，计算时间应尽可能少求解一个数学问题的数值方法一般有多种在保证所需精度的条件下，计算时间越少的数值方法越好对串行机而言，这相当于计算量越少越好1参考文献标注在右上角、小四号、新罗马字体、参考文献的编号要与后面参考文献一致，不能乱编 1.2 非线性方程简介1.2.1 非线性方程的背景非线性科学是近年来备受科学工作者和工程技术人员关注的科学，求非线性数学物理问题(包括常微、偏微边值问题，积分方程，积分微分方程等)、非线性力学、非线性优化、数理经济学等问题，又是非线性科学中最基本的问题，而上述问题
10、最终都归结为求非线性方程，因此研究非线性方程的求解方法有着十分重要的意义.1.2.2 非线性方程的研究内容科学计算中常常要求下列非线性方程组： (1.1)公式居中、用公式编辑器编辑；公式的编号按照章节编辑，编号右对齐其中为给定在维欧式空间中的区域上的实值函数.引进向量符号，令则(1.1)可写成 (1.2)这里表示定义在上取值于的非线性映像简记为若存在，使，则称为方程(1.2)的解2.非线性方程可分为两类：一类是多项式方程，这类方程可以定义为：，. 另一类是非多项式方程，它不能用多项式方程的形式表示，没有固定的形式.求解第一类多项式方程，现在已经有了比较成熟的理论和方法，现在比较常用的一种数值方
11、法是迭代法，能通过迭代次数的增加，从而越来越接近方程的解；求解第二类非多项式方程，是现在数学领域中的一个重点研究方向，一般来说，求解此类方程是采用随机搜索的办法.1.2.3 根的存在性定理设为实变量的实值函数，求方程的解，或者说求的根，一般分为两个步骤. 首先，利用数学分析工具确定其根的分布区域，然后在该区域中利用更精细的方法逐步逼近法求得其根.设为的根，即 如果存在正整数，使，且，则称为其重根.当时，称之为单根.下述介值定理是确定根存在区间的最简单方法.定理1：若在上连续，且，则在上至少有一个根.定理2：设在的邻域，连续可微且，则方程式在有解存在且唯一3.2 非线性方程的数值解法2.1 引言
12、本章重点介绍了二分法、牛顿迭代法和割线法，通过他们的原理及方法推理来说明这几种方法在求解非线性代数方程的重要性，以及为下面编写算法奠定基础.2.2 二分法2.2.1 二分法简介二分法是区间迭代法的一种.它是重复运用根的存在性定理，每次将区间压缩一半且其中一个区间至少包含一个根，逐步缩短区间，直至最终区间长度满足一定的精度要求为止.2.2.2 二分法的原理3 非线性方程的MATLAB实现3.1 二分法3.1.1 二分法的MATLAB程序function k，x，wuca，yx=erfen(a，b，abtol)% 和为闭区间的左右端点，为事先给定的精度，为二分法的次数，为精确根的近似值，而为函数在
13、的近似值处的函数值，而是使用次二分法得到的小区间长度的一半，即实根的近似值的绝对精度限.a(1)=a; b(1)=b;ya=fun(a(1); yb=fun(b(1); %程序中调用的fun.m，此函数所求的非线性方程.if ya* yb0， % 零点定理条件判定.disp (注意：ya*yb0，请重新调整区间端点a和b.)， returnendmax1=-1+ceil(log(b-a)- log(abtol)/ log(2); % ceil是向方向取整for k=1: max1+1a;ya=fun(a);b;yb=fun(b); x=(a+b)/2;yx=fun(x); wuca=abs(b-a)/2; k=k-1;k，a，b，x，wuca，ya，yb，yxif yx=0a=x; b=x;else if yb*yx0b=x;yb=yx;elsea=x; ya=yx;endif b-ax=-4:0.1:4;y=x.3-x +1; plot(x, y) grid，gtext(y=x3-x+1) 运行结果：图3.1 函数的图像图形的编号按照章节编写；图名在
