数学建模课后作业

数学实验与数学建模

第四章

1、有10个同类企业的生产性固定资产年平均价值和工业总产值资料如下：

（1）说明两变量之间的相关方向； （2）建立直线回归方程； （3）计算估计标准误差；

（4）估计生产性固定资产（自变量）

为1100万元时的总资产（因变量）

的可能值。

解：由表格易知：工业总产值是随

着生产性固定资产价值的增长而

增长的，而知之间存在正向相关性。

用MATLAB 可解：

运行结果为：b =395.5670 运行图为：

0.8958

stats =

1.0e+04 *

0.0001 0.0071 0.0000 1.6035 industry =

1.0e+03 *

1.0220 1.0680

1.1160 1.1663

1.2188 1.2736

construction =

1.0e+03 *

1.2190 0.3965

0.3965 0.3965

0.3965 0.3965

ans = 395.5670 0.8958

2、设某公司下属10个门市部有关资料如下： （1）、确定适宜

的 回归模型；

（2）、计算有关指

标，判断这三种经济

现象之间的紧密程

度。

（1） 设销售利润率（%）为y ，流通费用水平（%）为x2，职工平均销售额（万元）为x3

回归模型y=a1+a2*x1+a3*x2

利用MATLAB ：x1=[12.6 10.4 18.5 3.0 8.1 16.3 12.3 6.2 6.6 16.8; 6 5 8 1 4 7 6 3 3 7; 2.8 3.3 1.8 7.0 3.9 2.1 2.9 4.1 4.2 2.5]';

X = [ones(size(x1(:,1))),x1(:,2:3)]; Y = x1(:,1); [b,bint,r,rint,stats] = regress(Y ,X,0.05); b,bint,stats

运行结果为：b =-6.7691 2.9070 0.9578

bint = -15.7285 2.1902

2.0138

3.8003 -0.3676 2.2832

stats =0.9823 194.2113 0.0000 0.6002

所以，相关系数为0.9823 分析可得 y=-6.7691+2.9070*x1+0.9578*x2 （2） 利用MATLAB ：stepwise(X,Y ,[],0.05)

运行结果为：

此时，y 与x2线性关系紧密，而与x3的线性关系不是很密切，即使没有x3,y=-门市部编号 职工平均销售额（万元） 流通费用水平（%） 销售利润率（%）

1 6 2.8 12.6

2 5 3.

3 10.4

3 8 1.8 18.5

4 1 7.0 3.0

5 4 3.9 8.1

6 7 2.1 16.3

7 6 2.9 12.3

8 3 4.1 6.2

9 3 4.2 6.6 10 7 2.5 16.8

0.386+2.293*x2,相关系数为：0.975，已经很高了。

即，销售利润与职工平均销售额关系密切，销售利润与流通费用水平关系不是很密切.

第五章

1、比较5种品牌的合成木板的耐久性，对每个品牌取4个样品作摩擦实验测量磨损量，的以下数据：

品牌A 2.2 2.1 2.4 2.5 品牌B 2.2 2.3 2.4 2.6 品牌C 2.2 2.0 1.9 2.1 品牌D 2.4 2.7 2.6 2.7 品牌E 2.3 2.5 2.3 2.4

（1） 他们的耐久性有无明显差异？

（2） 有选择的做两品牌的比较，能得出什么结果? 运行代码为：

X=[2.2,2.1,2.4,2.5;2.2,2.3,2.4,2.6;2.2,2.0,1.9,2.1;2.4,2.7,2.6,2.7;2.3,2.5,2.3,2.4]; p=anova1(X)

运行结果为：

因为p=0.5737>0.05所以可认为五个品牌合成木板的耐久性无显著差异。

由图可知：

品牌A 的磨损量较小，品牌B 的磨损量较大。所以，品牌A 的质量最好。

2、将土质基本相同的一块耕地分成5块，每块又分成均等的4小块。在每块地内把4个品种的小麦分种在4小块内，每小块播种量相同，测得收获量如下：

考察地块和品种对小麦的收获量有无显著影响？并在必要时作进一步比较。

利用MATLAB求解：

function anova_2

fm1=[32.3 34.0 34.7 36.0 35.5;33.2 33.6 36.8 34.3 36.1;30.8 34.4 32.3 35.8 32.8;29.5 26.2 28.1 28.5 29.4;];

p=anova2(fm1,2); display(p);

运行结果为：

p = 0.7770 0.0121 0.9393

3、为了研究合成纤维收缩率和拉伸倍数对纤维弹性的影响进行了一些实验。收缩率取0,4,8,12,四个水平；拉伸倍数取460,520,580,640 四个水平，对二者的每个组合重复做两次试

（1）收缩率，拉伸倍数及其交互作用对弹性有无显著影响？

（2）使弹性达到最大的生产条件是什么？

利用MATLAB求解：

X=[71,72,75,77;73,73,73,75;73,76,78,74;75,74,77,74;76,79,74,74;73,77,75,73;75,73,70,69;73,72, 71,69]; p=anova2(X,2)

运行结果为：p = 0.1363 0.0000 0.0006

P1=0.136>0.05,则接受假设1，因素拉伸倍数对指标无显著影响；

p2=0.0000<0.01,则拒绝假设2，因素收缩率对指标的影响显著；

P3=0.0006<0.01，则拒绝假设12，拉因素伸倍数和因素收缩率对指标的交互影响非常显