全国大学生数学建模竞赛题目A题Word文档下载推荐.doc
- 文档编号:4640449
- 上传时间:2023-05-03
- 格式:DOC
- 页数:27
- 大小:666KB
全国大学生数学建模竞赛题目A题Word文档下载推荐.doc
《全国大学生数学建模竞赛题目A题Word文档下载推荐.doc》由会员分享,可在线阅读,更多相关《全国大学生数学建模竞赛题目A题Word文档下载推荐.doc(27页珍藏版)》请在冰点文库上搜索。
现进行成批生产,每批产量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)中,一当第i个零件的相对容差或第i个零件的容差:
确定以后,零件的生产成本就确定了,即这时不论零件的标定值在标定值取值范围中怎么取,生产成本都是不变的,为。
全部七个零件的相对容差只有种组合。
对所有108种情况可以逐次去计算;
实则上108种情况可以再简化。
从零件参数的标定值及不同容差等级的成本(元)如下表(符号/表示五此等级零件):
B等为
A等为
容差分为A、B、C三个等级,用与标定值的相对值表示,A等为,B等为,C等为.七个零件的参数标定值的标定值有一定的容许变化范围;
可见,A等(为)是不可能取到的。
用模拟的方法也可以确定,平均每个产品的生产成本一般在200-300之间,而损失费用为150上下。
如此,A等成本费用中200与500的费用太高,所增加的成本费用是损失费用无法弥补的,从而是不可能取到的。
于是容差等级总共不超过种。
考虑到100费用仍然也是难以弥补的,事实上,还应该是
种情况。
确定了一种容差方案后,上述问题实质上化为
(1.2)
,。
已经确定。
二、编程实现技术细节
1.容差等级、生产成本与标定值区间
R=[0.100.05];
%容差等级,最贵的一种舍去!
cost=[inf25%选定容差以后的生产成本费用
2050
50100
50inf
1025
inf25];
A=[0.0750.125
0.2250.375
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取定后,损失费用是个最小值优化问题。
其目标函数与约束条件是:
,已确定,而标定值是优化变量
已经可以计算。
其中已经没有零件参数的实际随机值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
相应的容差等级:
r=0.0500
0.0500
0.1000
零件总费用(成本+损失):
My=
4.2136e+005
每个零件成本与损失:
Cy=
275
Ly=
146.3603
程序求解运行时间:
T=
16.0051
(2).计算方法二:
用网格搜索方法
(3).计算方法三:
用随机模拟的方法
(4).计算方法四:
用模拟退火算法
(5).计算方法五:
用遗传算法
二.初始方案费用计算
1.原设计方案费用计算----用Matlab
见mydesign2
零件参数标定值:
x=
0.1000
0.3000
1.5000
16.0000
0.7500
r=0.0500
0.0500
3.0748e+003
零件成本与损失:
200
2.8748e+003
计算运行时间:
T=0.6351
2.原设计方案费用计算----用随机模拟
见mydesign3
用随机模拟法(n=10000次随机试验)开始求原设计方案费用。
请稍候---
x=0.1000
0.3000
1.5000
16.0000
0.7500
产品性能指标参数y分布的统计标准差:
sigmay=
0.1113
产品统计正品率:
p1=
0.1154
产品统计次品率:
p2=0.6242
产品统计费品率:
p3=0.2604
3.1678e+003
2.9678e+003
4.0799
模型检验与误差分析
对建模的结果,我们往往要作结果检验与敏感性分析和误差分析等,以提高论文的完整性和理论层次(可以参考上海交通大学与中国科大的获奖论文)。
1、敏感性分析
首先,在最优标定值点处,因为实际参数值
从而
即目标函数的梯度(绝对敏感性系数),除以函数均值,乘上每个参数的当前标定值,即为误差传递的敏感性系数。
从而第i个参数的相对敏感性系数为
又
grad_y=[24.1992-3.393211.5039-2.6496-1.3243-0.0623-0.9385]
或
敏感系数
1
2
3
4
5
6
7
估计值
1.2125
0.8501
0.8500
0.2124
1.0000
0.5000
其次,在最优标定值点处,固定六个参数,仅让其中一个参数改变,由此得到的函数值的改变情况比较图形为:
同样可见,函数y值的变化对参数x1,x3比较最为敏感,其次是x2,x4,x5,x7,而对x6是最不敏感的。
2、误差分析
(1)关于函数y的值的误差估计。
当标定值确定以后,由于零件实际的参数值带有随机性,因而函数y的值与均值也有一定随机性误差。
则对y取值的相对误差,也有用
来表示的。
由此可以估计出当标定值、容差等级确定后,产品的性能指数函数值有多大的随机性波动。
(2)关于函数y的值的非线性误差,以及所带来的对函数值y的分布误差的影响是不是会很显著的问题,可以见程序myhess。
当标定值达到允许的最大误差时,函数y的值的非线性误差为
=0.0338
其中是函数y在最优标定值点处的海森矩阵(即二阶导数矩阵)。
而非线性误差的相对值为
=2.26%。
所以,说经验函数值近似服从正态分布,其误差是很小的,即近似度是很高的。
模型改进
-----模型II建立
前面的优化问题模型
另外,用到质量损失函数
因而这使得该问题的目标函数高度非线性,以至于非线性规划问题怎样求解困难!
实际上,从另外一个角度来看,就是从产品的社会效用来说,性能指标差|y-y0|越接近于0.1的产品,其接近于正品,因而应该几乎没有什么损失;
另一方面,性能指标差|y-y0|越接近于0.3的产品,其接近于废品,因而损失应该几乎等同于废品,即它的损失费用应该接近于9000,而不是恒定在1000;
我们可以凭经验,合理地假设损失函数是y-y0的二次函数,即
由|y-y0|=0.1时,L(y)=1000,|y-y0|=0.3时,L(y)=9000可得K=105。
于是平均质量损失
其中,即标定值对应的产品质量性能指标的平均值,。
而
所以对应的修改优化模型是
(1.2)
,,,
这个优化问题是个确定性的,混合整数规划问题。
与分布函数没有关系。
求解稍简单一些。
求解结果见程序(modelII)。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 全国大学生 数学 建模 竞赛 题目