非线性方程组求解及matlab实现分解教学课件_第1页
非线性方程组求解及matlab实现分解教学课件_第2页
非线性方程组求解及matlab实现分解教学课件_第3页
非线性方程组求解及matlab实现分解教学课件_第4页
非线性方程组求解及matlab实现分解教学课件_第5页
已阅读5页,还剩49页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

非线性方程(组)求解■非线性方程(组)数值求解基本原理多项式求根函数—roots非线性方程求解函数一fzeo非线性方程组求解函数—fsolve复习与练习按以下要求编写一个函数计算A=y/x+sin(45°)×x的值,其中x>0时,y=x;x<0时,y=2;X=0时,返回错误信息(ⅹcann'tbezero)要求:1)主函数名称为exce1,X作为输如变量,A作为输出变量;2)主函数中包括一个子函数myun用于计算y的值。引言在94536kPa(9.33atm)、300.2K时,容器中充以2mo氮气,试求容器体积。已知此状态下氮气的P√-T关系符合范德华方程,其范德华常数为a=4.17atmL/mo2,b=0.0371L/mo数学模型:范德华方程变形可得关于V的非线性方程f()=(×④nb)-nRT=0非线性方程(组)在化学计算中的作用多组分混合溶液的沸点、饱和蒸气压计算流体在管道中阻力计算·多组分多平衡级分离操作模拟计算·平衡常数法求解化学平衡问题定态操作的全混流反应器的操作分析非线性方程非线性方程包括:高次代数方程、超越方程及其它们的组合与线性方程相比,非线性方程求解问题无论从理论上还是从计算公式上都要复杂得多对于高次代数方程,当次数>4时,则没有通解公式可用,对于超越方程既不知有几个根,也没有同样的求解方式。实际上,对于n3代数方程以及超越方程都采用数值方法求近似根。非线性方程数值求解原理逐步扫描法开始输入初值A,步长hX1<=A,Y<=fX1)X,<=X1+hY2<=八(X2)否W=fCx打印A,h,x」结束逐步扫描法效率较低,常用于求根的初始近似值逐步扫描法计算示例方程x2=0的正数解解所在区域△11.1+14[1.4,11,411,411.42[1.41,1.42]1,411.4111.4124131.4141.41[1.414,1.415]0.001二分法若函数x)在区间[a,b内单调连续,且f(a)f(b)<0,则在闭区间[ab]内必然存在方程f(x)=0的根x二分法的图形解释分法的MATLAB程序k=whileabs(b-a)>epsX=(ab)/2ifsign(f(x))==sign(f(b))elseendk=k+1;end二分法是一种可靠的算法,但计算速度较慢分法计算示例方程x2=0的正数解fa)xf(e)解的范臣l.50.5.25~1.5.25l.251.51.3751.425.425l.425~1.40.仍51.351.4251.401.400~1.421.400l.4251.4125~1.4250.0125.:251,41821.4l25~1.418250.006751.4251.4151.45371.4125~1,4153751.412541535L4139381,413938

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论