用matlab解析实际案例

Matlab大作业专业：东凌经济管理学院班级：小组成员：2012年5月1成员分工：**（组长）：第一题模型一建立，文档编写***：第一题模型二建立，文档编写***：第二题模型建立2第一题某小型超市出售某一品牌八宝粥，其需求量与消费者平均收入和商品价格密切相关，根据近期几个月每一个月的消费记录以及消费者收入市场调查，统计如下表。现在在一个地区新建一所同样的超市，出售同样一款八宝粥，该超市附近消费者平均收入为4000元，超市经理想知道八宝粥定价6元时，进多少货才会比较合适。需求量10075807050659010011060收入4000240048002000120016005200440052001200价格57668754392.1基本假设1）假设该品牌八宝粥的超市库存量与需求量一致，不存在多余库存。2）假设超市每个月就进一次货。3）假设超市之前调查的数据充分准确。4）假设在新超市，人群收入以及商品价格对需求量的影响与之前规律类似。5）假设该品牌八宝粥的需求量除了与消费者收入和商品价格有关，其他因素影响很小，可以忽略不计。2.2符号设定：消费者收入向量。)(incomeI：商品价格向量。)(priceP：商品需求向量。)(trequiremenR2.3模型建立2.3.1模型分析：因为有商品价格P和消费者收入I两个参数对商品需求量R产生影响，所以我们选择采用回归模型解决这个问题。模型一：多元二项式回归模型222211210PIPIRMatlab程序：I=[4000240048002000120016005200440052001200];P=[5766875439];R=[10075807050659010011060];f=[I'P'];rstool(f,R','purequadratic')程序运行结果以及图像：beta=110.53130.0366-26.5709-4.7229*10^(-6)1.8475rmse=4.5362residuals=5.2724-0.7162-4.5158-1.9390-3.33153.45663.4843-3.4452-0.09761.8320预测结果结果分析：由实验可得回归模型为：228475.1)6(^10*7229.4-5709.260366.05313.110PIPIR，因为剩余标准差为4.5362，说明此回归模型的显著性较好。根据回归模型，我们可以得到当消费者平均收入为4000，并且商品价格为6元时，商品需求量为88.4791，变化区间为[-16.3387,16.3387]。模型二：多元线性回归222211210PIPIRMatlab程序：I=[4000240048002000120016005200440052001200];P=[5766875439];R=[10075807050659010011060]';f=[ones(10,1)I'P'(I.^2)'(P.^2)'];[b,bint,r,rint,stats]=regress(R,f);程序结果：b=110.53130.0366-26.5709-4.7229*10^(-6)1.8475bint=57.2602163.80240.01020.0630-43.2247-9.9171-0.0000084748762-0.00000097099030.37453.3205stats=0.970240.66560.0005结果分析：由实验可得回归模型为：228475.1)6(^10*7229.4-5709.260366.05313.110PIPIR参数参数估计置信区间0110.5313[57.2602163.8024]10.0366[0.01020.0630]2-26.5709[-43.2247-9.9171]11-4.7229*10^(-6)[-0.0000084748762-0.0000009709903]221.8475[0.37453.3205]9702.02R6656.40F0005.0p因为相关系数2R充分接近1，并且F的概率p充分接近0，说明回归模型成立并且显著性水平高。当消费者平均收入为4000，并且商品价格为6元时，我们代入回归模型可以得到商品的需求量为88.4791。2.3.2结论与成果1）在新超市，经理每月进货量为88件商品时比较合理。2）通过建立回归模型，我们可以得到商品需求量与商品价格，消费者人均收入之间的关系。应用在实践中，有助于商家更好地管理商品库存，减少库存冗余，从而创造更多利润。3第二题首钢董事长张先生想在工厂搬迁后扩大新厂子的规模，投资钢熔炼和钢拉伸两个项目进行改造，提高钢材生产效率，为此，他特意请来了投资专家，专家估计，投资到钢熔炼技术上可以赚回70倍的投资款，而投资到钢拉伸技术中的可以赚回66倍，同时，项目技术的投资也会带来一定的风险，一般来说风险损失会随着总投资和单项投资的增加而增加。他准备投入5000万在钢熔炼和钢拉伸两个项目的改造上，他想知道如何分配投资资金才能让钢厂获得的收益最大，同时面临的风险最小。3.1基本假设1）假设此专家所作的估计充分正确。2）假设风险损失只与总投资和单项投资有关，其他因素影响很小，可忽略不计。3）假设除了全部5000元都用于钢熔炼和钢拉伸两个项目的改造外，不做任何其他投资与改变。4）假设5000元不一定全部用完。3.2参数设定：分析：因为有两个元素（钢熔炼与钢拉伸中分别投入资金）共同决定了收益值与风险值，所以进行以下设定:1）钢熔炼项目所投入改造资金为x12）钢拉伸项目所投入改造资金为x23）所得收益为f14）所面临风险损失为f23.3模型建立分析：通过题目要求提取出收益函数maxf1（x）=70x1+66x2通过上网查询资料将风险损失函数定义为f2=并根据假设对x1、x2进行约束由于所求问题条件为收益最大同时风险最小，所以收益与风险权重比为1：1。005000)(04.001.002.0)(min6670)(max212122122212211xxxxxxxxxfxxxf线性加权构造目标函数化为最小值问题3.4Matlab程序：首先编辑目标函数M文件f11.mfunctionf=f11(x)f=-0.5*(70*x(1)+66*x(2))+0.5*(0.02*x(1)^2+0.01*x(2)^2+0.04*(x(1)+x(2))^2)调用单目标规划求最小值的函数输入程序x0=[1000,1000];%初始猜测x1=1000,x2=1000;A=[11];b=5000;%满足A*x=b,即x1+x2=5000;lb=zeros(2,1);%x1=0,x2=0;2212221)(04.001.002.0xxxx)(5.0)(5.0max21xfxff)(5.0)(5.0)min(21xfxff[x,fval,exitflag]=fmincon(@f11,x0,A,b,[],[],lb,[])f1=70*x(1)+66*x(2)f2=0.02*x(1)^2+0.01*x(2)^2+0.04*(x(1)+x(2))^23.5程序运行结果x=307.1426414.2862fval=-1.2211e+004exitflag=5f1=4.8843e+004f2=2.4421e+004fmincon函数分析:[x,fval,exitflag,output]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub)用来找到具有初步估计的有约束多标量方程组的最小解。x（1），x（2）……x(n)是人为规定顺序的变量的统一表示法。参数中fun为目标函数，x0为变量的初始值，x为返回的满足要求的变量的值。A和b表示线性不等式约束,即满足A*x=b.Aeq，Beq表示线性等式约束，lb和ub分别为变量的下界和上界约束.返回值fval为目标函数。exitflag0表示优化结果收敛于解，exitflag=0表示优化超过了函数值的计算次数，exitflag0表示优化不收敛。3.6结果分析与成果展示:当给钢熔炼投入资金307元,给钢拉伸投入414元时,可获得最大收益约48843元,面临最小风险损失24421元。通过此模型的建立，我们可以在知道精确回报数字的前提下计算出任意N个项目的资金投入方案，以达到最大收益与最小风险损失共存的目的。有了这个辅助工具，我们能更加直观理性的衡量项目投入的资金对收益值和风险值的影响力。