決策支持系統(tǒng)作業(yè)(最終)

1、決策支持系統(tǒng)導(dǎo)論期末作業(yè)姓名：齊鵬 學(xué)號(hào)： 日期：2012年7月9日1、設(shè)某企業(yè)生產(chǎn)多種最終產(chǎn)品Y=（yij）,各種產(chǎn)品的單價(jià)為 Pi,它們的投入產(chǎn)出直接消耗系數(shù)為A=（aij）,企業(yè)的資源（煤、電力、勞力）的約束方程為BX<=>h （"<=>"表布<、=、>）, 其中，B=（bij）是資源消耗系數(shù)矩陣，X=（xi）是企業(yè)總產(chǎn)品向量，h是資源約束向量。為使企業(yè)凈產(chǎn)值最大，其目標(biāo)方程S=EPiyi -max,試安排生產(chǎn)計(jì)劃（求總產(chǎn)品X和最終產(chǎn)品Y）。 請(qǐng)?jiān)O(shè)計(jì)該企業(yè)的生產(chǎn)計(jì)劃決策支持系統(tǒng)，畫(huà)出DSS運(yùn)行結(jié)構(gòu)圖，并對(duì)總控程序、模型程序、數(shù)據(jù)庫(kù)進(jìn)

2、行結(jié)構(gòu)和功能說(shuō)明。提示：該決策支持系統(tǒng)需要利用3個(gè)模型（投入產(chǎn)出模型、線性規(guī)劃模型和報(bào)表模型（打印投入產(chǎn)出表）和兩個(gè)數(shù)據(jù)庫(kù)（投入產(chǎn)出數(shù)據(jù)庫(kù)和線性規(guī)劃數(shù)據(jù)庫(kù)）。在DSS總控程序中要詳細(xì)說(shuō)明何時(shí)調(diào)用哪個(gè)模型運(yùn)行，何時(shí)存取哪個(gè)數(shù)據(jù)庫(kù)中的數(shù)據(jù)，何時(shí)進(jìn)行數(shù)據(jù)計(jì)算。該DSS需要兩次調(diào)用投入產(chǎn)出模型：一次計(jì)算中間結(jié)果，一次計(jì)算最后結(jié)果。請(qǐng)注意，模型程序應(yīng)該是一個(gè)標(biāo)準(zhǔn)程序，在一定的參數(shù)控制下，可得到中間結(jié)果，也可得到最終結(jié)果。該模型程序既適合于該問(wèn)題的 DSS也適合于其他問(wèn)題的 DSS,不能是一個(gè)專用的 模型程序。（40分）一、模型1 .投入產(chǎn)出模型：可以確定W與七的關(guān)系，* =（三一5）匕,2 .線性規(guī)劃模

3、型：卞據(jù)約束方程BX<=>h與目標(biāo)方程 S=!2 Piyi -max可以計(jì)算出最優(yōu)的yi。3 .報(bào)表模型：根據(jù)最終產(chǎn)品 Y=（yij） , X=（xi）是企業(yè)總產(chǎn)品自動(dòng)生成報(bào)表。 二、數(shù)據(jù)庫(kù)1.投入產(chǎn)出數(shù)據(jù)庫(kù)字段名數(shù)據(jù)類型長(zhǎng)度是否可為空int8是int8是int8是int8是int8是int16是int16是int16是int16是int16是float16是float16是float16是float16是float16是2.線性規(guī)劃數(shù)據(jù)庫(kù)字段名數(shù)據(jù)類型長(zhǎng)度是否可為空f(shuō)loat8一 是float8是float8是float8是float8是h1float16是h2float16是h

4、3float16是float16是hifloat16是float16是float16是float16是float16是float16是三、DSS運(yùn)行結(jié)構(gòu)圖2、卡拉 OK廳唱歌的時(shí)候，是否要等待包間的問(wèn)題。域內(nèi)的實(shí)例:規(guī)I定如下屬性可用于描述該翎I(1)ers （其他地點(diǎn)）：附近是否有其他卡拉 OK廳口（等候條件）：供顧客等候的地方是否舒適(3) Weekend (周末): 若是周六或周日，則為真;（4） Conssumers （顧客）：店中有多少顧客（值為 Nor 座）；F ul調(diào)用線性規(guī)劃數(shù)據(jù)庫(kù):（5） Price （價(jià)格）：價(jià)格范圍（值為 Cheep （便宜），Middle士EXPenB=

5、（bij）星資源消耗系數(shù)（6） Raining （下雨）：外面是否在下雨;（7） Reservation （預(yù)約）：是否預(yù)約過(guò)；nsive （較貴）；WitEstimate （等候時(shí)間估計(jì)）:估計(jì)的等候時(shí)間（獻(xiàn)懶藥蛾型：根據(jù) X=(xi) , Y=(yij)編制 據(jù)0-10, 10-30, 3060>60,調(diào)四S融規(guī)劃模型屬性WCondWEnd單位為分鐘）。訓(xùn)練集見(jiàn)表：No0 4rJNoPr OthersX1Yes1顯示報(bào)表-ConsPriceSome0EX-1NoO調(diào)用線性規(guī)弧模型： 計(jì)算出最終產(chǎn)“Y=（yij）他一標(biāo)RW:解。Res : WEst WYeslT 0-10 1投入產(chǎn)出數(shù)

6、據(jù)庫(kù)： 將計(jì)算出的 X=（xi）, Y=（yij）填入表中。AllWaitYesX2Yes1No0No0Full11CH1No0No030-60-1NoX3No0Yes1No0Some0CH1No0No00-10 1YesX4Yes1No0Yes1Full1CH1Yes1No010-300YesX5Yes1No0Yes1Full1EX-1No0Yes1>60 -2NoX6No0Yes1No0Some0MID0Yes1Yes10-10 1YesX7No0Yes1No0None-1CH1Yes1No00-10 1NoX8No0No0No0Some0MID0Yes1Yes10-10 1YesX

7、9No0Yes1Yes1Full1CH1Yes1No0>60 -2NoX10Yes1Yes1Yes1Full1EX-1No0Yes110-300NoX11No0No0No0None-1CH1No0No00-10 1NoX12Yes1Yes1Yes1Full1CH1No0No030-60-1Yes要求：建立BP神經(jīng)網(wǎng)絡(luò)模型，并進(jìn)行容錯(cuò)性分析。（30分）、為訓(xùn)練集賦值屬性值名屬T1值屬性2值屬性3值屬性4值OthersAYes1No0WCondBYes1No0WEndCYes1No0ConsDFull1Some0None-1PriceECheap1Middle0Expensive-1Rain

8、FYes1No0ResGYes1No0WEstH0-10110-30030-60-1>60-2WillWaitIYes1No0在Matlab的命令窗口中運(yùn)行如下代碼：p=1 1 0 1 1 0 0 0 0 1 0 1;0 0 1 0 0 1 1 0 1 1 0 1;0 0 0 1 1 0 0 0 1 1 0 1;0 1 0 1 1 0 -1 0 1 1 -1 1;-1 1 1 1 -1 0 1 0 1 -1 1 1;0 0 0 1 0 1 1 1 1 0 0 0;1 0 0 0 1 1 0 1 0 1 0 0;1 -1 1 0 -2 1 1 1 -2 0 1 -1;t=1 0 1 1 0

9、 1 0 1 0 0 0 1;net=newff(minmax(p),15 1,'tansig','purelin','traincgb');二5;=300;=0.000001;net,tr二train(net,p,t);得到如下曲線圖: 容錯(cuò)性分析：實(shí)例輸入輸 出WillWait結(jié)果OthersWCondWEndConsPriceRainResWEstX11000-10111等X21001100-10不等X3010010011等X4101111001等X51011-101-20不等X6010001111等X7010-111010不等X80000

10、01111等X90111110-20不等X101111-10100不等X11000-110010不等X121111100-11等完成機(jī)器學(xué)習(xí)后，對(duì)樣本進(jìn)行改變條件輸入，有如下3種情況:（1）缺1個(gè)條件的情況（2）缺2個(gè)條件的情況（3）介于中間的情況（1）實(shí)例輸入輸出WillWait(1)輸出WillWait(0)結(jié)果OthsWConWEnConPriceRaiResWEstX11100-10110.84560.1526等（1）X21011100-10.07230.9281不等（0）X3010110010.86480.1352等（1）X4100111000.86340.1366等（1）X5111

11、1-101-20.23140.7686不等（0）X6011000110.91250.0875等（1）X7011-101010.34100.6595不等（0）X8011001110.97140.0286等（1）X90000110-20.49280.5012不止X10110010100.56810.4319不止X11111-110010.41230.5877不止X121100000-10.64280.3572不止3、編制旅行商路徑優(yōu)化問(wèn)題的遺傳算法程序，并計(jì)算一個(gè)實(shí)例。（30分）一，問(wèn)題描述：所謂旅行商問(wèn)題，即給定幾個(gè)城市，旅行商從中選擇一 條最短的路線，使他能夠訪問(wèn)到每個(gè)城市一次，然后返回起點(diǎn)O

12、,運(yùn)行結(jié)果:o 00 oooo三、代碼實(shí)現(xiàn)下面是程序的主要核心代碼：using System;using ;using System.ComponentModel;using System.Data;using System.Drawing;using System.Text;using ;using System.Diagnostics;namespace TSPGApublic partial classForml : Formprivate World _world = new World ();public Form1() InitializeComponent();private v

13、oid pictureBox1_MouseUp( object sender, MouseEventArgs e) Point point = new Point(e.X / (pictureBox1.Width / 200), e.Y / (pictureBox1.Height / 200);if (e.Button = MouseButtons .Left)_world.AddPointNeedToVisit(point);cityNum.Text = _(); else _world.StartPoint = point;ReDrawPB();private void Reset_Cli

14、ck( object sender, EventArgs e)_world.Reset();ReDrawPB();private void ReDrawPB()Bitmap bm = new Bitmap (pictureBox1.Width, pictureBox1.Height);using (Graphics g = Graphics.FromImage(bm) int x, y, radix;foreach (Point point in _world.PointsNeedToVisit) x = point.X * (pictureBox1.Width / 200);y = poin

15、t.Y * (pictureBox1.Height / 200);radix = 5;g.DrawEllipse( Pens.DarkCyan, new Rectangle(x - radix, y - radix, 2 * radix, 2 * radix);if (_world.StartPoint != Point.Empty) x = _world.StartPoint.X * (pictureBox1.Width / 200);y = _world.StartPoint.Y * (pictureBox1.Height / 200);radix = 5;g.DrawEllipse( P

16、ens.Red, new Rectangle(x - radix, y - radix, 2 * radix, 2* radix);pictureBox1.Image = bm;private void pictureBox1_MouseMove( object sender, MouseEventArgs e) Point point = new Point(e.X / (pictureBox1.Width / 200), e.Y / (pictureBox1.Height/ 200);infos1.Text = point.ToString();private void Start_Cli

17、ck( object sender, EventArgs e)_world.Population = Int32.Parse(population.Text);_world.Generation = Int32.Parse(generation.Text);_world.CrossOverRate = Double.Parse(crossOverRate.Text);_world.MutationRate = Double .Parse(mutationRate.Text);drawVisitRoute();/_world.Reset();private void drawVisitRoute

18、()ReDrawPB();Bitmap bm = (Bitmap ();using (Graphics g = Graphics.FromImage(bm)Point startPoint = _world.StartPoint;if (startPoint = Point.Empty)MessageBox.Show("You must define a start point." ); return;Point pointsNeedToVisit = _world.PointsNeedToVisit;Stopwatch sw = new Stopwatch();sw.St

19、art();int visitRoute = _world.GetVisitRoute();sw.Stop();infos3.Text = "Used: " + sw.ElapsedMilliseconds + "ms." ;List<Point> visitPoints = new List<Point>();visitPoints.Add(startPoint);foreach (int index in visitRoute)visitPoints.Add(pointsNeedToVisitindex);int td = 0;Point tp = visitPoints0;for (int i = 1; i < visitPoints.Count; i+)td += GetManhattanDis(tp, visitPointsi);tp = visitPointsi;totalDis.Text = "Total Distace: " + td;for (int i = 0; i