实验十五零件参数的设定

1、实验十五零件参数的设定【实验目的】1. 了解随机模拟法（即 Monte Carlo法）的基本原理。2. 学习随机模拟变量产生的基本方法，初步培养随机模拟的建模思想。3. 学习掌握MATLAB件中随机模拟的相关命令。【实验内容】一件产品由若干个零件组装而成，标志产品性能的某个参数取决于这些零件的参数。零件参数包括标定值和容差两部分。进行成批生产时，标定值表示一批零件该参数的平均值，容差则给出了参数偏离其标定值的容许范围。若将零件参数视为随机变量，则标定值代表期望值，在生产部门无特殊要求时，容差通常视为标准差的3倍。粒子分离器某参数（记作y）由7个零件的参数（记作洛，x2，x7） 决定，经验公式为

2、32/ r “， c cX4、0.56/X4、1.161 2.62 1 0.36（）（）",X1、/ X3 0.85x2x2y = 174.42 （）（） VX5 X2 X1IX6X7当各零件组装成成品时，如果产品参数偏离预先设定的目标值，就会造成 质量损失，偏离越大，损失就越大。y的目标值（记作y°）为1.50，当y偏离y 0.1时，产品为次品，质量损失为1000 （元）；当y偏离y。0.3时，产品为废品，质量损失为9000 （元）。给定某设计方案7个零件参数标定值及容差，如表1所示：容差分为A、B、C三个等级，用及标定值的 相对乘积值表示，A等为1% B等为 5% C等

3、为15%表1零件参数标定值和容差X1X2X3X4X5X6X7标定值0.20.30.10.11.5160.75容差BBBCCBB求每件产品的平均损失。【实验准备】在现实生活中，有大量问题由于模型中随机因素很多，很难用解析式模型来进行描述求解，这时就需要借助模拟的方法。 随机模拟法也叫MonteCarlo法，它是用计算机模拟随机现象，通过大量仿真实验，进行分析推 断，特别是对一些复杂的随机变量，不能从数学上得到它的概率分布，而 通过简单的随机模拟便可得到近似解答。象这类大容量的仿真实验，如果 用实物来做，需要大量人力物力且可能无法实现，但如果我们有了问题的 数学模型，用计算机模拟就轻而易举了。由于

4、Monte Carlo法计算量大，精度不是很高，因而适合一些用解析方法或常规数值方法难以解决问题的 低精度求解，或用于对一些计算结果的验证。1. 随机模拟的一些基本概念自然界发生的现象可分为两类，一类现象在一定条件下发生的结果是 完全可以预知的，称为必然现象。另一类现象发生的结果在事先是无法准 确预知的，称为 偶然现象或随机现象。下面两个试验都是随机现象：试验一：有10枚均匀硬币，随手抛在地上，有几枚正面向上？试验二：按身份证号码随意挑10个中国女子，他们的平均体重是多少？ 尽管随机现象的发生结果是不确定的，但还是有一定的规律可循：试验一 中正面向上的枚数一定是 010，5枚向上的可能性比8枚

5、向上的可能性 要大；试验二中平均体重基本在 40kg到70kg之间，且在45kg左右的可 能性比65kg左右的可能性要大。一个随机事件A发生的可能性的大小，用一个介于0及1之间的数表示，称为A的概率，记为P(A)。概率的意义在类似的现象大量重复发生时 会表现出来。比如，在试验一中若P (5枚向上)=0.25，那么意味着“若 把试验一做 100遍，大致有 25次左右出现 5枚向上的情况。 ”在随机现象中，变量的取值往往是不确定的，称为 随机变量 。描述随 机变量取各种值的概率函数称为 概率分布 。对于随机变量，通常主要关心 它的两个主要数字特征： 数学期望 用于描述随机变量的平均值， 方差 和标

6、 准差用于描述随机变量分布的差异程度。另外， 协方差和相关系数用于描 述两个随机变量的线性关联程度。 (数字特征的定义跟前面实验定义的一 致，且均能在概率统计的书籍中查找相关定义)随机变量的分布，根据其取值特点不同主要分为离散型和连续型两类。若用变量 表示试验一“正面向上次数”，其取值可能为0, 1, 2， 10(离散点集) ，则为离散型随机变量。典型的离散型分布有二项分布、 Poisson 分布等。若用变量 表示试验二中“平均体重”，其取值可能为30， 80中的任何值，则为连续型随机变量。典型的连续型分布有均匀分布、 正态分布、指数分布、2分布、t分布、F分布等。2、模拟随机数的产生为了产生

7、具有一定分布的随机数，一般采用一定的生成程序。首先要 有一个等概率密度随机数发生器，一般计算机上都有专门的程序，产生01 之间等概率密度分布的随机数， 使用时直接调用即可； 此01之间的 随机数进行一定的数字转换即可获得所要求的随机数，怎样进行数字转换 则视所要求的分布函数来定。假定将0，1 区间的均匀随机数记作 R，则a，b 区间的均匀随机数可按下述公式由 0， 1区间的均匀随机数产生：x = a + R ( b a )(1)逆转换法这是求概率分布的逆函数从而产生随机数的方法。因概率分布函数F(x)为定义在0，1 区间的单调递增函数，设 R为区间0，1 的均匀 随机变量，令F(x)= R，只

8、要求出逆函数x = F 1(R)，x即为具有概率分 布函数F(x)的随机数。组合法组合法是利用某些容易产生随机数数列的随机变量，通过组合得到所要求的随机变量的一种方法。例：产生泊松分布的模拟随机数列如果相继两个事件出现的间隔时间为负指数分布，则在某一时间间隔 内事件出现的次数服从泊松分布。根据此关系，可以用负指数分布的随机 变量来组合产生泊松分布的随机数序列。设yi , y2，yn为参数 的负指数分布的随机数序列，因为有1,cy InRj(2)所以将yi值按序累加，使得满足关系式：< 1 <(3)则求得的x就是参数 的泊松分布。的随机数。近似法这种方法一般用于随机变量的分布函数无法

9、求出的情形。此时可运用大数定理，当样本数量趋于无穷时，样本平均值趋向于总体平均值，它是数字特征随机模拟的理论根据。3. 及随机数相关的MATLA命令max min最大值，最小值sum各兀素和mea n均值cumsum各元素累计和media n中值prod各元素积std标准差cumprod各元素累计积cov协方差矩阵bar直方图corrcoef相关系数矩阵hist数据分组及直方数据分析函数 max min, mean, median, std , cov, sum, prod , cumprod等标准用法都是对列状数据进行的。bar (Y)作向量Y的直方图；bar(X，Y)作向量Y相对于X的直方

10、图；hist (X，k)将向量X中数据等距分为k组，并作出直方图，缺省值为 k = 10;有关它们更详细的内容可查阅帮助文件。随机数生成采用下面命令形式：R = rand( m , n ) 生成0, 1 区间上均匀分布的 m行n列随 机矩阵；R = randn( m , n )生成标准正态分布的 m行n列随机矩阵；R = randperm( N ) 生成1, 2,，N的一个随机排列；R = unidrnd( N , m , n )生成1，2，N的等概率 m行n列随机矩阵；R = unifrnd( a , b , m , n )生成a，b区间上均匀分布的m行n列随机矩阵；R = normrnd(

11、 mu , sigma , m , n )生成均值为 mu 标准差为sigma的m行n列正态分机随机数矩阵；R = binornd( k , p , m , n )生成参数为k，p的m行n列正态分机随机数矩阵，它模拟在k次重复试验中某事件(发生概率为p)出现的次数；R = mvnrnd( mu , sigma , m )生成 n 维正态分布数据，这里 mu为n维均值向量，sigma为n阶协方差矩阵(它必须是正定的)，R 为mx n矩阵，每行代表一个随机数。R = poissrnd (mu , m , n )生成均值为 mu的m行n列泊松分布的随机数矩阵；可以通过帮助文件查阅上述命令的详细内容。

12、【实验方法及步骤】1. MATLAB令的基本用法下面用几个例子来予以说明：>> data二13 76 356;11 89 278;10 86 302;8 92 362;15 69 311;14 83 299;1173 336;>> max(data)ans =1592 362>> mean( data)ans =11.7143 81.1429 320.5714>> sum(data)ans =825682244>> std(data)230 / 11>>>>>>>>ans =2.4300

13、8.6300 31.4211prod(data)ans =1.0e+017 *0.00000.00023.3805cov(data)%将三列看成三个随机变量ans =5.9048 -15.1190 -22.9762-15.1190 74.4762 -34.4286-22.9762 -34.4286 987.2857 corrcoef(data)%将三列看成三个随机变量ans =1.0000 -0.7210 -0.3009-0.72101.0000 -0.1270-0.3009 -0.12701.0000bar(data)% 作向量data的直方图2引例问题的分析求解在这个问题中，主要的困难是产

14、品的参数值 y 是一个随机变量，而由 于 y 及各零件参数间是一个复杂的函数关系，无法解析地得到 y 的概率分 布。本实验可以考虑采取随机模拟的方法计算。其基本思路是：用计算机 模拟工厂生产大量“产品” (如 1000 件)，计算产品的总损失，从而得到 每件产品的平均损失。对于大样本容量的随机变量，我们可以假设 7 个零件参数均服从正态 分布。根据题设里标定值和容差的定义，我们可以得到 7 个零件参数所对 应正态分布的均值及方差：x1 N(0.1 , (0.005/3)2), x2 N(0.3, (0.005)2) , x3 N(0.1,(0.005/3)2)2 2 2x4 N(0.1,(0.

15、005)2)， x5 N(1.5,(0.075)2)， x6 N(16,(0.8/3)2)2x7 N(0.75,(0.0125)2)下面在脚本文件 eg6_1.m 中产生 1000 个对零件 7 个参数的随机数， 通过随机模拟法求解零件平均损失的近似解。% 脚本 eg6_1.m 文件clear;% 清除内存变量mu=0.1,0.3,0.1,0.1,1.5,16,0.75; sigma=0.005/3,0.005,0.005/3,0.005,0.075,0.8/3,0.0125;for i=1:7 x(:,i)=normrnd(mu(i),sigma(i),1000,1);endp=(1-2.6

16、2*(1-0.36*(x(:,4)./x(:,2).A(-0.56).A1.5.*(x(:,4)./x(:,2).".16)./x(:,6)./x(:,7); q=(x(:,1)./x(:,5).*(x(:,3)./(x(:,2)-x(:,1).985;y=174.42*q.*p40.5;d=abs(y-1.5);%及目标值差的绝对值f二sum(9000*(d>0.3)+1000*(dv=0.3).*(d>0.1)/1000%求 零件的平均损失%注意此处使用的是数组的点乘、点除、和点幂运算。>> f =2948【结果分析】第一次运行脚本文件eg6_1.m时得到

17、的解为2948，是否每次运行结果都一致呢？很显然，每次运行的结果应该不同，并且有一定的差别，因为 我们是按计算机内部算法取1000个正态分布的随机模拟数，下表是连续10次运行的结果表2模拟1000对零件参数运行次数12345678910f (元)2897313330212894296728842873289629662918下面我们加大参数随机模拟的容量，提高两个数量级，取100000，同样我们取10次运行结果作成表2：表2模拟100000对零件参数运行次数123456789f (e+003)2.90852.92582.91522.89822.93102.91422.90832.91042.9

18、119这时，我们可以观察到，零件平均损失费用在2910附近波动，且波动辐度较小容量时小很多，此时我们可以确认所得的解是比较接近零件平 均损失的真实值。通过该实验也验证，随机模拟在很多实际问题的求解中 能够取得比较理想的效果。【练习及思考】1. 一个加油站服务员每天工作 8小时，工资为15元/天。要求加油的汽车按 =35辆/小时的泊松流到达。每个服务员分别为一辆汽车加油，又 每服务一辆汽车后，加油站盈利1.25元。设每辆汽车的加油时间为负指数分布，1/= 8min。当等待加油的汽车超过两辆时，后来的汽车就不再排队等待而离去。试用模拟方法确定该加油站合理的服务员人数。2. 已知零件C由零件A和零件B连接而成，已知A、B的长度均为随机变 量，具体数值如下表。试抽取 100个样本以计算C的平均长度。零件A的长度56789概率0.070.