




已阅读5页,还剩1页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
非线性规划建模实验一、二次规划标准型为:Min Z= 1/2XTHX+cTXs.t. AX=b VLBXVUB用 MATLAB 软件求解,其输入格式如下:1. x=quadprog(H,C,A,b);2. x=quadprog(H,C,A,b,Aeq,beq);3. x=quadprog(H,C,A,b,Aeq,beq,VLB,VUB);4. x=quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0);5. x=quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0,options);6. x,fval=quaprog(.);7. x,fval,exitflag=quaprog(.);8. x,fval,exitflag,output=quaprog(.);第一题: min f(x1,x2)=-2x1-6x2+x12-2x1x2+2x22s.t. x1+x22-x1+2x22x10, x20 1、写出标准形式为1 112 2 2 1 -min(,) 6Txxzx 1212 0xs.t.2、 输入命令:H=1 -1; -1 2; c=-2 ;-6;A=1 1; -1 2;b=2;2;Aeq=;beq=; VLB=0;0;VUB=;x,z=quadprog(H,c,A,b,Aeq,beq,VLB,VUB)3、运算结果为:x =0.6667 1.3333 z = -8.2222二、一般的非线性规划标准型为:min F(X)s.t AX=b G(X)Ceq(X)=0 VLB X VUB其中 X为 n维变元向量,G(X)与 Ceq(X)均为非线性函数组成的向量,其它变量的含义与线性规划、二次规划中相同.用 Matlab求解上述问题,基本步骤分三步:1. 首先建立 M文件 fun.m,定义目标函数 F(X):function f=fun(X);f=F(X);2. 若约束条件中有非线性约束: G(X) 或 Ceq(X)=0,则建立 M文件nonlcon.m定义函数 G(X)与 Ceq(X):function G,Ceq=nonlcon(X)G=.Ceq=.3. 建立主程序. 非线性规划求解的函数是 fmincon,命令的基本格式如下:(1) x=fmincon(fun,X0,A,b)(2) x=fmincon(fun,X0,A,b,Aeq,beq)(3) x=fmincon(fun,X0,A,b, Aeq,beq,VLB,VUB)(4) x=fmincon(fun,X0,A,b,Aeq,beq,VLB,VUB,nonlcon)(5)x=fmincon(fun,X0,A,b,Aeq,beq,VLB,VUB,nonlcon,options) (6) x,fval= fmincon(.)(7)x,fval,exitflag= fmincon(.)(8)x,fval,exitflag,output= fmincon(.)第二题: 编程求解 211minfxxs.t. 1、 写成标准形式: s.t. 0546321x20212minxf2、 先建立 M-文件 fun3.m:function f=fun3(x);f=-x(1)-2*x(2)+(1/2)*x(1)2+(1/2)*x(2)23、再建立主程序 youh2.m:x0=1;1;A=2 3 ;1 4; b=6;5;Aeq=;beq=;VLB=0;0; VUB=;0546310xx,fval=fmincon(fun3,x0,A,b,Aeq,beq,VLB,VUB)4、运算结果为:x = 0.7647 1.0588fval = -2.0294第三题: 求解 122122min()(441)xfexxs.t. x1+x2=01.5+x1x2 - x1 - x2=0-x1x2 10=01 先建立 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);2再建立 M 文件 mycon.m 定义非线性约束:function g,ceq=mycon(x)g=x(1)+x(2);1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10;3主程序 youh3.m 为 :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. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 三农信息人才培养计划及实施方案
- 多领域融合的智能科技项目开发协议
- 工程项目居间的合同
- 化工厂员工劳动协议
- 经营承包合同
- 劳务输出合作协议书
- 第2课 奖品数量好计算(教学设计)2024-2025学年五年级上册信息技术泰山版
- 新媒体主播合约协议书8篇
- Unit 8 Knowing the world Lesson 2 My home country 教学设计 2024-2025学年冀教版英语七年级上册
- 甘肃幼儿园塑胶施工方案
- 法律意见书(适用于股权投资)
- 新苏教版五年级科学上册活动手册答案
- 小区物业服务投标方案(技术标)
- 2025年高考作文专练(25道真题+审题立意+范文)- 2025年高考语文作文备考总复习
- 乡村研学规划方案
- 普洱市直属机关遴选笔试真题
- 2024-2030年中国电竞耳机行业市场发展分析及发展趋势与投资前景研究报告
- Unit1Myfamily单词解读(课件)Joinin外研剑桥英语五年级上册
- 中职汽修专业《汽车底盘构造与维修》说课稿
- 员工聘用合同范本(2024版)
- DL∕T 5161.6-2018 电气装置安装工程质量检验及评定规程 第6部分:接地装置施工质量检验
评论
0/150
提交评论