数学建模两个综合模型案例分析课件.ppt
《数学建模两个综合模型案例分析课件.ppt》由会员分享,可在线阅读,更多相关《数学建模两个综合模型案例分析课件.ppt(73页珍藏版)》请在三一办公上搜索。
1、2023/3/13,1,这时要考虑两方面因素:一是当零件组装成产品时,如果产品参数偏离预先设定的目标值,就会造成质量损失,偏离越大,损失越大;二是零件容差的大小决定了其制造成本,容差设计得越小,成本越高,试通过如下的具体问题给出一般的零件参数设计方法。粒子分离器某参数(记作y)由7个零件的参数(记作x1,x2,x7)决定,经验公式为,2023/3/13,2,y的目标值(记作y0)为1.50。当y偏离y00.1时,产品为次品,质量损失为1000(元);当y偏离y00.3时产品为废品,损失为9000(元)。零件参数的标定值有一定的容许变化范围;容差分为A、B、C三个等级,用与标定值的相对值表示,A
2、等为1%,B等为5%,C等为10%。7个零件参数标定值的容许范围及不同容差等级零件的成本(元)如表1(符号/表示无此等级零件)。,2023/3/13,3,表1,2023/3/13,4,现进行成批生产,每批产量1000个。在原设计中,7个零件参数的标定值为:x1=0.1,x2=0.3,x3=0.1,x4=0.1,x5=1.5,x6=16,x7=0.75;容差均取最便宜的等级。请你综合考虑y偏离y0造成的损失和零件成本,重新设计零件参数(包括标定值和容差),并与原设计比较,总费用降低了多少。,2023/3/13,5,二、模型的假设及符号说明(1)模型的假设 假设各零件的参数为随机变量,且是相互独立
3、的,它们都服从以零件标定值为均值,容差的三分之一为均方差的正态分布。(假设生产1000个批量产品时,每个零件参数只按一种容差等级进行生产,7个零件可以按不同的容差等级进行组合生产。,2023/3/13,6,考虑到题目所给条件,产品参数y偏离y00.1时,产品为次品,质量损失为1000元,y偏离y00.3时,产品为废品,损失为9000元,可假设产品参数y偏离目标值y0造成的单件产品质量损失函数L(y)与(y-y0)2成正比即(1)易见,2023/3/13,7,(2)符号说明 第i个零件的参数,它是随机变量;第i个零件的标定值,即xi的期望值;第i个零件的均方差;第i个零件的容差,即 第i个零件的
4、相对容差,即;,2023/3/13,8,y0:产品质量参数的目标值;y:产品质量参数,它是随机变量;y:y的均方差;ai:标定值xi0的取值下限(已知的)bi:标定值xi0的取值上限(已知的)x0=(x10,x20,x70)标定值向量t=(t1,t2,t7):相对容差向量,2023/3/13,9,三、问题的分析及数学模型(1)问题的分析 问题的目标函数是总费用函数最小,而总费用是由产品参数偏离目标值引起的质量损失费用和产品的成本费用两部分组成。一般来说,标定值设计不合理或容差设计得太大,会使产品参数远离目标值,造成质量损失;而容差设计得太小,又会增加零,2023/3/13,10,零件制造成本,
5、综合考虑产品质量损失费和成本费是问题的关键。由题目所给条件知,单件产品的零件成本取决于相对容差等级,第种零件的成本记作,于是单件产品的成本,即7个零件的总成本为(2)产品参数y由决定,,2023/3/13,11,记作y=f(x1,x2,x3,x4,x5,x6,x7),而xi是随机变量,且xiN(xi0,i2),所以产品参数y也是随机变量。进行成批生产时,平均每件产品的质量损失费用应该用损失函数L(y)的期望来度量,它取决于零件参数标定值x0和容差t,记为:,2023/3/13,12,由(2)式为了得到 的简单表达式,在 处对 作Taylor展开,并略去二阶及二阶以上项,有,2023/3/13,
6、13,其中 于是所以,2023/3/13,14,通过以上分析讨论可知,成批生产时平均每件产品的总费用可表示为(2)数学模型 问题的目标函数是总费用函数Z(x0,t)达到最小,且约束条件是标,2023/3/13,15,定值x0落在给定的容许范围内,所以该问题的数学模型为:(3),2023/3/13,16,四、模型的分析及求解 1.模型的分析 注意到在模型(3)中,决策变量x0是连续的,t是离散型的,而且t的取值共有23332=108种,对每种固定的t值,求解如下一系列子问题:(4),2023/3/13,17,得到最优解 和最优值,然后对108个t比较(5)得到全局最优解 和最优值。子问题(4)是
7、非线性规划问题,它的计算量是很大的,为了减少计算子问题的次数,在程序设计时,每解出一个子问题,由(5)式算出当前的最优值,对待求解的子问题,先判断 是否成立,若成立,则该子问题不必再解。,2023/3/13,18,2.模型的求解 按照上面对模型(3)的分析,下面设计在MATLAB中计算模型(3)的程序。为了方便起见,记,2023/3/13,19,分别建立函数文件v.m,u.m,fy.mfunction v=fun(x)v=1-0.36*x(2)0.56*x(4)(-0.56);function u=fun(x)u=1-2.62*v(x)1.5*(x(4)/x(2)1.16;function f
8、y=fun(x)fy=174.42*(x(1)/x(5)*(x(3)/(-x(1)+x(2)0.85*sqrt(u(x)/(x(6)*x(7);,2023/3/13,20,分别建立函数u(x)对x2,x4的偏导数函数文件uu2.m,uu4.mfunction uu2=fun(x)uu2=0.792288*v(x)0.5*x(2)(-1.6)*x(4)0.6+3.0392*v(x)1.5*x(2)(-2.16)*x(4)1.16;function uu4=fun(x)uu4=(-0.792288)*v(x)0.5*x(2)0.6*x(4)(-0.4)-3.0392*v(x)1.5*x(2)(-1
9、.16)*x(4)0.16;,2023/3/13,21,分别建立函数fy(x)对x1,x2,x7的偏导数函数文件y11.m,y21.m,y71.mfunction y11=fun(x)y11=fy(x)*(1/x(1)+0.85/(x(2)-x(1);function y21=fun(x)y21=fy(x)*(-0.85)/(x(2)-x(1)+0.5*(1/u(x)*uu2(x);function y31=fun(x)y31=fy(x)*(0.85/x(3);,2023/3/13,22,function y41=fun(x)y41=fy(x)*0.5*(1/u(x)*uu4(x);funct
10、ion y51=fun(x)y51=fy(x)*(-1/x(5);function y61=fun(x)y61=fy(x)*(-0.5/x(6);function y71=fun(x)y71=fy(x)*(-0.5/x(7);,2023/3/13,23,建立子问题(4)的目标函数文件ch721.mfunctionch721.g=ch721(x,p2,p3,p4,p6,p7)ch81=100000*(fy(x)-1.5)2+(100000/9)*(y11(x)*x(1)*0.05)2+(y21(x)*x(2)*p2)2+(y31(x)*x(3)*p3)2+(y41(x)*x(4)*p4)2+(y
11、51(x)*x(5)*0.1)2+(y61(x)*x(6)*p6)2+(y71(x)*x(7)*p7)2;g=-x(1);,2023/3/13,24,用MATLAB语言设计计算模型(7-8)的程序,可存为M文件ch722.mtic%启动计时器p2=0.1;p3=0.1;p4=0.1;p6=0.1;p7=0.05;%原设计方案的相对容差x0=0.1,0.3,0.1,0.1,1.5,16,0.75;%原设计方案的标定y0=fy(x0);%原设计方案产品的参数值,2023/3/13,25,Q0=ch81(x0,p2,p3,p4,p6,p7);%原设计方案的质量损失费Z0=Q0+200;%原设计方案的
12、总费用Zmin=Z0;%原设计方案的总费用作为优化初值t1(2)=0.05;%对各等级零件相对容差赋值t2(2)=0.05;t2(3)=0.1;t3(1)=0.01;t3(2)=0.05;t3(3)=0.1;,2023/3/13,26,t4(1)=0.01;t4(2)=0.05;t4(3)=0.1;t5(3)=0.1;t6(1)=0.01;t6(2)=0.05;t6(3)=0.1;t7(1)=0.01;t7(2)=0.05;c2(2)=50;c2(3)=20;c3(1)=200;c3(2)=50;c3(3)=20;c4(1)=500;c4(2)=100;c4(3)=50;c6(1)=100;c
13、6(2)=25;c6(3)=10;c7(1)=100;c7(2)=25;,2023/3/13,27,for i2=2:3for i3=1:3for i4=1:3for i6=1:3for i7=1:2C=75+c2(i2)+c3(i3)+c4(i4)+c6(i6)+c7(i7);%计算单件产品的成本p2=t2(i2);p3=t3(i3);p4=t4(i4);p6=t6(i6);p7=t7(i7);,2023/3/13,28,Zminif C=Zminx0=0.1,0.3,0.1,0.1,1.5,16,0.75;%优化的初始值v1b=0.075,0.225,0.075,0.075,1.125,1
14、2,0.5625%变量x的下界vub=0.125,0.375,0.125,0.125,1.875,20,0.935%变量x的上界options=;%优化中的参数均采用缺省值,2023/3/13,29,jacob=;%将空阵赋给雅可比阵x,options=constr(ch81,x0,options,vlb,vub,jacob,p2,p3,p4,p6,p7);%计算非线性规化问题(7)Z=C+cs97(x,p2,p3,p4,p6,p7);%当t值给定时的最小总费用if Z=ZminZmin=Z;%当前的最小总费用j2=i2;j3=i3;j4=i4;j6=i6;j7=i7;Xmin=x;,2023
15、/3/13,30,Cmin=C;Qmin=ch81(x,p2,p3,p4,p6,p7);y=fy(x)endendendendendendend,2023/3/13,31,T=0.05,t2(j2),t3(j3),t4(j4),0.1,t6(j6),t7(j7)%输出最优设计方案的相对容差Xmin%输出最优设计方案的标定值Cmin%输出最优设计方案的成本费Qmin%输出最优设计方案的质量损失费Zmin%输出最优设计方案的总费用,2023/3/13,32,y%输出最优设计方案产品的参数值 Q0%输出原设计方案的质量损失费 Z0%输出原设计方案的总费用 y0%输出原设计方案产品的参数值 toc%关
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 两个 综合 模型 案例 分析 课件
链接地址:https://www.31ppt.com/p-3498224.html