非线性方程-概念应用及解法
一:什么是线性方程
线性方程也称一次方程式,指变量均是一次幂的方程,其一般的形式是ax+by+…+cz+d=0,线性方程的本质是等式两边乘以任何相同的非零数,方程的本质都不受影响。因为在笛卡尔坐标系上任何一个一次方程的表示都是一条直线,组成一次方程的每个项必须是常数或者是一个常数和一个变量的乘积,且方程中必须包含一个变量,因为如果没有变量只有常数的式子是代数式而非方程式。
二:什么是非线性方程
非线性方程,就是因变量与自变量之间的关系不是线性的关系,这类方程很多,平方关系x^2、对数关系log(x)、指数关系2^x、三角函数关系sin(x)等等。例如e^x-con(x)=0就是一个非线性方程。与线性方程相比,无论是解的存在性,还是求解的计算公式,非线性方程都要复杂得多。
三:非线性方程的应用
随着科技技术的飞速发展,科学计算越来越显示出其重要性。科学计算的应用之广已遍及各行各业,例如气象资料的分析图像,飞机、汽车及轮船的外形设计,高科技研究等都离不开科学计算。因此经常需要求非线性方程 f(x)=0的根。
例如我们要求导弹的轨迹y,y又与导弹的速度v,空间位置p,加速度a,时间参数t等等有关,v又可以分为vx,vy,vz,p也可以分为px,py,pz,每个变量之间不能仅仅是简单的线性关系,如速度v和时间t之间就要考虑重力,空气阻力等等,往往最终列出来的就是一长串非线性方程组,但是这样的方程组的解又极具现实意义,往往科学工程经济学的问题都会转为对非线性方程组的求解。
四:非线性方程的解法
1.实根对分法
对分法又称二分法,设函数f(x)在[a,b]上连续,且f(a)f(b)<0,则f(x)在[a,b]上至少有一零点,这是微积分中的介值定理,也是使用对分法的前提条件。
注意:对分法算法简单,然而,若[a,b]上有几个零点时,只能算出其中一个零点;另一方面,即使[a,b]上有零点,也未必有f(a)f(b)<0,这就限制了对分法的使用范围;对分法只能计算方程的实根。
2.不动点迭代法
对给定的非线性方程f(x)=0,将它转换成等价形式:x=φ(x).给定初始值x0,构造迭代序列x(k+1)=φ(x(k)),k=1,2,…,如果迭代收敛
lim x(k+1)=lim φ(x(k))=α
有α=φ(α),则α就是方程f(x)的根。在计算中,当|x(k+1)−x(k)|小于给定精度控制量时,取x(k+1)为方程的根。
3.牛顿迭代法
对于非线性方程f(x)=0可构造多种迭代格式x(k+1)=φ(x(k)),牛顿迭代法是借助于对函数f(x)在x0处做泰勒展开,取其线性部分构造的一种迭代格式。
将f(x)=0在初始值x0做泰勒展开:
f(x)=f(x0)+f′(x0)/(x−x0)+f′′(x0)*(x−x0)^2/2!+…
取展开式的线性部分作为f(x)的近似,则
f(x0)+f′(x0)(x−x0)=0
假设f′(x0)≠0,则
x=x0−f(x0)/f′(x0)
因此得到牛顿迭代法的迭代格式:
x(k+1)=x(k)−f(x(k))/f′(x(k)) , k=1,2,..
4.弦截法
在牛顿迭代法中:x(k+1)=x(k)−f(x(k))/f′(x(k)),k=1,2,…
需要求解函数f(x)的一阶导数,通常比较耗时,因此可用差商代替导数:
f′(x(k))=(f(x(k))−f(x(k-1))) / (x(k)−x(k-1))
给定初始值x0和x1,那么弦截法的迭代格式为:
x(k+1) = x(k)−((x(k)−x(k-1)) * f(x(k)) / (f(x(k))−f(x(k-1)) , k=1,2,…
PS:以上内容只是对非线性方程的简单介绍,今后会不定期更新。
非线性方程-概念应用及解法相关推荐
- 计算机基础-二元一次方程组的概念及基本解法
二元一次方程组 二元一次方程组的概念及基本解法 二元一次方程组的概念及基本解法
- matlab多元方程整数解,matlab多元非线性方程组解法
流放之路一天赚几百是真的吗 你想太多了.如果是搬砖收入的话一天能有2E就不错了.当然脸正的时候也可以一次就赚上10E.5173上1E=3+. 用MATLAB解多元非线性方程组,CSS布局HTML小编今 ...
- 阿当姆斯matlab,数值计算方法与MATLAB应用
第1篇 数值计算的基本方法和概念 第1章 算法与误差 1.1 算法 1.2 误差 第2章 方程求解 2.1 引言 2.2 二分法 2.3 迭代法 2.4 牛顿法 2.5 弦截法 2.6 解非线性方程组 ...
- matlab求解helmholtz,MATLAB与科学计算(第2版)
第1章 安装及使用前的准备1.1 MATLAB 6.x简介1.1.1 21世纪的科学计算语言1.1.2 MATLAB的发展历史1.1.3 MATLAB 6.x的新特点1.1.4 MATLA ...
- matlab与科学计算 王沫然,MATLAB与科学计算(第3版) 王沫然著 电子工业出版社 9787121180521...
商品描述: 基本信息 书名:MATLAB与科学计算(第3版) 定价:49.80元 作者:王沫然 编著 出版社:电子工业出版社 出版日期:2012-10-01 ISBN:9787121180521 字数 ...
- matlab求解常微分方程——从原理到实践(代码详解)
目录 理论知识 一.概念 二.解法 matlab微分方程求解 一.解析解 1.1 解析解的存在 1.2 解析解的解法 1.3 实例 二.数值解 2.1 概述 2.2 优化措施 2.3 解法 2.4 检 ...
- matlab科学计算 书,MATLAB与科学计算(第3版)
<MATLAB与科学计算(第3版)> 第1章 安装及使用前的准备1 1.1 MATLAB简介1 1.1.1 21世纪的科学计算语言1 1.1.2 MATLAB的发展历史3 1.1.3 MA ...
- hihocoder 网络流二·最大流最小割定理
网络流二·最大流最小割定理 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi:在上一周的Hiho一下中我们初步讲解了网络流的概念以及常规解法,小Ho你还记得内容么? ...
- 手把手教你用WOS数据库解决英文文献
一.确定关键检索词. 1.在看文献之前,首先需要明确我们的研究方向,再进行主题分析."主题"即为文献或检索课题所论述的.研究的具体对象或问题,就是我们通常说的中心内容.而主题的确定 ...
- 控制理论中的几种稳定性
文章目录 1.李雅普诺夫稳定性 (Lyapunov Stable) 1.1 概念 1.2 解法 2.渐进稳定性 (Asympototically Stable) 2.1 概念 2.2 解法 3.指数稳 ...
最新文章
- cocos2d-js 3.0 RC0 监听返回键、菜单键、进入后台(home键)、恢复显示等事件
- Maven私服安装及配置——(十二)
- Itext导出pdf文件
- Java知识点:条件编译
- java nio与io_Java NIO和IO的区别(转)
- 记一次CPU飙升的问题分析解决思路(转)
- python中qt有哪些控件_使用PyQt5调用Qt程序,基础PythonQt控件的使用方法
- Programming: Principles and Practice Using C++
- 目标跟踪(1)基于OpenCV实现单目标跟踪
- FLASH知识及读写擦除规则
- python对象转json字符串,及json字符串的格式化
- 安装ruby-2.3.1
- 粤嵌开发板之手机WIFI摄像头
- 吐槽大会,加个好友,分享资源
- Qt安装包百度网盘下载分享
- Re-脱壳技术 脱壳实战(3): 脱壳及修复IAT
- android 普通蓝牙源码解读
- ipad2/3分辨率
- Android 红米Note5 刷机魔趣系统
- 定义并理解下列术语,说明它们之间的联系与区别
热门文章
- html设置单元格不可编辑状态,excel单元格不能编辑 如何取消Excel 2013中单元格的直接编辑状态...
- Windows-Exploit-Suggester
- xgboost的plot_importance绘图时出现的f0、f1、f2、f3、f4、f5等改为对应特征的字段名
- 计算机无法识别银行usb是,电脑不能识别USB设备(加密设备、网银U盾等)的解决方法...
- python 小海龟 教案_幼儿园小班健康活动《勇敢的小海龟》教案
- [2021时空AI白皮书]时空人工智能:关键技术
- [IOS APP]毕淑敏经典有声小说
- 计算机Word2010在线做题,Word2010试题练习题库(答案).doc
- 对一个浮点数保留两位小数
- LINUX下运行.sh文件出现Syntax error: end of file unexpected (expecting “then”)最方便解决方法