




版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領
文檔簡介
線性擬合方法第1頁,課件共54頁,創(chuàng)作于2023年2月第五章實驗數(shù)據(jù)及模型參數(shù)擬合方法第一節(jié)問題的提出第2頁,課件共54頁,創(chuàng)作于2023年2月第一節(jié)問題的提出在化工設計及化工模擬計算中,需要大量的物性參數(shù)及各種設備參數(shù)。這些參數(shù)有些可以通過計算得到,但大量的參數(shù)還是要通過實驗測量得到。實驗測量得到的常常是一組離散數(shù)據(jù)序列(xi,yi)。
如果數(shù)據(jù)序列(xi,yi)(為一般起見),i=1,2,…,m,含有不可避免的誤差(或稱“噪聲”,如圖5-1所示),或者無法同時滿足某特定的函數(shù)(如圖5-2所示),那么,只能要求所作逼近函數(shù)ψ(x)最優(yōu)地靠近樣點,即向量Q=(ψ(x1),ψ(x2),…,ψ(xm))T與Y=(y1,y2,…,ym)T的誤差或距離最小。按Q與Y之間誤差最小原則作為“最優(yōu)”標準構造的逼近函數(shù),稱為擬合函數(shù)。第3頁,課件共54頁,創(chuàng)作于2023年2月第一節(jié)問題的提出圖5-1含有噪聲的數(shù)據(jù)圖5-2無法同時滿足某特定函數(shù)的數(shù)據(jù)序列第4頁,課件共54頁,創(chuàng)作于2023年2月第一節(jié)問題的提出除了物性數(shù)據(jù)及設備參數(shù)需要利用數(shù)據(jù)擬合外,在化學化工中,許多模型也要利用數(shù)據(jù)擬合技術,求出最佳的模型和模型參數(shù)。如在某一反應工程實驗中,我們測得了如表5-1所示的實驗數(shù)據(jù)。
現(xiàn)在要確定在其他條件不變的情況下,轉(zhuǎn)化率y和溫度T的具體關系,現(xiàn)擬用兩種模型去擬合實驗數(shù)據(jù),兩種模型分別是
第5頁,課件共54頁,創(chuàng)作于2023年2月第一節(jié)問題的提出如何求取上述模型中的參數(shù),并判斷兩種模型的優(yōu)劣,是化學化工工作者經(jīng)常要碰到的問題,這個問題的求解將在本章下面的有關章節(jié)中進行詳細的講解。
第6頁,課件共54頁,創(chuàng)作于2023年2月第二節(jié)擬合的標準第7頁,課件共54頁,創(chuàng)作于2023年2月第二節(jié)擬合的標準
前面已經(jīng)提到按Q與Y之間誤差最小原則作為“最優(yōu)”標準構造的逼近函數(shù),稱為擬合函數(shù),而向量Q與Y之間的誤差或距離有各種不同的定義方法,一般有以下幾種。
(1)用各點誤差絕對值的和表示
(2)用各點誤差按絕對值的最大值表示
(3)用各點誤差的平方和表示
第8頁,課件共54頁,創(chuàng)作于2023年2月式中R稱為均方誤差。由于計算均方誤差的最小值的原則容易實現(xiàn)而被廣泛采用。按均方誤差達到極小構造擬合曲線的方法稱為最小二乘法。同時還有許多種其他的方法構造擬合曲線,感興趣的讀者可參閱有關教材。本章主要講述用最小二乘法構造擬合曲線。第二節(jié)擬合的標準
第9頁,課件共54頁,創(chuàng)作于2023年2月第二節(jié)擬合的標準__實例1
實驗測得二甲醇(DME)的飽和蒸氣壓和溫度的關系,見表5-2。
序號溫度℃蒸氣壓MPa1-23.70.1012-100.174300.2544100.3595200.4956300.6627400.880表5-2DME飽和蒸氣壓和溫度的關系由表5-2的數(shù)據(jù)觀測可得,DME的飽和蒸氣壓和溫度有正相關關系,如果以函數(shù)p=a+bt來擬合,則擬合函數(shù)是一條直線。通過計算均方誤差Q(a,b)最小值而確定直線方程。(見圖5-3)圖5-3DME飽和蒸汽壓和溫度之間的線性擬合
第10頁,課件共54頁,創(chuàng)作于2023年2月第二節(jié)擬合的標準__實例1
擬合得到直線方程為:
相關系數(shù)R為0.97296,平均絕對偏差SD為0.0707。
第11頁,課件共54頁,創(chuàng)作于2023年2月擬合的標準
——
實例2如果采用二次擬合,通過計算下述均方誤差
擬合得二次方程為
相關系數(shù)R為0.99972,平均絕對偏差SD為0.00815,具體擬合曲線見圖5-4。圖5-4DME飽和蒸氣壓和溫度之間的二次擬合
第12頁,課件共54頁,創(chuàng)作于2023年2月擬合的標準實例2比較圖5-3和圖5-4以及各自的相關系數(shù)和平均絕對偏差可知:對于DME飽和蒸汽壓和溫度之間的關系,在實驗溫度范圍內(nèi)用二次擬合曲線優(yōu)于線性擬合。二次擬合曲線具有局限性,由圖5-4觀察可知,當溫度低于-30℃時,飽和壓力有升高的趨勢,但在擬合的溫度范圍內(nèi),二次擬合的平均絕對偏差又小于一次擬合,故對物性數(shù)據(jù)進行擬合時,不僅要看在擬合條件下的擬合效果,還必須根據(jù)物性的具體性質(zhì),判斷在擬合條件之外的物性變化趨勢,以便使擬合公式在已做實驗點數(shù)據(jù)之外應用。第13頁,課件共54頁,創(chuàng)作于2023年2月第三節(jié)單變量擬合和多變量擬合
給定一組數(shù)據(jù)(xi,yi),i=1,2,…,m,作擬合直線p(x)=a+bx,均方誤差為
由數(shù)學知識可知,Q(a,b)的極小值需滿足:
整理得到擬合曲線滿足的方程:
5.3.1單變量擬合1.單變量擬合——
線性擬合第14頁,課件共54頁,創(chuàng)作于2023年2月1.單變量擬合——線性擬合該方程可用消元法或克萊姆方法解出方程(如下式所示)
第15頁,課件共54頁,創(chuàng)作于2023年2月第16頁,課件共54頁,創(chuàng)作于2023年2月單變量擬合——
線性擬合實例例:下表為實驗測得的某一物性和溫度之間的關系數(shù)據(jù),表中x為溫度數(shù)據(jù),y為物性數(shù)據(jù)。請用線性函數(shù)擬合溫度和物性之間的關系。
解:設擬合直線p(x)=a+bx,并計算得下表:
x7911131517192123252729y91215182124273033363942x313335373941434547
y454851545760636669
編號xyxyx212345…21Σ79111315…47567912151821…6981963108165234315…3243267334981121169225…220918389第17頁,課件共54頁,創(chuàng)作于2023年2月將數(shù)據(jù)代入法方程組(1-12)中,得到:解方程得:a=-1.5,b=1.5
擬合直線為:
相關系數(shù)R為1。單變量擬合——
線性擬合實例第18頁,課件共54頁,創(chuàng)作于2023年2月線性擬合VB清單PrivateSubCommand1_Click()Dimx(5),y(5),c,d,m,p,a,b,eerConstn=5Fori=1To5x(i)=InputBox(“x(i)=”)y(i)=InputBox(“y(i)=”)Print“x(i)=”;x(i)Print“y(i)=”;y(i)NextIc=0d=0m=0p=0
Fori=1To5c=c+x(i)d=d+x(i)^2m=m+y(i)p=p+x(i)*y(i)Nextia=(m*d-c*p)/(n*d-c^2)b=(n*p-c*m)/(n*d-c^2)‘參數(shù)計算
a=Int(a*1000+0.5)/1000b=Int(b*1000+0.5)/1000Text1.Text=Str(a)Text2.Text=Str(b)‘參數(shù)輸出Fori=1To5eer=eer+(a+b*x(i)-y(i))^2‘誤差計算
eer=Int(eer*100000+0.5)/100000Nextieer=eer/5Text3.Text=Str(eer)EndSub第19頁,課件共54頁,創(chuàng)作于2023年2月第20頁,課件共54頁,創(chuàng)作于2023年2月有關線性擬合變型問題例如要擬合y=a+b/x2,只需在數(shù)據(jù)輸入后增加一語句x(i)=1/x(i)^2,而在程序后面的誤差eer的計算中則不需要修改。第21頁,課件共54頁,創(chuàng)作于2023年2月2.單變量擬合二次擬合函數(shù)給定數(shù)據(jù)(xi,yi),i=1,2,…,m,用二次多項式函數(shù)擬合這組數(shù)據(jù)。
設,作出擬合函數(shù)與數(shù)據(jù)序列的均方誤差表達式由數(shù)學知識可知,Q(a0,a1
,a2)的極小值滿足:整理左式得二次多項式函數(shù)擬合的滿足條件方程(5-14):
(5-14)第22頁,課件共54頁,創(chuàng)作于2023年2月解此方程得到在均方誤差最小意義下的擬合函數(shù)p(x)。式(5-14)稱為多項式擬合的法方程,法方程的系數(shù)矩陣是對稱的。當擬合多項式n>5時,法方程的系數(shù)矩陣是病態(tài)的,在用通常的迭代方法求解線性方程時會發(fā)散,在計算中要采用一些特殊算法以保護解的準確性。關于線性方程的求解方法,已在第三章中介紹。2.單變量擬合二次擬合函數(shù)第23頁,課件共54頁,創(chuàng)作于2023年2月3.二次擬合函數(shù)的拓展和一次擬合一樣,二次擬合也可以有多種變型,例如套用上面的公式,可以得到關于求解此擬合函數(shù)的法方程(5-15)。值得注意的是在此法方程的構建過程中,進行了變量的代換。首先是擬合函數(shù)中變量的代換:。
P(x)=a0+a1x3+a2x5
(5-15)第24頁,課件共54頁,創(chuàng)作于2023年2月其次是法方程的代換:將相應擬合函數(shù)中的代換引入法方程中。同時應注意法方程中x的4次冪是由兩個2次冪相乘得到,x的3次冪是由一個2次冪和一個1次冪相乘得到,而2次冪就是變量本身,而非兩個1次冪相乘得到。這個概念至關重要,在以后的二次擬合的各類變型中,均需利用這個概念,千萬不要用常規(guī)的思路去進行代入計算。
3.二次擬合函數(shù)的拓展如果我們需要求解是下面的擬合函數(shù):
第25頁,課件共54頁,創(chuàng)作于2023年2月參照上面的方法,我們很容易得到求解該擬合函數(shù)的法方程3.二次擬合函數(shù)的拓展第26頁,課件共54頁,創(chuàng)作于2023年2月4.二次擬合實例請用二次多項式函數(shù)擬合下面這組數(shù)據(jù)。解:設
并計算得下表序號xyxyX2X2yX3X41-34-12936-27812-22-448-8163-13-313-114000000051-1-11-11162-2-44-881673-5-159-452781Σ01-3928-70196第27頁,課件共54頁,創(chuàng)作于2023年2月4.二次擬合實例將上面數(shù)據(jù)代入式
(5-14),相應的法方程為解方程得:
a0=0.66667,a1=-1.39286,a2=-0.13095所以:-3-2-10123-6-4-2024y=0.66667-1.39286x-0.13095x2yYX圖
5-6擬合曲線與數(shù)據(jù)序列
第28頁,課件共54頁,創(chuàng)作于2023年2月二次擬合--VB程序清單PrivateSubCommand1_Click()Dimn,mAsIntegern=InputBox("n=方程次數(shù)")m=InputBox("m=實驗次數(shù)")ReDimx(m),y(m),a0(n+1),a1(n+1),aa(m,n+1)ReDimqq(n+1,m),pp(n+1,n+1),b(n+1),g(n+1),tt(n+1,n+1)'ForI=1Tom:'讀入數(shù)據(jù)'READx(I),y(I)'NextI'Data-3,4,-2,2,-1,3,0,0,1,-1,2,-2,3,-5Fori=1Tomx(i)=InputBox("x(i)")y(i)=InputBox("y(i)")Nexti第29頁,課件共54頁,創(chuàng)作于2023年2月二次擬合--VB程序清單omiga=InputBox("omiga=松弛因子")Forj=1Ton+1'為計算法方程中的系數(shù)做準備
Fori=1Tomaa(i,j)=x(i)^(j-1)NextiNextjForj=1TomFori=1Ton+1qq(i,j)=aa(j,i)NextiNextj'計算法方程中的右邊項Fori=1Ton+1b(i)=0Forj=1Tomb(i)=b(i)+aa(j,i)*y(j)NextjNexti'開始計算法方程中的右邊項的系數(shù)Fori=1Ton+1Forj=1Ton+1pp(i,j)=0Fork=1Tompp(i,j)=qq(i,k)*aa(k,j)+pp(i,j)NextkNextjNexti第30頁,課件共54頁,創(chuàng)作于2023年2月二次擬合--VB程序清單'開始用松弛迭代法求解法方程中的變量Fori=1Ton+1a0(i)=0a1(i)=0.11NextiFori=1Ton+1g(i)=b(i)/pp(i,i)Forj=1Ton+1Ifj=iThentt(i,j)=0Elsett(i,j)=-pp(i,j)/pp(i,i)EndIfNextjNexti50eer=0Fori=1Ton+1eer=eer+Abs(a0(i)-a1(i))NextiIfeer<0.00001ThenGoTo100ElseFori=1Ton+1a0(i)=a1(i)Nexti第31頁,課件共54頁,創(chuàng)作于2023年2月二次擬合--VB程序清單'關鍵的迭代計算公式Fori=1Ton+1s=g(i)Forj=1Ton+1s=s+tt(i,j)*a0(j)Nextja1(i)=omiga*s+(1-omiga)*a0(i)'Print"a1("&i&")="&a1(i)NextiGoTo50EndIf'打打印結(jié)果100Fori=1Ton+1Print"a("&i-1&")="&a1(i)NextiEndSub第32頁,課件共54頁,創(chuàng)作于2023年2月3.2多變量的曲線擬合
前面介紹的曲線擬合方法只涉及單變量函數(shù)的曲線擬合,但實際在化工實驗數(shù)據(jù)處理及模型參數(shù)擬合時,通常會碰到多變量的參數(shù)擬合問題。一個典型的例子是傳熱實驗中努塞爾數(shù)、雷諾數(shù)及普朗特數(shù)之間的擬合問題:
(1-16)
根據(jù)若干組實驗測得的數(shù)據(jù),如何求出式(1-16)中的參數(shù)c1、c2、c3,這是一個有2個變量的參數(shù)擬合問題,為不失一般性,我們把它表達成以下形式。給定數(shù)據(jù)序列用一次多項式函數(shù)擬合這組數(shù)據(jù)。
設,作出擬合函數(shù)與數(shù)據(jù)序列的均方誤差
(1-17)
(x1i,x2i,yi),i=1,2,3…,m
第33頁,課件共54頁,創(chuàng)作于2023年2月多變量的曲線擬合由多元函數(shù)的極值原理,Q(a0
,a1,a2
)的極小值滿足
整理得多變量一次多項式函數(shù)擬合的法方程
(1-18)第34頁,課件共54頁,創(chuàng)作于2023年2月通過求解方程(1-18)就可以得到多變量函數(shù)線性擬合時的參數(shù),由于方程(1-16)不是線性方程,我們可以通過對方程(1-16)兩邊同取對數(shù),就可以得到以下線性方程只要作如下變量代換:并將實驗數(shù)據(jù)代入法方程(1-18)就可以求出方程(1-16)中的系數(shù)。對于變量數(shù)多于2個,并且擬合曲線模型是非線性型時,可參照本節(jié)的方法,推導得到法方程,通過對法方程的求解就可以求得各種擬合曲線參數(shù)。靈活運用上面介紹的方法,可以解決大部分實驗數(shù)據(jù)及模型參數(shù)的擬合問題。多變量的曲線擬合第35頁,課件共54頁,創(chuàng)作于2023年2月多變量的曲線擬合—VB程序清單第36頁,課件共54頁,創(chuàng)作于2023年2月多變量的曲線擬合實例根據(jù)某傳熱實驗測得如下數(shù)據(jù),請用方程(1-16)的形式擬合實驗曲線。
解:利用上面的VB程序,將數(shù)據(jù)依次輸入,就可以得到方程(1-16)中的三個參數(shù)
C1=0.023C2=0.8C3=0.3
則式(1-16)就變成了常見的光滑管傳熱方程第37頁,課件共54頁,創(chuàng)作于2023年2月多變量的曲線擬合實例值得注意的是程序中對c2(1)的處理,不是直接將計算結(jié)果顯示出來,而是進行指數(shù)運算后才顯示出來。這是由于我們在進行擬合計算的時候,對方程(1-16)進行了對數(shù)運算。如果擬合方程的形式和方程(1-16)不同,則需對上面提供的程序作適當修改。例如以下兩個自變量的擬合函數(shù)第38頁,課件共54頁,創(chuàng)作于2023年2月第四節(jié)解矛盾方程組第39頁,課件共54頁,創(chuàng)作于2023年2月第四節(jié)解矛盾方程組本節(jié)中將用最小二乘法求解線性矛盾方程的方法來構造擬合函數(shù),并將其推廣至任意次和任意多個變量的擬合函數(shù),為在化學化工中實驗數(shù)據(jù)處理及模型參數(shù)擬合提供更為一般性的方法。給定數(shù)據(jù)序列(xi,yi),i=1,2,…,m
,做擬合直線p(x)=a0
+a1x,如果要直線p(x)過這些點,那么就有
p(xi)=a0+a1xi
=yi,i=1,2,…,m,即寫成矩陣形式為第40頁,課件共54頁,創(chuàng)作于2023年2月上述方程組中有2個未知量m個方程(m>>2)。一般地,將含有n個未知量m個方程的線性方程組其一般形式為第四節(jié)解矛盾方程組寫成矩陣形為
第41頁,課件共54頁,創(chuàng)作于2023年2月一般情況下,當方程數(shù)m多于變量數(shù)n
,且m個方程之間線性無關,則方程組無解,這時方程組稱為矛盾方程組。方程組在一般意義下無解,也即無法找到n個變量同時滿足m個方程。這種情況和擬合曲線無法同時滿足所有的實驗數(shù)據(jù)點相仿,故可以通過求解均方誤差極小意義下矛盾方程的解來獲取擬合曲線。由數(shù)學知識還可證明:方程組ATAX=ATY的解就是矛盾方程組AX=Y
在最小二乘法意義下的解。這樣我們只要通過求解ATAX=ATY就可以得到矛盾方程組的解,進而得到各種擬合曲線,為擬合曲線的求解提供了另一種方法。第四節(jié)解矛盾方程組min║AX-B║22第42頁,課件共54頁,創(chuàng)作于2023年2月例如,擬合直線p(x)=a0+a1x的矛盾方程組ATAX=ATY的形式如下:
化簡得到與式(1-12)相同的法方程第四節(jié)解矛盾方程組第43頁,課件共54頁,創(chuàng)作于2023年2月這里需要注意的是變量X和系數(shù)(a0,
a1)之間的相互轉(zhuǎn)換關系。即對于n次多項式曲線擬合,要計算Q(a0,a1,…,an)
的極小值問題,這與解矛盾方程組第四節(jié)解矛盾方程組與求
的極小問題是一回事。或第44頁,課件共54頁,創(chuàng)作于2023年2月在這里
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年綿陽市游仙區(qū)定向招聘社區(qū)專職工作人員考試筆試試題(含答案)
- 2025年法務管理專業(yè)考試試題及答案
- 2025年城市規(guī)劃與設計考試試題及答案
- 汽修廠焊接火花防護措施執(zhí)行規(guī)定
- 汽修廠漆面流掛起泡處理規(guī)范制度
- 安全在我心中演講稿15篇
- 品質(zhì)組長年度工作總結(jié)
- 危險化學品運輸事故應急預案 8篇
- 廣東汕尾圖書館招聘試題帶答案分析2024年
- 安全施工考試試題及答案
- 房顫患者用藥護理
- 2025年全國統(tǒng)一高考數(shù)學試卷(全國二卷)
- 玻璃幕墻工程技術規(guī)范 JGJ 102-2003 知識培訓
- 2024年04月高等教育自學考試《00067財務管理學》試題
- 2025年中考英語考前沖刺卷(北京卷)(解析版)
- 農(nóng)村房產(chǎn)抵押協(xié)議書
- 2025年中國主控操作鍵盤市場現(xiàn)狀分析及前景預測報告
- 境外旅游包車協(xié)議書
- 分揀物流考試試題及答案
- 抱負與使命主題范文引領-2024-2025學年高一語文單元寫作深度指導(統(tǒng)編版必修下冊)
- 簽訂屠宰協(xié)議書范本
評論
0/150
提交評論