第6章面向方程模擬法N

第6章面向方程模擬法第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）第一節(jié)原理EquationOrientedMethod聯(lián)立方程模擬法面向方程法將整個(gè)系統(tǒng)聯(lián)合建模Subflow(S1,ALFA,S2,S3,S4)S2=S1+S4’以上混合器模型方程

S3=S2*ALFAS4=S2*(1–ALFA)’以上

分割器模型方程EndSub整個(gè)系統(tǒng)方程組，整體求解

序貫法模：每個(gè)單元為獨(dú)立子程序?qū)ν馓峁┙Y(jié)果需要迭代循環(huán)物流的去哪里？過程系統(tǒng)的模型方程模型方程：模型方程決策變量狀態(tài)變量①各種平衡方程②單元聯(lián)結(jié)方程③設(shè)計(jì)規(guī)定方程④物性方程⑤現(xiàn)象方程模型變量決策變量所有單元模塊的設(shè)備參數(shù)進(jìn)料流股變量狀態(tài)參數(shù)所有中間變量產(chǎn)品流股變量?jī)?nèi)部變量等模擬結(jié)構(gòu)循環(huán)物流方程設(shè)計(jì)規(guī)定方程僅提供方程多!!!特殊處理物性系統(tǒng)的處理占用70~80%的計(jì)算量處理方法：1）不提供方程組，僅提供子程序方程數(shù)量大幅減少，但需經(jīng)常調(diào)用物性子程序，效果有限2）提供k,h

的方程組，其它物性由子程序調(diào)用

k,h最常用，其它較少用到。效果很好物性：焓、相平衡常數(shù)、密度、粘度、導(dǎo)熱系數(shù)等第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）第二節(jié)面向方程模擬法與序貫?zāi)M法的比較模擬方法結(jié)構(gòu)特征序貫?zāi)K模擬法面向方程模擬法1.序貫?zāi)K法優(yōu)點(diǎn)：

①與實(shí)際過程的直觀聯(lián)系強(qiáng)，軟件的建立、維護(hù)和擴(kuò)充很方便；

②易于通用化；

③需要的計(jì)算機(jī)內(nèi)存較??；

④易于診斷出錯(cuò)位置；缺點(diǎn)：

循環(huán)物流及設(shè)計(jì)問題：慢2.面向方程法優(yōu)點(diǎn)

①便于實(shí)際用戶的提出的質(zhì)量性能要求。有設(shè)計(jì)規(guī)定的系統(tǒng)，僅增加幾個(gè)方程

②用空間換取時(shí)間。相同的過程單元出現(xiàn)多次，將占用更多的存儲(chǔ)空間，自由度的選擇也可能各不相同

③以空間換取求解的數(shù)值穩(wěn)定性。有較多回流高度交互作用流程相對(duì)容易收斂2.面向方程法缺點(diǎn)

①要求提供較好的初值

②難于自動(dòng)產(chǎn)生流程的模型方程組，難以通用化

③不容易診斷出現(xiàn)錯(cuò)誤的地方。序貫法與聯(lián)立方程法對(duì)比內(nèi)容序貫法聯(lián)立方程法占用存儲(chǔ)空間小大迭代循環(huán)圈多少計(jì)算效率低高指定設(shè)計(jì)規(guī)定不靈活靈活初值要求低高錯(cuò)誤診斷易難編程、調(diào)試較易較難聯(lián)立模塊法序貫?zāi)K法建立流程簡(jiǎn)單，處理循環(huán)流與設(shè)計(jì)規(guī)定需迭代聯(lián)立方程法同時(shí)求解，但建模困難聯(lián)立模塊法又稱雙層法，將過程系統(tǒng)的近似模型方程與單元模塊交替求解聯(lián)立模塊法兼有序貫?zāi)K法和面向方程法的優(yōu)點(diǎn)。既能使用序貫?zāi)K法積累的大量模塊，又能將最費(fèi)計(jì)算時(shí)間的流程收斂和設(shè)計(jì)約束收斂等迭代循環(huán)合并處理，通過聯(lián)立求解達(dá)到同時(shí)收斂斷裂變量不迭代簡(jiǎn)化模型的產(chǎn)生??？聯(lián)立模塊法特點(diǎn)以簡(jiǎn)化模型聯(lián)立求解取代序貫法的回路迭代計(jì)算和設(shè)計(jì)規(guī)定計(jì)算可利用原有的序貫法資源方程數(shù)量較聯(lián)立方程法要少得多，求解難度大大降低難點(diǎn)：簡(jiǎn)化方程的產(chǎn)生及其適用范圍過程系統(tǒng)穩(wěn)態(tài)模擬三種方法的比較方法優(yōu)點(diǎn)缺點(diǎn)代表軟件系統(tǒng)序貫?zāi)K法與工程師直觀經(jīng)驗(yàn)一致，便于學(xué)習(xí)使用；易于通用化，已積累了豐富的單元模塊；需要計(jì)算機(jī)內(nèi)存較小；有錯(cuò)誤易于診斷檢查；再循環(huán)引起的收斂迭代很費(fèi)機(jī)時(shí)；進(jìn)行設(shè)計(jì)型計(jì)算時(shí)，很費(fèi)機(jī)時(shí)；不宜用于最優(yōu)化計(jì)算；PROCESS(美)CONCEPT(英)CAPES(日)ASPEN(美)FLOWTRAN(美)面向方程法解算快；模擬型計(jì)算與設(shè)計(jì)型計(jì)算一樣；適合最優(yōu)化計(jì)算，效率高；便于與動(dòng)態(tài)模擬聯(lián)合實(shí)現(xiàn)；要求給定較好的初值，否則可能得不到解；計(jì)算失敗后診斷錯(cuò)誤所在困難；形成通用化程序有困難有，故使用不便；難以繼承已有的單元操作模塊。ASCEND-Ⅱ(美)SPEEDUP(英)（雙層法）聯(lián)立模塊法可以利用前人開發(fā)的單元操作模塊；可以避免序貫?zāi)K法中的循環(huán)流迭代；比較容易實(shí)現(xiàn)通用。將嚴(yán)格模型做成簡(jiǎn)化模型時(shí)，需要花費(fèi)機(jī)時(shí)；用簡(jiǎn)化模型來尋求優(yōu)化時(shí)，其解與嚴(yán)格與嚴(yán)格模型優(yōu)化解是否一致，有爭(zhēng)論。TISFLO（德）FLOWPACK-Ⅱ(英)第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）第三節(jié)大型稀疏非線性方程組的降維解法系統(tǒng)方程多，整體求解難度很多處理方法：方程組分解過程系統(tǒng)方程組的特點(diǎn)： 方程數(shù)多，變量多每個(gè)方程包含的變量不多每個(gè)變量存在的方程(出現(xiàn)次數(shù))不多稀疏方程組稀疏比將方程組轉(zhuǎn)化為有向圖，再應(yīng)用系統(tǒng)分解的手段分解成可順序求解的子方程組第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）

一、輸出變量的指定Hall各異條件？方程組有解的必要條件！輸出變量？介紹方程組決策變量影響時(shí)學(xué)過！將f2(x,y,z)=0改寫為y=f2’(x,z)，y則稱為f2的輸出變量N維方程組可指定N個(gè)不同的輸出變量失?。2出現(xiàn)最少，應(yīng)先考慮f3包含的變量最少，也應(yīng)優(yōu)先考慮！成功！包含變量少的方程和出現(xiàn)次數(shù)少的變量應(yīng)先考慮輸出變量的指定方法：1)選事件矩陣中元素最少的行和其與元素最少的列的交點(diǎn)處元素對(duì)應(yīng)的變量作為優(yōu)先指定的輸出變量，然后從事件矩陣中刪去相應(yīng)的行和列。2)重復(fù)上述過程∶①若矩陣的所有行和列被刪除，則指定完畢。②若有行或列無法刪除，則表示與剩余列對(duì)應(yīng)的變量不存在于與剩余的行對(duì)應(yīng)的方程中。“Steward”通路∶a)從矩陣的不飽和行(無輸出變量的行)或不飽和列(無輸出變量的列)的某一非輸出變量開始，垂直找到與該元素在一列或一行的輸出變量。再轉(zhuǎn)90，找到另一個(gè)非輸出變量，再垂直找到輸出變量，直到找到不飽和列或行中的某一非零元素。b)此軌跡就是“Steward”通路。將此通路上的輸出變量與非輸出變量互換，即可增加一個(gè)輸出變量。f1f2f3指定輸出變量的作用？將方程組轉(zhuǎn)化為有向圖！可應(yīng)用過程系統(tǒng)分解的手段對(duì)方程系統(tǒng)進(jìn)行分解可及矩陣法索引矩陣法Steward通路搜索法Sargent-Westerberg搜索法（圖解法）例×不飽和行不飽和列×從不飽和行的任一非零開始將出現(xiàn)死循環(huán)將Steward通路上的輸出變量與非輸出變量互換可增加一個(gè)輸出變量！Steward通路例第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）

二、回路搜索法分解方程組回路搜索法分解方程組步驟:1）寫出方程組事件矩陣，指定輸出變量

2）如果每個(gè)方程都指定了一個(gè)不同的輸出變量，則方程組滿足Hall各異條件，轉(zhuǎn)6；否則轉(zhuǎn)3。

3）從不飽和行的一個(gè)非零開始搜索Steward通路。如果能從不飽和行走到不飽和列，轉(zhuǎn)4；否則轉(zhuǎn)5。4）不能增加輸出變量數(shù)目，方程組無解，停止。

5）將Steward通路上的輸出變量與非輸出變量類型互換，則可增加一個(gè)輸出變量。如果此時(shí)輸出變量的數(shù)目與方程數(shù)目相等，轉(zhuǎn)6；否則轉(zhuǎn)3）

6）將方程以節(jié)點(diǎn)形式排成一行，將每個(gè)方程與其輸出變量出現(xiàn)的方程節(jié)點(diǎn)以有向弧相連，即形成了方程組的有向圖。

7）可以用可及矩陣法、索引矩陣法、Steward通路搜索法和圖解法等分解方程組。不飽和行不飽和列0從不飽和行開始找Steward通路不飽和列將Steward通路上的輸出變量與非輸出變量對(duì)換指定輸出變量：最終的輸出變量指定：f1繪制有向圖：f2f3f4f5f1x4f1f2x3f2f3x2f3f4x1f4f5x5f5f1f4f5f2f3搜索回路：f2f3f1f2f5f2

，

f2f5構(gòu)成回路，將f2f5作為組合節(jié)點(diǎn)，改寫有向圖f1f2f5f3f4f2f5f1f4f5內(nèi)部弧內(nèi)部弧f1f2f3f4f5搜索回路：f1f2f5，無路可走！刪除f2f5，計(jì)入次序表，改寫有向圖f1f3f4次序表f2f5f2f5f2f5f1f3，無路可走！刪除f3，計(jì)入次序表，改寫有向圖f3f3f1f4f1，構(gòu)成回路，形成最后一個(gè)組合節(jié)點(diǎn)，刪除f1f4，計(jì)入次序表f1f4計(jì)算順序第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）

三、不可分解方程組的斷裂降維解法大規(guī)模方程組的處理：分解成小規(guī)模方程組，順序求解問題：子方程組規(guī)模仍很大整個(gè)方程組本身不可分解？必須聯(lián)立求解！？斷裂降維求解方程組的斷裂求解步驟如下∶

1)選輸出變量，分解方程組；

2)在分解后的子方程組中選擇包含變量數(shù)最少的方程(=k)中的k-1變量作為斷裂變量；選擇k-1個(gè)包含變量數(shù)最多的方程作為驗(yàn)算方程；

3)降維后用回路搜索法進(jìn)一步分解；

4)給斷裂變量賦初值，求解降維后的方程組；

5)驗(yàn)算斷裂變量直到收斂。例整個(gè)方程組不可分解斷裂降維求解––迭代求解迭代求解––要求最少的斷裂變量！––選變量數(shù)最少的方程斷裂斷裂x1，即給定x1初值x3=f4’(x1)斷裂變量，賦初值斷裂1個(gè)變量，可排除一個(gè)方程，2個(gè)變量例斷裂變量狀態(tài)變量驗(yàn)算方程驗(yàn)算方程應(yīng)選擇含變量數(shù)最多的方程方程數(shù)比變量數(shù)多1個(gè)指定輸出變量分解方程組f1f3f5f6求解順序：f5f6f1f3計(jì)算步驟：假設(shè)x1初值由f4計(jì)算出x3聯(lián)合求解f1

f3

，得出x2,x4將x1,x2,x3,x4代入f5

f6聯(lián)合求解得x5,x6驗(yàn)算f2是否滿足，若滿足則結(jié)束，否則轉(zhuǎn)6)重新假設(shè)x1初值轉(zhuǎn)2)四、n<m型稀疏非方程組的處理選決策變量獲得的方程組的結(jié)構(gòu)方程組進(jìn)行分解或斷裂降維按順序聯(lián)立求解獲迭代求解分解后的子方程組五、n>m型稀疏非方程組的處理①利用最小二乘法取得一組妥協(xié)解。②求解其中的m個(gè)方程第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）第四節(jié)稀疏矩陣的壓縮存儲(chǔ)n=1000，N=2000，=0.2%全部存儲(chǔ)，占用n2即106存儲(chǔ)單位其中998000個(gè)為零一次高斯消去運(yùn)算需1/3n3次加法和乘法運(yùn)算，約3.33億次以1微秒/次計(jì)，耗機(jī)時(shí)333秒以上采用稀疏矩陣存儲(chǔ)法只存儲(chǔ)非零元素，只對(duì)非零元素進(jìn)行運(yùn)算需占用10000個(gè)存儲(chǔ)單位，進(jìn)行20000次運(yùn)算需20毫秒第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針非零元素屬性∶l

行標(biāo)；l

列標(biāo)；l

指向本行中下一個(gè)非零元素的指針；l

指向本列中下一個(gè)非零元素的指針；l

指向本行中前一個(gè)非零元素的指針；l

指向本列中前一個(gè)非零元素的指針；第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）

二、靜態(tài)存儲(chǔ)（列表法）一、靜態(tài)存儲(chǔ)(列表法)需三個(gè)數(shù)組：LA矩陣元素值；LB非零矩陣元素的列標(biāo)；LC每行的第一個(gè)非零元素在LA中的存儲(chǔ)位置；特點(diǎn)：按行或按列順序存儲(chǔ)例LA-346-11233LB12431434LC1457

元素aij訪問步驟∶1)找i行的第一個(gè)非零元素的存儲(chǔ)位置，K1=LC(i);2)找i

行的最后一個(gè)非零元素的存儲(chǔ)位置，K2=LC(i+1)1;3)比較該行的非零元素下標(biāo)，找出aij。①令k=K1，若LB(k)=j，則說明已找到aij；aij=LA(k)；②否則，令k=k+1，返回①；最多從d1做到d2，若找不到，則aij值即為零例:訪問a34(i=3,j=4)LA-346-11233LB12431434LC1457

1）找第3行第一個(gè)非零位置：

LC(3)=5(d1)2）找第3行最后一個(gè)非零位置：

LC(3+1)–1=7–1=6(d2)3）對(duì)比列標(biāo)：

LB(5)=1≠j4）對(duì)比本行下一非零的列標(biāo)：

LB(5+1)=4=j5）已找到aij

，取元素值：

LA(5+1)=2第一節(jié)原理第二節(jié)面向方程模擬法與序貫?zāi)M法的比較第三節(jié)大型稀疏非線性方程組的降維解法

一、輸出變量的指定

二、回路搜索法分解方程組

三、不可分解方程組的斷裂降維解法第四節(jié)稀疏矩陣的壓縮存儲(chǔ)

一、稀疏矩陣壓縮存儲(chǔ)的信息鏈指針

二、靜態(tài)存儲(chǔ)（列表法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）

三、動(dòng)態(tài)存儲(chǔ)（三地址法）1)Bending-Hutchison法采用屬性①和②及非零元素值構(gòu)成鏈表包括五列∶①序號(hào)∶可作為存儲(chǔ)數(shù)組的下標(biāo)數(shù)值，不占用存儲(chǔ)空間；②變量號(hào)(列標(biāo))；③方程號(hào)(行標(biāo))；④非零元素的數(shù)值；⑤非零元素的狀態(tài)屬性。2)三地址法（按列存儲(chǔ)）由屬性①、②和④及非零元素值構(gòu)成：每列第一個(gè)非零元素在存儲(chǔ)鏈中的位置

──列首址BC[1，n]；2)信息條INF[1∶3(N+M)]；信息條由三個(gè)單元構(gòu)成∶①非零元素值的行號(hào)；②非零元素的數(shù)值；③該列的下一個(gè)非零元素存儲(chǔ)位置。3)第一個(gè)空白區(qū)的位置─空白區(qū)首址S[L]；123d1d2d3例123456789101112131415161718192021222324列首址BC：信息條INF：空白區(qū)首址S：13行號(hào)元素值下一非零位置本列結(jié)束開始下列1-34350072-1224301320最后一列去空白區(qū)190結(jié)束按列存儲(chǔ)位置動(dòng)態(tài)！1016(1)找aij①計(jì)算第j列的第一個(gè)非零元素存儲(chǔ)位置∶d=BC(j)；②找j列中的第一個(gè)非零元素所在行∶k=INF(d)；③若k=i，則aij=INF(d+1)；

ki，令d=INF(d+2)，轉(zhuǎn)②；例:找a431234567891011121314151617181920212223243203502-12201-34101643013071列首址BC：信息條INF：空白區(qū)首址S：19第3列第一個(gè)非零對(duì)比行號(hào)不等于4去下行對(duì)比行號(hào)找到啦!取值i=4j=3(2)插入aij，i1<i<i2①找到d1，d2;(j列中與aij相鄰的上下兩個(gè)非零位置)②令L=S=d3；③S=INF(L+2);(將d4作為第一個(gè)空白區(qū))④INF(L)=i

(在d3保存行號(hào))INF(L+1)=aij(保存元素值)INF(L+2)=d2(將ai,j與ai2,j相連)⑤INF(d1+2)=L

(將ai1,j與ai,j相連)Sd3…d1d1+1d1