第六讲 求非线性方程和方程组的近似解.pptVIP

第六讲 非线性方程和方程组的近似解 一、实例与数学模型 数学模型 二、非线性方程和方程组的Matlab求解 【例3】 【例5】 补充:最优化问题的求解 一、实例与数学模型 模型的建立 二、最优化问题Matlab求解 例1 【例2】 【例3】 三、实例求解 （二）有约束最优化问题的求解 按照静态优化问题的结构是否线性分为线性规划和非线性规划．线性规划的特征是目标函数和约束条件中的函数都是决策变量的线性函数，并且约束是必不可少的（否则不存在有实际意义的解）， 1、线性规划问题求解 线性规划问题数学描述为 matlab中优化工具箱的求解函数 在MATLAB优化工具箱中，求解线性规划的函数为linprog（）函数，该函数的调用格式如下 [1]若没有不等式：AX ≤ b 存在，则令A [ ]，b [ ]. [2] 若没有等式约束: Aeq · X beq, 则令Aeq [ ], beq [ ]. [3]其中X0表示初始点 注意： 2、非线性规划问题的求解 非线性规划的一般数学描述为 为求解方便，约束条件可以进一步细分，这原规划可写为 该细分模型的求解命令如下 * * * 补充：部分优化问题的求解 问题描述 在一条20m宽的道路两侧，安装了两只2kW的路灯，它们离地面的高度分别为5m。在漆黑的夜晚，当两只路灯开启时，两只路灯连线的路面上最暗的点和最亮的点在哪里？如果路灯的高度可以在3m和9m之间变化，如何使路面上最暗点的亮度最大？ 路灯照明问题 建立如图所示的坐标系，即路面的宽度为s，每只路灯的功率分别是P，高度是h 。设两只路灯连线的路面上某点Q的坐标为（x，0），其中，假设两个光源都可以看成点光源，并记两个光源到点Q的距离分别为 r1和r2，从光源到点Q的光线与水平面的夹角分别为和，两个光源在点Q的照度分别为I1 和I2 ，则 x y h h P P 0 s x,0 r1 r2 Q 其中k是量纲单位决定的比例系数，不妨记为k 1，且 得到点Q的照度为 问题1，求路面上最暗点和最亮点的问题化为求C（x）的最小值点与最大值点 先计算C（x）的驻点，C（x）的阶导函数为 另 得 由数分的知识我们知道，只要能求出上述方程的根，我们就能求出照度函数的驻点，从而得到最小与最大值点。 问题2 ：把照度函数看出是关于变量x和h 的二元函数 ，即 那么该问题就转化成求二元函数的极值问题，分别对两个变量求偏到可得以下方程组 1、图解法 （1）一元方程的图解法 用ezplot（）函数可以绘制出给定的隐函数f x 0曲线，再绘出直线y 0，读出两条曲线的交点的横坐标即为方程的解。 ezplot函数的调用格式， ezplot fun %隐函数表达式 【例1】 （2）二元方程组的图解法 二元方程也可以用图解法求解，用ezplot（）函数将第一个方程对应的曲线绘制出来，再在同一个坐标系下绘制出第二个方程对应的曲线，得出曲线后就可以通过读取交点坐标的方式得出联立方程的根。 【例2】 2、准解析解 非线性方程中的多项式方程或者多项式方程组，根据代数中的根与系数的关系，对于次数低的我们可以求出它的解析解，对于次数高的我们可以得到高精度的数值也叫准解析解。 在matlab符号工具箱中，solve（）函数就是求解这一类问题，以下为该函数的调用格式 3、数值解法 （1）fzero 函数 fzero函数用于求单变量方程的根，所采用的算法主要是二分法、割线法等的混合方法。fzero至少要输入两个参数：函数和迭代初始值（或者有根区间），其最简单调用格式如下 S fzero eqn,x0 其中，eqn是指方程左端的函数，如果函数复杂可以通过编写函数M文件输入，如果函数形式简单，直接用inline函数输入方程左端的函数。 2 fsolve 函数 fsolve函数一般用于非线性方程组的求解（当然也可以用于方程的求解，但效果一般不如fzero函数），其调用格式如下 OPT optimset; OPT.TolX 1e-10 [x,Y,c,d] fsolve f,[-1;0],OPT 求解函数 3 roots 函数 roots 函数专门用于求单变量代数方程求根，即当f（x）为一元多项式时的求解，其调用格式如下 r roots p %p为多项式的系数向量（按降幂排列），输出r为f（x） 0的全部根（包括复根） 三、实例的求解 1、求路面上最暗点与最亮点 Q的照度函数 照度函数的导数 将所给的实际数据，P 2,h 5 , s 20 代入上面两个式子，在[0,20]上用matlab画出这两个函数的图形，找出照度函数的驻点大约在那些数值附近，然后以这些值为初始值，用fzero 函数求出这些驻点的横坐标，然后再求出相应的函数值，比较这些函数值