资源描述
单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,工艺系统热变形对加工精度的影响,工艺系统热变形对加工精度的影响,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,工艺系统热变形对加工精度的影响,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,工艺系统热变形对加工精度的影响,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,工艺系统热变形对加工精度的影响,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,二级,三级,四级,五级,*,专题,4,工艺系统热变形对加工精度,的影响及误差统计分析,目录,一、概述,二、工件热变形对加工精度的影响,三、刀具热变形对加工精度的影响,四、机床热变形对加工精度的影响,五、减少工艺系统热变形对加工精度影响的措施,六、误差统计分析,附录:误差分析的,MATLAB,实现,一、概 论,概念,工艺系统的热源,工艺系统的热平衡和温度场,1,、概念,:,在机械加工过程中,工艺系统会受到各种热的影响而产生变形,一般称为,热变形,。,这种变形将破坏刀具与工件的正确几何关系和运动关系,造成工件的加工误差。,热变形对加工精度影响比较大,特别是在精密加工和大件加工中,热变形所引起的加工误差通常会占到工件加工总误差的,40,70,。,一、概 论,*,工艺系统热变形不仅影响加工精度,而且还影响加工效率,因为为减少受热变形对加工精度的影响,通常需要:,预热机床以获得热平衡;,或降低切削用量以减少切削热和摩擦热;,或粗加工后停机以待热量散发后再进行精加工;,或增加工序(使粗、精加工分开)等等。,目前,无论在理论上还是在实践上都有许多问题尚待研究解决。,一、概 论,2,、工艺系统的热源,热的传递方式,热总是由高温处向低温处传递的。热的传递方式有三种:,传导,对流,辐射,一、概 论,一、概 论,热传导,一、概 论,热对流,一、概 论,热辐射,引起工艺系统变形的热源可分为,:,内部热源,主要指切削热和摩擦热,它们产生于工艺系统内部,其热量主要是以热传导的形式传递的。,外部热源,主要是指工艺系统外部的、以对流传热为主要形式的环境温度(它与气温变化、通风、空气对流和周围环境等有关)和各种辐射热(包括由阳光、照明、暖气设备等发出的辐射热)。,一、概 论,内部热源,:,切削热,切削热是切削加工过程中最主要的热源,对工件加工精度的影响最为直接,在切削(磨削)过程中,消耗于切削层的弹、塑性变形能及刀具、工件和切屑之间摩擦的机械能,绝大部分都转变成了切削热。,一、概 论,影响切削热传导的因素,切削热,Q,(,J,),的大小计算,式中,Pz -,主切削力(,N,),v-,切削速度(,m,min,),t-,切削时间(,min,),影响切削热传导的主要因素:工件、刀具、夹具、机床等的材料的,导热性能,,以及周围介质的情况。,若工件材料热导率大,则由切屑和工件传导的切削热较多;,若刀具材料热导率大,则从刀具传出的切削热也会较多。,各类切削加工方法热传导情况,车削,切屑所带走的热量最多,可达,50,80,(切削速度越高,切屑带走热量占总切削热百分比就越大),传给工件热量次之(约,30,),传给刀具的热量很少,一般不超过,5,铣、刨削加工,传给工件热量在总切削热的,30,以下,钻削和卧式镗孔,大量的切屑滞留在孔中,传给工件的热量就比车削时要高,磨削时磨屑小,带走的热量也少(约为,4,),大部分热量(,84,左右)传入工件,使磨削表面的温度高达,800,1000,左右,磨削热既影响工件的加工精度,又影响工件的表面质量。,内部热源,:,摩擦热,工艺系统中的摩擦热,主要是机床和液压系统中运动部件产生的,如电动机、轴承、齿轮、丝杠副、导轨副、离合器、液压泵、阀等各运动部分产生的摩擦热。,摩擦热比切削热少,但摩擦热在工艺系统中是,局部发热,引起局部温升和变形,,破坏了系统原有的几何精度,对加工精度也会带来严重影响。,一、概 论,外部热源,如:,需要昼夜连续加工时,昼夜温度不同,引起工艺系统的热变形就不一样,从而影响了加工精度,照明灯光、加热器等对机床的热辐射,往往是局部的,日光对机床的照射不仅是局部的,而且不同时间的辐射热量和照射位置也不同,因而会引起机床各部分不同的温升和变形,这在大型、精密加工时尤其不能忽视。,一、概 论,3,、工艺系统的热平衡和温度场,一、概 论,工艺系统受各种热源的影响,其温度会逐渐升高。同时,它们也通过各种传热方式向周围散发热量。,热 平 衡,当单位时间内传入和散发的热量相等时,工艺系统达到了热平衡状态。,而工艺系统的热变形也就达到某种程度的稳定。,一、概 论,温度场,稳态温度场,热平衡,物体中各点的温度,分布称为,温度场,T=,f,(,x,y,z,t),当物体未达热平衡时,各点温度不仅是坐标位置的函数,也是时间的函数。这种温度场称为,不稳态温度场,物体达到热平衡后,各点温度将不再随时间而变化,只是其坐标位置的函数。这种温度场称为,稳态温度场,不稳态温度场,一、概 论,机床在开始工作的一段时间内,其温度场处于不稳定状态,其精度也是很不稳定的,工作一定时间后,温度才逐渐趋于稳定,其精度也比较稳定。,因此,,精密加工应在热平衡状态下进行。,在生产中,必须注意:,二、工件热变形对加工的影响,(一)概述,(二)工件比较均匀受热,(三)工件不均匀受热,(一)概述,在工艺系统热变形中,机床热变形最为复杂,工件、刀具的热变形相对来说要简单一些,使工件产生热变形的热源,主要是,切削热,。,对于精密零件,周围,环境温度,和局部受到日光等,外部热源,的辐射热也不容忽视。,二、工件热变形对加工的影响,二、工件热变形对加工的影响,(,二,),工件均匀受热,对于一些形状简单、对称的零件,如轴、套筒等,加工时(如车削、磨削)切削热能较均匀地传入工件,工件热变形量可按下式估算:,L=Lt,式中,工件材料的热膨胀系数,单位为,1/,;,L,工件在热变形方向的尺寸,单位为,mm,;,t,工件温升,单位为。,均匀热变形对工件加工精度要求的影响,加工盘类和长度较短的销轴、套类零件时,由于走刀行程很短,可以忽略,车削较长工件时,由于温升逐渐增加,工件直径随之逐渐胀大,因而车刀的背吃刀量将随走刀而逐渐增大,工件冷却收缩后外圆表面就会产生圆柱度误差,当工件以两顶尖定位,工件受热伸长时,如果顶尖不能轴向位移,则工件受顶尖的压力将产生弯曲变形,对加工精度产生影响。宜采用弹性或液压尾顶尖。,二、工件热变形对加工的影响,一般,工件热变形在精加工中影响比较严重,特别是长度长而精度要求很高的零件,如:磨削丝杠。若丝杠长度为,2m,,每磨一次其温度相对于机床母丝杠就升高约,3,,则丝杠的伸长量为,0.07mm,,而,6,级丝杠的螺距累积误差在全长上不允许超过,0,02mm,,由此可见热变形的严重性。,二、工件热变形对加工的影响,(三)工件不均匀受热,二、工件热变形对加工的影响,在,刨削、铣削、磨削,加工,平面,时,,工件单面受热,,上下平面间产生,温差,,导致工件向上凸起,凸起部分被工具切去,加工完毕冷却后,加工表面就产生了中凹,造成了几何形状误差。,减少受热引起的误差对加工精度的影响的措施,对于大型精密板类零件,工件单面受热引起的误差对加工精度的影响是很严重的。为了减少这一误差,通常采取的措施是:,在切削时使用充分的切削液以减少切削表面的温升;,也可采用误差补偿的方法,即在装夹工件时使工件上表面产生微凹的夹紧变形,以此来补偿切削时工件单面受热而拱起的误差。,二、工件热变形对加工的影响,三、刀具热变形对加工精度的影响,刀具热变形主要是由切削热引起的。通常传入刀具的热量并不太多,但由于热量集中在切削部分,以及刀体小,热容量小,故仍会有很高的温升。,例如车削时,高速钢车刀的工作表面温度可达,700,800,,而硬质合金刀刃可达,1000,以上。,(一)概述,(二)刀具不同切削状态下的变形,连续切削时,,刀具的热变形在切削初始阶段增加很快,随后变得较缓慢。经过不长的时间后(约,10,20min,)便趋于热平衡状态。此后,热变形变化量就非常小。刀具总的热变形量可达,0.03,0.05mm,。,三、刀具热变形对加工精度的影响,间断切削时,,由于刀具有短暂的冷却时间,故其热变形曲线具有热胀冷缩双重特性,且总的变形量比连续切削时要小一些,J,最后趋于稳定在,范围内变动。,当切削停止后,刀具温度立即下降,开始冷却较快,以后逐渐减慢。,三、刀具热变形对加工精度的影响,加工大型零件,刀具热变形往往造成几何形状误差。如车长轴时,可能由于刀具热伸长而产生锥度(尾座处的直径比主轴箱附近的直径大),为了减小刀具的热变形,应合理选择切削用量和刀具几何参数,并给以充分冷却和润滑,以减少切削热,降低切削温度。,三、刀具热变形对加工精度的影响,在生产中,应注意,:,四、机床热变形对加工精度的影响,1,、概念,机床在工作过程中,受到内外热源的影响,各部分的温度将逐渐升高。由于各部件的热源不同、分布不均匀,以及机床结构的复杂性,因此不仅各部件的温升不同,而且同一部件不同位置的温升也不相同,形成不均匀的温度场,使机床各部件之间的相互位置发生变化,破坏了机床原有的几何精度而造成加工误差。,2,、机床热态几何精度,机床空运转时,各运动部件产生的摩擦热基本不变。运转一段时间之后,各部件传入的热量和散失的热量基本相等,即达到热平衡状态,变形趋于稳定。机床达到热平衡状态时的几何精度称为热态几何精度。,在机床达到热平衡状态之前,机床几何精度变化不定,对加工精度的影响也变化不定。因此,精密加工应在机床处于热平衡之后进行。,四、机床热变形对加工精度的影响,不同机床类型的热变形对加工精度的影响也不相同。,车、铣、钻、镗类机床,主轴箱中的齿轮、轴承摩擦发热,润滑油发热是其主要热源使主轴箱及与之相连部分如床身或立柱的温度升高而产生较大变形。,例如车床主轴发热使主轴箱在垂直面内和水平面内发生偏移和倾斜。,四、机床热变形对加工精度的影响,四、机床热变形对加工精度的影响,a),车床受热变形形态,b,)温升与热变形,四、机床热变形对加工精度的影响,四、机床热变形对加工精度的影响,立式铣床受热变形形态,四、机床热变形对加工精度的影响,外圆磨床受热变形形态,四、机床热变形对加工精度的影响,导轨磨床受热变形形态,对于不仅在水平方向上装有刀具,在垂直方向和其它方向上也都可能装有刀具的自动车床、转塔车床,其主轴热位移,无论在垂直方向还是在水平方向,都会造成较大的加工误差。,因此在分析机床热变形对加工精度影响时,还应注意分析热位移方向与误差敏感方向的相对角位置关系。对于处在误差敏感方向的热变形,需要特别注意控制。,四、机床热变形对加工精度的影响,龙门刨床、导轨磨床等大型机床,它们的床身较长,如导轨面与间稍有温差,就会产生较大的弯曲变形,故床身热变形是影响加工精度的主要因素。,四、机床热变形对加工精度的影响,大型机床热变形,五、减少工艺系统热变形对加工精度影响的措施,(一)减少热源的发热和隔离热源,(二)均衡温度场,(三)采用合理的机床部件结构及装配基准,(四)加速达到热平衡状态,(五)控制环境温度,立式平面磨床采用热空气加热温升较低的立柱后壁,以均衡立柱前后壁的温升,减小立柱的向后倾斜。图中热空气从电动机风扇排出,通过特设的软管引向立柱的后壁空间。采取这种措施后,磨削平面的平面度误差可降到未采取措施前的,1,3,1,4,。,五、减少工艺系统热变形对加工精度影响的措施,(三)采用合理的机床部件结构及装配基准,l,采用热对称结构,在变速箱中,将轴、轴承、传动齿轮等对称布置,可使箱壁温升均匀,箱体变形减小。,五、减少工艺系统热变形对加工精度影响的措施,加工中心机床,在热源影响下,单立柱结构会产生相当大的扭曲变形,而双立柱结构由于左右对称,仅产生垂直方向的热位移,很容易通过调整的方法予以补偿。,五、减少工艺系统热变形对加工精度影响的措施,2,合理选择机床零部件的装配基准,图,3-52,表示了车床主轴箱在床身上的两种不同定位方式。由于主轴部件是车床主轴箱的主要热源,故在图,a,中,主轴轴心线相对于装配基准,H,而言,主要在,Z,方向产生热位移,对加工精度影响较小。而在图,b,中,方向,Y,的受热变形直接影响刀具与工件的法向相对位置,故造成的加工误差较大。,五、减少工艺系统热变形对加工精度影响的措施,(四)加速达到热平衡状态,对于精密机床特别是大型机床,达到热平衡的时间较长。为了缩短这个时间,可以在加工前,使机床作高速空运转,或在机床的适当部位设置控制热源,人为地给机床加热,使机床较快地达到热平衡状态,然后进行加工。,五、减少工艺系统热变形对加工精度影响的措施,(五)控制环境温度,精密机床应安装在恒温车间,其恒温精度一般控制在士,l,以内,精密级为土,0.5,、恒温室平均温度一般为,20,,冬季可取,17,,夏季取,23,。,五、减少工艺系统热变形对加工精度影响的措施,实例分析,实例分析,六、加工误差的综合分析,一、加工误差的性质及分类,加工误差,随机误差,系统误差,常值误差,变值误差,连续加工一批工件,误差大小和方向保持不变。,连续加工一批工件,误差大小和方向有规律变化。,连续加工一批工件,误差大小和方向无规律变化,但,具有一定的统计规律。,误差种类,常值系统误差,变值系统误差,随机误差,特点,与加工顺序无关;,预先可以估计;,较易完全消除;,不会引起工件尺寸波动;,不会影响尺寸分布曲线的形状,与加工顺序有关;,预先可以估计;,较难完全消除;,会造成工件尺寸的增大或缩小;,影响尺寸分布曲线的形状,预先不能估计;,不能完全消除,只能减小;,工件尺寸忽大忽小,造成一批工件的尺寸分散,六、加工误差的综合分析,随机性误差:,可采用统计分析法,缩小它们的变动范围。,变值系统性误差:,查明其大小和方向随时间变化的规律后,采用,自动连续补偿或自动周期补偿,的方法消除。,常值系统性误差:,查明其大小和方向后,通过,调整,消除。,不同性质误差的解决途径,六、加工误差的综合分析,二、加工误差的统计分析法,加工误差的统计分析法,指以生产现场观察、检测所得的结果为基础,运用数理统计的方法进行归纳、分析和判断,找出产生误差的原因,从而采取相应的措施。,误差的统计分析常用,2,类方法:,实际分布图(直方图),理论分布曲线(正态分布曲线),分布图分析法,点图分析法,六、加工误差的综合分析,分布图分析法,1,、实际分布图(直方图),加工一批工件,由于各种误差的存在,加工尺寸的实际数值是各不相同的,这种现象称为,尺寸分散,。,以工件尺寸为横坐标,以频数或频率为纵坐标,,,即可作出该工序工件加工尺寸的实际分布图,直方图,。,直方图:,可以判断生产过程是否稳定,估计生产过程的加工质量及产生废品的可能性。,六、加工误差的综合分析,图例,零件尺寸直方图,六、加工误差的综合分析,直方图的作法与步骤,1)收集数据,抽取一个样本,样本容量一般取100件左右,测量各零件的尺寸,并找出,x,max,和,x,min,。,2),分组,组数过多,分布图会被频数的随即波动所歪曲;组数太少,分布特征将被掩盖。,3)确定组距及分组组界,组距:,h=(x,max,-x,min,)/(k,-,1,),第一组上界值:,s,1,=,x,min,+,h/2,第一组下界值:,x,1,=,x,min,-,h/2,4)统计频数分布,将各组的尺寸频数、频率和频率密度填入表中。,5)绘制直方图,按表列数据以频率为纵坐标,组距为横坐标画出直方图。,六、加工误差的综合分析,取在一次调整下加工出来的轴件,200,个,经测量,得到最大轴颈为,15.145mm,,最小轴颈为,15.015mm,,统计每组的工件数,结果如下表,画出实际分布曲线。,组号,尺寸间隔,频数,频率,频率密度,组号,尺寸间隔,频数,频率,频率密度,1,15.01-15.02,2,0.010,1.0,8,15.08-15.09,58,0.290,29.0,2,15.02-15.03,4,0.020,2.0,9,15.09-15.10,26,0.130,13.0,3,15.03-15.04,5,0.025,2.5,10,15.10-15.11,18,0.090,9.0,4,15.04-15.05,7,0.035,3.5,11,15.11-15.12,8,0.040,4.0,5,15.05-15.06,10,0.050,5.0,12,15.12-15.13,6,0.030,3.0,6,15.06-15.07,20,0.100,10.0,13,15.13-15.14,5,0.025,2.5,7,15.07-15.08,28,0.140,14.0,14,15.14-15.15,3,0.015,1.5,工件频数分布表,六、加工误差的综合分析,六、加工误差的综合分析,2、理论分布图,正态分布曲线,在分析工件的加工误差时,通常用,正态分布曲线,代替实际分布曲线,可使问题的研究大大简化。,大量实践经验表明,在用调整法加工时,当所取工件数量足够多,且无任何优势误差因素的影响,则所得一批工件尺寸的实际分布曲线便非常接近正态分布曲线。,六、加工误差的综合分析,正态分布曲线方程式为:,尺寸分布概率密度,工序的标准偏差,工件平均尺寸(分散中心),六、加工误差的综合分析,1,)直线以,x=,左右对称;,2,)对,正负偏差相等;,3,)分布曲线与横坐标所围成的面积包括了全部零件数,(l00,),,故其面积等于,1,;取正态分布曲线的分布范围为土,3,(占,99.73,)。,-,+,六、加工误差的综合分析,六、加工误差的综合分析,正态分布曲线的特征参数 和,算术平均值,是确定曲线位置的参数。,决定一批工件尺寸,分散中心的坐标位置。改变时,整个曲线沿,轴平移,,但曲线形状不变。,值主要由常值系统误差确定。,(2),工序标准偏差,决定了分布曲线的形状和分散范围。,值减小时曲线形状陡峭,尺寸分散的范围小;,值增大时曲线形状平坦,尺寸分散的范围大。,值主要由随机误差和变值系统误差决定。,六、加工误差的综合分析,1),判断加工方法是否合适,工序能力系数:,判断工艺能力能否满足加工精度要求。,3,、分布曲线法的应用,六、加工误差的综合分析,工序能力系数,工序等级,说明,Cp1.67,特级,工艺能力高,可以允许有异常波动,不一定经济,1.67Cp1.33,一级,工艺能力足够,可以允许有一定的异常波动,1.33Cp1.00,二级,工艺能力勉强,必须密切注意,1.00Cp0.67,三级,工艺能力不足,可能出现少量不合格品,0.67Cp,四级,工艺能力很差,必须加以改进,六、加工误差的综合分析,2),估算工序加工的合格率及废品率,令,则有:,(z),为图中阴影线部分的面积。,对于不同,z,值的,(z),,可由表查出。,(z),标准正态分布,Q,废品率,0.5-,(Z),六、加工误差的综合分析,六、加工误差的综合分析,六、加工误差的综合分析,计算简单,能及时提供主动控制信息,可用于稳定过程、也可用于不稳定过程。,分布图分析法缺点,1)不能反应误差的变化趋势,很难把随机性误差与变值系统误差区分开来;,2,)不能及时提供工艺过程精度的信息,,事后分析,。,点图分析法,六、加工误差的综合分析,点图分析法,六、加工误差的综合分析,1,、个值点图,依次测量每件尺寸,记入以零件号为横坐标,以尺寸为纵坐标的图表中,,能较清楚地揭示出加工过程中误差的性质及其变化趋势。,瞬时尺寸分散中心,变值系统误差。,六、加工误差的综合分析,2,、均值-极差点图,采用顺序小样本(46),由小样本均值点图和极差点图组成,横坐标为小样本组序号。小样本组一般为2030。,反映了系统性误差、随机误差及其变化趋势。,定期测小样本尺寸;,计算均值和极差,R,;,确定中心线和上下控制线,定期描点。,六、加工误差的综合分析,图例,六、加工误差的综合分析,均值点图反映了质量指标分布中心(系统误差)的变化,极差点图反映了质量指标分布范围(随机误差)的变化,六、加工误差的综合分析,点子无明显规律性。,生产过程稳定的标志:,没有点子超出控制线;,大部分点子在中线附近波动,小部分点子在控制线附近,;,六、加工误差的综合分析,44,20,46,32,20,40,52,33,40,25,43,38,40,41,30,36,49,51,38,34,22,46,38,30,42,38,27,49,45,45,38,32,45,48,28,36,52,32,42,38,40,42,38,52,38,36,37,43,28,45,36,50,46,33,30,40,44,34,42,47,22,28,34,30,36,32,35,22,40,35,36,42,46,42,50,40,36,20,16,53,32,46,20,28,46,28,54,18,32,33,26,45,47,36,38,30,49,18,38,38,表中数据为实测尺寸与基本尺寸之差。单位,um,实例分析,1,79,参考答案,80,组号,组界,um,中心值,x1,频数,频率(),频率密度,um(%),1,13.5-18.5,16,3,3,0.6,2,18.5-23.5,21,7,7,1.4,3,23.5-28.5,26,8,8,1.6,4,28.5-33.5,31,13,13,2.6,5,33.5-38.5,36,26,26,5.2,6,38.5-43.5,41,16,16,3.2,7,43.5-48.5,46,16,16,3.2,8,48.5-53.5,51,10,10,2,9,53.5-58.5,56,1,1,0.2,3),整理频数分布表,参考答案,81,4,)根据数据画直方图,参考答案,82,参考答案,83,参考答案,84,实例分析,2,参考答案,86,87,参考答案,88,一、加工误差的性质及分类,二、加工误差的统计分析法,1,、分布图分析法,2,、点图分析法,小结,六、加工误差的综合分析,随机变量,名称,MATLAB,密度函数,随机变量,名称,MATLAB,密度函数,Beta,分布,betapdf,标准正态分布,normpdf,二项分布,binopdf,泊松分布,poisspdf,卡方分布,chi2pdf,瑞利分布,raylpdf,指数分布,exppdf,T,分布,tpdf,F,分布,fpdf,均匀分布,unifpdf,伽马分布,gampdf,Weibull,分布,weibpdf,几何分布,geopdf,非中心,F,分布,ncfpdf,超几何分布,hygepdf,非中心,T,分布,nctpdf,对数正态分布,lognpdf,非中心卡方布,ncx2pdf,如果将上述命令中的后缀,pdf,分别改为,cdf,,,inv,,,rnd,,,stat,就得到相应的随机变量的分布函数、分位数、随机数的生成以及均值与方差,.,一,.,随机变量与分布,附录:加工误差分析的,MATLAB,实现,我们利用,MATLAB,中的函数来计算正态分布的分布函数、概率密度函数值、做出密度函数曲线、分位数,.,在,MATLAB,中计算上述分布函数的命令为:,P=NORMCDF(x,mu,sigma),Default values for MU and SIGMA are 0 and 1 respectively.,例,1,已知 试求:,解:,normcdf(3,2,0.5)=0.9772;,normcdf(2,2,0.5)-normcdf(1,2,0.5)=0.4772,1.,计算分布函数与概率密度函数值,:,2.,做出密度函数曲线、求分位数,已知,X,的均值和标准差及概率,p=PXx,求,x,的命令为,:,X=NORMINV(P,MU,SIGMA),例,2.XN(1,0.04),pXx=0.6827,求,x,解:,x=norminv(0.6827,1,0.2)=,1.0951,NORMSPEC(a,b,MU,SIGMA),用于做出随机变量在区间,a,b,上的正态密度曲线,例,3.,若,XN(2,4),,作出,X,在,-1,,,3,上的曲线,解:,normspec(-1,3,2,2),图形如右图所示,(,图,7.1),二,.,数据特征,设 是取自总体,X,的一个简单随机样本,在,n,次抽样以后得到样本的一组观测值 我们通过对数据的分析研究可以得到总体,X,的有关信息,在,MATLAB,中有专门的函数分析数据特征,如下表所示,.,位置特征,MATLAB,函数,变异特征,MATLAB,函数,算术平均,mean,极差,range,中位数,median,方差,var,切尾平均,trimmean,标准差,std,几何平均,geomean,四分位极差,iqr,调和平均,harmmean,平均绝对偏差,mad,例,4.,已知数据:,459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680 926 653 164 487 734 608 428 1153 593 844 527 552 513 781 474 388 824 538 862 659 775 859 755 649 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638 699 634 555 570 84 416 606 1062 484 120 447 654 564 339 280 246 687 539 790 581 621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851,计算该数据特征,.,解:,a=,459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680,926 653 164 487 734 608 428 1153 593 844 527 552 513 781 474 388 824 538 862 659,775 859 755 649 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638,699 634 555 570 84 416 606 1062 484 120 447 654 564 339 280 246 687 539 790 581,621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851;,b=a(:);%,将矩阵变成数列,T=mean(b),median(b),trimmean(b,10),geomean(b),harmmean(b),range(b),var(b),std(b),iqr(b),mad(b),位置特征,计算结果,变异特征,计算结果,算术平均,600,极差,1069,中位数,599.5,方差,38663.03,切尾平均,600.64,标准差,196.629,几何平均,559.68,四分位极差,243.5,调和平均,499.06,平均绝对偏差,150.86,计算结果如下:,例,5.,已知数据:,1,,,1,,,1,,,1,,,1,,,1,,,100,;计算其数据特征,由此你有何发现?,解:,x=1,1,1,1,1,1,100,;,y=,mean(x),median(x),geomean(x),harmmean(x),trimmean(x,25);,range(x),var(x),std(x),iqr(x),mad(x),五,.,直方图与概率纸检验函数,1.,直方图:,为了直观地了解随机变量的分布特征,如对称性,峰值等,直方图是广泛使用的方法,.,格式:,hist(,data,k,).,说明:,data,是原始数据,该命令将区间,(min(,data,),,,max(,data,),分成,k,等份,并描绘出频数直方图,,k,的缺省值为,10.,如果需要事先给出小区间,则将区间的中点存放在向量,nb,然后用命令:,n,x,=hist(,data,),,,或,n,x,=hist(,data,k,),,,或,n,x,=hist(,data,nb,),,,其中,,n,返回,k,个小区间的频数,,x,返回小区间的中点,.,2.MATLAB,中的概率纸检验函数,.,格式:,normplot(,data,)%,如果数据,data,服从正态分布,则做出的图形基本上都位于一条直线上,.,weibplot(,data,)%,如果数据,data,服从威布尔分布,则做出的图形基本上都位于一条直线上,,例,8.,作出例,4,数据的直方图,该数据服从正态分布还是威布尔分布?,解:输入命令,hist(,a,(:),30),得到图,7.2,从图,7.2,所示的直方图发现数据比较接近于正态分布,我们用命令,normplot(a(:),进行检验,.,从图,7.3,中可见数据点基本上都位于直线上,故可认为该数据服从正态分布,由于已经计算出该数据的均值为,600,,标准差为,196.629,,所以数据服从,正态分布的检验:,1.,大样本,h=jbtest(x),h=0,接受正态分布,,h=1,拒绝正态分布,2.,小样本,h=lillietest(x),h=0,接受正态分布,,h=1,拒绝正态分布,The End,此课件下载可自行编辑修改,供参考!,感谢您的支持,我们努力做得更好!,
展开阅读全文