非线性方程(组)求解.pptVIP

* 非线性方程（组）求解 一、实验目的 1.熟悉用matlab软件求解非线性方程（组）的命令函数. 2.掌握求非线性近似解的常用数值方法——迭代法. 3.了解分叉与混沌现象. 二、实验内容 1.用matlab求方程的解. 2.用matlab求方程组的解. 3.二分法 4.迭代法：牛顿切线迭代法，弦截法 5.了解分叉与混沌现象. 1.用matlab软件求方程的解 Matlab软件求方程f(x)=0近似解的命令是fzero,具体用法为： （1）建立函数：f=inline(‘表达式’) （2）求函数零点：c=fzero(f,[a,b]) c=fzero(f,x0) %求函数在区间内的零点 %求函数f在x0附近的零点 roots(c) 示例1 求五次曲线y=x5与直线y=1-x交点的横坐标. [0,1]内必有实根 [0,1]内必有唯一实根 f=inline(‘x^5+x-1’); c=fzero(f,[0,1]) c = 0.75487766624669 c=roots([1,0,0,0,1,-1]) c = -0.87743883312335 + 0.74486176661974i -0.87743883312335 - 0.74486176661974i 0.50000000000000 + 0.86602540378444i 0.50000000000000 - 0.86602540378444i 0.75487766624669 c=fzero(f,0.8) c = 0.75487766624669 示例2 求螺旋线 与平面 的交点. f=inline(4*cos(x)+4*sin(x)+0.5*x-2); c=fzero(f,[0,pi]) c = 2.19606411841355 x=0:0.1:8*pi f1=4*cos(x)+4*sin(x)+0.5*x-2; f2=0; plot(x,f1,x,f2,‘go’),grid 2.用matlab求方程组的解 Matlab软件求上述非线性方程组的数值解命令是： [x,fval]=fsolve(fun,x0) 示例3 求方程组 的近似解. [x,fval]=fsolve(fun,x0) function f=group(x) f=[sin(x(1))+x(2) +x(3)^2*exp(x(1))-4; x(1)+x(2)*x(3); x(1)*x(2)*x(3)]; [x,fval]=fsolve(‘group’,[1,1,1]) x = Columns 1 through 2 0.00046164326023 -0.00023082156767 Column 3 1.99948065120099 fval = 1.0e-006 * -0.26357352345130 0.12000180119623 -0.21305910169336 x = Columns 1 through 2 0.00046164326023 -0.00023082156767 Column 3 1.99948065120099 fval = 1.0e-006 * -0.26357352345130 0.12000180119623 -0.21305910169336 f=inline([sin(x(1))+x(2)+ x(3)^2*exp(x(1))-4; x(1)+x(2)*x(3); x(1)*x(2)*x(3)]); [x,fval]=fsolve(f,[1,1,1]) fzero命令虽然简单，但从科学的态度出发，我们应该充分了解 “隐藏”在该命令背后的原理、算法以及实现过程. 因此，掌握求解一般非线性方程近似根的常用数值方法，熟悉 相关原理等对从事科学研究是十分必要的. 基本理论知识 1.零点定理 2.压缩映射原理（证明参看工科数学分析上册第77页） 设函数f(x)定义在[a,b]上，若满足以下两个条件： 则当L1时，f(x)在[a,b]上有唯一的不动点. 即f(x)=0有唯一的根. 零点定理 则 3.二分法求方程的近似解 开始 精度m,初始值a, b f(x) break m=(a+b)/2 a=m b=m f (m)=0 ? f(a)f(m)＞0 |a-b|＜0.0005 结束 输出所求的近似根m m=(a+b)/2 是 否 否 是 否 示例4 用二分法求方程 的近似解. 在[0,1]上函数至少有一个零点. a=0;b=1; f=inli