第五讲MATLAB程序设计及其调试讲解.ppt

脚本文件first.m 程序调试 1 程序调试概述 一般来说，应用程序的错误有两类，一类是语法错误，另一类是运行时的错误。语法错误包括词法或文法的错误，例如函数名的拼写错、表达式书写错等。 程序运行时的错误是指程序的运行结果有错误，这类错误也称为程序逻辑错误。 3 调试命令 除了采用调试器调试程序外，MATLAB还提供了一些命令用于程序调试。 命令的功能和调试器菜单命令类似，具体使用方法请读者查询MATLAB帮助文档。 非线性优化 例4 计算下面函数在区间(0，1)内的最小值. 例6 求 的最小值点. 解：X=fminsearch(2*x(1)^3+4*x(1)*x(2)^3-10*x(1)*x(2)+x(2)^2, [0,0]) 结果为 X = 1.0016 0.8335 或在MATLAB编辑器中建立函数文件. function f=myfun(x) f=2*x(1)^3+4*x(1)*x(2)^3-10*x(1)*x(2)+x(2)^2; 保存为myfun.m，在命令窗口键入 X=fminsearch (myfun, [0,0]) 或 X=fminsearch(@myfun, [0,0]) 结果为： X = 1.0016 0.8335 有约束的多元函数最小值 非线性有约束的多元函数的标准形式为： sub.to 其中：x、b、beq、lb、ub是向量，A、Aeq为矩阵，C(x)、Ceq(x)是返回向量的函数，f(x)为目标函数，f(x)、C(x)、Ceq(x)可以是非线性函数. 在MATLAB5.x中，它的求解由函数constr实现. 函数 fmincon 格式 x = fmincon(fun,x0,A,b) x = fmincon(fun,x0,A,b,Aeq,beq) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options) [x,fval] = fmincon(…) [x,fval,exitflag] = fmincon(…) [x,fval,exitflag,output] = fmincon(…) [x,fval,exitflag,output,lambda] = fmincon(…) [x,fval,exitflag,output,lambda,grad] = fmincon(…) [x,fval,exitflag,output,lambda,grad,hessian] = fmincon(…) 参数说明： fun为目标函数，它可用前面的方法定义； nonlcon的作用是通过接受的向量x来计算非线性不等约束和等式约束分别在x处的估计C和Ceq，通过指定函数柄来使用，如： x = fmincon(@myfun,x0,A,b,Aeq,beq,lb,ub,@mycon)，先建立非线性约束函数，并保存为mycon.m：function [C,Ceq] = mycon(x) C = … % 计算x处的非线性不等约束的函数值. Ceq = … % 计算x处的非线性等式约束的函数值. lambda是Lagrange乘子，它体现哪一个约束有效. output输出优化信息； grad表示目标函数在x处的梯度； hessian表示目标函数在x处的Hessian值. 例7求下面问题在初始点（0，1）处的最优解 min sub.to 解：约束条件的标准形式为： sub.to 先在MATLAB编辑器中建立非线性约束函数文件： function [c, ceq]=mycon (x) c=(x(1)-1)^2-x(2); ceq=[ ]; %无等式约束. 然后，在命令窗口键入如下命令或建立M文件： fun=x(1)^2+x(2)^2-x(1)*x(2)-2*x(1)-5*x(2); %目标函数. x0=[0 1]; A=[-2 3];