本公開涉及電磁場數(shù)值計算領(lǐng)域，特別涉及一種對于時域電磁場進行仿真計算的方法與裝置。
背景技術(shù)：
：：時域電磁場仿真計算能夠提供時域信息，便于分析時域脈沖、頻域?qū)拵У碾姶彭憫?，并且能夠?qū)r變媒質(zhì)進行處理，有其獨特的應用優(yōu)勢。常見的時域電磁場數(shù)值計算方法有時域有限差分(finite-differencetime-domain，fdtd)、時域有限元(finite-elementtime-domain，fetd)、時域有限體積(finite-volumetime-domain，fvtd)、時域不連續(xù)加略金(discontinuousgalerkintime-domain，dgtd)等方法。時域電磁場的仿真計算按時間步迭代的方式分為顯式迭代和隱式迭代。其中顯式迭代由于易于實現(xiàn)，計算量小而獲得了廣泛應用。一些時域顯式迭代計算方法具有：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的特點。如dgtd方法，在邊界處通過數(shù)值流與相鄰單元邊界的電磁場產(chǎn)生能量交換。相比頻域電磁場計算，時域電磁場計算多了一個時間維度，所以未知量更多，計算量更大。為了提高時域電磁場的計算效率，如果當前時間歩的電磁場分布滿足一定要求，則可以對計算方法進行相應的近似，從而減少計算量，加快求解速度。但在決定可否選用近似算法時，當前時間的電磁場分布是未知的，只有前一時間步的電磁場是已知的。如果在不知道當前時間電磁場分布的情況下，盲目采用近似算法，將帶來較大誤差，并有可能導致算法不穩(wěn)定。在文獻【3】中，我們提出一種方法，通過cfl條件和因果性來對下一時間步的電磁場分布做出預測。該方法基于cfl條件和因果性，因此適用性強，適用于一切受cfl和因果性限制的顯示時域迭代方法。但是對于一些特殊方法，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響，如dgtd方法。對這些方法而言，使用文獻【3】中的方法判斷效率低，且有一些本來可以簡化的單元被錯誤判斷，進一步影響效率。參考文獻：【1】j.-m.jin,thefiniteelementmethodinelectromagnetics,thirdedition.hoboken.newjersey:johnwiley&sonsinc,2014.【2】y.j.zhangetal.,“anefficienthybridfinite-elementanalysisofmultipleviassharingthesameanti-padinanarbitrarilyshapedparallel-platepair,”ieeetrans.microw.theorytech.,vol.63,no.3,pp.883–890,mar.2015.【3】w.mai,j.hu,p.li,andh.zhao,“anefficientandstable2-d/3-dhybriddiscontinuousgalerkintime-domainanalysiswithadaptivecriterionforarbitrarilyshapedantipadsindispersiveparallel-platepair,”ieeetrans.microw.theorytech.,vol.pp,no.99,pp.1–11,2017技術(shù)實現(xiàn)要素：為了用更高效的方法，使用前一時間步的電磁場值來對當前時間歩的電磁場值做出預測，進而判斷當前時間步可否采用近似算法，從而提高效率，本公開實施例提供了一種顯式時間迭代時域電磁場計算方法和裝置。所述技術(shù)方案如下：一種在時域電磁場數(shù)值計算中的預測方法，該時域電磁場數(shù)值計算對時間采用顯式迭代，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響，其中：由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測。其有益效果包括：利用前一時間歩的電磁場值直接對當前時間歩的電磁場值做出預測，計算量小。一種在時域電磁場數(shù)值計算中的預測裝置，該時域電磁場數(shù)值計算對時間采用顯式迭代，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響，其中：由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測。其有益效果包括：利用前一時間歩的電磁場值直接對當前時間歩的電磁場值做出預測，計算量小。一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代時域電磁場計算方法，其特征在與，所述方法包括：將待求網(wǎng)格相鄰邊界上的電磁場值組合，形成一個與待求網(wǎng)格形狀相同，位置相同的鏡像網(wǎng)格；將每個待求網(wǎng)格和相應鏡像網(wǎng)格前一時間步的電磁場值與判斷條件比較，如果都滿足判斷條件，則該網(wǎng)格屬于區(qū)域二，否則屬于區(qū)域一；對區(qū)域一中的所有待求網(wǎng)格采用算法一，對區(qū)域二中的所有待求網(wǎng)格采用算法二；所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其有益效果包括：計算量小、判斷準確、精度和穩(wěn)定性好。在一些實施方式中，其中，所述判斷條件為：待求網(wǎng)格中的電、磁場值基本上為零；所述算法二為：無操作，或?qū)⒋缶W(wǎng)格電、磁場值置零。其有益效果包括：計算量小、判斷準確、精度和穩(wěn)定性好。在一些實施方式中，其中，所述判斷條件為：待求網(wǎng)格單元中的電場只主要存在于一個方向，稱為主要電場方向，且沿該方向電場值基本不變，磁場只主要存在于與主要電場方向正交的平面上；所述算法一采用棱柱網(wǎng)格，所述算法二采用多邊形網(wǎng)格；其中多邊形網(wǎng)格是棱柱網(wǎng)格的直截面，與棱柱網(wǎng)格一一對應；每層棱柱網(wǎng)格對應一層多邊形網(wǎng)格；所述算法一為：遍歷待求網(wǎng)格單元的相鄰網(wǎng)格單元，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元都屬于區(qū)域一，則對待求網(wǎng)格單元直接采用三維算法求解，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元有網(wǎng)格單元屬于區(qū)域二，則將該區(qū)域二網(wǎng)格單元前一時間步的二維網(wǎng)格電磁場值轉(zhuǎn)化為三維網(wǎng)格電磁場值之后再對待求網(wǎng)格單元采用三維算法求解；其中三維算法在任一時間步任一網(wǎng)格單元的待求電磁場值只由前一時間步本網(wǎng)格單元和相鄰邊界單元的電磁場值決定；所述算法二為：遍歷待求網(wǎng)格單元的相鄰網(wǎng)格單元，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元都屬于區(qū)域二，則對待求網(wǎng)格單元直接采用二維算法求解，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元有網(wǎng)格單元屬于區(qū)域一，則將該區(qū)域一網(wǎng)格單元前一時間步的三維網(wǎng)格電磁場值轉(zhuǎn)化為二維網(wǎng)格電磁場值之后再對待求網(wǎng)格單元采用二維算法求解；其中二維算法在任一時間步任一網(wǎng)格單元的待求電磁場值只由前一時間步本網(wǎng)格單元和相鄰邊界單元的電磁場值決定。其有益效果包括：計算量小、便于并行計算、判斷準確、精度和穩(wěn)定性好。在一些實施方式中，其中，所述棱柱網(wǎng)格為三棱柱網(wǎng)格；所述多邊形網(wǎng)格為三角形網(wǎng)格；所述三維算法和二維算法采用時域不連續(xù)迦略金算法求解。其有益效果包括：計算量小、便于并行計算、判斷準確、精度和穩(wěn)定性好。一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代時域電磁場計算裝置，其中，所述裝置包括：預處理模塊：將待求網(wǎng)格相鄰邊界上的電磁場值組合，形成一個與待求網(wǎng)格形狀相同，位置相同的鏡像網(wǎng)格；判斷模塊：將每個待求網(wǎng)格和相應鏡像網(wǎng)格前一時間步的電磁場值與判斷條件比較，如果都滿足判斷條件，則該網(wǎng)格屬于區(qū)域二，否則屬于區(qū)域一；計算模塊：對區(qū)域一中的所有網(wǎng)格采用算法一，對區(qū)域二中的所有網(wǎng)格采用算法二；所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其有益效果包括：計算量小、判斷準確、精度和穩(wěn)定性好。在一些實施方式中，其中，所述判斷條件為：待求網(wǎng)格中的電、磁場值基本上為零；所述算法二為：無操作，或?qū)⒋缶W(wǎng)格電、磁場值置零。其有益效果包括：計算量小、判斷準確、精度和穩(wěn)定性好。在一些實施方式中，其中，所述判斷條件為：待求網(wǎng)格單元中的電場只主要存在于一個方向，稱為主要電場方向，且沿該方向電場值基本不變，磁場只主要存在于與主要電場方向正交的平面上；所述算法一采用棱柱網(wǎng)格，所述算法二采用多邊形網(wǎng)格；其中多邊形網(wǎng)格是棱柱網(wǎng)格的直截面，與棱柱網(wǎng)格一一對應；每層棱柱網(wǎng)格對應一層多邊形網(wǎng)格；所述算法一為：遍歷待求網(wǎng)格單元的相鄰網(wǎng)格單元，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元都屬于區(qū)域一，則對待求網(wǎng)格單元直接采用三維算法求解，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元有網(wǎng)格單元屬于區(qū)域二，則將該區(qū)域二網(wǎng)格單元前一時間步的二維網(wǎng)格電磁場值轉(zhuǎn)化為三維網(wǎng)格電磁場值之后再對待求網(wǎng)格單元采用三維算法求解；其中三維算法在任一時間步任一網(wǎng)格單元的待求電磁場值只由前一時間步本網(wǎng)格單元和相鄰邊界單元的電磁場值決定；所述算法二為：遍歷待求網(wǎng)格單元的相鄰網(wǎng)格單元，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元都屬于區(qū)域二，則對待求網(wǎng)格單元直接采用二維算法求解，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元有網(wǎng)格單元屬于區(qū)域一，則將該區(qū)域一網(wǎng)格單元前一時間步的三維網(wǎng)格電磁場值轉(zhuǎn)化為二維網(wǎng)格電磁場值之后再對待求網(wǎng)格單元采用二維算法求解；其中二維算法在任一時間步任一網(wǎng)格單元的待求電磁場值只由前一時間步本網(wǎng)格單元和相鄰邊界單元的電磁場值決定。其有益效果包括：計算量小、便于并行計算、判斷準確、精度和穩(wěn)定性好。在一些實施方式中，其中，所述棱柱網(wǎng)格為三棱柱網(wǎng)格；所述多邊形網(wǎng)格為三角形網(wǎng)格；所述三維算法和二維算法采用時域不連續(xù)迦略金算法求解。其有益效果包括：計算量小、便于并行計算、判斷準確、精度和穩(wěn)定性好。附圖說明圖1是本公開實施例2提供的一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代時域電磁場計算方法的流程圖。圖2是本公開實施例3進行101預處理的示意圖，其中201是四面體網(wǎng)格的待求網(wǎng)格和待求網(wǎng)格的相鄰邊界的電磁場值，202是由四面體待求網(wǎng)格的相鄰邊界的電磁場值組成的鏡像網(wǎng)格；203是三棱柱網(wǎng)格的待求網(wǎng)格和待求網(wǎng)格的響鈴邊界的電磁場值，204是三棱柱待求網(wǎng)格的相鄰邊界的電磁場值組成的鏡像網(wǎng)格。圖3是本公開實施例3的計算區(qū)域示意圖，其中301是整個計算區(qū)域及網(wǎng)格劃分，302是t＝tn-1時刻經(jīng)過與101預處理，將待求網(wǎng)格和鏡像網(wǎng)格前一時間歩的電磁場值與102判斷條件比較后不都滿足的區(qū)域一示意圖，即簡化后的計算區(qū)域示意圖。圖4是本公開實施例4計算的多導孔平行板結(jié)構(gòu)示意圖，其中401是俯視圖，402是平視圖，403是隔離盤，404是金屬導孔。圖5是本公開實施例4提供的一種基于三棱柱網(wǎng)格和對應下頂面三角形網(wǎng)格的一階棱邊基函數(shù)位置示意圖，其中，501為三棱柱一階棱邊基函數(shù)，502為三角形一階棱邊基函數(shù)，503為簡化的三角形一階棱邊基函數(shù)。圖6是本公開實施例4提供的一種基于長方體網(wǎng)格和對應下頂面矩形網(wǎng)格的一階棱邊基函數(shù)位置示意圖，其中，601為長方體一階棱邊基函數(shù)，602為矩形一階棱邊基函數(shù)，603為簡化的矩形一階棱邊基函數(shù)。圖7是在本公開實施例4提供的，在區(qū)域二，由同一三角形網(wǎng)格拉伸而成的兩層三棱柱網(wǎng)格的三棱柱一階棱邊基函數(shù)，其中701為兩層三棱柱網(wǎng)格的三棱柱一階棱邊基函數(shù)，702為簡化的單層三棱柱網(wǎng)格的三棱柱一階棱邊基函數(shù)。圖8是本公開實施例5提供的一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算裝置的框圖。圖9是本公開又一實施例提供的一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算裝置的框圖。具體實施方式為使本公開的目的、技術(shù)方案和優(yōu)點更加清楚，下面將結(jié)合附圖對本公開實施方式作進一步詳細描述。實施例1一種時域電磁場數(shù)值計算中的預測方法，該時域電磁場數(shù)值計算對時間采用顯式迭代，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響，其特征在于：由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測。其中時域電磁場數(shù)值計算包括但不限于dgtd等數(shù)值計算方法，本公開實施例對此不做具體限定。統(tǒng)一特點是，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響其中對時間采用顯式迭代包括但不限于二階蛙跳、四階龍格庫塔等迭代方式，本公開實施例對此不做具體限定。其中由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測的具體方法，本公開實施例給出以下優(yōu)選實施例，并不用以限制本公開。例如：前一時間歩待求網(wǎng)格與待求網(wǎng)格的相鄰邊界的電磁場值都小于一判據(jù)時，即認為在前一時間歩，待求網(wǎng)格與待求網(wǎng)格的相鄰邊界的電磁場值基本上為零，故可以預測當前時間歩待求網(wǎng)格的電磁場值也為零。其中判據(jù)可以由經(jīng)驗給出，也可由控制誤差的算法給出，本公開實施例對此不做具體限定。本實施例提供的上述方法，由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測。本實施例提供的上述方法適用廣泛，計算量少。實施例2參見圖1，本實施例提供了一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算方法，包括：101將待求網(wǎng)格相鄰邊界上的電磁場值組合，形成一個與待求網(wǎng)格形狀相同，位置相同的鏡像網(wǎng)格；102將每個待求網(wǎng)格和相應鏡像網(wǎng)格前一時間步的電磁場值與判斷條件比較，如果滿足判斷條件，則該網(wǎng)格屬于區(qū)域二，否則屬于區(qū)域一；103對區(qū)域一中的所有網(wǎng)格采用算法一，對區(qū)域二中的所有網(wǎng)格采用算法二；所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其中顯式迭代包括但不限于二階蛙跳、四階龍格庫塔等迭代方式，本公開實施例對此不做具體限定。在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算方法有很多種，本實施例不作具體限定，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。其中網(wǎng)格可以是任意形狀，本公開實施例對此不做具體限定。其中對每個網(wǎng)格的處理可以采用逐一處理的方式，也可采用并行處理的方式，本公開實施例對此不做具體限定。其中101待求網(wǎng)格和鏡像網(wǎng)格可以是任意形狀，基函數(shù)可以是各種形式，本實施例分別以四面體網(wǎng)格、三棱柱網(wǎng)格，棱邊即函數(shù)為例說明，并不用以限制本發(fā)明。如圖201和203所示待求網(wǎng)格i，其電磁場值為ei，與其相鄰的其它網(wǎng)格的邊界上的電磁場值為e′,e″。將這些邊界上的電磁場值組合起來形成鏡像網(wǎng)格ei*。鏡像網(wǎng)格ei*的組合方式有多種，如圖202和204，本公開實施例給出兩種優(yōu)選例，如ei*＝e′+e″或ei*＝(e′+e″)/2，并不用以限制本發(fā)明。其中103中的算法一是當前時間歩、待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，包括但不限于dgtd等數(shù)值計算方法，本公開實施例對此不做具體限定，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。其中103中的算法二是當待求網(wǎng)格當前時間步的電磁場值滿足102中判斷條件時方可采用的一種近似算法。算法二與判斷條件對應，例如：例一：102中判斷條件為待求網(wǎng)格最大電磁場值小于一判據(jù)，即認為待求網(wǎng)格中的電磁場值基本上為零。則對應的103中算法二為無操作或電磁場值置零；例二：102中判斷條件為網(wǎng)格中的電場只主要存在于一個方向，稱作主要電場方向，且沿該方向電場強值基本不變，磁場只主要存在于與主要電場方向正交的平面上，具體地，電場在兩個正交方向的最大值小于一判據(jù)，磁場在另一個正交方向的最大值小于一判據(jù)。則認為網(wǎng)格單元中的電磁場值可以用二維網(wǎng)格準確擬合，對應的103中的算法二為二維算法。其中二維算法包括但不限于dgtd等數(shù)值計算方法，本公開實施例對此不做具體限定，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。相反，如果待求網(wǎng)格當前時間歩的電磁場值不滿足102中的判斷條件，例如：待求網(wǎng)格中有電磁場值大于經(jīng)驗判據(jù)，即不能認為待求網(wǎng)格電磁場值基本上為零，此時采用算法二，無操作或電磁場值置零，將帶來不可控的誤差，從而影響精度和穩(wěn)定性；待求網(wǎng)格中電磁場值不滿足：電場只主要存在于一個方向，稱作主要電場方向，且沿該方向電場強值基本不變，磁場只主要存在于與主要電場方向正交的平面上。此時待求網(wǎng)格中電磁場值不能用二維網(wǎng)格準確擬合，采用算法二，二維算法，將帶來不可控誤差，從而影響精度和穩(wěn)定性。通過上述優(yōu)選實施例，本專業(yè)相關(guān)技術(shù)人員應該了解，103中的算法二是當待求位置當前時間步的電磁場值滿足102中判斷條件時方可采用的一種近似算法?；诒竟_的基本思想，本專業(yè)相關(guān)技術(shù)人員可以對判斷條件和相應的算法二進行修改或替換，而不脫離本公開的精神范圍。上述的兩種判斷條件和對應的算法二為優(yōu)選實施例，并不用以限制本公開。本實施例提供的上述方法，利用在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的性質(zhì)，由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測，從而判定可否采用近似算法。在實現(xiàn)計算量少，判斷準確的效果同時，保證了精度與算法穩(wěn)定性。實施例3參見圖301，為一塊三維計算區(qū)域及網(wǎng)格劃分，本實施例以此結(jié)構(gòu)作為一具體實施例，該結(jié)構(gòu)并不用以限制本公開。本實施例提供了一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算方法，流程如圖1，包括：101將待求網(wǎng)格相鄰邊界上的電磁場值組合，形成一個與待求網(wǎng)格形狀相同，位置相同的鏡像網(wǎng)格；102將每個待求網(wǎng)格和相應鏡像網(wǎng)格前一時間步的電磁場值與判斷條件比較，如果滿足判斷條件，則該網(wǎng)格屬于區(qū)域二，否則屬于區(qū)域一；103對區(qū)域一中的所有網(wǎng)格采用算法一，對區(qū)域二中的所有網(wǎng)格采用算法二；所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其中顯式迭代包括但不限于二階蛙跳、四階龍格庫塔等迭代方式，本公開實施例采用四階龍格庫塔迭代方式作為一優(yōu)選實施例，并不用以限制本公開。在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的計算方法有多種，本公開實施例以四面體網(wǎng)格的節(jié)點基dgtd為例，并不用以限制本公開，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。其中對每個網(wǎng)格的處理可以采用逐一處理的方式，也可采用并行處理的方式，本公開實施例采用并行處理的方式，并不用以限制本公開。其中103中的算法一是當前時間歩、待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，包括但不限于dgtd等數(shù)值計算方法，本公開實施例采用dgtd方法，并不用以限制本公開，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。其中103中的算法二是當待求網(wǎng)格當前時間步的電磁場值滿足102中判斷條件時方可采用的一種近似算法。其中算法二與判斷條件對應。本實施例采用：102中的判斷條件為網(wǎng)格單元最大電磁場值小于一判據(jù)，即認為網(wǎng)格單元內(nèi)的電磁場值基本上為零。則對應的103中的算法二為無操作或置零，本實施例采用置零，并不用以限制本公開。以上所述的經(jīng)驗判據(jù)可由經(jīng)驗取得，本公開實施例采用1e-8為一優(yōu)選實施例，也可采用該位置歷史最大值的1e-3，并不用以限制本公開。上述的判斷條件和對應的算法二為一優(yōu)選實施例，并不用以限制本公開。圖302為t＝tn-1時刻的由電磁場分布經(jīng)過101預處理，將待求網(wǎng)格與鏡像網(wǎng)格電磁場值進行102判斷比較后不都滿足判斷條件的網(wǎng)格單元的示意圖，這些不都滿足判斷條件的網(wǎng)格單元即為區(qū)域一。通過上述方法預測在t＝tn時刻，電磁場只主要存在于圖302所示的網(wǎng)格單元。故在t＝tn時刻的計算區(qū)域為圖302中的區(qū)域一，而將其它網(wǎng)格單元的電磁場值簡單置零。相比圖301的整個計算區(qū)域，經(jīng)過本實施例方法的圖302計算區(qū)域減少。本實施例提供的上述方法，利用在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的性質(zhì)，由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測，從而判定可否采用近似算法。在實現(xiàn)計算量少、判定準確的效果同時，保證了精度與算法穩(wěn)定性。實施例4參見圖4，為一多導孔平行板結(jié)構(gòu)，本公開實施例以此結(jié)構(gòu)為一實施例，該結(jié)構(gòu)并不用作限制本公開。本實施例提供了一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算方法，流程如圖1，包括：101將待求網(wǎng)格相鄰邊界上的電磁場值組合，形成一個與待求網(wǎng)格形狀相同，位置相同的鏡像網(wǎng)格；102將每個待求網(wǎng)格和相應鏡像網(wǎng)格前一時間步的電磁場值與判斷條件比較，如果滿足判斷條件，則該網(wǎng)格屬于區(qū)域二，否則屬于區(qū)域一；103對區(qū)域一中的所有網(wǎng)格采用算法一，對區(qū)域二中的所有網(wǎng)格采用算法二；所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其中顯式迭代包括但不限于二階蛙跳、四階龍格庫塔等迭代方式，本公開實施例采用四階龍格庫塔迭代方式作為一優(yōu)選實施例，并不用以限制本公開。其中對每個網(wǎng)格的處理可以采用逐一處理的方式，也可采用并行處理的方式，本公開實施例采用并行處理的方式，并不用以限制本公開。所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其中103中的算法一采用棱柱網(wǎng)格，算法二采用多邊形網(wǎng)格。根據(jù)待求結(jié)構(gòu)生成三維棱柱網(wǎng)格及二維多邊形網(wǎng)格；其中多邊形網(wǎng)格是棱柱網(wǎng)格的直截面，與棱柱網(wǎng)格一一對應；每層棱柱網(wǎng)格對應一層多邊形網(wǎng)格。本公開實施例提供一種基于三棱柱和對應的三角形網(wǎng)格的一階棱邊基函數(shù)作為較佳實施例，并不用以限制本公開，參見圖5：其中l(wèi)i(i＝1:6)為定義在【1】中的第i個頂點的節(jié)點基函數(shù)，ζ沿側(cè)棱從1線性變化到0。圖4中，501為三棱柱網(wǎng)格及其一階棱邊基函數(shù)的位置，502為對應的下底面三角形網(wǎng)格及其一階棱邊基函數(shù)位置。基函數(shù)ni前的系數(shù)ei或hi分別表示該基函數(shù)用于表征電場或磁場。由于在區(qū)域二中，電場主要只存在于與側(cè)棱平行的方向，磁場主要只存在于與側(cè)棱正交的方向，所以n7，n8，n9的磁場系數(shù)h7、h8、h9為零，n1，n2，n3的電場系數(shù)e1、e2、e3為零，由此得到503，簡化的三角形一階矢量基函數(shù)。本公開實施例還提供了一種基于長方體網(wǎng)格和對應下底面矩形網(wǎng)格的一階棱邊基函數(shù)作為又一較佳實施例，并不用以限制本公開，參見圖6：其中分別為長方體網(wǎng)格在x、y、z方向的長度，分別為長方體網(wǎng)格中心的x、y、z方向坐標。圖6中，601為長方體網(wǎng)格及其一階棱邊基函數(shù)的位置，602為對應的下底面矩形網(wǎng)格及其一階棱邊基函數(shù)位置。基函數(shù)ni前的系數(shù)ei或hi分別表示該基函數(shù)用于表征電場或磁場。由于在區(qū)域二中，電場主要只存在于與側(cè)棱平行的方向，磁場主要只存在于與側(cè)棱正交的方向，所以n9、n10、n11、n12的磁場系數(shù)h9、h10、h11、h12為零，n1、n2、n3、n4的電場系數(shù)e1、e2、e3、e4為零，由此得到603，簡化的矩形一階矢量基函數(shù)。其中103中的算法一是當前時間歩、待求網(wǎng)格單元有任意電磁場值都可采用的一種普適算法，包括但不限于dgtd等數(shù)值計算方法，本公開實施例采用三棱柱dgtd方法，并不用以限制本公開，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。其中102判斷條件為網(wǎng)格單元中的電場只主要存在于一個方向，稱作主要電場方向，且沿該方向電場強值基本不變，磁場只主要存在于與主要電場方向正交的平面上。具體的，本公開實施例給出兩種優(yōu)選實施例，并不用以限制本公開：其中為波阻抗，α為依經(jīng)驗設(shè)置的數(shù)值零，本公開實施例將α設(shè)為1e-8作為一種較佳實施例，并不用以限制本公開?；颍浩渲笑聻榻?jīng)驗值，本公開實施例取β＝1％，并不用以限制本公開。其中103中的的算法一為遍歷待求網(wǎng)格單元的相鄰網(wǎng)格單元，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元都屬于區(qū)域一，則對待求網(wǎng)格單元直接采用三維算法求解，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元有網(wǎng)格單元屬于區(qū)域二，則將該區(qū)域二網(wǎng)格單元前一時間步的二維網(wǎng)格電磁場值轉(zhuǎn)化為三維網(wǎng)格電磁場值之后再對待求網(wǎng)格單元采用三維算法求解；其中將二維網(wǎng)格電磁場值轉(zhuǎn)化為三維網(wǎng)格電磁場值采用如下方式：對三棱柱和對應三角形網(wǎng)格：e3d1:6＝0；e3d7:9＝e2d7:9；h3d1:3和h3d4:6＝h2d1:3；h3d7:9＝0；對長方體和對應矩形網(wǎng)格：e3d1:8＝0；e3d9:12＝e2d9:12；h3d1:4和h3d5:8＝h2d1:4；h3d9:12＝0；其中ei和hi和分別為第i個電場和磁場基函數(shù)前的系數(shù)，上角標3d、2d分別表示該系數(shù)用于三維網(wǎng)格和二維網(wǎng)格。其中三維算法，本實施例采用時域不連續(xù)迦略金方式，并不用以限制本公開，其統(tǒng)一特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。時域不連續(xù)迦略金方式可以采用包括但不限于迎風流、中心流等各種數(shù)值流方式來實現(xiàn)。本公開實施例采用棱邊基函數(shù)、迎風流作為較優(yōu)實施例，并不用以限制本公開。本公開實施例提供了一種對棱邊基函數(shù)、迎風流的三維算法：第i個棱柱網(wǎng)格單元求解所需的矩陣方程為：其中方程右端的分別表示與i相鄰的所有網(wǎng)格單元j，在前一時間步的電場、磁場；可以注意到，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。方程中的各個子矩陣表達式為：其中z和y分別為特性阻抗和特性導納，為電場基函數(shù)，為磁場基函數(shù)。將前述的棱邊基函數(shù)代入上述矩陣方程并求解，即完成算法一。其中103中的算法二為：遍歷待求網(wǎng)格單元的相鄰網(wǎng)格單元，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元都屬于區(qū)域二，則對待求網(wǎng)格單元直接采用二維算法求解，如果待求網(wǎng)格單元的相鄰網(wǎng)格單元有網(wǎng)格單元屬于區(qū)域一，則將該區(qū)域一網(wǎng)格單元前一時間步的三維網(wǎng)格電磁場值轉(zhuǎn)化為二維網(wǎng)格電磁場值之后再對待求網(wǎng)格單元采用二維算法求解；其中將三維網(wǎng)格電磁場值轉(zhuǎn)化為二維網(wǎng)格電磁場值采用如下方式：對三棱柱和對應三角形網(wǎng)格：e2d7:9＝e3d7:9；h2d1:3＝h3d4:6或h3d1:3；對長方體和對應矩形網(wǎng)格：e2d9:12＝e3d9:12；h2d1:4＝h3d5:8或h3d1:4；其中ei和hi和分別為第i個電場和磁場基函數(shù)前的系數(shù)，上角標3d、2d分別表示該系數(shù)用于三維網(wǎng)格和二維網(wǎng)格。其中二維算法，本實施例采用時域不連續(xù)迦略金方式，并不用以限制本公開，其統(tǒng)一特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。時域不連續(xù)迦略金方式可以采用包括但不限于迎風流、中心流等各種數(shù)值流方式來實現(xiàn)。本公開實施例采用棱邊基函數(shù)、迎風流作為較優(yōu)實施例，并不用以限制本公開。本公開實施例提供了一種對棱邊基函數(shù)、迎風流的二維算法：第i個棱柱網(wǎng)格單元求解所需的矩陣方程為：其中方程右端的ej、hj分別表示與i相鄰的所有網(wǎng)格單元j，在前一時間步的電場磁場；可以注意到，在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。方程中的子矩陣表達式與本實施例中三維處理方式中的表達式相同。盡管二維處理與三維處理矩陣方程和表達式形式相同，但是對邊界積分分別采用二維曲線積分和三維曲面積分，所以得到矩陣是不相同的。將上述的棱邊基函數(shù)代入上述矩陣方程并求解，即完成二維處理?？蛇x的，本公開實施例提出一種可選步驟，在二維處理中，如果網(wǎng)格中存在由同一多面體拉伸而成的多層棱柱網(wǎng)格單元，且所述各層棱柱網(wǎng)格單元都屬于區(qū)域二，則將該多層網(wǎng)格簡化為單層網(wǎng)格，對該同一多面體網(wǎng)格單元采用二維處理。具體的，參見圖7，例如有兩層三棱柱網(wǎng)格單元由同一三角形網(wǎng)格單元拉伸而成，且都屬于區(qū)域二，如701，則將其簡化為一單層網(wǎng)格處理，如702。上述二維網(wǎng)格電磁場值轉(zhuǎn)化為三維網(wǎng)格電磁場值步驟改為：對由同一三角形拉伸而成的多層三棱柱網(wǎng)格中的每個三棱柱網(wǎng)格，e3d1:6＝0；e3d7:9＝e2d7:9；h3d1:3和h3d4:6＝h2d1:3；h3d7:9＝0；三維網(wǎng)格電磁場值轉(zhuǎn)化為二維網(wǎng)格電磁場值步驟改為：從同一三角形拉伸而成的多層三棱柱網(wǎng)格中任選一個三棱柱網(wǎng)格，使得e2d7:9＝e3d7:9；h2d1:3＝h3d4:6或h3d1:3；通過此可選步驟，可以進一步減少待求未知量。上述二維算法、三維算法盡管只采用了時域不連續(xù)迦略金方式為例進行說明，但本領(lǐng)域?qū)I(yè)技術(shù)人員應該注意到，在分別進行二維處理和三維處理，并使區(qū)域二與區(qū)域一聯(lián)接起來的過程中，本公開實施例只利用了時域不連續(xù)迦略金方法的以下性質(zhì)，即：其中三維算法采用棱柱網(wǎng)格，二維算法采用多邊形網(wǎng)格。三維算法和二維算法在任一時間步任一網(wǎng)格單元的待求電磁場值只由前一時間步本網(wǎng)格單元和相鄰邊界單元的電磁場值決定。所以根據(jù)本實施例足以合理推斷，具有上述性質(zhì)的包括但不限于時域不連續(xù)迦略金方式的其它算法同樣適用于本公開所提出的方法。本公開實施例采用時域不連續(xù)迦略金方式作為一種較佳實施例，并不用以限制本公開。本實施例提供的上述方法，通過將三維結(jié)構(gòu)分為區(qū)域一和區(qū)域二分別處理計算，在區(qū)域二利用了電磁場分布的物理特性以減少未知量。本實施例提供的上述方法，利用在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的性質(zhì)，由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測，從而判定可否采用近似算法。在實現(xiàn)計算量少、判斷準確的效果同時，保證了精度與算法穩(wěn)定性。實施例6參見圖8，本實施例提供了一種在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的顯式迭代的時域有限電磁場計算裝置，包括：801預處理模塊：將待求網(wǎng)格相鄰邊界上的電磁場值組合，形成一個與待求網(wǎng)格形狀相同，位置相同的鏡像網(wǎng)格；802判斷模塊：將每個待求網(wǎng)格和相應鏡像網(wǎng)格前一時間步的電磁場值與判斷條件比較，如果都滿足判斷條件，則該網(wǎng)格屬于區(qū)域二，否則屬于區(qū)域一；803計算模塊：對區(qū)域一中的所有網(wǎng)格采用算法一，對區(qū)域二中的所有網(wǎng)格采用算法二；所述算法一是待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，算法二是當待求網(wǎng)格當前時間步的電磁場值滿足所述判斷條件時方可采用的一種近似算法。其中顯式迭代包括但不限于二階蛙跳、四階龍格庫塔等迭代方式，本公開實施例對此不做具體限定。其中在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的方法有多種，包括但不限于dgtd方法，本公開實施例對此不做具體限定。其中801待求網(wǎng)格和鏡像網(wǎng)格可以是任意形狀，基函數(shù)可以是各種形式，本實施例分別以四面體網(wǎng)格、三棱柱網(wǎng)格，棱邊即函數(shù)為例說明，并不用以限制本發(fā)明。如圖201和203所示待求網(wǎng)格i，其電磁場值為ei，與其相鄰的其它網(wǎng)格的邊界上的電磁場值為e′,e″。將這些邊界上的電磁場值組合起來形成鏡像網(wǎng)格ei*。鏡像網(wǎng)格ei*的組合方式有多種，如圖202和204，本公開實施例給出兩種優(yōu)選例，如ei*＝e′+e″或ei*＝(e′+e″)/2，并不用以限制本發(fā)明。其中803中的算法一是當前時間歩、待求網(wǎng)格有任意電磁場值都可采用的一種普適算法，包括但不限于dgtd等數(shù)值計算方法，本公開實施例對此不做具體限定，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。其中803中的算法二是當待求網(wǎng)格當前時間步的電磁場值滿足102中判斷條件時方可采用的一種近似算法。算法二與判斷條件對應，例如：例一：802中判斷條件為待求網(wǎng)格最大電磁場值小于一判據(jù)，即認為待求網(wǎng)格中的電磁場值基本上為零。則對應的803中算法二為無操作或電磁場值置零；例二：802中判斷條件為網(wǎng)格中的電場只主要存在于一個方向，稱作主要電場方向，且沿該方向電場強值基本不變，磁場只主要存在于與主要電場方向正交的平面上，具體地，電場在兩個正交方向的最大值小于一判據(jù)，磁場在另一個正交方向的最大值小于一判據(jù)。則認為網(wǎng)格單元中的電磁場值可以用二維網(wǎng)格準確擬合，對應的803中的算法二為二維算法。其中二維算法包括但不限于dgtd等數(shù)值計算方法，本公開實施例對此不做具體限定，統(tǒng)一的特征是：在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響。相反，如果待求網(wǎng)格當前時間歩的電磁場值不滿足802中的判斷條件，例如：待求網(wǎng)格中有電磁場值大于經(jīng)驗判據(jù)，即不能認為待求網(wǎng)格電磁場值基本上為零，此時采用算法二，無操作或電磁場值置零，將帶來不可控的誤差，從而影響精度和穩(wěn)定性；待求網(wǎng)格中電磁場值不滿足：電場只主要存在于一個方向，稱作主要電場方向，且沿該方向電場強值基本不變，磁場只主要存在于與主要電場方向正交的平面上。此時待求網(wǎng)格中電磁場值不能用二維網(wǎng)格準確擬合，采用算法二，二維算法，將帶來不可控誤差，從而影響精度和穩(wěn)定性。通過上述優(yōu)選實施例，本專業(yè)相關(guān)技術(shù)人員應該了解，803中的算法二是當待求位置當前時間步的電磁場值滿足802中判斷條件時方可采用的一種近似算法?；诒竟_的基本思想，本專業(yè)相關(guān)技術(shù)人員可以對判斷條件和相應的算法二進行修改或替換，而不脫離本公開的精神范圍。上述的兩種判斷條件和對應的算法二為優(yōu)選實施例，并不用以限制本公開。本實施例提供的上述裝置，利用在前一時間歩，只有待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值會對當前時間歩的待求網(wǎng)格電磁場值產(chǎn)生影響的性質(zhì)，由前一時間歩待求網(wǎng)格和待求網(wǎng)格相鄰邊界的電磁場值，對當前時間歩待求網(wǎng)格的電磁場值做出預測，從而判定可否采用近似算法。在實現(xiàn)適用廣泛、計算量少的效果同時，保證了一定的精度與算法穩(wěn)定性。本領(lǐng)域的技術(shù)人員應明白，本公開的實施例可提供為方法、裝置(設(shè)備)、或計算機程序產(chǎn)品。因此，本公開可采用完全硬件實施例、完全軟件實施例、或結(jié)合軟件和硬件方面的實施例的形式。而且，本公開可采用在一個或多個其中包含有計算機可用程序代碼的計算機可用存儲介質(zhì)(包括但不限于磁盤存儲器、制度光盤、光學存儲器等)上實施的計算機程序產(chǎn)品的形式。本公開是參照根據(jù)本公開實施例的方法、裝置(設(shè)備)和計算機程序指令實現(xiàn)流程圖和/或方框圖來描述的。應理解可由計算機程序指令實現(xiàn)流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合?？商峁┻@些計算機程序指令到通用計算機、專用計算機、嵌入式處理機或其它可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn)生一個機器，使得通過計算機或其它可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實現(xiàn)在流程圖的一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的裝置。這些計算機程序指令也可存儲在能引導計算機或其它可編程數(shù)據(jù)處理設(shè)備以特定方式工作的計算機可讀存儲器中，使得存儲在該計算機可讀存儲器中的指令產(chǎn)生包括指令裝置的制造品，該指令裝置實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能。這些計算機程序指令也可裝載到計算機或其它可編程數(shù)據(jù)處理設(shè)備上，使得在計算機或其它可編程設(shè)備上執(zhí)行一系列操作步驟以產(chǎn)生計算機實現(xiàn)的處理，從而在計算機或其他可編程設(shè)備上執(zhí)行的指令提供用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的步驟。盡管已描述了本公開的優(yōu)選實施例，但本領(lǐng)域內(nèi)的技術(shù)人員一旦得知了基本創(chuàng)造性概念，則可對這些實施例做出另外的變更和修改。所以，所附權(quán)利要求意欲解釋包括優(yōu)選實施例以及落入本公開范圍的所有變更和修改。顯然，本領(lǐng)域的技術(shù)人員可以對本公開進行各種改動和變型而不脫離本公開的精神范圍。這樣，倘若本公開的這些修改和變型屬于本公開權(quán)利要求及其等同技術(shù)的范圍之內(nèi)，本公開也意圖包含這些改動和變型在內(nèi)。當前第1頁12當前第1頁12