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

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

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

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

4、oat16是h3float16是float16是hifloat16是float16是float16是float16是float16是float16是三、DSS運行結(jié)構(gòu)圖投入產(chǎn)出數(shù)據(jù)庫:開始jO潘幽為直接消耗系數(shù)為皿 調(diào)用投入產(chǎn)出模一-咽用掾慮去串模納KT唱歌的時候，是否要等待包間的問題。規(guī)定婦下屬性可調(diào)m拱甦醇窿塾內(nèi)的實例:調(diào)用投入產(chǎn)出模型線性規(guī)劃數(shù)據(jù)庫廠(1) Others（其他地點）：附近是否有其惜榭 瞅麴瀾潮 加（2） WaitCond （等候條件）：供顧客等候用地方是否舒適;調(diào)備31讖ekend（周末）：若是周六或周日，口陰嵬報表模型: 網(wǎng)為真；調(diào)用;線ticonssumers （顧客

5、）：店中有多少顧得些）或Full （滿座）;（值為 None （僅入型:（5） Price （價格）小格范圍（值為Cheep便宜）上M由dle （中等） 投入廣出數(shù)據(jù)庫：,Some(,Expensive將計算出的(6)Raining （下雨）：外面是否在下雨;（7） Reservation （預(yù)約）：是否預(yù)約過；（8） WaitEstimate （等候時間估計）：估計的等候時間（值為 010, 1030, 3060, >60,單位為分鐘）訓(xùn)練集見表：實例屬性目標(biāo)WillW aitOthe rsWCondWEndConsPric eRainResWEstX1Yes1No0No0Some0E

6、X-1No0Yes10-101YesX2Yes1No0No0Full11CH1No0No030-60 -1NoX3No0Yes1No0Some0CH1No0No00-101YesX4Yes1No0Yes1Full1CH1Yes1No010-30 0YesX5Yes1No0Yes1Full1EX-1No0Yes1>60-2NoX6No0Yes1No0Some0MID0Yes1Yes10-101YesX7No0Yes1No0None-1CH1Yes1No00-101NoX8No0No0No0SomeMID0Yes1Yes10-10Yes01X9No0Yes1Yes1Full1CH1Yes1N

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

8、tH0-10110-3030-60-1>60-20WillWaitIYes1No0在Matlab的命令窗口中運行如下代碼: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 1 0 1 0 0 0 1;net=new

9、ff(minmax(p),15 1,'tansig','purelin','traincgb');=5;=300;=0.000001;net,tr=train(net,p,t);得到如下曲線圖：容錯性分析：實例輸入輸出WillW結(jié)果OtheWCondWEndConsPricRainResWEstrseaitX11000-10111等X21001100-10不等X3010010011等X4101111001等X51011-101-20不等X6010001111等X7010-111010不等X8000001111等X90111110-20不等X101

10、111-10100不等X11000-110010不等X121111100-11等完成機器學(xué)習(xí)后，對樣本進(jìn)行改變條件輸入，有如下 3種情況:(1)缺1個條件的情況(2)缺2個條件的情況(3)介于中間的情況(1)實例輸入輸出WillWait (1)輸出WillWait(0)結(jié)果Oth sWConWEnCo nPri ceRaiRe sWEstX11100-10110.84560.1526等（1）X21011100-10.07230.9281不等（0）X3010110010.86480.1352等（1）X4100111000.86340.1366等（1）X51111-101-20.23140.768

11、6不等（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)化問題的遺傳算法程序，并計算一個實例。（30分）一，問題描述：所謂旅行商問題，即給定幾個城市，旅行商從中選擇一條最短的路線，使他能夠訪問到每個城市一次，然后返回起OO O 00。0二，運行

12、結(jié)果:三、代碼實現(xiàn)F面是程序的主要核心代碼:using System;using ;using System.ComponentModel;using System.Data;using System.Drawing;using System.Text;using ;using System.Diagnostics;namespaceTSPGA1GUQOO oooo_J-public partial class Forml : Formprivate World _world = new World();public Form1()(InitializeComponent();private

13、void 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_Cl

14、ick( 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 / 20

15、0);y = point.Y * (pictureBoxl.Height / 200);radix = 5;g.DrawEllipse(PensDarkCyan, new Rectangle (x -radix, y - radix, 2 * radix, 2 * radix);if (_world.StartPoint !=Point .Empty)x = _world.StartPoint.X * (pictureBoxl.Width /200);y= _world.StartPoint.Y * (pictureBoxl.Height /200);radix = 5;g.DrawEllip

16、se(PensRed, newRectangle (x - radix,y - radix, 2 * radix, 2 * radix);)pictureBoxl.Image = bm;)private void pictureBox1_MouseMove( object sender, MouseEventArgs e) (Point point = new Point (e.X / (pictureBoxl.Width / 200), e.Y / (pictureBoxl.Height / 200);infosl.Text = point.ToString();)_world.Popula

17、tion =private void Start_Click( object sender, EventArgs e)Int32 .Parse(population.Text);_world.Generation =Int32 .Parse(generation.Text);_world.CrossOverRateDouble.Parse(crossOverRate.Text);_world.MutationRate =Double.Parse(mutationRate.Text);drawVisitRoute();_world.Reset(); private void drawVisitR

18、oute()ReDrawPB();Bitmap bm = ( Bitmap();using ( Graphics g = Graphics .FromImage(bm) Point startPoint = _world.StartPoint;if (startPoint = Point .Empty) MessageBo)Show( "Youmust define a startpoint." );return ;Point pointsNeedToVisit_world.PointsNeedToVisit;Stopwatch sw = new Stopwatch ();

19、sw.Start();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 < vi