1、论文原创性声明本人郑重声明: 所呈交的论文,是本人在老师的指导下,独立进行研究工作所取得的成果。除文中已经注明引用的内容外,本论文不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。数学论文作者签名:日期:2011 年 1 月 1 日I模糊规划摘 要本文应用了一种求解多目标模糊线性规划题的解法就是应用单纯形法先求解相应的普通线性规划问题,然后通过模糊约束集与模糊目标集的隶属度比较,获得两个集合交集的最优隶属度,将此最优隶属度代入最优单纯形表中,即可求得到相应多目标模糊线性规划的解。用此方法,
2、可以很好的解决本文中水资源在农业规划中的应用。并且在文中对不同的解模糊规划的方法进行了比较。关键词:单纯形法;多目标模糊线性规划;隶属度1目录模糊规划 I摘 要 I前 言 11 背景知识 22 模糊规划 122.1 问题提法 22.2 模糊极值 32.2.1 有界函数 f的极值和模糊极值 32.2.2 普通限制下,目标函数 f的极值与模糊极值 32.2.3 模糊限制下,目标函数 的模糊极值 32.3 模糊线性规划 42.4 多目标线性规划与模糊线性规划 43 实际问题求解 53.1 实际 问题 63.2 解题步骤 64 各种方法比较 7参考文献 7附 录 9致 谢 .101前 言中国是世界人口
3、大国,用水量也居世界前列。虽然中国有长江与黄河,但中国仍然是一个重度取缺水的国家,例如中国西北部的部分城市,长年降水量不足1厘米,水井经常处于缺水的状态。我国水资源总量不少, 但人均占有量不足世界平均值的1/ 4, 且时空分配不均。总的来说南多北少, 而开发利用程度南小北大。为了缓解这种状况,有了南水北调工程,但这并没有完全解决中国缺水这一现状。随着人口增长、经济发展, 将有更多的城市和地区严重缺水, 因此, 水资源短缺对社会、经济活动制约的形势日趋严峻, 尤其是解决华北、西北和沿海一些城市供水不足的问题, 已成为政府亟待解决的一个战略问题。改革开放20 年间, 我国农业得到长足发展, 长期困
4、扰我国的粮食紧缺状况, 由于北方粮食生产充分发展, 使历史形成的“南粮北调”局面得到扭转, 成功地以占世界7% 的耕地面积, 养活了世界22%的人口, 举世瞩目, 堪称典范。北方旱区为了发展生产, 克服水障碍, 大力发展水利, 拦截径流, 开发地下水, 强化水、热资源的调控力度,有力地促进了北方粮食生产的发展, 这是当今“北粮南运”的基础所在。然而这种粮食生产格局的变化, 是以过量地消耗本已短缺的北方地表水和地下水资源为代价的。可以说我国农业在取得长足发展的同时,也在农业生产上浪费了大量的水资源。只有解决了水资源的浪费问题,才能在根本上有所缓解中国缺水的情形。本文就是在保证农民最大利益的基础上
5、解决如何合理的利用水资源(本文中合理用水不是减少用水,是必要的用水;只要不是浪费水,而能满足用水量,才叫合理用水)来合理安排水资源在农业生产方面的应用。农业生产中合理利用水资源的方式有很多种,例如将灌溉改为滴灌、减少植物因蒸发而失去的水量。但由于水平有限,本文只考虑在已知农作物在不同方式下用水量以及利益问题来改变农作物的种植量从而达到节约水资源的利用。由于农作物使用的水资源是一个模糊量,即只有一个大概的范围,因此为了解决本文相应的问题,我们采用模糊规划的知识进行求解。21 背景知识数学规划问题可归结为给定一个目标函数 与某种限制条件,寻求在接受限制条件f下尽可能取得最优的目标。如果限制条件可用
6、经典的集合 表示,就是普通规划问题,A如果限制条件是不确定的,需要用模糊集 来表示,就是模糊规划问题。这种模糊规划,A就是在模糊集 的限制条件下来求目标函数 的模糊极值。A f2 模糊规划 12.1 问题提法普通规划问题:给定目标函数 : , 为有界函数,给定 X 的子集 作为f()XYRf A限制集,要求在 中选择元素 ,使得 取得极大值,这是一个要求在 上的条件极A0x0x值问题。经常遇到的限制集是模糊的。模糊规划问题: 给定目标函数 : , 为有界函数,给定 X 的模糊子集f()XYRf作为限制集,要求选择元素 ,使得 对 的隶属度及 对于目标函数值 都尽可A 0x0A0x0()fx能地
7、达到高水平。这是一个在 的限制下的条件模糊极值问题。32.2 模糊极值2.2.1 有界函数 的极值和模糊极值f定义:设 : ( 有界),令 (表示模糊函数集合) ,fXYRfiMX9,()nf()supiffxx称 为 的无条件模糊优越集,并称 为 的无条件模糊极大值。fM ()fYF9f,()f ffxyM当 = 时, =1;当 = 时, =0;当1()fxma()Xfx1fx2()min()X2fMx时, 。因此 反映了在模糊的意义下, 的优越程度。2f1f2ffx又当 = 时, ;当 = 时,1yax()Xf11()f ffxyM2yin()xXf;当 时, ;因此, 反映22()0f
8、ffxyM()X()0TMTMy了在模糊意义下, 对 的模糊极大值的隶属程度。2.2.2 普通限制下,目标函数 的极值与模糊极值f定义:设 : ( 有界函数),令 。令 ,称 为 在 上fXYRf AX9ffAMfA(条件)模糊极大值。 ()()()f f fffxyfxyAx= () ,ffM()in()supffxfx反映了接受 的限制,又反映了在在整个目标函数中所处地位。()fAyA42.2.3 模糊限制下,目标函数 的模糊极值f定义:设目标函数 : ( 有界函数)对模糊限制 。令 ,称fXYRf AX9ffAM为 在 上(条件)模糊优越集,称 为 在 上(条件)模糊极大值,这里fA (
9、)fA22()()in,sup()()(), ().ff ff f fffxyfxyfxMAAMx 由于 说明当 时,既要求 对限制集 的隶属程度达到或()()f f f A超过 水平。因此 既反映 接受 的限制的程度,又反映 达到理想目标的程度。fAx ()fx关于目标函数 的模糊极小值问题,可以转化为 的模糊极大值问题。f f2.3 模糊线性规划在经典规划论中,线性规划师最基本的而且运用较广泛的一个分支,其数学模型是目标函数: 12,nScxcx2112210,0,nmmnaaxbxx求一组变量 的值满足约束条件并且使得目标函数 取得最小值(或最大)值。2n S2.4 多目标线性规划与模糊
10、线性规划设有一个多目标线性规划问题5(2.4.1)11212212max0,0,nrrrnzcxcxz (2.4.2)1212120,0,nmmnxaxbaxx约 束 :记 ,()ijmnAa, , , ,ijrnCc12,)Tmbb 12(,)Tmxx 12(,)Tmzz则上述问题可简写为 a,0.zCAxb由于目标函数不止一个,要想同时使得这些目标函数都达到极大值往往是不可能的,这就需要提出某种折衷方案,使各目标都相对的“极大” ,这就需要要用模糊数学的方法进行处理。设 是(2.4.1) 中目标函数 在约束条件(2.4.2)下所得到的最大值*iz1niijzcx.对每个目标,根据其重要性给
11、出伸缩指标 , 愈小说明目标 俞(1,2)ir id(0)i iz重要,相应构造一个模糊目标集 ,其隶属函数为fM*1* *1 110,()(), nijiiin nniiij iijiiijii iijiicxzdMxgcxzcxzddcxz 记 ,称为多目标线性规划问题的综合模糊目标,记12rM为满足约束方程的可行解集(它是 的普通子集) 。,0DxAb nR63 实际问题求解3.1 实际问题某地区农民王某在收获了玉米后准备开始冬小麦的种植。现在市场上冬小麦的品种种类繁多,其中有“京冬 8 号”、 “京旺十号”、 “农大 3214”、 “北农 66”、 “京 411”、 “豫教 5 号”、
12、“周麦 18 号 “等,每种小麦都有自己的优势与劣势,相应的各种小麦的价格也不尽相同。考虑到自己买小麦种子所带的资金,该农民最终选择了优质品种农大 135 和节水品种农大 3488。该农民选择的农大 135 的特征特性是:冬性,长芒,红粒。幼苗健壮,叶片宽厚。株高 85 厘米左右,抗倒性较强。抗寒、抗病性较好,成熟期与京 411 相近。产量要素比较协调,亩穗数 3540 万,穗粒数 2830 粒,千粒重 5055 克,产量水平在 450 公斤左右。品质较优,做饺子透明度好,有咬劲,不混汤;做面条色泽白,有咬劲。而农大 3488 的特征特性是:冬性,中熟,成熟比京冬 8 号晚 12 天。幼苗半匍
13、匐,分蘖力强,成穗率较高。株高 80 厘米左右,长芒,白粒,穗粒数 30 粒,千粒重 42 克。茎杆粗壮,抗倒伏。抗寒性中等。中感条锈病、高感叶锈和白粉病。耐旱性强,抗干热风。丰产性好,节水栽培下亩产 350 公斤左右。现在已知该农民家中拥有 12 亩土地,他准备全部用来种植小麦。当地政府为了响应国家的号召,节约用水,于是对部分土地使用了滴灌式灌溉来节约用水。该农民王某恰有两亩土地采用了这种滴管式灌水。在这种土地上,浇灌小麦所需的用水量大量减少。王某根据以往经验,熟知了这两种小麦的用水情况及收益,大致如下 :2土地类型小麦品种 全生长期浇水量(百方/亩)开花期浇水量(百方/亩)小麦单产(吨/亩
14、)小麦净产值(元/亩)农大 135 9.7 1.65 0.45 800普通土地 农大 3488 7.23 1.4 0.35 580农大 135 7.55 0 0.44 685滴灌式土地 农大 3488 5.89 0 0.34 5757其中滴灌式土地由于前期投入的资金较多,导致小麦净产值有所偏低,但一段时期之后,成本得到回收,小麦净产值会增加(净产值=小麦单产(吨/ 亩)X 小麦价格(元/ 吨) )。当地政府为了节约用水,限制了当地农民的用水总量,即全生长时期平均每亩用水8.4 百方,开花期每亩用水 1.3 百方。农民王某现在对种植的小麦有如下打算:小麦产量最多,同时盈利最多。针对上述情形,该农
15、民应该怎么安排种植情形才能达到他的要求?3.2 解题步骤将上述实际问题转化为模糊规划问题如下:用 表示农大 135 与农大 3488 的总产量,用 表示该农民种植小麦的总收益,则转1z 2z化为如下方程:设: 表示普通土地种植的农大 135 的亩数1x表示普通土地种植的农大 3488 的亩数2表示滴灌式土地种植的农大 135 的亩数3x表示滴灌式土地种植的农大 3488 的亩数4根据原问题列出规划方程:(1) 土地限制下的约束: 12340,x(2) 水资源限制下的约束: 12349.7.7.5.8910.,6546,xx最后根据问题构造目标函数(1) 小麦产量最大的目标函数: 112340.
16、45+.0.+.,zxx(2) 小麦最大收益的目标函数: 2123486857,z8综上所述,就可以得到这一问题的多目标模糊线性规划:max112340.45+.0.+.,zxx286857约束:(3.2.1)123412342340,.65.5.6,977.8910.,0,xxx对目标 采用单纯形法 :1z234xx1234aa0.450.00.4530.401 112.897.235.898.7967.55640056400T .310.171673.9541 221005.43.67293.54.02.06 6. .0.4.1.8.4.416.35773.954710.2 25.43.6
17、729354.3026 .06. .90.48.710215.43.672.93.5460 所以最优解, , , , ,最大值 ,此时 .15.4x2.3x40*14.92z28364z对目标 ,利用单纯形法:2z9最优解 , , , ,最大值 ,此时 .19.8x20x342x*28491.z148z同时考虑两个目标,合理的方案应该使 我们取伸缩指标 ,1.l217.lmaxS约束123412340.6515.609778921.8049250.18767xxx,12340,0,xx利用单纯形法得最优解 *16.3,.5,.6,1.,0.76此时 ,489z24z4 各种方法比较上一节中解线
18、性规划用的是单纯形法,但运算过于繁琐,一步错,全部错。单纯形法解普通线性规划问题要求极高的耐心与极高的细心,而且还要把握好接的过程。可以说,单纯形法一种计算量大、步骤繁琐的方法。现在我们对该种方法进行改进,改用相应的软件进行求解。4.1 Lingo 求解法 相应程序见附录 1对线性规划方程(3.2.1) ,运用 lingo 软件可得最优解: , , ,15.43x2.732x,最大值 ,此时 。40x*14.92z28364.z10对线性规划方程(3.2.2) ,运用 lingo 软件可得最优解: , , , 19.8x20x3,最大值 ,此时 。42x*28491.z14.8z类似单纯线性方
19、法可得方程(3.2.3)对线性规划方程(3.2.3) ,运用 lingo 软件可得最优解:,此时 , 。*12346.3,.5,0.6,1.,0.76xx14.89z251.4z参考文献请按下列例子的格式输入“专著、论文集、学位论文、报告”:序号 主要责任者. 文献题名文献类型标识. 出版地:出版者,出版年. 起止页码(任选).1. 罗承忠. 模糊集引论.北京:北京师范大学出版社,2005 年 8 月.2. 国家发改委. 2010 小麦最低收购价执行预案.中国期刊全文数据库 2010.3. 张筑生. 微分半动力系统的不变集D. 北京:北京大学数学系数学研究所, 1983.4. 冯西桥. 核反应
20、堆压力管道与压力容器的 LBB 分析R. 北京:清华大学核能技术设计研究院,1997.请按下列例子的格式输入“期刊文章”:序号 主要责任者. 文献题名J. 刊名,年,卷(期):起止页码.5. 何龄修. 读顾城南明史 J. 中国史研究,1998,(3):167-173.6. 金显贺,王昌长,王忠东,等. 一种用于在线检测局部放电的数字滤波技术J. 清华大学学报(自然科学版),1993,33(4):62-67.请按下列例子的格式输入“论文集中的析出文献”:序号 析出文献主要责任者. 析出文献题名A. 原文献主要责任者(任选) . 原文献题名C. 出版地:出版者,出版年. 析出文献起止页码.7. 钟
21、文发. 非线性规划在可燃毒物配置中的应用A. 赵玮. 运筹学的理论与应用中国运筹学会第五届大会论文集C. 西安:西安电子科技大学出版社,1996. 468-471.请按下列例子的格式输入“报纸文章”:序号 主要责任者. 文献题名N. 报纸名,出版日期 (版次).8. 谢希德. 创造学习的新思路N. 人民日报,1998-12-25(10).11请按下列例子的格式输入“国际、国家标准”:序号 标准编号,标准名称S.9. GB/T 16159-1996,汉语拼音正词法基本规则S.请按下列例子的格式输入“专利”:序号专利所有者. 专利题名P. 专利国别:专利号,出版日期.10. 姜锡洲. 一种温热外敷
22、药制备方案 P. 中国专利:881056073,1989-07-26.请按下列例子的格式输入“电子文献”:序号 主要责任者. 电子文献题名电子文献及载体类型标识. 电子文献的出处或可获得地址,发表或更新日期/引用日期(任选).11. 王明亮. 关于中国学术期刊标准化数据库系统工程的进展EB/OL. http:/ 万锦坤. 中国大学学报论文文摘 (1983-1993). 英文版DB/CD. 北京:中国大百科全书出版社,1996.请按下列例子的格式输入“各种未定义类型的文献”:序号 主要责任者. 文献题名Z.出版地:出版者,出版年.12附 录附录 1Lingo 程序:(1) 的求法zModel:m
23、ax=0.45*x1+0.35*x2+0.44*x3+0.34*x4;x1=0;x2=0;x3=0;x4=0;x1+x2=0;x2=0;x3=0;x4=0;x1+x2=0;x2=0;x3=0;x4=0;y=0;0.45*x1+0.35*x2+0.44*x3+0.34*x4-0.112*y=4.81;800*x1+580*x2+685*x3+575*x4-127.26*y=8364.6;x1+x2=10;x3+x4=2;9.7*x1+7.23*x2+7.55*x3+5.89*x4+2*y=102.8;1.65*x1+1.4*x2+0.1*y=15.7;y=1;end运行结果Global opti
24、mal solution found at iteration: 14Objective value: 0.758022114Variable Value Reduced CostY 0.7580221 0.000000X1 6.630214 0.000000X2 3.345960 0.000000X3 0.6021600 0.000000X4 1.397840 0.000000Row Slack or Surplus Dual Price1 0.7580221 1.0000002 6.630214 0.0000003 3.345960 0.0000004 0.6021600 0.000000
25、5 1.397840 0.0000006 0.7580221 0.0000007 0.000000 -2.1738948 0.000000 -0.1902234E-029 0.2382551E-01 0.00000010 0.000000 0.319124711 0.000000 0.257009112 0.000000 0.4273311E-0213 0.2419779 0.000000附录 215致 谢由衷感谢 老师,这一学期以来对我们的悉心的教导。在他风趣,负责,循循善诱的教学课堂中,使我们领略到了模糊数学的魅力。 老师上课的风格深深的感染了我们,老师上课所举的风趣幽默的例子,常常在让我们开怀大笑的同时也理解了所教授的概念与知识,使我们感受到数学的魅力不仅仅局限于追求精确的解,模糊的结果也是数学的另一种美。这种别样的数学思维方式,对我们的学习和社会生活都产生了深刻的影响,凡事无绝对,模糊的数学,很美!作者名 2011 年 1 月