DSS实验报告例子

1、 本科实验报告课 程 名 称： 决策支持系统 课 程 编 号： 07010192 学 生 姓 名： 学 号： 学 院： 信息科学技术学院 系： 数 学 系 专 业： 信息管理与信息系统 指 导 教 师： 谭 满 春 教 师 单 位： 数 学 系 开 课 时 间：2014 2015学年度第 1学期暨南大学教务处2014 年 12 月 1 日 决策支持系统 课程实验项目目录学生姓名：林培纯 学号：2012051378序号实验项目编号实验项目名称*实验项目类型成绩指导教师107010192901基金使用模型与方法验证性谭满春207010192902按揭购房的利率设计验证性谭满春30701019290

2、3生产计划中的产量问题综合性谭满春407010192904最佳广告编排方案设计验证性谭满春507010192905供应与选址决策问题综合性谭满春607010192906商品的需求量预测与决策综合性谭满春707010192907决策支持系统分析与设计综合性谭满春89101112131415161718*实验项目类型：演示性、验证性、综合性、设计性实验。*此表由学生按顺序填写。暨南大学本科实验报告专用纸课程名称 决策支持系统 成绩评定 实验项目名称 基金的使用模型与计划 指导教师 谭满春 实验项目编号07010192901实验项目类型验证性实验地点南海楼209学生姓名 学号 学院 信息科学技术学院

3、 系 数学系 专业 信息管理与信息系统实验时间 2014 年 11 月 3日 上午11月3日下午 温度 湿度【实验目的】1.介绍与线性方程组有关的基本概念。2.了解线性方程组的消去法、迭代法等基本求解方法。3.学习MATLAB软件中有关线性方程组运算的命令。【实验内容】某校基金会有一笔数额为M元的基金，打算将其存入银行，当前银行存款及各期的利率见下表，取款政策参考银行的现行政策。校基金会计划在n年内每年用部分本息奖励优秀师生，要求每年的奖金额大致相同，且在n年末仍保留原基金数额。校基金会希望获得最佳的基金使用计划，以提高每年的奖金额。请你帮助校基金会在上述情况下设计基金存款使用方案，并对M=5

4、000万元，n=10年给出具体结果。银行存款随税后年利率%活期0.792半年期1.664一年期1.800二年期1.944三年期2.160五年期2.304【实验方法与步骤】1问题的分析问题本身含有一些不确定的因素，比如说基金到位的时间，每年奖学金发放的日期，银行利率的变动情况等。为使问题简化，先做如下假设：假设1：该笔资金于年底一次性到位，自下年起每年年底一次性发放的奖金，每年发放的奖金额尽可能的相同;假设2: 银行存款利率执行现行利率标准，且在n年内不发生变化。 将总额M分成11份，为用作奖学金的钱，这些钱经过存入银行加息的过程，到第i年取出用作第i年的奖学金钱；M11为用作奖池的钱，即经过1

5、0年的银行加暨南大学本科实验报告专用纸(附页)过程，又变成原来的总额M。设占总额的比例为，M10与M11加起来占总额比例为x10，则有对于存款方案，可知，能一次存时间长点，则利息就更多，所以应该尽量选足够多的五年期，其次再到三年期，再到二年期，最后才考虑一年期，而半年期和活期可以忽略。所以对于用作第i年奖学金钱的Mi，它所存的五年期的次数为m(i)=fix(i/5)，三年期的次数为k(i)=fix(rem(i,5)/3)，二年期的次数为l(i)=fix(rem(rem(i,5),3)/2)，一年期的次数为t(i)=i-5m(i)-3k(i)-2l(i)。由于每年的奖学金额度都是相同的，所以我们

6、有：这样就联立了一个十元线性方程组，将后面九个式子中的M消去，可得增广矩阵为用MATLAB编入代码计算：d1=1+1.8/100;d2=1+2*1.944/100;d3=1+3*2.16/100;d5=1+5*2.304/100;n=10;M=5000;for j=1:na(1,j)=1;endfor i=2:na(i,1)=1+1.8/100;endfor i=2:n %构造系数矩阵a(i,j)m(i)=fix(i/5);k(i)=fix(rem(i,5)/3); l(i)=fix(rem(rem(i,5),3)/2); t(i)=i-5*m(i)-3*k(i)-2*l(i); a(i,i)

7、=-d5m(i)*d3k(i)*d2l(i)*d1t(i);endb(1,:)=1;b(n,:)=-1;x=ab; %用常数向量除以系数矩阵的方法求解线性方程组for i=1:neachyear(i)=M*x(i);endrichaward=d1*M*x(1)；最终得到结果：richaward = 109.8169【结果分析】每年大约发放109.8万元的奖学金，同时在M=5000万元，n=10年的情形下不同年限基金的存款方式如下表（单位：万元）：i12345678910Mxi107.9105.7103.1101.398.596.794.892.590.84108.7用于第10年发放奖金和剩余基

8、金的总额为4108.7万元，经过10年的利息累积后，将增长为大约5109.8万元。109.8万元用于当年发放，剩余5000万。我们可以看到，基金的使用计划具有周期性，当n年到期，便可按原方案进入下一周期；如果利率或政策有变，我们只需在下一周期开始前，对利率等做些改变或引入其他参数，整体基金投资方案仍可以沿用。 教师评语暨南大学本科实验报告专用纸课程名称 决策支持系统 成绩评定 实验项目名称 按揭购房的利率设计 指导教师 谭满春 实验项目编号07010192902实验项目类型验证性实验地点南海楼209学生姓名 学号 学院 信息科学技术学院 系 数学系 专业 信息管理与信息系统实验时间 2014

9、年 11 月 8日 上午11月8日下午 温度 湿度【实验目的】1. 介绍与非线性方程组有关的基本概念。2. 了解非线性方程组的图解法、二分法、迭代法等基本求解方法3. 学习MATLAB的有关命令【实验内容】计算一下贷款年利率。建筑面积总价30%首付70%按揭月还款85.98m236万10.8万30年1436元【实验方法与步骤】设xk为第k个月的欠款数，a为月还款数，r为月利率。可得那么有 根据a=0.1436，x0=25.2，x360=0得到这是一个月利率r的高次代数方程，从中解得r，年利率R=12r。相应的MATLAB代码： r=fzero(25.2*(1+x)360-(1+x)360-1)

10、/x*0.1436,0.0198/12) R=12*rR =0.0553得到年利率为5.53%。【结果分析】本例是公积金和商业性贷款的平均利率，所得解与当时的实寄情况相符。教师评语暨南大学本科实验报告专用纸课程名称 决策支持系统 成绩评定 实验项目名称 生产计划中的产量问题 指导教师 谭满春 实验项目编号07010192903实验项目类型综合性实验地点南海楼209学生姓名 学号 学院 信息科学技术学院 系 数学系 专业 信息管理与信息系统实验时间 2014 年 11 月 12日上午11月12日下午 温度 湿度【实验目的】1介绍无约束最优化方法的一些基本概念。2了解几种常见的无约束优化问题的求解

11、方法，如迭代算法、最速下降法（梯度法）、牛顿法（Newton）、拟牛顿法。3学习掌握用MATLAB优化工具箱中的命令来求解无约束优化问题。【实验内容】某公司生产一种产品有甲、乙两个品牌，试讨论产销平衡下的最大利润。所谓产销平衡指公司的产量等于市场上的销量。利润既取决于销量和单件价格，也依赖于产量和单件成本。按照市场规律，甲种品牌的价格固然会随其销量的增长而降低；同时乙品牌销量的增长也会使甲的价格有稍微下降，根据实际情况，可以确定价格与销量成线性关系，即3002.350.09乙的价格遵循同样的规律，有4800.142.98甲品牌的成本会随着其产量的增长而降低，按实际情况可假设为负指数关系，即有3

12、8116乙品牌的成本遵循同样的规律，有94145试确定甲、乙两种品牌的产量，使公司获得的总利润最大。【实验方法与步骤】1MATLAB优化工具箱中解无约束问题的命令和参数options的基本用法下面用例题来予以说明例1求解 5，12输入命令： x=fmin( exp(x) - 5*x , 1 , 2 ) , f=exp( x ) - 5*xx = 1.6094 f = -3.0472同时，该问题可以用fminbnd求解，得到的解答一样输入命令： x , fval=fminbnd( exp( x ) - 5*x , 1 , 2 )x = 1.6094fval = -3.0472例2用拟牛顿法的DF

13、P公式求解，初值为（1，2），要求精度为7，给出函数计算次数以及函数值首先建立M文件，文件名取函数名fun.mfunction f = fun( x )f = 4 * x( 1 )2 + x( 2 )2 - x( 1 )3 * x( 2 )输入命令： x0=1, 2; opt( 2 )=1e-6;%设置自变量要求的精度 opt( 3 )=1e-6;%设置函数所要求的精度 opt( 6 )=1;%搜索算法用拟牛顿法的DFP公式求解 x, opt=fminu( fun , x0, opt);%返回参数options的值给向量opt可得到结果x = 1.0e-008 * -0.3493 0.7566

14、参数options的值返回给向量opt，全部元素略。 n=opt( 10 );%函数计算次数赋值给nn = 35 y=opt( 8 );%在极值点处函数值赋值给yy = 1.0606e-0162问题的分析与求解由题设可知，甲品牌产品单件获利为，乙品牌产品单件获利为，由产销平衡原理，所有产品的销量即为产量，则甲、乙两种产品总获利为（）（）容易看出，原问题实际上转化为求二元函数的极大值，为用MATLAB优化工具箱中的fminunc求解，需将其转化为求函数的最小值。为确定初始值，先忽略成本，并令价格中项的较小系数0.09和中项较小的系数0.14等于零（因为它们对价格的作用比较微弱，暂时可忽略不计），

15、则确定初值问题转化为求（3002.35）（4802.98）的极值，很容易可以求得63.83，80.54，我们用它作为原问题的初始解。首先建立M文件，文件名取函数名fun1.mfunction y = fun1( x )p1 = 300 - 2.35 * x( 1 ) - 0.09 * x( 2 );q1 = 38 * exp( - 0.023 * x( 1 ) ) + 116;p2 = 480 - 0.14 * x( 1 ) - 2.98 * x( 2 );q2 = 94 * exp( - 0.018 * x( 2 ) )+145;y = - ( p1 - q1 ) * x( 1 ) - (

16、p2 - q2 ) * x( 2 )输入命令： x0=63.83; 80.54; x, fval=fminunc( fun1 , x0);可得到结果：x = 35.8482 54.7380fval = -1.0015e+004甲种品牌产量为35.8482，乙种品牌产量为54.7380，最大利润和为1.0015e+004。用fminu可得到同样结果，同时根据参数options可观察到函数采用拟牛顿的DFP公式并迭代计算了20次。教师评语暨南大学本科实验报告专用纸课程名称 决策支持系统 成绩评定 实验项目名称 最佳广告编排方案设计 指导教师 谭满春 实验项目编号07010192904实验项目类型验

17、证性实验地点南海楼209学生姓名 学号 学院 信息科学技术学院 系 数学系 专业 信息管理与信息系统实验时间 2014 年 11 月 19日上午11月19日下午 温度 湿度【实验目的】1了解线性规划问题及其可行解、基本解、最优解的概念。2通过对实际应用问题的分析，初步掌握建立线性规划模型的基本步骤和方法。3学习掌握MATLAB软件求解有关线性规划的命令。【实验内容】一家广告公司想在电视、广播上做公司的宣传广告，其目的是争取尽可能多地招徕顾客。下表是公司进行市场调研的结果：电视网络媒体杂志白天最佳时段每次做广告费用（千元）45862512受每次广告影响的顾客数（千人）350880430180受每

18、次广告影响的女顾客数（千人）260450160100这家公司希望总广告费用不超过750（千元），同时还要求：（1）受广告影响的妇女超过200万；（2）电视广告的费用不超过450（千元）；（3）电视广告白天至少播出4次，最佳时段至少播出2次；（4）通过网络媒体、杂志做的广告要重复5到8次。【实验方法与步骤】建立线性规划模型有三个基本步骤：第一步，找出待定的未知变量（决策变量），并用代数符号来表示它们；第二步，找出问题的所有限制或约束条件，写出未知变量的线性方程或线性不等式；第三步，找到模型的目标，写成决策变量的线性函数，以便求其最大或最小值。1问题的分析与模型的建立首先，确定决策变量，要求如何安

19、排白天电视、最佳时段电视、网络媒体、杂志广告的次数，用符号表示，分别设定为，；其次，确定所有的约束条件，广告总费用不超过750（千元），则有45862512750受广告影响的女顾客数不少于200万，则有2604501601002000电视广告费用不超过450（千元），且白天至少播4次，最佳时段至少播出2次，则有4586450 ，4 ，2由于网络媒体和杂志广告要重复5到8次，则有58 ，58最后，确定问题的目标函数，由题意知确定广告编排方案，使得受各种广告影响的潜在顾客总数：350880430180最多。故该问题完整的线性规划模型如下：35088043018045862512750 260450

20、1601002000458600450000800084 ，2，5，52MATLAB计算机求解用MATLAB求解的程序代码： c=-350 -880 -430 -180;%取将目标函数标准化 a=45 86 25 12; -260 -450 -160 -100; 45 86 0 0; 0 0 1 0; 0 0 0 1; b=750; -2000; 450; 8; 8; lb=4; 2; 5; 5; x , Fval=linprog(c, a, b, , , lb, )%无等式约束条件和的上界，取表缺省Optimization terminated successfully.x = 4.0000

21、 3.1395 8.00008.0000Fval =-9.0428e+003【结果分析】引例问题的目标函数是求受广告影响的最多顾客人数，而MATLAB命令linprog针对线性规划模型（5）求最小值，那么我们取，将目标函数化成标准形式，在求得的最小值后，我们即可得到的最大值，根据约束条件，受广告影响的最多潜在顾客人数为9042800人。在这里，用命令lp可以求得相同的结果。教师评语暨南大学本科实验报告专用纸课程名称 决策支持系统 成绩评定 实验项目名称 供应与选址决策问题 指导教师 谭满春 实验项目编号07010192905实验项目类型综合性实验地点南海楼209学生姓名 学号 学院 信息科学技

22、术学院 系 数学系 专业 信息管理与信息系统实验时间 2014 年 11 月 23日上午11月23日下午 温度 湿度【实验目的】1了解非线性规划问题的基本概念和求解方法。2通过对应用问题的分析、建模、求解，加深对非线性规划理论的理解。3学习掌握MATLAB有关非线性规划求解的命令。【实验内容】某建筑公司有6个建筑工地要开工，每个工地的位置（用平面坐标，表示，距离单位：千米）及水泥日用量（吨）由表1给出。目前有两个废旧的料场位于（5,1），（2,7）处，现需要重新建设、两个料场，日储存量各为20吨。试问料场、各设在何处，并且如何制定每天的原料供应计划，使得从、两料场分别向各工地运送水泥总的吨千米

23、数最小。表1工地的位置（，）及各工地水泥日用量位置及用量123456横坐标1.258.750.55.7537.25纵坐标1.250.754.7556.57.75（吨）3547611【实验方法与步骤】1问题的分析记工地的位置为（，），水泥日用量为，1，6；料场位置为（，），日储量为，1，2；从料场向工地的运送量为。这个优化问题的目标函数（总吨千米数）可表为（6）各工地的日用量必须满足，所以有，1，6（7）各料场的运送量不能超过日储量，所以，1，2（8）则该问题的决策变量为料场位置，和、两料场往各工地的运送量，问题归结为在约束条件（7）、（8）及决策变量为非负的情况下求料场位置（，）和运送量使（6

24、）的总吨千米数最小。由于目标函数对，是非线性的，所以在求新建料场位置和用料时是非线性规划模型。2MATLAB计算机求解首先定义非线性规划的M文件函数：function f,g=liaoch(x)a=1.25, 8.75, 0.5, 5.75, 3, 7.25;b=1.25, 0.75, 4.75, 5, 6.5, 7.75;d=3, 5, 4, 7, 6, 11;e=20, 20;f1=0;%f1是料场A到各工地的吨千米总数，其中x(1)至x(6)为料场A往各工地的运送量, (x(13), x(14)为A的位置for i=1:6 s(i)=sqrt(x(13)-a(i)2+(x(14)-b(i

25、)2); f1=s(i)*x(i)+f1;endf2=0;%f2是料场B到各工地的吨千米总数，其中x(7)至x(12)为料场B往各工地的运送量, (x(15), x(16)为B的位置for i=7:12 s(i)=sqrt(x(15)-a(i-6)2+(x(16)-b(i-6)2); f2=s(i)*x(i)+f2; endf=f1+f2;for i=1:6 g(i)=x(i)+x(i+6)-d(i);%各工地用量必须满足endg(7)=sum(x(1:6)-e(1);%各料场运送量不超过日储量g(8)=sum(x(7:12)-e(2);然后，在MATLAB命令框里输入求解命令： x0=zer

26、os(1,12) 5 1 2 7;%取零为往各工地运送量的初值，取废弃料场位置为新料场的初值 vlb=zeros(1,16);%求解下界为零 op(13)=6;op(14)=2000;%确定等式约定的数目和命令求解的最大迭代次数 x,op=constr(liaoch,x0,op,vl), y=op(8)将结果作成列表的形式为工地123456新料场位置354710（5.6959，4.9284）0000511（7.2500，7.7500）最小的总吨千米数为89.8835。【结果分析】在命令框输入命令 text(1.25,1.25,+3); text(8.75,0.75,+5); text(0.5,

27、4.75,+4) text(5.75,5,+7); text(3,6.5,+6); text(7.25,7.75,+11) text(5.6959,4.9284,A) text(7.2500,7.7500,B)图9.1图9.1画出了工地、新料场的位置（+为工地，旁边的数字为用量，A、B分别表示新料场的位置，可以看出，新料场应建在两个用量最大的工地旁边，这个结果预先估计到了吗？所求得的解是全局最优解吗？实际上改变决策变量的初始值，那么给各工地的供应量和新料场的位置都要发生变化，不妨试试。【教师评语】暨南大学本科实验报告专用纸课程名称 决策支持系统 成绩评定 实验项目名称 商品的需求量预测与决策

28、指导教师 谭满春 实验项目编号07010192906实验项目类型综合性实验地点南海楼209学生姓名 学号 学院 信息科学技术学院 系 数学系 专业 信息管理与信息系统实验时间 2014 年 12 月 1日上午12月1日下午 温度 湿度【实验目的】1了解回归分析的基本原理和方法。2学习用回归分析的方法解决问题，初步掌握对变量进行预测和控制。3学习掌握用MATLAB命令求解回归分析问题。【实验内容】现有某种商品的需求量、消费者的平均收入、商品价格的统计数据如表1所示，试用所提供的数据预测消费者平均收入为1000、商品价格为6时的商品需求量。需求量10075807050659010011060收入1

29、0006001200500300400130011001300300价格5766875439【实验方法与步骤】1引例问题的分析求解由问题提供的数据，我们可以初步判断，商品的需求量与消费者的平均收入和商品价格之间存在某种相关关系，具体的函数关系式我们还不清楚。输入三组数据，我们先独立分析商品需求量与消费者平均收入，商品需求量与价格之间存在何种关系： x1=1000 600 1200 500 300 400 1300 1100 1300 300;%消费者的平均收入 x2=5 7 6 6 8 7 5 4 3 9;%商品价格 y=100 75 80 70 50 65 90 100 110 60;%商品

30、的需求量 plot(x1,y,+)%以消费者的平均收入和商品的需求量所对应的离散点作图 plot(x2,y,+)% 以商品的价格和商品的需求量所对应的离散点作图由上面两图我们看到商品的需求量随着消费者平均收入增加呈线性递增的趋势，而随着商品的价格增加呈线性递减趋势，这样我们可初步判断商品需求量与消费者平均收入和商品价格之间存在某种线性相关的关系。接下来用多元线性回归来进行分析检验： x=ones(10,1) x1 x2; b,bint,r,rint,stats=regress(y,x)b = 111.6918 0.0143 -7.1882bint = 56.0503 167.3334-0.0120 0.0406-13.2306 -1.1458stats = 0.8944 29.6533 0.0004可知回归系数111.6918，0.0143，-7.1882，它们的置信区间为bint，均包含了回归系数的估计值，stats第一个分量为0.8944，第三个分量p0.00040.05，拒绝H0，说明回归方程系数不为0，线性回归方程模型111.69180.01437.1882（12）成立。继续对残差进行分析，作残差图： rcoplot(r