数学建模讲座之数学建模案例分析.pdf
《数学建模讲座之数学建模案例分析.pdf》由会员分享,可在线阅读,更多相关《数学建模讲座之数学建模案例分析.pdf(239页珍藏版)》请在咨信网上搜索。
1、简要提纲伟tw斶介2.LINDO/LINGO软件的简单使用简介3.建模与求解实例(结合软件使用)4.LINGO软件基本使用方法2014-12-10数学建模1.优化软件简介2014-12-10数学建模I优化模型和优化软件的重要意义(最)优化:在一定条件下,寻求使目标最大(小)的决策最优化是工程技术、经济管理、科学研究、社会生活中经常遇到的问题,如:结构设计资源分配 生产计划 运输方案 解决优化问题的手段经验积累,主观判断作试验,比优劣建立数学模型(优化模型),求最优策略(决策)CUMCM赛题:约一半以上与优化有关,需用软件求解2014-12-10数学建模(最)优化理论是运筹学的基本内容运筹学(O
2、R:Operations/Operational Research)1 OR/MS/管理科学(MS:Management Science)DS 决策科学(DS:Decision Science)优化(Opti m i zati on),规划(Program m i ng)无约束优化线性规划2014-12-非线性规划整数规划组合优化多目标规划目标规划网络优化模优化问题的一般形式约束条件1问题三要素:决策变量;目标函数;约束条件min/(x)一s.t.(x)=0,i=1,mg/(x)4 0,j=1,,I策变G。况”可行解(满足约束)与可行域(可行解的集合)最优解(取到最小/大值的可行解)2014-
3、12-10数学建模无约束优化:最优解的分类和条件|给定一个函数/(X),寻找X*使得/(X*)最小,即 Min/(%)其中X=(占.2,,匕尸金况X局部最优解口 必要条件 充分条件 W*)=(4,-,4)r=oV/(x*)=0,V2/(x*)0 Hessian 阵全局最优解最优解在可行域边界上取得时不能用无约束优化方法求解2014-12-10数学建模约束优化的 简单分类数学规划mm f(%)s.t.a,(x)=0,7=1,mgy(x)0,j=1,,/Q c 况连续优化离散优化 线性规划(LP)目标和约束均为线性函数 非线性规划(NLP)目标或约束中存在非线性函数/二次规划(QP)目标为二次函数
4、、约束为线性 整数规划(IP)决策变量(全部或部分)为整数/整数线性规划(ILP),整数非线性规划(INLP)/纯整数规划(PIP),混合整数规划(MIP),一般整数规划,04(整数)规划2014-12-10 数学建模常用优化软件1.LINDO/LIN8软件2.MATLAB优化工具箱3.EXCEL软件的优化功能4.SAS(统计分析)软件的优化功能5.其他2014-12-10数学建模2014-12-10数学建模J LINDO公司软件产品简要介绍 对哥(Ch icago)大学的Linus Sch rage教授于 1980 年前后开发,后来成立LINDO系统公司(LINDO Systems Inc.
5、),网址:LINDO:Linear INteractive and Discrete Optimizer(V6.1)LINGO:Linear INteractive General Optimizer(V9.0)LINDO API:LINDO Application Programming Interface(V3.0)Wh ats Best!:(Spread sh eet e.g.EXCEL)(V8.0)演示(试用)版、学生版、高级版、超级版、工业版、扩展版(求解问题规模和选件不同)2014-12-10 数学建模2014-12-10数学建模2014-12-10数学建模】建模时需要注意的几个基
6、本问题TI尽量使用实数优化,减少整数约束和整数变量2、尽量使用光滑优化,减少非光滑约束的个数如:尽量少使用绝对值、符号函数、多个变量求 最大/最小值、四舍五入、取整函数等3、尽量使用线性模型,减少非线性约束和非线性变量的个数(如工。5改为xv5y)4、合理设定变量上下界,尽可能给出变量初始值5、模型中使用的参数数量级要适当(如小于103)2014-12-10数学建模3.LINDO/LINGO软件的使用简介2014-12-10数学建模需要掌握的几个重要方面LINDCM-正确阅读求解报告(尤其要掌握敏感性分析)2、LINGO:掌握集合 的应用;正确阅读求解报告;正确理解求解状态窗口;学会设置基本的
7、求解选项;掌握与外部文件的基本接口方法2014-12-10 数学建模例1加工奶制品的生产计划 Ihs I牛奶或12小时8小时2斤A14公斤A2获利24元/公斤获利16元/公斤每天:50桶牛奶 时间480小时 至多加工100公斤A制订生产计划,使每天获利最大 35元可买到1桶牛奶,买吗?若买,每天最多买多少?可聘用临时工人,付出的工资最多是每小时几元?A1的获利增加到30元/公斤,应否改变生产计划?2014-12-10数学建模1桶 牛奶12小时3公斤A获利24元/公斤获利16元/公斤或、时4公斤A2天50桶牛奶时间480小时至多加工100公斤A1决策变量巧桶牛奶生产A1整桶牛奶生产A2目标函数获
8、利24X3均 获利16X4x2每天获利 Max z-72x+64x?-L 乙原料供应 约束条件劳动时间 加工能力 非负约束 2014-12-1G西+/V 5012再+8x2 4803匹 0数学建模性划型D 线规模缶模型求解IOBJECTIVE FUNCTION VALUEMil kOUt x1)3360.000VARIABLE VALUEXI 20.000000X2 30.000000REDUCED COST0.0000000.000000max 72x1+64x2st2)xl+x2503)12x1+8x24804)3xl100endROW SLACK OR SURPLUS DUAL PRIC
9、ES2)0.0000003)0.0000004)40.00000048.0000002.0000000.000000DO RANGE(SENSITIVITY)ANALYSIS?NoNO.ITERATIONS=220桶牛奶生产A0 30桶生产A2,利润3360元。2014-12-10数学建模模型求解OBJECTIVE FUNCTION VALUE1)3360.000VARIABLE VALUEXI 20.000000X2 30.000000ROW2)0.0000003)0.0000004)40.000000NO.ITERATIONS=REDUCED COST0.0000000.000000SLA
10、CK OR SURPLUS DUAL PRICESred uced cost 值表 示当该非基变量 增加一个单位时(其他非基变量 保持不变)目标 函数减少的量(对 max型问题)48.0000002.0000000.000000为了使该非基变 量变成基变量,目标函数中对应 系数应增加的量2014-12-10数学建模结果解释max 72x1+64x2st-2)xl+x2503)12x1+8x24804)3xl100OBJECTIVE FUNCTION VALUE1)3360.000VARIABLE VALUEXI 20.000000X2 30.000000REDUCED COST0.000000
11、0.0000005三种资源ROW原料无剩余 2)时间无剩余 3)加工能力剩余40 4)SLACK OR SURPLUS DUAL PRICES0.000000 48.0000000.000000 2.00000040.000000 0.000000资源”剩余为零的约束为紧约束(有效约束)2014-12-10数学建模结果解释最3下“资源”增 加rf位时“双益”的影子神寤OBJECTIVE FUNCTION VALUE1)3360.000VARIABLE VALUEXI 20.000000X2 30.000000REDUCED COST0.0000000.000000ROW 原料增1单位,利润增4
12、8 2)时间加1单位,利润增2 3)能力增减不影响利润 4)SLACK OR SURPLUS DUAL PRICES0.0000000.00000040.00000048.0000002.0000000.000000 35元可买到1桶牛奶,要买吗?35(或V 55)号与“=(或“V=35)功能相.同一2.变量与系数间可有空格(甚至回车),但无运算符3.变量名以字母开头,不能超过8个字符4.变量名不区分大小写(包括LINDO中的关键字)5.目标函数所在行是第一行,第二行起为约束条件6.行号(行名)自动产生或人为定义。行名以“)”结 束7.行中注有!”符号的后面部分为注释。如:!Its Comme
13、nt.8.在模型的任何地方都可以用“TITLE”对模型命名(最多72个字符),如:数学建模使用LINDO的一些注意事项y 荐量不附tn也仕一 t 约 尿余什 则力顺1最达式中不接受括号“()”和逗号“等任何符号,例:400(X1+X2)需写为400X1+400X211.表达式应化简,如2X1+3X2-4X1应写成2X1+3X212.缺省假定所有变量非负;可在模型的“END”语句 后用“FREE name”将变量name的非负假定取消13,可在“END”后用“SUB”或“SLB”设定变量上 下界例如:“sub x1 1 0”的作用等价于“x1 v=1 0”但用“SUB”和“SLB”表示的上下界约
14、束不计入 模型的约束,也不能给出其松紧判断和敏感性分 机 814-12-10 数学建模i 二次规划(QP)问题4后。可求解二次规划(QP)问题,但输入方式 较复杂,因为在LINDO中不许出现非线性表达式 需要为每一个实际约束增加一个对偶变量(LAGRANGE乘子),在实际约束前增加有关变 量的一阶最优条件,转化为互补问题“END”后面使用QCP命令指明实际约束开始的 行号,然后才能求解 建议总是用UNGO解QP注意对QP和I P:敏感性分析意义不大2014-12-10数学建模状态窗口(LINDO Solver StatusLINDO Solver StatusX 当前状态:已达最优解Optim
15、izer StatusS tatus:Iterations:Infeasibility0 blec We;Besl IP:IP B oundB ranches:Elapsed Tme:Optimal180949493.5100:00:00llpdate Interval:迭代次数:18次 约束不满足的“量”(不 是“约束个数”):0 当前的目标值:94 最好的整数解:94 整数规划的界:93.5 分枝数:1 所用时间:0.00秒(太 快了,还不到0.005秒)刷新本界面的间隔:1(秒)数学建模2t 不二区/Options.Nonzero Limit:非零系数的个数上限;Iteration Li
16、mit:最大迭代步数;Initial Contraint Tol:约束的初始误差上限;Final Contraint Tol:约束的最后误差上限;Entering Var Tol:进基变量的REDUCED COST的误差限;Pivot Size Tol:旋转元的误差限I_ Preprocess:预处理(生成割平面);Preferred Branch:在先的分枝方式:“Default”(缺省方式)、“Up”(向上取整优先)、“Down”(向下取整优先);IP Optimality Tol:IP最优值允许的误 差上限(一个百分数,如5%叫005);IP Obj ective Hurd le:IP
17、目标函数的篱 笆值,即只寻找比这个值更优最优解(如当知道当前模型的某个整数可行解 时,就可以设置这个值);IP Var Fixing Tol:固定一个整数变量 取值所依据的一个上限(如果一个整数 变量的判别数(REDUCED COST)的 值很大,超过该上限,则以后求解中把 该整数变量固定下来)。2014-12-10数学建模Rannrt/QtatiaticaROWS=5 VARS=4 INTEGER VARS=2(0=0/1)QCP=4NONZEROS=19 CONSTRAINT NONZ=12(6=+-1)DENSITY=0.760SMALLEST AND LARGEST ELEMENTS
18、IN ABSOLUTE VALUE=0.300000 277.000OBJ=MIN,NO.:2 0 2,GUBS=0SINGLE COLS=0 REDUNDANT COLS=0第一行:模型有5行(约束4行),4个变量,两个整数变量(没有 0-1变量),从第4行开始是二次规划的实际约束。第二行:非零系数19个,约束中非零系数12个(其中6个为1或-1),模型密度为0.760(密度=非零系数/行数*(变量数+l)o第三行的意思:按绝对值看,系数最小、最大分别为0.3和277。第四行的意思:模型目标为极小体;a工等于、等于一大壬等于约 束分别看2、0、2不;广义上界约束(GUBS)不急过1个;变量上
19、界约束(VUBS)不少于0个。所谓GUBS,是指一组不 含有相同变量的约束;所谓VUBS,是指一个蕴涵变量上界的约 束,如从约束Xl+X2-X3=0可以看出,若X3=0,贝!|X1=0,X2=0(因为有非负限制),因此Xl+X2-X3=0是一个VUBS约束。第五行的意思:只含1个变量的约束个数=0个;冗余的列数=0个2014-12-10数学建模(LINDO行命令、命令脚本文件WINDOWS环境下行命令的意义不大批处理:可以采用命令脚本(行命令序列)Mil kO2JpkExample 演示BatO1.txt IffiFILE/TAKE COMMANDS(F11)命令调入 2014-12-10Mi
20、l k03.l pkSAVE行命令 FILE/SAVE命令必须是以LINDO PACKED形式(压缩)保存的文件数学建模 LINGO软件简介/NGO模型的优点方包含了 LINDO的全部功能提供了灵活的编程语言(矩阵生成器)LINGO模型的构成:5个段 目标与约束段 集合段(SETS ENDSETS)数据段(DATA ENDDATA)初始段(INIT ENDINIT)计算段(CALC ENDCALC)-LINGO9.0乙UI4T4TU 奴子廷吴LINGO模型一例:选址问口七某公司有6个建筑工地,位置坐标为Q)(单位:公里),水泥日用量4(单位:吨)123456假设:料场a1.258.750.55
21、.7537.25b1.250.754.7556.57.75和工地之间d3547611有直线道路1)现有2料场,位于A(5 J),B(2,7),记(Xj,yj)j=l,2,日储量R各有20吨。目标:制定每天的供应计划,即从A,B两料场分别 向各工地运送多少吨水泥,使总的吨公里数最小。2014-12-10数学建模2 6l、i 2 2 1/2决策变量:.mm ZZ引a厂田)+(y厂与)(料场j到工症的 2=1”运量)12维 s.t.Z Cy=dj,i=线性规划模型用例中数 据计算,最优解为j=l6Z c,j eji=l1123456%(料场A)350701a(料场B)0040610数学建模总吨公里数
22、为13622014-12-10选址问题:NLP2)改建两个新料场,需要确定新料场位置(琴”)和 运量.,在其它条件不变下使总吨公里数最小。2 6min V V Cy(xy-6ZZ-)2+(j/y-%)21/2 决策变量:勺占 为,(与”)76维2s.t.7 g=dj,i=1.6jJ 非线性规划模型片16 Z c产泞 j=2i=l2014-12-10 数学建模IT LINGO Model-locationMODEL:Tit l e Lo ca t io n Pr o b l em;se t s:d ema n d/1.6/:a,b,d;sup p l y/1.2/:xzyze;l in k(d
23、ema n d,sup p l y):c;en d set s d a t a:LINGO模型的构成:4个段集合段(SETS ENDSETS)!l o ca t io n s f o r t he d em&n d(需求点的位置);a=1.25,8.75,0.5,5.75,3,7.25;b=l.25,0.75,4.75,5,6.5,7.75;!qua n t it ies o f t he d ema n d a n d sup p l y(供需量);*3,5,4,7,6,11;e=20,20;en d d a t a数据段(DATA ENDDATA)in it:LP:移到数据段初始段(INI
24、T ENDINIT)en d in it!Ob ject ive f un ct io n(目标);OBJ min=sum(l in k(izj):c(i,j)*(x(j)-a(i)A2+(y(j)-b(i)A2)A(1/2);!d ema n d co n st r a in t s(需求约束);o r(d ema n d(i):DEMAND_CON sum(sup p l y(j):c(izj)=d(i););!sup p l y co n st r a in t s(供应约束);o r(sup p l y(i):SUPPLY_CON sum(d ema n d(j)=e(i););o r
25、(sup p l y:f r ee(X);Of r ee(Y););局部最优:89.8835(吨公里)目标与 约束段二12014-12-10数学建模I 集合的类型那同型邓 郡ist)setname/member list/memberjist/7):attributejist?SETS:CITIES/A1,A2,A3,B1,B2/;ROADS(CITIES,CITIES)/Al,Bl A1,B2 A2.B1 A3,B2/:D;ENDSETSSETS:STUDENTS/S1.S8/;PAIRS(STUDENTS,STUDENTS)|&2#GT#&1:BENEFIT,MATCH;ENDSETS集合
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 讲座 案例 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【曲****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【曲****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。