决策支持系统作业最终)

决策支持系统导论期末作业姓名：齐鹏学号：日期：2012年7月9日1、设某企业生产多种最终产品Y=(yij)，各种产品的单价为Pi，它们的投入产出直接消耗系数为A=(aij)，企业的资源（煤、电力、劳力）的约束方程为BX<=>h （“<=>”表示<、=、>），其中，B=(bij)是资源消耗系数矩阵，X=(xi)是企业总产品向量，h是资源约束向量。

为使企业净产值最大，其目标方程S=∑Piyi→max，试安排生产计划（求总产品X和最终产品Y）。

请设计该企业的生产计划决策支持系统，画出DSS运行结构图，并对总控程序、模型程序、数据库进行结构和功能说明。

提示：该决策支持系统需要利用3个模型（投入产出模型、线性规划模型和报表模型（打印投入产出表））和两个数据库（投入产出数据库和线性规划数据库）。

在DSS总控程序中要详细说明何时调用哪个模型运行，何时存取哪个数据库中的数据，何时进行数据计算。

该DSS需要两次调用投入产出模型：一次计算中间结果，一次计算最后结果。

请注意，模型程序应该是一个标准程序，在一定的参数控制下，可得到中间结果，也可得到最终结果。

该模型程序既适合于该问题的DSS，也适合于其他问题的DSS,不能是一个专用的模型程序。

（40分）一、模型1.投入产出模型：可以确定与的关系，。

2.线性规划模型：根据约束方程BX<=>h与目标方程S=∑Piyi→max可以计算出最优的yi。

3.报表模型：根据最终产品Y=(yij)，X=(xi)是企业总产品自动生成报表。

二、数据库1.投入产出数据库2.线性规划数据库三、DSS 运行结构图 ：附近是否有其他卡拉OK 厅； 线性规划数据库： B=(bij)是资源消耗，根据公式计算（6）Raining（下雨）：外面是否在下雨；（7）Reservation（预约）：是否预约过；（8）WaitEstimate（等候时间估计）：估计的等候时间（值为0—10，10—30，30—60，>60，单位为分钟）。

训练集见表：要求：建立BP神经网络模型，并进行容错性分析。

（30分）一、为训练集赋值在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=newff(minmax(p),[15 1],{'tansig','purelin'},'traincgb'); =5;=300;=0.000001;[net,tr]=train(net,p,t);得到如下曲线图：容错性分析：完成机器学习后，对样本进行改变条件输入，有如下3种情况：（1）缺1个条件的情况（2）缺2个条件的情况（3）介于中间的情况（1）3、编制旅行商路径优化问题的遗传算法程序，并计算一个实例。

（30分）一，问题描述：所谓旅行商问题，即给定几个城市，旅行商从中选择一条最短的路线，使他能够访问到每个城市一次，然后返回起点。

二，运行结果：三、代码实现下面是程序的主要核心代码：using System;using ;using ponentModel;using System.Data;using System.Drawing;using System.Text;using ;using System.Diagnostics;namespace TSPGA{public partial class Form1 : Form{private World _world = new World();public Form1(){InitializeComponent();}private 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_Click(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 = 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.Width / 200);y = _world.StartPoint.Y * (pictureBox1.Height / 200);radix = 5;g.DrawEllipse(Pens.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_Click(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(){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.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(pointsNeedToVisit[index]);}int td = 0;Point tp = visitPoints[0];for (int i = 1; i < visitPoints.Count; i++){td += GetManhattanDis(tp, visitPoints[i]);tp = visitPoints[i];}totalDis.Text = "Total Distace: " + td;for (int i = 0; i < visitPoints.Count; i++){Point p = new Point();p.X = visitPoints[i].X * pictureBox1.Width / 200; p.Y = visitPoints[i].Y * pictureBox1.Height / 200; visitPoints[i] = p;}g.DrawLines(Pens.SeaShell, visitPoints.ToArray());}pictureBox1.Image = bm;}private void random_Click(object sender, EventArgs e){_world.Reset();Random rnd = new Random();int c = Int32.Parse(cityNum.Text);_world.StartPoint = new Point(rnd.Next(0, 200), rnd.Next(0, 200));for (int i = 0; i < c; i++){_world.AddPointNeedToVisit(new Point(rnd.Next(0, 200), rnd.Next(0, 200)));}ReDrawPB();}private static int GetManhattanDis(Point p1, Point p2){return Math.Abs(p1.X - p2.X) + Math.Abs(p1.Y - p2.Y); }}}。