2012数学建模5.ppt

第五章 数值计算 引言 数值计算是MATLAB最强大的功能之一，其矩阵的计算能力强、编程简单、计算代码经过精心地优化，计算效率很高。在数值计算方面，涉及的数学学科面广，例如矩阵代数计算、微积分计算、快速傅立叶变换、常微分方程计算、偏微分方程计算、数据拟合等等。MATLAB提供了非常丰富的解决这类问题的命令。在这一章，我们将主要以实例来介绍解决问题的方法，并假定读者具有线性代数、数学分析、计算方法等背景知识。 计算结果为： Optimization terminated successfully. x = 0.6667 1.3333 最优解 fval = -8.2222 最优值 exitflag =1 output = iterations: 3 algorithm: medium-scale: active-set firstorderopt: [] cgiterations: [] lambda = lower: [2x1 double] upper: [2x1 double] eqlin: [0x1 double] ineqlin: [3x1 double] 5.5.3 非线性有约束最优化问题 对于非线性函数的求最优化问题情况较为复杂，没有一个固定的方法。非线性函数的求最优化模型为： 其中f(x)，c(x)，ceq(x)，为非线形函数。 求非线性函数的求最优化问题命令的语法为： 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 = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options,P1,P2, ...) [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为非线性函数。 【例5.5.5】求以下非线性函数的最小值及最优解。搜索初始值为[10，10，10] 解：这里 x0 = [10; 10; 10]; % Starting guess at the solution A=[-1,-2,-2;1,2,2] b=[0;72] [x,fval] = fmincon(f,x0,A,b) 计算结果为： x = 24.00 12.00 12.0 fval = -3456 5.5.4 多目标线性规划 在实际问题中目标函数可能有多个，例如一个大型的工程其最终的目标是多方面的，既要考虑工程的成本尽量少，又要考虑企业投产带来的生态和社会效应等多方面的因素。这就是一个多目标的规划问题。 多目标规划命令的语法为： x = fgoalattain(fun,x0,goal,weight) x = fgoalattain(fun,x0,goal,weight,A,b) x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq) x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub) x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon) x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq, lb,ub,nonlcon,options) x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq, lb,ub,nonlcon,options,P1,P2,...) [x,fval] = fgoa