决策支持系统作业最终

1、决策支持系统导论 期末作业 姓名：齐鹏 学号： 日期： 2012年 7月 9日 1、设某企业生产多种最终产品丫=（yij），各种产品的单价为Pi,它们的投入产出直接消耗系 数为A=（aij），企业的资源（煤、电力、劳力）的约束方程为BXh （“ ”表示）， 其中，B=（bij）是资源消耗系数矩阵，X=（xi）是企业总产品向量，h是资源约束向量。为使企 业净产值最大，其目标方程S=EPiyi-max,试安排生产计划（求总产品X和最终产品丫）。 请设计该企业的生产计划决策支持系统， 画出DSS运行结构图，并对总控程序、模型程序、 数据库进行结构和功能说明。 提示：该决策支持系统需要利用 3 个模型

2、（投入产出模型、线性规划模型和报表模型（打 印投入产出表）和两个数据库（投入产出数据库和线性规划数据库）。在DSS总控程序中 要详细说明何时调用哪个模型运行，何时存取哪个数据库中的数据，何时进行数据计算。 该DSS需要两次调用投入产出模型：一次计算中间结果，一次计算最后结果。 请注意，模型程序应该是一个标准程序，在一定的参数控制下，可得到中间结果，也可得 到最终结果。该模型程序既适合于该问题的 DSS,也适合于其他问题的DSS不能是一个专 用的模型程序。（40 分， 、模型 1. 投入产出模型：可以确定鞫I 与的关系,X - 5 -力)1 。 2. 线性规划模型：根据约束方程 BXh与目标方程

3、S=E Piyi max可以计算出最优的 yi。 3. 报表模型：根据最终产品丫=(yij)，X=(xi)是企业总产品自动生成报表。 二、数据库 1.投入产出数据库 字段名 数据类型 长度 是否可为空 int 8 是 int 8 是 int 8 是 int 8 是 int 8 是 int 16 是 int 16 是 int 16 是 int 16 是 int 16 是 float 16 是 float 16 是 float 16 是 float 16 是 float 16 是 2. 线性规划数据库 字段名 数据类型 长度 是否可为空 float 8 是 float 8 是 float 8 是

4、float 8 是 float 8 是 hl float 16 是 h2 float 16 是 h3 float 16 是 float 16 是 hi float 16 是 float 16 是 float 16 是 float 16 是 float 16 是 float 16 是 、DSS运行结构图 投入产出数据库: 开始投入产出数据库: 调用投入产出模型调用投入产出模型 2、考虑去卡拉OK厅唱歌的时候，是否要等待包间的问题 调用内产实出模莫型 产如出模型生可用于描述该领 （1） Others （其他地点）：附近是否有其他卡拉 0K厅; 调用报表模型: 调用报表模型Cond （等候条件）：供

5、顾客等候的地方是否舒适； （3） Weekend （周末）： 若是 根据 X=(xi), 丫=(yij)编制 I （4） Conssumers（顾客）:店中有多少顾客（值为投入入产出没数据库Some （一些）或Full （满 显示报表 座）； (5) Price （价格）：价格范围（值为Cheep （便宜），Middle （中等），Expensive （较贵）; (6) Raining （下雨）：外面是否在下雨; (7) Reservation （预约）：是否预约过; （8） WaitEstimate（等候时间估计）：估计的等候时间（值为010, 10 30, 30 60, 60， 单位为分钟

6、）。 训练集见表: 实例 属性 目 标 WillWai t Other WCon WEnd Cons Price Rai n Res WEst s d X1 Yes1 NoO NoO Some O EX-1 NoO Yes1 O-1O 1 Yes X2 Yes1 NoO NoO Fullll CH1 NoO NoO 3O-6O -1 No X3 NoO Yes1 NoO Some O CH1 NoO NoO O-1O 1 Yes X4 Yes1 NoO Yes1 Fulll CH1 Yes1 NoO 1O-3O O Yes X5 Yes1 NoO Yes1 Fulll EX-1 NoO Yes

7、1 6O -2 No X6 NoO Yes1 NoO Some O MIDO Yes1 Yes1 O-1O 1 Yes X7 NoO Yes1 NoO None- 1 CH1 Yes1 NoO O-1O 1 No X8 NoO NoO NoO Some O MIDO Yes1 Yes1 O-1O 1 Yes X9 NoO Yes1 Yes1 Fulll CH1 Yes1 NoO 60 -2 No X10 Yes1 Yes1 Yes1 Fulll EX-1 NoO Yes1 1O-3O O No X11 NoO NoO NoO None- 1 CH1 NoO NoO O-1O 1 No X12

8、Yes1 Yes1 Yes1 Fulll CH1 NoO NoO 3O-6O -1 Yes 要求：建立BP神经网络模型，并进行容错性分析。（30 分） 一、 为训练集赋值 属性 值名 属性 1 值 属性2 值 属性3 值 属性 4 值 Others A Yes 1 No 0 WCo nd B Yes 1 No 0 WEnd C Yes 1 No 0 Cons D Full 1 Some 0 None -1 Price E Chea P 1 Middl e 0 Expe nsiv e -1 Rai n F Yes 1 No 0 Res G Yes 1 No 0 WEst H 0-10 1 10-

9、30 0 30-60 -1 60 -2 WillWait I Yes 1 No 0 在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；

10、 net=n ewff(mi nm ax(p),15 1,ta nsig,pureli n,tra in cgb); =5; =300; =0.000001; n et,tr=trai n(n et,p,t); 得到如下曲线图： 容错性分析： 实例 输入 输出 WillWai t 结果 Other s WCon d WEnd Cons Price Rai n Res WEst X1 1 0 0 0 -1 0 1 1 1 等 X2 1 0 0 1 1 0 0 -1 0 不等 X3 0 1 0 0 1 0 0 1 1 等 X4 1 0 1 1 1 1 0 0 1 等 X5 1 0 1 1 -1 0

11、 1 -2 0 不等 X6 0 1 0 0 0 1 1 1 1 等 X7 0 1 0 -1 1 1 0 1 0 不等 X8 0 0 0 0 0 1 1 1 1 等 X9 0 1 1 1 1 1 0 -2 0 不等 X10 1 1 1 1 -1 0 1 0 0 不等 X11 0 0 0 -1 1 0 0 1 0 不等 X12 1 1 1 1 1 0 0 -1 1 等 完成机器学习后，对样本进行改变条件输入，有如下3种情况: （1）缺1个条件的情况 （2）缺2个条件的情况 （3）介于中间的情况 （1） 实例 输入 输出 输出 结果 Oths WCon WEn Co Price Rai Res WE

12、s WillWait WillWait n t （1） （0） X1 1 1 0 0 -1 0 1 1 0.8456 0.1526 等 (1) X2 1 0 1 1 1 0 0 -1 0.0723 0.9281 不 等 （0） X3 0 1 0 1 1 0 0 1 0.8648 0.1352 等 (1) X4 1 0 0 1 1 1 0 0 0.8634 0.1366 等 (1) X5 1 1 1 1 -1 0 1 -2 0.2314 0.7686 不 等 （0） X6 0 1 1 0 0 0 1 1 0.9125 0.0875 等 (1) X7 0 1 1 -1 0 1 0 1 0.3410

13、 0.6595 不 等 （0） X8 0 1 1 0 0 1 1 1 0.9714 0.0286 等 (1) X9 0 0 0 0 1 1 0 -2 0.4928 0.5012 不- 一疋 X10 1 1 0 0 1 0 1 0 0.5681 0.4319 不- 一疋 X11 1 1 1 -1 1 0 0 1 0.4123 0.5877 不- 一疋 X12 1 1 0 0 0 0 0 -1 0.6428 0.3572 不- 一疋 3、编制旅行商路径优化问题的遗传算法程序，并计算一个实例 (30 分) ，问题描述：所谓旅行商问题，即给定几个城市，旅行商从中选择一 条最短的路线，使他能够访问到每个

14、城市一次，然后返回起 点。 ，运行结果: O O O OO o OgO O OOQqqo o o so 0 o O OOoVP o O )O o o o OO 008 0 三、代码实现 下面是程序的主要核心代码: using System; using ; using System.Comp onen tModel; using System.Data; using System.Drawing; using System.Text; using ; using System.Diagnostics; namespaceTSPGA public partial classForm1 : For

15、m private World _world = new World(); public Form1() InitializeComponent(); private void pictureBox1_MouseUpo( bject sender,MouseEventArgse) Point point = new Point(e.X / (pictureBox1.Width / 200), e.Y / (pictureBox1.Height / 200); if (e.Button = MouseButtons.Left) _world.AddPointNeedToVisit(point);

16、 cityNum.Text = _(); else _world.StartPoint = point; ReDrawPB(); private void Reset_Clicko( bject sender,EventArgs e) _world.Reset(); ReDrawPB(); private void ReDrawPB() Bitmap bm = new Bitmap(pictureBox1.Width, pictureBox1.Height); using (Graphicsg = Graphics.FromImage(bm) int x, y, radix; foreach

17、(Point point in _world.PointsNeedToVisit) x = point.X * (pictureBox1.Width / 200); y = point.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.W

18、idth / 200); y = _world.StartPoint.Y * (pictureBox1.Height / 200); radix = 5; radix); g.DrawEllipse(Pens.Red,new Rectangle(x - radix, y - radix, 2 * radix, 2 * pictureBox1.Image = bm; private void pictureBox1_MouseMoveo( bject sender,MouseEventArgse) Point point = new Point(e.X / (pictureBox1.Width

19、/ 200), e.Y / (pictureBox1.Height / 200); infos1.Text = point.ToString(); private void Start_Click(object sender,EventArgs e) _world.Population = Int32.Parse(population.Text); _world.Generation = Int32.Parse(generation.Text); _world.CrossOverRate =Double.Parse(crossOverRate.Text); _world.MutationRat

20、e = Double.Parse(mutationRate.Text); drawVisitRoute(); /_world.Reset(); private void drawVisitRoute() ReDrawPB(); Bitmap bm = (Bitmap(); using (Graphicsg = Graphics.FromImage(bm) Point startPoint = _world.StartPoint; if (startPoint = Point.Empty) MessageBo.xShow(You must define a start point.); retu

21、rn; Point pointsNeedToVisit = _world.PointsNeedToVisit; Stopwatchsw = new Stopwatch(); sw.Start(); int visitRoute = _world.GetVisitRoute(); sw.Stop(); infos3.Text = Used: + sw.ElapsedMilliseconds +ms.; List visitPoints = new List(); 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 = v