(计算机在材料科学和 与工程中的应用)第二章材料科学和 与工过程中的数据处理.ppt

假设随机变量y与p个自变量 之间存在着线性相关关系，实际样本量为n,其第i次观测值为 则其n次观测值可写为如下形式： 其中 是未知参数， 是p个可以精确测量并可控制的一般变量, 是随机误差, (2) (3) 假定 是相互独立且服从同一正态分布N(0, ) 的随机变量 。 2.1 数据处理的基本理论—2.1.3、多元线性回归 若将方程组(3)用矩阵表示，则有 (4) 式中: 建立多元线性回归方程 (5) 2.1 数据处理的基本理论—2.1.3、多元线性回归 来描述多元线性模型 (6) (7) (8) 2.1 数据处理的基本理论—2.1.3、多元线性回归 与一元线性回归分析相同，其基本思想是根据最小二乘原理，求解 使全部观测值 与回归值 由于残差平方和 的残差平方和达到最小值。 (9) 是 的非负二次式，所以它的最小值一定存在。 根据极值原理，当Q取得极值时， 应满足 (10) 2.1 数据处理的基本理论—2.1.3、多元线性回归 由（9）式，即满足 (11) 2.1 数据处理的基本理论—2.1.3、多元线性回归 上式(11)称为正规方程组。它可以化为以下形式 如果用A表示上述方程组的系数矩阵可以看出A是对称矩阵。则有 (12) (13) 2.1 数据处理的基本理论—2.1.3、多元线性回归 ?式中X是多元线性回归模型中数据的结构矩阵， (13)式右端常数项也可用矩阵D来表示? 是结构矩阵X的转置矩阵。 即 Ab=D 或 (14) (15) (17) 如果A满秩（即A的行列式 则由(13)式和(14)式得 的最小二乘估计为 ）那么A的逆矩阵A-1存在， (16) ?b就是多元线性回归方程的回归系数。 为了计算方便往往并不先求 而是通过解线性方程组来求b ，再求b， (14)是一个有p+1个未知量的线性方程组，它的第一个方程可化为 2.1 数据处理的基本理论—2.1.3、多元线性回归 式中 (18) 将（17）式代入（12）式中的其余各方程，得 (19) 2.1 数据处理的基本理论—2.1.3、多元线性回归 其中 (20) 将方程组（19）式用矩阵表示，则有 Lb=F 其中 ? 于是 (21) (22) (23) 2.1 数据处理的基本理论—2.1.3、多元线性回归 2.1 数据处理的基本理论—2.1.3、多元线性回归 例：表观变形激活能计算 高温塑性变形最显著的特点之一便是变形速度受热激活过程控制。表观变形激活能表示原子跃迁所需克服的能垒大小，是反映塑性变形难易程度的重要物理参量。变形温度、应变速率对流动应力的影响可用Arrhenius方程表示： 为应变速率(s-1)； 为表观变形激活能(kJ·mol-1)； 为流动应力(MPa)； 为绝对变形温度(K)； 为普适气体常数(8.3145 J·mol-1·K-1)； 为应变速率敏感性指数。 2.1 数据处理的基本理论—2.1.3、多元线性回归 Ti-6Al-4V合金在不同温度下压缩变形时的应力-应变曲线 2.1 数据处理的基本理论—2.1.3、多元线性回归 2.1 数据处理的基本理论—2.1.3、多元线性回归 用多元线性回归求解 ? ? ? ? ? ? ? 利用多元线性回归一次性求解系数a、b1、b2即可得Q值！ 2.2 Excel和Origin软件在材料科学与工程中的数据处理应用 Excel应用举例 研究某种材料在某种腐蚀液中腐蚀时间与腐蚀量的关系，实验数据如下表所示： 例2.1、一元(直线)回归分析应用实例 方法一： 1、在A,B列输入x，y值； 2、用鼠标选中x，y数据区； 3、用鼠标指到“图表向导”工具按钮,选择x，y散点图，单击下一步； 4、在弹出的对话框中进行相应选择； 5、在产生的散点图上单击鼠标右键,在弹出的菜单上选择“添加趋势线”，进行相应选择即可得到所需结果。 2.2 Excel和Origin软件在材料科学与工程中的数据处理应用 2.2 Excel和Origin软件在材料科学与工程中的数据处理应用 2.2 Excel和Origin软件在材料科学与工程中的数据处理应用 2.2 Excel和Origin软件在材料科学与工程中的数据处理应用 方法二（以2013版Office-Excel为例） 1、在A,B列输入x,y值 2、用鼠标点击文件-选项--加载项- 3.再次用鼠标 3、点击数据-数据分析-回归 4.在“回归”分析对话框中进行相关选择后，即可得到所要结果。相关选择： 输入区域选择 y值区域,x值区域 输出区域选择 置信度选择 95% 2.2 Excel和Origin软件在材料科学与工程中的数据处理应用 2.2 Excel和Origin