原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
function [G,Ceq]=mycon1(x) G(1)=x(1)*x(1)+x(2)*x(2)+x(3)*x(3)-100 G(2)=60-x(1)*x(1)+10*x(3)*x(3) Ceq(1)=x(1)+x(2)*x(2)+x(3)- 80 Ceq(2)=x(1)^3+x(2)*x(2)+x(3)- 80 注意： [1] fmincon函数提供了大型优化算法和中型优化算法。默认时，若在fun函数中提供了梯度（options参数的GradObj设置为’on’），并且只有上下界存在或只有等式约束，fmincon函数将选择大型算法。当既有等式约束又有梯度约束时，使用中型算法。 [2] fmincon函数可能会给出局部最优解，这与初值X0的选取有关。 1、写成标准形式： 2x1+3x2 6 s.t x1+4x2 5 x1,x2 0 例 2、先建立M-文件 fun3.m: function f=fun3(x); f=-x(1)-2*x(2)+(1/2)*x(1)^2+(1/2)*x(2)^2 3、再建立主程序youh2.m： x0=[1;1]; A=[2 3 ;1 4]; b=[6;5]; Aeq=[];beq=[]; VLB=[0;0]; VUB=[]; [x,fval]=fmincon(fun3,x0,A,b,Aeq,beq,VLB,VUB) 4、运算结果为： x = 0.7647 1.0588 fval = -2.0294 1．先建立M文件 fun4.m,定义目标函数: function f=fun4(x); f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); x1+x2=0 s.t. 1.5+x1x2 - x1 - x2 0 -x1x2 –10 0 例 2．再建立M文件mycon.m定义非线性约束： function [g,ceq]=mycon(x) g=[1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10]; ceq=[]; 3．主程序为: x0=[-1;1]; A=[];b=[]; Aeq=[1 1];beq=0; vlb=[];vub=[]; [x,fval]=fmincon(fun4,x0,A,b,Aeq,beq,vlb,vub,mycon) 4. 运算结果为： x = -1.2250 1.2250 fval = 1.8951 习题：求解以下非线性规划 非负条件下线性最小二乘 lsqnonneg 适合如下模型： 注意：约束只有非负约束 yonkey_mwh@163.com 2015.5.15 语法： x =lsqnonneg(c,d) x =lsqnonneg(c,d,x0) x =lsqnonneg(c,d,x0,options) yonkey_mwh@163.com 2015.5.15 语法： x = lsqlin(C,d,A,b) x = lsqlin(C,d,A,b,Aeq,beq) x = lsqlin(C,d,A,b,Aeq,beq,lb,ub) x = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0) x = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0,options) [x,resnorm] = lsqlin(...) [x,resnorm,residual] = lsqlin(...) [x,resnorm,residual,exitflag] = lsqlin(...) [x,resnorm,residual,exitflag,output] = lsql
文档评论（0）