上海交通大学计算方法课件(宋宝瑞)CH.5

1第六章函数的最优逼近与拟合1线性赋范空间中的逼近问题11函数逼近与函数空间逼近的思想和方法渗透于几乎所有的学科，其中包括自然科学和人文科学中的学科。逼近论既是一门研究函数的各类逼近性质的学科，属于函数论的范畴，同时又是计算数学和科学工程计算诸多数值方法的理论基础和方法的依据。本章讨论科学计算中基于逼近论的一些函数逼近方法。函数逼近方法与函数插值方法相类似，它也是在某一函数类中求函数，使它与被逼近函数之间满足一定的近似条件。在插值方法中，这个近似条件是在插值结点上使插值函数与被插值函数的函数值对应相等（包括各阶导数值对应相等）；在逼近方法中，这个近似条件用逼近函数与被逼近函数之间的某种距离来表达。数学上常在各种集合中引入某些确定关系，称为赋予集合以某种空间结构，并将这样的集合称为空间。例如，在线性代数中将所有实N维向量组成的集合，按向量加法及向量与数的乘法构成实数域上的线性空间，记作，称为N维向量空间。类似地，对次数不超过N（N为正整数）的实R系数多项式全体，按通常多项式加法及数与函数乘法构成数域R上的线性空间，用表示，称为多项式空间，又如所有定义在区间上的连续函数集NP合，按函数加法和数与函数乘法构成R上的线性空间，记作，称为,BAC连续函数空间。定义11设集合S是数域P上的线性空间，如果存在SXN,12不全为零的数使得PAN,1（11）02NXAX则称是线性相关的；若（11）式只对成立，NX,1021NA则称是线性无关的。如果S包含一个由N个元素组成的极大线性无关组，则称S是N维的，如果对任意自然数N，S中存在N个线性无关的元素，则称S是无穷维的。显然是无穷维的，但对于和0，可以,BAC,FXCABXF用有限维空间中的元素逼近，使误差，这就NPXPMPF是著名的WEIERSTRASS逼近定理。定理11设，则0，使得,BAFXNPXPFN在上一致成立。,BA此定理可在数学分析书中找到证明。1912年BERNSTEIN构造了一个多项式（12）KNKNKNXFXFB10并证明在上一致成立，若在上MLIMN,XF1,0阶可导，则还有。这也就从理论上给出定理11XFMN的构造性证明。称为F的N次BERNSTEIN多项式。但由于XFB收敛于很慢，因此实际计算的近似值时，很少用这种方法。XFN3连续函数还可用其他函数集合逼近。一般地，可用一组在XF上线性无关的函数集合的线性组合来逼近，,BACNIIX0,BACXF即用01SPAN,CAB（13）XN逼近。于是函数逼近问题就是对，在线性子空间F,FX中找某一元素,使在某种意义下最小。换句话说，F也就是找一组系数，使（13）式中的成为在10NAXF中的最佳逼近元。12赋范线性空间中的最佳逼近既然是在线性子空间中寻找某一函数的逼近，必须引入一个度量的概念来衡量逼近的好坏，即衡量误差的大小。这对于不仅要|FXP计算函数在个别点上的值，而且要考虑函数在区间上的整体性质时，特别有意义。定义12设为线性空间，如果对每一向量，都有一实数与X之对应，把这一实数记为，并且这一对应具有下列性质|X1），等号成立当且仅当，（的零向量）正性|0X02），正齐性|3）三角不等式|YY则称为X的范数，一个线性空间，如果其中定义了满足上述三条公理的范数，我们称之为赋范线性空间，记为X；注范数的概念是很广泛的，只要满足上述3条，都可作为范数。以4前我们引入的向量和矩阵的范数都是相应空间的范数。同一线性空间可以赋予不同的范数。定义13设为赋范线性空间，其范数为，若序列，K0NK，使F0LIMFN则称序列依范数收敛于F，记作0N|LINFA现在运用赋范线性空间的概念来讲解函数的逼近问题。逼近论的两个基本问题1给定赋范线性空间，以及的一个真子空间。对于如KFK果在中存在这样的元素使得对所有的有SSFF成立，则称为F在范数意义下在中的最佳逼近元。立即就会出现S|这样的问题（1）最佳逼近元是否存在是否唯一（2）当最佳逼近元唯一存在时，如何构造最佳逼近元。2设有的一系列子空间，若对每一K12NK个子空间，问题1中的最佳逼近元存在，是否有，以及JJSLIMNSF收敛的速度。我将对于不同的范数定义，逐一研究这些问题。2最佳一致逼近给出线性空间，取其范数和子空间为,BAC5，（21）MAXFFB1,NSPXNP根据定理11，问题2的回答是肯定的。考虑问题1对中任一元素F的最佳逼近问题称为最佳一致逼近（通常又称为,BACCHEBYSHEV意义下的逼近）。定理21如果F在上连续，则在集合中存在一个元素，是F,BA的最佳一致逼近元。证略。4最佳平方逼近大家已经学过三维欧氏空间，知道在欧氏空间里一个向量的长度，两个向量的夹角，向量到子空间的投影等是什么意思，在这一节里，我们要把函数逼近问题与欧氏空间联系起来，研究欧氏范数下的最佳逼近。41线性内积空间定义41设为实线性空间，如果对每一对向量，都有一YX,实数与之对应，把这一实数记为，并且这一对应具有下列性质,YX1），,XY2），3），,ZZ64），等号成立当且仅当，（的零向量）。,0X0X则称为X和Y的内积，一个线性空间，如果其中定义了满足上述四条公理的内积，我们称之为内积空间。例41设，对，定义内积,BAC,FXGCAB，易验证它满足公理14）。DXGFFBA,例42令（N维实向量空间），对其中向量R，TN,21TNYY,21定义内积NXXYX21,易知它满足公理14）。令，可以证明满足第二章中关于范数的公理，即,NN是一种范数，称为由内积导出的范数，通常记作，在不发生混淆X2X的情况下也可简记为。X定义42设为内积空间，若，则称X与Y正YX,0,Y交，记作42线性内积空间的最佳逼近设，是线性内积空间的N1个线性无关的元素，子01N集，在中寻求对的某一元素F的最佳逼近，,SPAS即对有SFSF7定理41是集合中对F的最佳逼近元素，其充要条NIICS0件是与所有正交。F,1J证明先证充分性，设S是中任一元素，考虑FSF由于与所有正交，从而也与正交，因此FSJ,2FSFFSF22SF从而就是最佳逼近元素，充分性得证。S再证必要性，假设是中元素，与，中某一元素1SFS10N不正交，记K，则KK且0,1KSF记KS12易知，现估计的范数2SF,1122KKSFFF,21SSF28从而推得2SF21F即不是最佳逼近元，必要性得证。1S推论最佳逼近元素如果存在，必定是唯一的。事实上，如果在子集中有两元素都是F的最佳逼近，则由定理41必有，J0，1，N,21JJSFSF于是和都与正交，于是有210212212121,SFSFSS这就表示。下面，我们证明最佳逼近元素是存在的并将其构造出来，由定理41最佳逼近元，必须满足JNJCS0，K0，1，N（41）0,NJKJF式（41）即，K0，1，N（42）,0KKJNJFC从而由下列方程组所决定9（43）,1001110100NNNNFCCF方程组（43）通常称为法方程组。现在我们研究方程组的系数行列式（44）0010110101,NNNG称为关于的GRAM行列式。,10NN定理42子集的元素，线性相关的充要条件是它们1F2NF的GRAM行列式等于零。,1NFG证明先证必要性。由于，线性相关，因此有一组不全为零1FNF的数，使得12N021KFF上式分别对，作内积有1FKF0,21221211KKKKKFFFFFF这是关于，的齐次线性方程组，由于，不全为112K零，故推出其系数行列式必为零，从而100,21KFFG再证充分性，假设,研究下列方程组0,21221211KKKKKFFFFFF容易明白它必有非零解，设为。令KKFFF21易知，另一方面，对上式分别用作内积可得F,1，0,1F,2F0FK从而有,21FFFK由于不全为零，故线性相关。定理的充分性得证。K,21,1现在我们回到方程组（43）的讨论，由于是的线性无N,10关元素，故,10NG从而（43）的解存在且唯一。也就是说，最佳逼近元素是存在的并可由（43）构造出来。下面再给出误差估计式，记为最佳逼近误差，SF则有11,22SFFSF,FF,S（45）,10FCFFCFN至此，关于内积空间的最佳逼近元素的特征，存在唯一性，构造及误差估计已全部论述清楚了。同一列在不同的K中，其收敛性可能不同NS例1,022,0,1NXNSXN220,103NLSD而,CX43函数的最佳平方逼近设A,B为有限或无限的区间，定义在它上面的函数，如果具有下X12列性质1）,0XBA2）0,DC,0CDC3）积分存在，N0，1，。XBAN则称其为A,B上的权函数。对于在A,B上给定的函数，引入内积DXGFXGFBA,22这里为权函数，总假设积分X（46）DXFBA2是存在的。我们研究满足（46）存在的函数全体组成的内积空间，并选子集。在子集上寻找一函数,10NSP为中某一函数FX的最佳逼近，是指对于0XCXSJNJ，都有（47）DSFBA2DXSFBA2（47）式的意义就是误差的平方在积分意义下达到极NJJXCXF0小，因此对于这种逼近就称为函数的最佳平方逼近，或称为最小二乘逼近，由于它是一个特殊的线性内积空间，因此最佳逼近的存在性，唯一性，解13的构造，误差估计等等已由42节全部回答了。我们在这里指出，（47）也可从另外的观点来得解，我们知道，上任一函数可写成，而积分NJJXC0DXCFNJJBA20是关于的二次多元函数，记N,10（48）DXXFCINJJBAN2010,在子集寻找对F的最佳平方逼近函数，就是寻找函数的极,10NCI小值。其必要条件是，K0，1，NKCI从而有0NBBKIIKAAIXXFDX写成内积符号就是，K0，1，N0,NKJKJCF上式恰巧就是（43）式。当然，从函数极小值的讨论去构造最佳逼近函数，其存在性，唯一性等都要建立一套理论来回答这些问题，但是由于我们已建立了一套最佳逼近理论，所以在解决具体问题时运用求极小值的办法，常常会带来演算的方便。,10NCI14例43在空间给定元素，子集由1，X的线性1,4CXF组合构成（线性多项式空间），试求在中的最佳平方逼近（取）。1X解由于，故0X1，143,DT3215,01041TD62,14DT另外又有1407,2FT143,80FTD设最佳逼近元为，据（43）列出法方程式01AXY80316423570A解得135,270即线性函数8XY15为在子集中的最佳平方逼近函数。XY从例43可以看出，平方逼近算法简单。例44取，1X，K0，1，N，在中求的N次KKX1,0CXFNPF最佳平方逼近多项式NXAAS10此时，J,K0，1，N10,KJJKXDJ则法方程组（43）的系数矩阵为（49）1213211,1NNNXGN这是一个HILBERT矩阵，前面我们已经知道当N较大时，系数矩阵（49）是高度病态的，求解时舍入误差很大，这是很不利于计算的，为避免这种情况，需要引入正交基的概念。44正交基如果是两两互相正交的（称为的一组正N,10N,10交基），法方程组的系数矩阵为对角阵，（43）成为16,11010NNFCFC从而直接可得，I0，1，N（410）,IIFC这样就不存在方程组病态的问题了。如果假设的范数均为1（称为的一组标N,10N,10准正交基），那么进一步有，I0，1，N（411）,ICF在这种情况下，1,NIISF而（45）式的误差估计式有（412）220|NIIFFC上式左端恒正，故又有BESSEL不等式（413）20|NIIF如果对所有的I，都存在，称为F的广义FOURIER展开，I0IIC称为广义FOURIER系数。IC17下面用正交化过程来证明正交基的存在性。定理43任何N维空间都存在正交基证明按照N维空间的定义，存在一组基，由这组基，通过NF,1正交化过程，可以构造出两两正交的基。NE,1令，在所在“平面”上找，也就是找具有下列形式的1FE1,F12FE选择数，使，即使0,120,12F由此得,/,12EF设两两正交而异于零的向量已经构造出来，向量要求1KEKE具有形式121EFEKKKK选择系数使与向量正交，即满足11,1,E01FKK,210,11KKEEF由于两两正交，故上述等式简化为1,KE，,11FKK，022EE18。0,11KKEEF由此得,/,11FKK222E,/,111KKF最后我们利用线性无关的条件来证明，注意到NF,20KE是和的线性组合，又可表示为和的线性KEF1,KE1KE1F12,组合，依此类推最后可将表示为K11KKFF的系数为1，而线性无关，所以，证毕KFKF,0E以上的正交化方法通常称为SCHMIDT正交化过程。5正交多项式若首项系数的N次多项式，满足0NAXGN（J,K0，1，）KJADXGXKKJBA0则称多项式序列在A,B上带权正交，并称是,10XXGNA,B上带权的N次正交多项式。X一般来说，当权及区间A,B给定后，从序列就可用,1244节所述的SCHMIDT正交化过程构造出正交多项式。用上述方法只能19一个接一个地构造出正交多项式，在使用上有所不便，利用多项式的某些性质，我们可以得到一些更直接，更方便的方法来构造正交多项式，较重要的有下列几类51勒让德（LEGENDRE）多项式当区间为1,1，权函数时，由正交化所得的多1X,2X项式就作为LEGENDRE多项式并用表示，这一类,0PN正交多项式有如下的简单表达式，N0，1，（51）21,20NNXDXP由于是2N次多项式，求N阶导数后得X12011AXANN于是得到首项的系数，显然最高系数为1的LEGENDRE多NX2AN项式为（52）122NNNXDXPLEGENDRE多项式有下述几个重要性质性质1正交性（53）NMDXPMN120120证明令，则（K0，1，N1）。NX120K设是在区间1,1上有N阶连续可微的函数，由分部积分知XQDXQDPNN211XN211DN1下面分两种情况讨论（1）若是次数小于N的多项式，则，故得XQ0XQN，当01DXPMM（2）若,22NNNXX，XQNN于是12212DXDXPNNN122NN由于，故10210234COSDXDXN1212NPN于是（53）得证。21性质2奇偶性（54）1XPXNN由于是偶次多项式，经过偶次求导仍为偶次多项式，2经过奇次求导则为奇次多项式，故N为偶数时为偶函数，N为奇数时为奇函数，于是（54）式成立。性质3递推关系考虑N1次多项式它可表示为110XPAXAPXN两边乘以，并从1到1积分，得PKDDKKN12当时，次数小于等于N1，上式左端积分为0，故得2NX，当KN时为奇函数，左端积分仍为0。故，于是0KA2PNNA11XPAXNN其中12422111NDXNAN11333NNP从而得到以下的递推公式（N1，2，）1N12XNNXXP（55）22由，利用（55）式就可推出10XPX，2/1322/353XXP，8/0544，765X，16/24P性质4在所有最高次项系数为1的N次多项式中，LEGENDRE多项式在1,1上与零的平方误差最小。XPN设是任意一个最高次项系数为1的N次多项式，它可表示为Q0XPAXKKNN于是DXQNN,21,0NKKNPAP当且仅当时等号才成立，即当时平110AXN方误差最小。性质5在区间（1，1）内有N个不同的实零点。XPN前几个LEGENDRE多项式的图形如下2352切比雪夫（CHEBYSHEV）多项式（31）ARCOSXNXTN如此定义的是X的多项式吗试看N0和N1的情形10XTCOSAR1注意到如下的三角恒等式COS21NN令，则有XARCOS，N1，2，（32）211XTTNN不难看出确实是N次多项式且与N的奇偶性相同，即随XTNN为奇数或偶数而成为上的奇函数或偶函数。,24显然的最高次项的系数为，即XTNNX12N（33）21的低次多项式称为N次CHEBYSHEV多项式。N前9个CHEBYSHEV多项式如下表表31132160251876483184214877243552330XXXTXXTXT显然（34）1MAX1N并且在点列COS0,KKN110X1X上，以正负交错的符号取到它的绝对值的最大值1。XTN25（36）KKNXT1根据定理22，首项系数为1的多项式N21XN为所有系数为1的N次多项式类中唯一的，在1,1上与零偏差最小的NX多项式。的N个零点全位于（1，1）内，且都是单重的，它们是。T，（37）2COSKX1,KN从几何上看，如果将以原点为圆心，以1为半径的上半圆周分成2N等分，再把圆周上所有奇分点位往X轴上投影，则恰好得到点列（37）。此外，实际计算中时常要求用，的线性组合表示，N0T1N其公式为2021NKKNXX这里规定，N18的结果见表32。20T表3226825631287461502643812164053742352043200TTXXTXTX它还是一类重要的正交多项式。CHEBYSHEV多项式有很多重要性质性质1正交性CHEBYSHEV多项式在区间1,1上带权正交，XTN21/XX且（57）0,2,12NMDXMN事实上，令，则，于是COSXSIN271200,COS02,NMNMTXDND性质2递推关系，N1，2，21XTXTNN，0这只要由三角恒等式COSCOSCOS1N令即得。X性质3只含X的偶次幂，只含X的奇次幂，这性质2TK12TK由递推关系直接得到。性质4在区间（1,1）上有N个零点，NNKK21COSK1，2，N。前几个TCHEBYSHEV多项式的图形如下2853无穷区间上的正交多项式正交多项式的正交区间A,B也可以是无界区域，当然此时权函数必须保证，N1，2，。XBNAXD1拉盖尔（LAGUERRE）多项式在区间上带权的正交多项,0XE式称为LAGUERRE多项式，其表达式为NXXNDLE它也具有正交性质200,XNMNEDX和递推关系，10L29，N1，2，。2112XLXNXLN2埃尔米特（HERMITE）多项式在区间上带权的正,2XE交多项式称为HERMITE多项式，其表达式为122XNXNEDH它满足正交关系NMXENNMX,2020并有递推关系，10HN1，2，。21XNXN54一般正交多项式的几个重要性质性质1在空间中，首一的N次正交多项式在所有2,LABNPX最高次项系数为1的N次多项式中与零的偏差最小。设是任意一个最高次项系数为1的N次多项式，它可表示为XQN0NNKPXAPX于是120,NNNKKXXP（勾股定理）30当且仅当时等号才成立，即当0110NAA时平方误差最小。XPQNN推论在所有最高次项系数为1的N次多项式中，LEGENDRE多项式在1,1上与零的平方误差最小。N性质2首一的CHEBYSHEV多项式在所有最高次项系数为1NTX的N次多项式中在1,1上与零的一致偏差最小。定理51是A,B上带权的N次正交多项式的充分必要条NPXX件是是N次多项式并且K0,1,N10BKNAD证明是容易的，从略。定理52设是A,B上带权的N次正交多项式，则NPXX的零点全位于A,B内，并且都是单重的。NPX证明设，如果，不妨设，令1NKKXX1XA，据定理51，2NKQXNP0BNAXPQDX但是上式中的被积函数为，积分21,AE,B应当大于0，矛盾。同理可证所有的零点不可能大于B，出现一对共轭复数或为二重的。316离散情况的最佳平方逼近在4论述了内积空间的最佳逼近理论。现在我们讨论离散情况的最佳平方逼近。对于N维欧氏空间中任二向量，NXX21NYY21其内积定义为，IITY1,相应的范数为2,XX给定欧氏空间的一个子集，其中为线性无关1LSPANIX组。子集对某一向量Y的最佳逼近，是指在中寻找一向量使它对的任意向量X都满足不等式YX如果，也即子集是一个N维欧氏空间，那么由线性代数的知识，NL恒有唯一的一组常数使C,21,1NIIXY此时，因此就是Y的最佳平方逼近，下面我们讨论0XYLN的情况。32由4定理41，与所有正交，即XYLX,21，K1，2，L0,K记，ILIXC112IINIX那么有,1KKILIYC从而由下列方程组所决定IC,21（61）,212221211LLLLLLXYCCXCCC其中NKJIJTIJIX1,1,NJIIKYXY（61）可写成33，,2121YXCXTJLLTJTJJJ1，2，L容易看出，上式又可写成,2122121YXCXXTLLLTL引入阶矩阵LN（62）12121212,LLNNLXXAX于是上式可写成（63）YACTT其中为L阶方阵，分别为L，N维向量，AT,，LC21NYY21容易看出，（63）式的系数矩阵是对称的。由于为线ATLX,21性无关组，我们可进一步推出矩阵是正定的。事实上，对于L维欧氏34空间的任一向量G，考虑的二次型有AT,0,G其中等号仅当时才成立即有，或即0T0,1221TILITLLXGG由是线性无关组推得，即。从而，当时LX,210I0G有0，,GAT也就是说矩阵是正定矩阵。由于矩阵是正定的，所以可建立求解ATT（63）的特殊方法，例如平方根法，乔列斯基（CHOLESKY）法。若两两正交，则矩阵成为对角形，即,21LIXT,2221LTXXA从而（63）的解变为，I1,LNKNKIITIIXYXYC1212/7数据拟合的最小二乘法3571问题的引入在工程实践和科学实验中，量与量之间的关系表现为1）确定性关系如电学中著名的欧姆定律就是确定性的关系。用V表示电压，R表示电阻，I表示电流，欧姆定律指出RIV有的确定性关系是由微分方程或积分方程来描述的。例如，阻尼振动中，位移X与时间的确定性关系由微分方程022XDTTX所决定，其中为固有频率，为阻尼系数。02）非确定性关系由于因素的复杂性或其它原因，变量之间找不到完全确定的关系。例如纱的回潮率与原棉含水量之间；钢水含碳量与冶炼时间；鱼的活动与海水温度；台风登陆路径与沿海各地的风向、气温、湿度之间等等，这些量之间，既存在密切关系，又不能由一个（或几个）变量的数值精确地求出另一个变量的值。但是通过人的实践，通过仪器，我们获得了大量的实验数据，这些大量的偶然现象，始终是服从内部隐藏着的规律的。现在我们又回到数值逼近范围内来谈这个问题，为了叙述方便，先讨论两个变量X，Y的情况。也就是说，通过观测变量X，Y积累了一组资料，I1，2，N，一般地说N都比较大。我们的任务是从积累得,I到的实验数据，I1，2，N，寻求一近似函数去逼近,IY。由于观测数据都带有观测误差，数组数目又较大，对于这类问,IYX题运用插值函数去描述Y往往是不适当的。我们以下面的例子来说明建36立近似函数的一种办法，即最小二乘法。X例71合成纤维抽丝工段，第一导丝盘的速度对丝的质量是很重要的参数，现发现它和电流周波有重要关系，由生产记录得到的数据如下表周波X492500493490490495498499502502第一导丝盘速度（Y）167170168166167168169170170171表4今要研究Y与X的关系，通常的步骤如下1）先用一坐标纸，将描于图上（图8）。,IY图712）凭视觉约略知道在一条直线的两侧附近，于是猜想Y与X,IYX37近似地成直线关系，BXAY上面直线关系式称为数学模型。在第I次观测数据中，与实测值有误IYI差，I1，2，N，IIIXY将它们平方后加起来得到总误差,2121BAIIINIINI我们当然希望，数学模型（主观猜想）应尽量接近客观实际，即总误差越小越好，也就是选取A,B使IA,B最小。问题又回到了内积空间21INI的最佳逼近。定义向量Y、X、E分别为，NY21NX21N维欧氏空间的一个子集SPANX，E对Y的最佳平方逼近就是选取A,B使向量BA的范数达到极小，即2XEYYT（72）MIN,21BAIXYIINI（72）与（71）的提法完全一致。从问题的来源看，确定使误差平38方和达到最小的方法称为最小二乘法，高斯对于天文观测数据的处理就运用了这个方法。依照我们建立的理论，A,B由下列方程组所决定NINIIIIIYAXBX112,解之有,1212NINIIIIIIXYYXA1212NINIIIIXYYB运用表4的数据求得A004，B0339,即3904XY必须指出，这里所说的欧氏空间最佳逼近，并不是说是YBXAY的最佳数学模型。因此就存在着另一个问题，上述数学模型是否符合客观实际呢也就是说，如何检查数学模型的质量呢当然，最根本的办法是拿到生产中去考验，但当观测数据积累多了以后，就能够建立一套数学理论去检验数学模型的准确性，偶然的现象，隐藏着必然的规律，概率统39计的课程里将介绍这个问题。例72某航空售票点（PVGIAD）机票销售加价和预期销售量之间的关系加价额度X（元）销售量P（张/周）加价额度X（元）销售量P（张/周）505822554755325052100592755012555300471506232535175623502720055对P和X的关系用三次多项式进行拟合，用最小二乘法可解得320268079043765813XXX40最小二乘法模型中的非线性函数例73已知及拟合这批数据的非线性数学模,0,1IXYN型（A,B为待定参数）BYAE1如何将非线性模型线性化2写出线性化模型中待定系数的法方程。3设数据,0,1IXYN如下I0123I201012100740004500IY2002712169073950449341求出拟合上述数据的非线性拟合函数。1、设则模型成为YLNBXAYLN2、设BBAA法方程式为XNIIINIIIYXBX1123、498150680712146973190682131001BABA最后拟合的非线性函数950XYXE10226E0042C,RENORMLSQNONLINFUN1117,X0LOCALMINIMUMFOUND0123Y2010121007400045000698101906030110798542C2007605008RENORM67226E005例74出钢时所用盛钢水的钢包，由于钢水对耐火材料的浸蚀，容积不断增大，我们希望找出使用次数与增大的容积之间的关系。试验数据如下表使用次数增大容积使用次数增大容积使用次数增大容积IXIYIXIYIXIY264271000121060382089931310804958999914106059501010491510906970111059161076解将数据标在坐标纸上（参见图72），我们看到，开始时浸蚀速度快，然后逐渐减弱，显然钢包容积不会无穷增加，于是可以想象它有一条平行于X轴的渐近线，根据这些特点我们选取数据拟合的曲线为双曲线。43图72假设选择的数学模型为,XBAY1令,1,XY于是上式变为BAY解得结果如下A00823,B01312从而,132081XY即132081XY44若将曲线拟合的双曲线模型改成指数形式,XBAEY将上式两边取对数LNX,L,1,LBAY则有XY求得107,4582BA从而,691EA最后求得78107XEY怎样比较数学模型的好坏呢双曲线模型的误差为,132081XY指数形式模式模型的误差为6791072XEY针对建立的模型比较实测值与拟合值的误差如下表45实测值拟合值拟合值误差误差实测值拟合值拟合值误差误差IY双曲模型指数模型1I2IIY双曲模型指数模型1I2I64267616702034102821049104791045100110039820793480650266013510591061210557002200339588687884708930733106010725106460125004695092129353028801471080108231072300230077970959997050101000510601090810788030801881000989699650104003510901098310845008300559931013110165020102351076110491089602890136999103221032203320333,1928026034112NII4365822212NII由于较小，所以我们选择指数模型作为钢包