全国大学生数学建模竞赛题目A题.docx

全国大学生数学建模竞赛题目A题.docx

《全国大学生数学建模竞赛题目A题.docx》由会员分享，可在线阅读，更多相关《全国大学生数学建模竞赛题目A题.docx（25页珍藏版）》请在冰豆网上搜索。

全国大学生数学建模竞赛题目A题

1997年全国大学生数学建模竞赛题目

A题零件的参数设计

一件产品由若干零件组装而成，标志产品性能的某个参数取决于这些零件的参数。

零件参数包括标定值和容差两部分。

进行成批生产时，标定值表示一批零件该参数的平均值，容差则给出了参数偏离其标定值的容许范围。

若将零件参数视为随机变量，则标定值代表期望值，在生产部门无特殊要求时，容差通常规定为均方差的3倍。

进行零件参数设计，就是要确定其标定值和容差。

这时要考虑两方面因素：

一是当各零件组装成产品时，如果产品参数偏离预先设定的目标值，就会造成质量损失，偏离越大，损失越大；二是零件容差的大小决定了其制造成本，容差设计得越小，成本越高。

试通过如下的具体问题给出一般的零件参数设计方法。

粒子分离器某参数（记作）由七个零件的参数（记作）决定，经验公式为：

的目标值（记作）为1.50。

当偏离时，产品为次品，质量损失为1000（元）；当偏离时，产品为废品，质量损失为9000（元）；

零件参数的标定值有一定的容许变化范围；容差分为A、B、C三个等级，用与标定值的相对值表示，A等为，B等为，C等为.七个零件的参数标定值的容许范围，及不同容差等级的成本（元）如下表（符号/表示五此等级零件）：

标定值容许范围

C等

B等

A等

[0.075,0.125]

/

25

/

[0.225,0.375]

20

50

/

[0.075,0.125]

20

50

200

[0.075,0.125]

50

100

500

[1.125,1.875]

50

/

/

[12,20]

10

25

100

[0.5625,0.935]

/

25

100

现进行成批生产，每批产量1000个。

在原设计中，七个零件参数标定值为，，，，，，；容差均取最便宜的等级。

请你综合考虑偏离造成的损失和零件成本，重新设计零件参数（包括标定值和容差），并与原设计比较，总费用降低了多少?

B题截断切割

某些工业部门（如贵重石材加工等）采用截断切割的加工方式。

这里“截断切割”是指将物体沿某个切割平面分成两部分。

从一个长方体中加工出一个已知尺寸、位置预定的长方体（这两个长方体的对应表面是平行的），通常要经过6次截断切割。

设水平切割单位面积的费用是垂直切割单位面积费用的r倍，且当先后两次垂直切割的平面（不管它们之间是否穿插水平切割）不平行时，因调整刀具需额外费用e。

试为这些部门设计一种安排各面加工次序（称“切割方式”）的方法，使加工费用最少。

（由工艺要求，与水平工作台接触的长方体底面是事先指定的）详细要求如下：

1）需考虑的不同切割方式的总数。

2）给出上述问题的数学模型和求解方法。

3）试对某部门用的如下准则作出评价：

每次选择一个加工费用最少的待切割面进行切割。

4）对于e=0的情形有无简明的优化准则。

5）用以下实例验证你的方法：

待加工长方体和成品长方体的长、宽、高分别为10、14.5、19和3、2、4，二者左侧面、正面、底面之间的距离分别为6、7、9（单位均为厘米）。

垂直切割费用为每平方厘米1元，r和e的数据有以下4组：

a.r=1,e=0;b.r=1.5,e=0;c.r=8,e=0;d.r=1.5;2<=e<=15.对最后一组数据应给出所有最优解，并进行讨论。

模型建立

一．符号说明

：

产品零件参数的标定值；：

第i个零件的标定值；

：

第i个零件的标定值取值的上、下界；

：

产品零件参数的实际值；

：

产品性能参数的目标值；

：

产品性能参数的经验公式；即

：

产品性能参数的平均值；

：

产品性能参数的实际值；

：

第i个零件的相对容差（绝对值）；

：

第i个零件的容差，；

：

第i个零件参数取第j个容差等级时所需成本，；

第1，2，3个容差分别表示C，B，A等级；

：

变量，；,如果第i个零件参数取第j个容差等级时取值1，否则为零；

：

第i个零件参数（的实际值）的均方差；

：

产品质量性能参数（的实际值）的均方差；

：

产品的生产成本；

：

产品的损失费用；

：

产品的生产成本与损失费用之总和。

二．关于零件参数的假设

在第i个零件取定其标定值为后，由于在加工过程中存在许随机因素，刀具磨损，测量的误差等，因此，由中心极限定理知零件参数的实际值可看成是服从正态分布的随机变量，即

并且设七个零件的加工过程是相互独立的。

概率统计学告诉我们，如果某个随机变量服从正态分布

则由“3”原则，有

所以，当要求第i个零件取定其标定值为，第i个零件的容差为，则意味着

于是

即当第i个零件的标定值，容差等级（C，B，A，既第i个零件的相对容差））确定后，第i个零件参数的实际值所服从的正态分布就完全确定了。

三．关于产品参数的分布

当七个零件的标定值，容差等级（C，B，A，即每个零件的相对容差））确定后，为确定产品性能参数的分布规律与求产品质量的损失费用，现在讨论的分布情况。

由于产品性能参数的经验公式

非常复杂，直接得出的精确分布有困难。

对此有两种办法：

①用随机模拟的方法

在零件参数标定值的允许范围内任意取一组值和任意一组容差等级，产生n组相互独立的正态分布随机数其中七个零件参数中的每个都是服从正态分布

，

的随机数。

由经验公式得到产品性能参数的n组（n很大）对应的样本值

，

画出直方图，并用检验法检验是否服从正态分布，并确定的分布的有关特征数字。

例如，取n=10000时,的直方图与是否服从正态分布的检验

具体见ytest1和ytest2。

②用理论方法推导

在零件参数标定值的允许范围内任意取一组值和一组容差等级，记，表示产品性能参数的平均值。

产生零件参数的实际值，由经验公式得到产品性能参数的一组对应的值

因为

，

所以

即

因为，服从相互独立的正态分布，所以也服从正态分布，并且对上式两边取方差，得

于是产品性能参数也近似服从正态分布，即

至于近似服从正态分布的误差有多大？

我们后面再讨论。

三．目标函数

由原问题要求，建模的目标函数应包括产品零件的的生产成本与质量损失费用两部分。

①生产成本

当零件的标定值与零件的容差确定以后，零件的的生产成本就完全确定了，为

其中是第i个零件参数取第j个容差等级时所需成本，；，第1，2，3个容差分别表示C，B，A等级；而为变量，如果第i个零件参数取第j个容差等级时取值1，否则为零，并且

，

还要注意到对每个，对容差等级来说，不一定都能取到；

②质量损失费用

质量损失函数为

。

由上一段，服从正态分布，即

其密度函数为

所以（）

正品的概率；

次品的概率；

废品的概率。

再由标准正态分布分布函数与一般正态分布分布函数分布函数之间的关系：

有

同理

和

从而1个产品的平均质量损失费用

其中为标准正态分布分布函数。

于是整个优化问题模型为

（1.1）

s.t.

，，

，

并且其中

。

注记：

模型中的函数是标准正态分布分布的分布函数。

严格地说，这是带变量上下界约束与部分整数变量限制的、目标函数高度非线性的非线性规划问题。

模型求解

一．模型I计算思想

（1.1）

s.t.

，，，

并且

其中

。

模型中的函数是标准正态分布分布的分布函数。

这个带变量上下界约束与整数变量的、目标函数高度非线性的非线性规划问题怎样求解呢？

注意到（1.1）中，一当第i个零件的相对容差或第i个零件的容差：

确定以后，零件的生产成本就确定了，即这时不论零件的标定值在标定值取值范围中怎么取，生产成本都是不变的，为。

全部七个零件的相对容差只有种组合。

对所有108种情况可以逐次去计算；

实则上108种情况可以再简化。

从零件参数的标定值及不同容差等级的成本（元）如下表（符号/表示五此等级零件）：

标定值容许范围

C等

B等为

A等为

[0.075,0.125]

/

25

/

[0.225,0.375]

20

50

/

[0.075,0.125]

20

50

200

[0.075,0.125]

50

100

500

[1.125,1.875]

50

/

/

[12,20]

10

25

100

[0.5625,0.935]

/

25

100

容差分为A、B、C三个等级，用与标定值的相对值表示，A等为，B等为，C等为.七个零件的参数标定值的标定值有一定的容许变化范围；

可见，A等（为）是不可能取到的。

用模拟的方法也可以确定，平均每个产品的生产成本一般在200-300之间，而损失费用为150上下。

如此，A等成本费用中200与500的费用太高，所增加的成本费用是损失费用无法弥补的，从而是不可能取到的。

于是容差等级总共不超过种。

考虑到100费用仍然也是难以弥补的，事实上，还应该是

种情况。

确定了一种容差方案后，上述问题实质上化为

（1.2）

s.t.

其中

，。

已经确定。

二、编程实现技术细节

1.容差等级、生产成本与标定值区间

R=[0.100.05];%容差等级，最贵的一种舍去！

cost=[inf25%选定容差以后的生产成本费用

2050

2050

50100

50inf

1025

inf25];

A=[0.0750.125

0.2250.375

0.0750.125

0.0750.125

0.1251.875

1220

0.56250.9375];%标定值区间

2.怎么选定容差等级

用i1，i2，i3，i4，i5，i6，i7表示选定的容差等级，值都在1，2中选择，并且i1=2，i5=1，i7=2。

于是全部容差等级的循环是4重循环。

从而选定的第k个零件的容差等级为

r（ik）=R（ik）

而对应的成本费用为

cost（k,ik）

所以，这个产品的生产成本为

Cy=cost（1,i1）+cost（2,i2）+cost（3,i3）+cost（4,i4）+

cost（5,i5）+cost（6,i6）+cost（7,i7）

3.损失费用的最小化计算

当标定值x与容差等级r取定后，损失费用是个最小值优化问题。

其目标函数与约束条件是：

（1.2）

s.t.

其中

，已确定，而标定值是优化变量

，。

已经可以计算。

其中已经没有零件参数的实际随机值x了。

目标函数中的是正态分布的分布函数，Matlab中有现成的函数normcdf可以调用（比较normpdf）。

4.多元函数的偏导数（或梯度）怎么计算

多元函数的偏导数（或梯度）与海森矩阵怎么计算，Matlab中有现成的函数jacobian可以调用，具体见objgrad.m。

5.当时的计算程序编写方法

（1）．计算方法一：

用Matlab中函数Fmincon来计算

见mydesign1

求解结果：

零件参数标定值：

x=0.0750

0.3750

0.1106

0.1200

1.1309

12.0100

0.7977

相应