最小二乘法的线性拟合.ppt

第4章 工程数据处理及数据库技术 4.2 线图的程序化处理 我们知道用线图来表示函数关系是一种很常用的方法，这种方法比较直观，并且能表现出函数的变化趋势。在机械设计资料中，有很多参数间的函数关系是通过线图来表示的。这些线图在对数坐标中，一般是直线或折线；在普通的直角坐标中，大多是曲线图。传统设计时，根据线图查得的参数直接代入公式，而在CAD作业中，目前还不能对线图直接编程，因此需对线图进行相应的处理，以实现对参数图的存贮和自动检索。 常用的处理方法有三种： （1）线图所表示的各参数之间本来就有计算公式，只是由于计算公式复杂．为了便于手工计算将公式绘成线图，以供设计时查用。对于这类线图处理的方法为：找到线图原有公式，将公式编写成程序。这是最精确的程序化处理方法，但难以找到。 （2） 线图所表示的各参数之间没有或找不到计算公式。将线图离散化为数表，然后按数表的处理方法进行处理。 （3）用曲线拟合的方法求出线图的拟合公式 （又称经验公式），再将公式编写成程序。 4.2.1线图的数表化处理 所谓线图数表化处理是将线图离散转化为一张数表，然后按数表的处理方法进行处理。 右图为渐开线齿轮的 一种齿形系数曲线图 在CAD作业中，若能找到所给参数线图曲线的计算公式，则应选择该公式直接计算出所需的参数值，以确保所获得参数的精度。对于直线线图(直角坐标或对数坐标)或是以折线构成的折线图，则应把它们先转化为代数表达式，然后通过计算来获得较准确的参数值，这就是下面介绍的线图公式化处理。 4.2.2 直线图的公式化处理 1、直角坐标直线图的公式化处理 （a）直齿轮 （b）斜齿轮 2、对数坐标直线图的公式化处理 对数坐标中的直线方程可写为： 注意：一般程序语言中，只有lnx （自然对数）无十进制对数lgx ，所以编程时，要进行换底运算。 3、区域图的公式化处理 （1）按区域图的中线取值 找出区域中线的表达式， 为此设齿面最小硬度为HB0 。 最高硬度为HBl ，SH0表示 最小硬度对应的极限应力中 值，SHl表示最大硬度对应 的极限应力中值，由此根据 直线的两点式方程可以写出 极限应力中线的表达式为： （2）按区域图的位置取值 在确定材料极限应力时只能取中值，不尽合理。为了使设计者能根据所用材料的不同性能，按实际情况在区域图内取不同的值，为此，增加两个参数，一个是极限应力的幅度值SF，另一个是表示极限应力在区域图中的位置参量ST。ST＝1时表示取上限值，ST＝0时表示取中值，ST＝-1时表示取下限值。此时极限应力的计算式变为： ST在+1~-1之间取不同的值时，就可以获得区域图 中任意位置上的极限应力值。 4.3 数表与线图的公式化处理 前面介绍的数表与线图的程序化处理方法，这种方法虽然解决了数表和线图在CAD作业中的存储和检索问题，但还存在下述一些缺点： 1)占用大量计算机内存。数表和线图的程序化处理，要将数表中的全部数据编进计算程序中，实现数据的自动检索。当数表很庞大时，所占内存很大。一般情况下，一个设计计算程序常常需要使用多个数表，则所占内存更加庞大，严重时甚至会影响程序的正常运行。 2)效率低，占机时间长。通常设计所使用到的仅是数表中的一小部分数据，有时甚至只是其中的一、二个。但数表程序化处理对数表中的每个数据，无论在当时的计算程序中 是否被用到，都必须顺序地将全部数据读入内存。 检索时，一般又得顺序地从头检索至所需的那个 数据为止。 4.3.1 曲线拟合 数表程序化处理一般只适用于数表较小(数据量较小)、计算程序使用数表个数不多的情况。对于比较大型的计算程序，常常需使用很多的数表，数据量很大，在这种情况下数表的处理就要采用其它的方法。其中一种方法就是本节所要介绍的曲线拟合。 此方法的基本原则是构造一曲线函数y＝P(x)作为函数f（x）的近似表达式，所构造曲线函数并不严格通过所有节点，而是尽可能反映函数的数据变化趋势，比较符合实际规律。曲线拟合的方法有很多，这里我们只介绍最小二乘法（线性、多项式、 指数曲线）的数据拟合。 线性拟合 在直角坐标系中点大致呈线性分布，