最新用matlab解析实际案例Word格式文档下载.docx
《最新用matlab解析实际案例Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《最新用matlab解析实际案例Word格式文档下载.docx（7页珍藏版）》请在冰豆网上搜索。
80
70
50
65
90
110
60
收入
4000
2400
4800
2000
1200
1600
5200
4400
价格
5
7
6
8
4
3
9
2.1基本假设
1）假设该品牌八宝粥的超市库存量与需求量一致，不存在多余库存。
2）假设超市每个月就进一次货。
3）假设超市之前调查的数据充分准确。
4）假设在新超市，人群收入以及商品价格对需求量的影响与之前规律类似。
5）假设该品牌八宝粥的需求量除了与消费者收入和商品价格有关，其他因素影响很小，可以忽略不计。
2.2符号设定
I（income）：
消费者收入向量。
P（price）：
商品价格向量。
R（requirement）：
商品需求向量
2.3模型建立
2.3.1模型分析：
因为有商品价格P和消费者收入I两个参数对商品需求量R产生影响，所以我们选择采用回归模型解决这个问题。
模型一：
多元二项式回归模型
R01I2P11I222P2
Matlab程序：
I=[4000240048002000120016005200440052001200];
P=[5766875439];
R=[10075807050659010011060];
f=[I'
P'
];
rstool（f,R'
'
purequadratic'
）
程序运行结果以及图像：
beta=
110.5313
0.0366
-26.5709
-4.7229*10^（-6）
1.8475
rmse=
4.5362
residuals=
5.2724
-0.7162
-4.5158
-1.9390
-3.3315
3.4566
3.4843
-3.4452
-0.0976
1.8320
结果分析：
由实验可得回归模型为：
R110.53130.0366I26.5709P-4.7229*10^（6）I21.8475P2，因为剩余标准差为4.5362，说明此回归模型的显著性较好。
根据回归模型，我们可以得到当消费者平均收入为4000，并且商品价格为6元时，商品需求量为88.4791，变化区间为[-16.3387,16.3387]。
模型二：
多元线性回归
RIPI2P2
R01I2P11I22P
R=[10075807050659010011060]'
;
f=[ones（10,1）I'
（I.^2）'
（P.^2）'
[b,bint,r,rint,stats]=regress（R,f）;
程序结果：
b=
bint=
57.2602163.8024
0.01020.0630
-43.2247-9.9171
-0.0000084748762-0.0000009709903
0.37453.3205
stats=
0.970240.66560.0005
R110.53130.0366I26.5709P-4.7229*10^（6）I21.8475P2
参数
参数估计
置信区间
[57.2602163.8024]
1
[0.01020.0630]
2
[-43.2247-9.9171]
11
[-0.0000084748762-0.0000009709903]
22
[0.37453.3205]
R20.9702F40.6656p0.0005
因为相关系数R2充分接近1，并且F的概率p充分接近0，说明回归模型成立并且显著性水平高。
当消费者平均收入为4000，并且商品价格为6元时，我们代入回归模型可以得到商品的需求量为88.4791。
2.3.2结论与成果
1）在新超市，经理每月进货量为88件商品时比较合理。
2）通过建立回归模型，我们可以得到商品需求量与商品价格，消费者人均收入之间的关系。
应用在实践中，有助于商家更好地管理商品库存，减少库存冗余，从而创造更多利润。
4第二题
首钢董事长张先生想在工厂搬迁后扩大新厂子的规模，投资钢熔炼和钢拉伸两个项目进行改造，提高钢材生产效率，为此，他特意请来了投资专家，专家估计，投资到钢熔炼技术上可以赚回70倍的投资款，而投资到钢拉伸技术中的可以赚回66倍，同时，项目技术的投资也会带来一定的风险，一般来说风险损失会随着总投资和单项投资的增加而增加。
他准备投入5000万在钢熔炼和钢拉伸两个项目的改造上，他想知道如何分配投资资金才能让钢厂获得的收益最大，同时面临的风险最小。
4.1基本假设
1）假设此专家所作的估计充分正确。
2）假设风险损失只与总投资和单项投资有关，其他因素影响很小，可忽略不计
3）
4）假设除了全部5000元都用于钢熔炼和钢拉伸两个项目的改造外，不做任何其他投资与改变。
5）假设5000元不一定全部用完。
4.2参数设定：
分析：
因为有两个元素（钢熔炼与钢拉伸中分别投入资金）共同决定了收益值与风险值，所以进行以下设定:
1）钢熔炼项目所投入改造资金为x1
2）
3）钢拉伸项目所投入改造资金为x2
4）
5）所得收益为f1
6）
7）所面临风险损失为f2
4.3模型建立
通过题目要求提取出收益函数
maxf1（x）=70x1+66x2通过上网查询资料将风险损失函数定义为f2=0.02x120.01x220.04（x1x2）2并根据假设对x1、x2进行约束由于所求问题条件为收益最大同时风险最小，所以收益与风险权重比为1：
maxf1（x）70x166x2
minf2（x）0.02x120.01x220.04（x1x2）2x1x25000
x10
x20
线性加权构造目标函数
maxf0.5f1（x）0.5f2（x）
化为最小值问题
min（f）0.5f1（x）0.5f2（x）
4.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];
%
A=[11];
b=5000;
%lb=zeros（2,1）;
%x
初始猜测x1=1000,x2=1000;
满足A*x<
=b,即x1+x2<
=5000;
1>
=0,x2>
=0;
[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））^2
4.5程序运行结果
x=307.1426414.2862fval=-1.2211e+004exitflag=5f1=4.8843e+004
f2=2.4421e+004
fmincon函数分析:
[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为目标函数。
exitflag>
0表示优化结果收敛于解，exitflag=0表示优化超过了函数值的计算次数，exitflag<
0表示优化不收敛。
4.6结果分析与成果展示:
当给钢熔炼投入资金307元,给钢拉伸投入414元时,可获得最大收益约48843元,面临最小风险损失24421元。
通过此模型的建立，我们可以在知道精确回报数字的前提下计算出任意N个项目的资金投入方案，以达到最大收益与最小风险损失共存的目的。
有了这个辅助工具，我们能更加直观理性的衡量项目投入的资金对收益值和风险值的影响力。
升级会员 