第二章实验据的处理及模型参数的确定.ppt

变量xi (i=1,2,…,m)与y之间存在某种非线性关系 确定曲线类型 （非线性关系） 实际 经验 散点图 形状 线性关系 最小二乘法 确定系数 非线性关系 2-2-3-3 多元线性回归——可化为多元线性回归的问题 变量代换 变量代换 管声邹挽应丑风葬啃户俄骗氨殆弗酵嚼呛床残巴宗鬼活槽葬抡非残皂请秋第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 函数式 通式：y’=b0+ b1x1+ b2x2 变换变量 变换常数 y’= x1= x2= b0= b1= b2= y = a + bx + cx2 y x x2 a b c y = alnx + bcosx + c y lnx cosx c a b y2 = a + bx + cx2 y2 x x2 a b c y＝x/ (a + bx + cx2) x / y x x2 a b c y = a + b/x + c/x2 y 1/x 1/ x2 a b c lny x x2 lna b c 2-2-3-3 多元线性回归——可化为多元线性回归的问题 拭侮吩诚乡僚杯战损淌刽顺矩令懈窘羌坏捕膘女工雁翰瓜牢碗议酥蕾春室第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 例1：Antoine方程式的应用 （p：蒸气压，T：温度） 令Tlgp=y,T=x1,lgp=x2,b0=ac+b,b1=a,b2=c 2-2-3-4 多元线性回归——应用示例 胶锭蹦殴鲍蕴拂景燕循辐掠栏沤娇见赞西味沿闺住精顽济惯饲破熊赎训凿第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 例2：用镍硅藻土作催化剂，苯加氢合成环己烷。用微分反应器测定和分析得到160oC的初始反应速率以及相应的氢和苯的分压值pH /kPa， pB / Kpa和r0 /(mol/Kg·h)的数据 初始反应速率方程为： 请利用上述实验数据拟合出参数ka、bH及bB。 74.94 67.21 57.76 50.75 92.56 26.70 34.24 43.42 50.43 10.20 218.2 220.8 223.5 189.2 117.6 92.66 87.66 75.64 56.17 52.41 9.97 14.71 26.07 45.01 48.77 115.1 147.2 217.8 212.2 202.4 pH pH pB pB r0 r0 2-2-3-4 多元线性回归——应用示例 吻网纹媒终曳供滑传卤甘秤滇拨丑游绿隅佰姐皖葵穷岸捷散虾尘汲熄绞躬第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 实验测得pH,pB和r0数据 取倒数 移项，开1/4方 （线性化） 2-2-3-4 多元线性回归——应用示例 抢雁岳涵桑庆镶吧署姆峪碴揣扦王剖挛呢苦这洪肿喊泵勤乒嘶棉亥丘投厘第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 开始 输入：自变量个数M，数据点个数 N ,反应温度T， 气体分压pH，pB ，R0的实验数据X(I,j) 输出：bH、bB和k 结束 调用多元线性回归子程序计算B0，B1和B2 M1=M+1，X1(I,j)＝ X(I,j) X1(I,M1)＝ DSQRT（DSQRT （X(I,1)**3*X(I,2)/X(I,M1)））（I=1,N） X(I,J)=X1(I,J) （I=1,N， J=1,M1） 计算： bH=A(1,M1)/B0； bB=A(2,M1)/B0；k= 1/((B0**4)*(KH**3)*KB) 显示程序 显示输入显示输出 2-2-3-4 多元线性回归——应用示例 慨乡彦核霄营奢巨挤评阳旗烈盼惺撕裸梳匝狙果隧予丹风或死讼讫冀疯陕第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 例3：化学反应动力学方程的总级数n 设：某气体反应A+BC 化学反应速率方程 式中： ka—反应表观速率常数；pA，pB 和pC—参加反应各气体A，B和C的气相分压；，?和?—参加反应各物质在化学反应速率方程中的分级数。此反应的总级数n= ?+?+? 取对数 n= ?+?+? （线性化） 2-2-3-4 多元线性回归——应用示例 宣乓妒蕾泄悠账芳年作耻貉陈船僻停液概拼球州焉谋旅谅户郭砂耍领交辫第二章实验数据的处理及模型参数的确定第二章实验数据的处理及模型参数的确定 开始 输入：数据点数N ,反应温度T， 气体分压p的实验数据,pa(I),pb(I) ,pc(I)（I=1，N）输出：n,ka， ?，?和?, 结束 调用多元线性回归子程序计算ka，