文档介绍：该【数学建模4-非线性规划模型求解 (2) 】是由【莫比乌斯】上传分享，文档一共【4】页，该文档可以免费在线阅读，需要了解更多关于【数学建模4-非线性规划模型求解 (2) 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。湖南第一师范学院数学系实验报告
姓名:
学号:
专业:
数学与应用数学
班级:
12级
课程名称:
线性规划与数学建模
实验名称:
非线性规划模型的Matlab求解
实验类型:
基础实验
实验室名称:
数学建模实验室
实验地点:
实A302
实验时间:
2015年6月25日
指导教师:
曾
成绩评定:
一、实验目的与要求:
掌握非线性规划问题的求解方法。
熟练掌握函数fmincon,fminbnd,fminsearch,fminunc等的用法。
会利用matlab优化工具箱求解简单的非线性规划问题。
二、实验环境(实验器材、环境要求):
计算机
Matlab软件
三、实验内容(实验原理、任务等):
1、求解下列非线性规划问题:
2、(供应与选址问题)某公司有6个建筑工地要开工,每个工地的位置(用平面坐标系a,b表示,距离单位:千米)及水泥日用量d(吨)由下表给出。目前有两个临时料场位于A(5,1),B(2,7),日储量各有20吨。假设从料场到工地之间均有直线道路相连。
(1)试制定每天的供应计划,即从A,B两料场分别向各工地运送多少吨水泥,使总的吨千米数最小。
四、实验具体步骤:
1、
1),定义目标函数:
functionf=fun4(x);
f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1);
2):
function[g,ceq]=mycon(x)
g=[+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10];
ceq=[];
3):
x0=[-1;1];
A=[];b=[];
Aeq=[11];beq=[0];
vlb=[];vub=[];
[x,fval]=fmincon('fun4',x0,A,b,Aeq,beq,vlb,vub,'mycon')
4)结果为:
x=-
fval=
2、建立数学模型:
min=*x11+*x21+*x31+*x41+*x51+*x61+*x12+*x22+*x32+*x42+*x52+*x62;
+x12=3;
x21+x22=5;
x31+x32=4;
x41+x42=7;
x51+x52=6;
x61+x62=11;
x11+x21+x31+x41+x51+x61<=20;
x12+x22+x32+x42+x52+x62<=20;
得到最优解:x11=3,x12=5,x14=7,x16=1,x23=4,x25=6,x26=10,其余的为0,。
建立M文件:
f=[,,,,,,,,,,,];
A=[1,1,1,1,1,1,0,0,0,0,0,0;0,0,0,0,0,0,1,1,1,1,1,1];
B=[20;20]
Aeq=[1,0,0,0,0,0,1,0,0,0,0,0;0,1,0,0,0,0,0,1,0,0,0,0;0,0,1,0,0,0,0,0,1,0,0,0;0,0,0,1,0,0,0,0,0,1,0,0;0,0,0,0,1,0,0,0,0,0,1,0;0,0,0,0,0,1,0,0,0,0,0,1];
beq=[3,5,4,7,6,11];
[xopt,fopt]=linprog(f,A,b,Aeq,beq)
B=2020
Exiting:Oneormoreoftheresiduals,dualitygap,ortotalrelativeerror
hasstalled:
thedualappearstobeinfeasible(andtheprimalunbounded).
(Theprimalresidual<TolFun=-008.)
xopt=+009*
-
-
-
-
-
-
fopt=-+010
五、实验心得(质疑、建议):
第一题把Aeq=[11]写成了Aeq=[1,1****惯性以逗号分隔,然而并非所有都是逗号做分隔号的,所以在对不同类型程序语法进行编写时,要注意结合其自身特点,不能被惯性思维影响。第二题建立模型由一定难度,需要认真分析,用软件编程时也挺麻烦,不知道lingo软件是否也能解出此题,二者相比哪个操作更便捷。