(19)國家知識產(chǎn)權(quán)局(12)發(fā)明專利39-7號國家電網(wǎng)有限公司國網(wǎng)信息通信產(chǎn)業(yè)集團有限公司謝可邱鎮(zhèn)劉勁松范維張亮張彬劉佳鑫廖逍王興濤宋純賀董雪情段方維崔世界有限公司21100專利代理師何學軍李叢GO6F30/23(2020.0審查員張曉娜一種基于Krylov子空間的變壓器溫度場模本發(fā)明公開了一種基于Krylov子空間的變較，誤差≥0.01%驗證降階模型是否符合誤差要求2步驟2:基于步驟1中所述簡化仿真物理模型，選用Maxwell3D分析模塊構(gòu)建電磁場仿真模H(s)=-C(A-s。E)?1B-C(A-s?E)?1E)(A-s?E)?1B(sT表示對矩陣求轉(zhuǎn)置；步驟43.根據(jù)ArnoldiWcRnx9,其中q<<n,如下式：colspan{V}=K,(A?-S?E?)?1E?;算法構(gòu)造出他們各自的標準列正交矩陣V,3模型的傳遞函數(shù)保持原系統(tǒng)的前r?+r?階距，模型由n階變成q階；步驟44.驗證模型是否符合要求，包括：將降階后的方程形變得到狀態(tài)空間方程，將矩陣輸入到MATLABsimulink平臺中的狀態(tài)空間方程模塊，記錄實驗數(shù)據(jù)，并與仿真前的全階模型的計算數(shù)據(jù)進行比較，如果誤差在0.01%以內(nèi)，則符合要求。2.根據(jù)權(quán)利要求1所述的一種基于Krylov子空間的變壓器溫度場模型降階方法，其特征是：所述根據(jù)磁場不變準則將繞組等效成圓環(huán)，具體是根據(jù)變壓器容量、電壓參數(shù)，推算出變壓器的幾何尺寸，建立變壓器仿真模型，根據(jù)磁場不變準則，對繞組進行簡化，用圓環(huán)代替繞組。3.根據(jù)權(quán)利要求1所述的一種基于Krylov子空間的變壓器溫度場模型降階方法，其特征是：所述構(gòu)建的電磁場仿真模型為電磁-流體多物理場耦合仿真模型，包括如下步驟：步驟(1)添加鐵芯、繞組參數(shù)、磁導率，構(gòu)建外部等效電路作為激勵源，選擇渦流求解步驟(2)將電磁場仿真模型以及求解結(jié)果導入到流體場中，在流體場中構(gòu)建流體區(qū)域命名出入口以及邊界，對電磁-流體多物理場耦合仿真模型進行分部網(wǎng)格拋分，獲得良好的步驟(3)在流體場檢查電磁場仿真模型是否符合標準，網(wǎng)格質(zhì)量是否良好，并選擇能量場和湍流場作為求解模型，把繞組作為流體場的激勵源，建立材料庫設(shè)定繞組以及變壓器油的參數(shù)；步驟(5)對電磁-流體多物理場耦合仿真模型進行求解。4.根據(jù)權(quán)利要求1所述的一種基于Krylov子空間的變壓器溫度場模型降階方法，其特征是：所述構(gòu)建的電磁場仿真模型為電磁-流體多物理場耦合仿真模型，是在ANSYS軟件中步驟a.電磁場仿真獲取電壓器損耗的實現(xiàn)；步驟b.電磁場--流體場獲取電壓器溫度場分布的實現(xiàn)。5.根據(jù)權(quán)利要求4所述的一種基于Krylov子空間的變壓器溫度場模型降階方法，其特②選取渦流場求解器，根據(jù)設(shè)計的線圈匝數(shù)，在圓環(huán)截面上添加繞組，并添加外電路作為仿真模型的激勵源；③對模型進行仿真，獲得模型的損耗分布。6.根據(jù)權(quán)利要求4所述的一種基于Krylov子空間的變壓器溫度場模型降階方法，其特征是：所述電磁場--流體場獲取電壓器溫度場分布的實現(xiàn)，包括：①將電磁場的模型以及求解結(jié)果導入到流體場中，在流體場中構(gòu)建流體區(qū)域，命名出入口以及邊界，對模型進行分部網(wǎng)格拋分，獲得②在流體場檢查模型是否符合標準，網(wǎng)格質(zhì)量是否良好，并選擇能量場和湍流場作為4③變壓器內(nèi)部的熱量傳遞主要是通過熱傳導方式來進行的，通過電磁場求解出來的渦流損耗，選取繞組作為流體場的激勵源，建立材料庫對繞組以及變壓器油的參數(shù)；⑤對模型進行求解。7.根據(jù)權(quán)利要求5所述的一種基于Krylov子空間的變壓器溫度場模型降階方法，其特8.一種計算機設(shè)備，包括存儲介質(zhì)、處理器及存儲在存儲介質(zhì)上并可在處理器上運行的計算機程序，其特征在于，所述處理器執(zhí)行所述計算機程序時實現(xiàn)權(quán)利要求1-7中任一權(quán)利要求所述的一種基于Krylov子空間的變壓器溫度場模型降階方法的步驟。9.一種計算機存儲介質(zhì)，其特征是：所述計算機存儲介質(zhì)上存有計算機程序，所述計算機程序被處理器執(zhí)行時實現(xiàn)權(quán)利要求1-7中任一權(quán)利要求所述的一種基于Krylov子空間的變壓器溫度場模型降階方法的步驟。5一種基于Krylov子空間的變壓器溫度場模型降階方法技術(shù)領(lǐng)域[0001]本發(fā)明屬于油浸式三相變壓器溫度場仿真技術(shù)領(lǐng)域，尤其涉及一種基于Krylov子空間的變壓器溫度場模型降階方法，更具體的是一種基于Krylov子空間的變壓器溫度場數(shù)字孿生模型降階方法。背景技術(shù)[0002]數(shù)字孿生技術(shù)從原本的航天航空領(lǐng)域逐步發(fā)展到各個制造業(yè)當中，并在智能制造領(lǐng)域表現(xiàn)出良好的應用前景。軍、民等多領(lǐng)域?qū)?shù)字孿生技術(shù)越來越重視，并且開始逐步探索其核心技術(shù)和開發(fā)潛力，數(shù)字孿生技術(shù)作為橋梁連接著物理世界和虛擬世界，以復雜的物理仿真，實時的數(shù)據(jù)共享與分析，數(shù)據(jù)處理等為關(guān)鍵技術(shù)，構(gòu)建物理世界與虛擬世界的數(shù)字孿生體，并實時顯現(xiàn)物理世界的物理狀態(tài).,數(shù)字孿生體能夠高效、精確地模擬物理實體，分析物理實體的狀態(tài)，對故障提前預警，并輔助操作人員的決策制定。[0003]變壓器作為變電站的重要組成部分，其熱特性直接影響設(shè)備能否安全可靠運行。變壓器運行過程中，繞組和鐵芯在漏磁場的作用下均會產(chǎn)生渦流損耗，繞組內(nèi)部電流與線圈電阻相互作用產(chǎn)生歐姆損耗，這些損耗均轉(zhuǎn)化為熱量，并通過變壓器油將熱量傳遞至外部環(huán)境。隨著變壓器的等級和容量的增大，其損耗與溫度也會逐漸增長。變壓器溫度升高，會加速絕緣的老化，降低變壓器壽命，故應當提高對繞組熱點溫度測量的準確性，保證其安全穩(wěn)定運行。當前變壓器測溫主要在壁面，油溫。對繞組熱點測溫的方法有在繞組內(nèi)部埋光纖傳感器，但無法準確定位熱點位置。因此，對變壓器進行溫度場仿真可以準確計算各點的溫度值。目前物理場仿真雖然能準確的模擬處變壓器各部分的溫度，但是仿真耗時長，計算量大，不適用于構(gòu)建變壓器數(shù)字孿生體。[0004]例如：現(xiàn)有專利號為2021108220457,名稱為一種換流變壓器溫度場模型構(gòu)建方法，該申請采用動態(tài)模態(tài)分解法對溫度場搭建降階模型，同樣建立在已經(jīng)計算完畢的全階模型基礎(chǔ)上，但進一步需根據(jù)計算獲得的離散時刻溫度數(shù)據(jù)樣本組成瞬像矩陣，構(gòu)建瞬像矩陣對離散時刻點的選取要求具備相關(guān)性，且樣本數(shù)量也直接影響計算精度。技術(shù)的ARNOLDI模型降階方法，該文章雖然也是采用Krylov子空間與Arnoldi算法實現(xiàn)模型降階，但該文章提出傳統(tǒng)Arnoldi算法在復雜動力系統(tǒng)降階中存在無法一次性獲得穩(wěn)定的降階系統(tǒng)的問題。[0006]因此，針對現(xiàn)有技術(shù)中存在的不足之處，對物理場構(gòu)建降階模型，在滿足誤差的范圍內(nèi)通過犧牲部分計算精度來減少計算量，大幅度減少計算時間，可以解決目前物理場計算存在的缺陷。發(fā)明內(nèi)容[0007]針對上述現(xiàn)有技術(shù)中存在的不足之處，本發(fā)明提供了一種基于Krylov子空間的變壓器溫度場模型降階方法。其目的是為了實現(xiàn)通過降階模型加快仿真效率，減少仿真耗時，6更好的適用于數(shù)字孿生體的搭建的發(fā)明目的。[0008]本發(fā)明為實現(xiàn)上述目的所采用的技術(shù)方案是：[0009]一種基于Krylov子空間的變壓器溫度場模型降階方法，包括以下步驟：[0010]步驟1:建立油浸式三相變壓器仿真物理模型，后續(xù)物理場仿真均基于所述簡化仿真物理模型完成計算，具體是根據(jù)變壓器實際物理尺寸，構(gòu)建1:1仿真模型；由于實際變壓器內(nèi)部結(jié)構(gòu)復雜，而本發(fā)明涉及的場仿真僅與鐵芯繞組構(gòu)件有關(guān)，因此對模型進行簡化，根據(jù)磁場不變準則將繞組等效成圓環(huán)；[0011]步驟2:基于步驟1中所述簡化仿真物理模型，選用Maxwell3D分析模塊構(gòu)建電磁場[0012]步驟3:將獲得的變壓器損耗分布作為流體場的激勵，建立流體場仿真模型，實現(xiàn)電磁場與流體場耦合分析，獲得油浸式三相變壓器的溫度場分布特性；基于變壓器電磁-流體場仿真計算獲得溫度分布，將計算結(jié)果導出至MATLABSimulink平臺完成后續(xù)溫度場降階模型搭建；[0013]步驟4:基于Krylov子空間的方法對變壓器溫度場搭建降階模型進行降階并驗證。[0014]更進一步的，所述根據(jù)磁場不變準則將繞組等效成圓環(huán)，具體是根據(jù)變壓器容量、電壓參數(shù)，推算出變壓器的幾何尺寸，建立變壓器仿真模型，根據(jù)磁場不變準則，對繞組進[0015]更進一步的，所述構(gòu)建的電磁場仿真模型為電磁-流體多物理場耦合仿真模型，包括如下步驟：解器，通過電路與磁場耦合仿真，獲得變壓器的損耗分布；[0017]步驟(2)將電磁場仿真模型以及求解結(jié)果導入到流體場中，在流體場中構(gòu)建流體區(qū)域命名出入口以及邊界，對電磁-流體多物理場耦合仿真模型進行分部網(wǎng)格拋分，獲得良好的網(wǎng)格；[0018]步驟(3)在流體場檢查電磁場仿真模型是否符合標準，網(wǎng)格質(zhì)量是否良好，并選擇能量場和湍流場作為求解模型，把繞組作為流體場的激勵源，建立材料庫設(shè)定繞組以及變壓器油的參數(shù)；[0020]步驟(5)對電磁-流體多物理場耦合仿真模型進行求解。[0021]更進一步的，所述基于Krylov子空間的方法對獲得的油浸式三相變壓器溫度場搭建降階模型，在MATLABsimulink平臺計算降階模型，與電磁-流體場仿真結(jié)果進行對比，驗[0022]步驟41.利用Fluent中的功能命令HBMAT,使Fluent以Harwell-boeing格式輸出整熱容矩陣，保存在輸出文件中；在Matlab中利用讀取文件命令，從輸出文件中得到熱傳導矩陣和載荷列陣及熱容矩陣，從而得到溫度場穩(wěn)態(tài)方程，如式(1)所示：7為對狀態(tài)變量求一階導數(shù)，上述即狀態(tài)方程。[0025]對系統(tǒng)進行拉氏變換，獲得傳遞函數(shù)H()=C(sE-A)1B,在s?處對傳遞函數(shù)進行Taylor展開，如式(2)所示：[0027]上式中：由第二項開始依次定義C((A-s。E)?1E)"(A-s,E)?1B=M,n=1,2,...,i為系統(tǒng)的第n階矩，H(s)為傳遞函數(shù)在s。處的泰勒級數(shù)展開式；[0028]步驟42.一個r維Krylov子空間K由1個正定矩陣A和1個向量b組成，即由一組基向量組成定義子空間表達式如下：K(A,b)=span{b,Ab,A?1b};對于歸一化系統(tǒng)構(gòu)造以下兩個子空間，如式(3)所示：矩陣，T表示對矩陣求轉(zhuǎn)置；[0031]步驟43.根據(jù)Arnoldi算法構(gòu)造出他們各自的標準列正交矩陣V,WcRnx9,其中q<<n,如式(4)所示：矩陣，得到原始系統(tǒng)的降階模型，如式(5)所示：是狀態(tài)變量，(t)∈R"是輸出變量，u(t)為輸入向量，E,A,B為降階后的輸入矩陣；V和W滿秩降階模型的傳遞函數(shù)保持原系統(tǒng)的前r?+r?階距，模型由n階變成q階；[0036]步驟44.驗證模型是否符合要求，包括以下步驟：[0037]將降階后的方程形變得到狀態(tài)空間方程，將矩陣輸入到MATLABsimulink平臺中的狀態(tài)空間方程模塊，記錄實驗數(shù)據(jù)，并與仿真前的全階模型的計算數(shù)據(jù)進行比較，如果誤差在0.01%以內(nèi)，則符合要求。[0038]更進一步的，所述構(gòu)建的電磁場仿真模型為電磁-流體多物理場耦合仿真模型，是在ANSYS軟件中建立電磁-流體多物理場耦合分析、設(shè)置材料，邊界條件、求解域；包括：[0039]步驟a.電磁場仿真獲取電壓器損耗的實現(xiàn)；[0040]步驟b.電磁場--流體場獲取電壓器溫度場分布的實現(xiàn)。[0041]更進一步的，所述電磁場仿真獲取電壓器損耗的實現(xiàn)，包括：[0042]①對模型鐵心，繞組設(shè)置參數(shù)，設(shè)置材料屬性；[0043]②選取渦流場求解器，根據(jù)設(shè)計的線圈匝數(shù)，在圓環(huán)截面上添加繞組，并添加外電8路作為仿真模型的激勵源；[0044]③對模型進行仿真，獲得模型的損耗分布。[0046]①將電磁場的模型以及求解結(jié)果導入到流體場中，在流體場中構(gòu)建流體區(qū)域，命名出入口以及邊界，對模型進行分部網(wǎng)格拋分，獲得良好的網(wǎng)格；[0047]②在流體場檢查模型是否符合標準，網(wǎng)格質(zhì)量是否良好，并選擇能量場和湍流場作為求解模型；[0048]③變壓器內(nèi)部的熱量傳遞主要是通過熱傳導方式來進行的，通過電磁場求解出來的渦流損耗，選取繞組作為流體場的激勵源，建立材料庫對繞組以及變壓器油的參數(shù)；[0050]⑤對模型進行求解。[0052]一種計算機設(shè)備，包括存儲介質(zhì)、處理器及存儲在存儲介質(zhì)上并可在處理器上運行的計算機程序，所述處理器執(zhí)行所述計算機程序時實現(xiàn)任一所述的一種基于Krylov子空間的變壓器溫度場模型降階方法的步驟。[0053]一種計算機存儲介質(zhì)，所述計算機存儲介質(zhì)上存有計算機程序，所述計算機程序被處理器執(zhí)行時實現(xiàn)任一所述的一種基于Krylov子空間的變壓器溫度場模型降階方法的步驟。[0054]本發(fā)明具有以下有益效果及優(yōu)點：[0055]本發(fā)明以繞組為研究主體，通過電磁-流體多物理場耦合分析，獲得變壓器繞組的溫度場分布。模型降階的思想是將物理模型的大規(guī)模狀態(tài)空間投影到以一組基矢量特征的低維空間。針對上述系統(tǒng)，本發(fā)明采用Krylov子空間的方法對系統(tǒng)進行降階，從Fluent的結(jié)果文件中獲取系統(tǒng)的整體矩陣的剛度、質(zhì)量、以及阻尼矩陣，基于上述矩陣可以構(gòu)造系統(tǒng)的狀態(tài)方程。因此，對系統(tǒng)降階的問題就可以轉(zhuǎn)化成對狀態(tài)方程降階。針對此方程的降階，保持降階前后兩個系統(tǒng)的傳遞函數(shù)的矩最大程度匹配，而Arnoldi算法則很好地解決了直接矩匹配方式存在數(shù)值不穩(wěn)定性的問題，更好的實現(xiàn)矩的匹配。最后，通過Arnoldi算法建立上述兩個子空間的標準列正交矩陣，從而得到本狀態(tài)方程的降階方程，實現(xiàn)模型的降階。[0056]本發(fā)明運用該降階模型計算電磁-流體耦合場仿真，其計算效率顯著提升，適用于變壓器溫度場的快速仿真計算，進而構(gòu)建變壓器數(shù)字孿生模型。數(shù)字孿生技術(shù)是實現(xiàn)電力裝備行業(yè)數(shù)字化轉(zhuǎn)型的重要手段，基于Krylov子空間的變壓器溫度場模型降階方法，可以實現(xiàn)物理場快速仿真計算，進一步構(gòu)建變壓器數(shù)字孿生模型。[0057]本發(fā)明通過降階模型加快仿真效率，減少仿真耗時，本降階模型能更好的適用于數(shù)字孿生體的搭建。附圖說明[0058]本發(fā)明的上述和/或附加的方面和優(yōu)點從結(jié)合下面附圖對實施例的描述中將變得[0059]圖1是本發(fā)明基于Krylov子空間的三相變壓器溫度場降階模型的構(gòu)建流程圖；[0060]圖2是本發(fā)明三相變壓器損耗分布圖；9[0061]圖3是本發(fā)明電磁-流體場的實現(xiàn)圖；[0062]圖4是本發(fā)明溫度場分布圖。具體實施方式[0063]為了能夠更清楚地理解本發(fā)明的上述目的、特征和優(yōu)點，下面將結(jié)合附圖和具體實施方式對本發(fā)明進行進一步的詳細描述。需要說明的是，在不沖突的情況下，本申請的實施例及實施例中的特征可以相互組合。[0064]在下面的描述中闡述了很多具體細節(jié)以便于充分理解本發(fā)明，但是，本發(fā)明還可以采用其他不同于在此描述的其他方式來實施，因此，本發(fā)明的保護范圍并不受下面公開的具體實施例的限制。[0065]下面參照圖1-圖4描述本發(fā)明一些實施例的技術(shù)方案。[0067]本發(fā)明提供了一個實施例，是一種基于Krylov子空間的變壓器溫度場模型降階方法，在傳統(tǒng)的有限元分析仿真的基礎(chǔ)上，對系統(tǒng)的模型進行降階處理，實現(xiàn)更快速的相應，更適合于數(shù)字孿生體對模型的要求。[0068]如圖1所示，圖1是本發(fā)明基于Krylov子空間的三相變壓器溫度場降階模型的構(gòu)建流程圖。[0069]本發(fā)明提供一種基于Krylov子空間的變壓器溫度場模型降階方法。具體包括以下[0070]步驟1:建立油浸式三相變壓器仿真物理模型，后續(xù)物理場仿真均基于所述簡化仿真物理模型完成計算，具體是根據(jù)變壓器實際物理尺寸，構(gòu)建1:1仿真模型；由于實際變壓器內(nèi)部結(jié)構(gòu)復雜，而本發(fā)明涉及的場仿真僅與鐵芯繞組構(gòu)件有關(guān)，因此對模型進行簡化，根據(jù)磁場不變準則將繞組等效成圓環(huán)。[0071]所述根據(jù)磁場不變準則將繞組等效成圓環(huán)，具體是根據(jù)變壓器容量、電壓參數(shù)，推算出變壓器的幾何尺寸，建立變壓器仿真模型，根據(jù)磁場不變準則，對繞組進行簡化，用圓環(huán)代替繞組。[0072]步驟2:基于步驟1中所述簡化仿真物理模型，選用Maxwell3D分析模塊構(gòu)建電磁場仿真模型，采用渦流求解器，獲得變壓器損耗分布。如圖2所示，圖2是本發(fā)明三相變壓器損耗分布圖。[0073]所構(gòu)建的電磁場仿真模型為電磁-流體多物理場耦合仿真模型，包括如下步驟：選擇渦流求解器，通過電路與磁場耦合仿真，獲得變壓器的損耗分布。[0075]步驟(2)將電磁場仿真模型以及求解結(jié)果導入到流體場中，并且在流體場中構(gòu)建流體區(qū)域命名出入口以及邊界，對電磁-流體多物理場耦合仿真模型進行分部網(wǎng)格拋分，獲得良好的網(wǎng)格。[0076]步驟(3)在流體場檢查電磁場仿真模型是否符合標準，網(wǎng)格質(zhì)量是否良好，并選擇能量場和湍流場作為求解模型，忽略鐵芯，把繞組作為流體場的激勵源，建立材料庫設(shè)定繞流體場的實現(xiàn)圖?；谧儔浩麟姶?流體場仿真計算獲得溫度分布，將計算結(jié)果導出至MATLABSimulink平臺完成后續(xù)溫度場降階模型搭[0080]步驟4:基于Krylov子空間的方法對獲得的油浸式三相變壓器溫度場搭建降階模利用讀取文件命令，可以從Cond_File.dat文件得到熱傳導矩陣和載荷列陣，從Cap_-s?E)?1B(s-so)-……(2)決直接矩匹配方式存在數(shù)值不穩(wěn)定性的問題。[0089]步驟42.一個r維Krylov子空間K由1個正定矩陣A和1個向量b組成，同時也被認為降階模型的傳遞函數(shù)保持原系統(tǒng)的前r?+r?階距，模型由n階變成q階。[0098]將降階后的方程形變得到狀態(tài)空間方程，將矩陣輸入到MATLABsimulink平臺中的狀態(tài)空間方程模塊，記錄實驗數(shù)據(jù)，并與仿真前的全階模型的計算數(shù)據(jù)進行比較，如果誤差在0.01%以內(nèi)，則符合要求。[0100]本發(fā)明提供了一個實施例，是一種基于Krylov子空間的變壓器溫度場模型降階方[0101]步驟1.建立油浸式三相變壓器仿真簡化模型，根據(jù)磁場不變準則將繞組等效成圓[0102]建立油浸式三相變壓器仿真簡化模型，是根據(jù)變壓器容量、電壓等參數(shù)，可以推算出變壓器的幾何尺寸，建立變壓器仿真模型，根據(jù)磁場不變準則，對繞組進行簡化，即用圓環(huán)代替繞組。解域等；[0104]在ANSYS軟件中建立電磁-流體多物理場耦合分析，ANSYSworkbench有限元分析軟真流程更加緊密地組合起來，通過簡單的步驟即可完成復雜的多物理場分析流程，具體包[0105]步驟a.電磁場仿真獲取電壓器損耗的實現(xiàn)。[0107]②選取渦流場求解器，根據(jù)設(shè)計的線圈匝數(shù)，在圓環(huán)截面上添加繞組，并添加外電路作為仿真模型的激勵源。[0109]步驟b.電磁場--流體場獲取電壓器溫度場分布的實現(xiàn)。[0110]①將電磁場的模型以及求解結(jié)果導入到流體場中，并且在流體場中構(gòu)建流體區(qū)[0111]②在流體場檢查模型是否符合標準，網(wǎng)格質(zhì)量是否良好，并選擇能量場和湍流場作為求解模型。[0112]③變壓器內(nèi)部的熱量傳遞主要是通過熱傳導方式來進行的，通過電磁場求解出來的渦流損耗，選取繞組作為流體場的激勵源，建立材料庫對繞組以及變壓器油的參數(shù)，如傳[0115]步驟3.基于Krylov子空間的方法對變壓器溫度場搭建降階模型，對系統(tǒng)模型降[0117]Fluent流體流動分析模塊中的功能命令HBMAT,可以使Fluent以Harwell-boeing文件格式輸出整體矩陣，且為稀疏矩陣，通過編程還原為滿矩陣。再運用命令流分別從ANSYS的.FULL文件中提取熱傳導矩陣和熱容矩陣，并保存在相應的輸出文件Cond_到熱傳導矩陣和載荷列陣，從Cap_File.dat文件得到熱容矩陣.從而得到溫度場的穩(wěn)態(tài)方[0127]可以根據(jù)Arnoldi算法構(gòu)造出他們各自的標準列正交矩陣V,WcR"x,其中q<<colspan{V}=K,(A?-s?E?)?1E?;colspan{V}=K,(A?-s?E?)?1E?;是狀態(tài)變量，(t)∈R"是輸出變量，u(t)為輸入向量，E,A,B為降階后的輸入矩陣。V和W滿秩降階模型的傳遞函數(shù)保持原系統(tǒng)的前r?+r?階距，模型由n階變成q階。[0132]步驟4.對于降階模型進行驗證，并提出降階模型的驗證方法，通過分析誤差判定降階模型是否合理。[0133]降階模型的驗證方法具體是將降階后的方程形變得到狀態(tài)空間方程，將矩陣輸入到MATLABsimulink平臺中的狀態(tài)空間方程模塊，記錄實驗數(shù)據(jù)，并與仿真前的全階模型的計算數(shù)據(jù)進行比較，如果誤差在0.01%以內(nèi)，則符合要求。[0135]基于同一發(fā)明構(gòu)思，本發(fā)明實施例還提供了一種計算機設(shè)備，包括存儲介質(zhì)、處理器及存儲在存儲介質(zhì)上并可在處理器上運行的計算機程序。所述處理器執(zhí)行所述計算機程序時實現(xiàn)實施例1或2所述的任意一種基于Krylov子空間的變壓器溫度場模型降階方法的[0137]基于同一發(fā)明構(gòu)思，本發(fā)明實施例還提供了一種計算機存儲介質(zhì)，所述計算機存