技術(shù)特征：

1.一種基于粒子密度不均勻性的迭代蒙卡全局權(quán)窗參數(shù)生成方法，其特征在于包括以下步驟：

(1)生成初始權(quán)窗網(wǎng)格：

(11)建立覆蓋整個模型空間的網(wǎng)格；

(12)將所有網(wǎng)格的權(quán)窗參數(shù)都初始化為0；

(2)在粒子輸運計算過程中，統(tǒng)計進入每個網(wǎng)格的粒子數(shù)量N_i、總粒子權(quán)重W_i和每個網(wǎng)格對粒子密度均勻性的貢獻，以及所有網(wǎng)格的平均粒子數(shù)量Na，i指第i個網(wǎng)格；在每個粒子計算過程中需要：

(21)記錄粒子的徑跡信息，即粒子所穿過的所有網(wǎng)格的編號；

(22)當粒子第一次進入網(wǎng)格i后，如果粒子權(quán)重w大于當前所有網(wǎng)格權(quán)窗參數(shù)的最小值，將進入網(wǎng)格i的粒子數(shù)量N_i增加1，進入網(wǎng)格i的總粒子權(quán)重W_i增加當前粒子的權(quán)重w；

(23)將粒子每一步輸運中，終點所在的網(wǎng)格的粒子數(shù)量N_i和所有網(wǎng)格的平均粒子數(shù)量Na做比較，如果N_i小于Na，則該粒子徑跡上的所有網(wǎng)格的貢獻增加，增加的貢獻即為粒子進入一個網(wǎng)格i對粒子密度均勻性的貢獻；

(3)計算權(quán)窗參數(shù)

(31)計算每個網(wǎng)格的重要性，

計算每個網(wǎng)格的重要性Imp_i，公式中C_i為粒子進入網(wǎng)格i對粒子密度均勻性的總貢獻，W_i為進入網(wǎng)格i的粒子總權(quán)重；

(32)將每個網(wǎng)格的重要性根據(jù)源所在網(wǎng)格的平均重要性進行歸一化；

(33)將每個網(wǎng)格的權(quán)窗參數(shù)設(shè)為歸一化后的網(wǎng)格重要性的倒數(shù)的一半；

(4)迭代優(yōu)化

(41)在權(quán)窗生成過程中，將粒子輸運分為若干個周期，根據(jù)前一個周期的網(wǎng)格重要性生成新的權(quán)窗參數(shù)，并應(yīng)用于當前周期的輸運，同時每個周期生成的權(quán)窗參數(shù)都通過文本進行輸出；

(42)在每個周期粒子輸運開始時，將每個網(wǎng)格的貢獻和進入權(quán)重都乘以一個折扣因子f；

(43)迭代自動進行，在達到了用戶設(shè)定的時間上限后停止計算，當前的權(quán)窗參數(shù)即為最終生成的全局權(quán)窗參數(shù)。

2.根據(jù)權(quán)利要求1所述的基于粒子密度不均勻性的迭代蒙卡全局權(quán)窗參數(shù)生成方法，其特征在于：所述步驟(31)中粒子進入網(wǎng)格i對粒子密度均勻性的貢獻計算方法如下：對粒子密度均勻性的貢獻的初始值為1.0，并根據(jù)網(wǎng)格粒子密度權(quán)重和粒子輸運模擬時間權(quán)重進行加權(quán)，粒子密度權(quán)重F_density定義為：

公式中N_a為所有網(wǎng)格平均粒子數(shù)量，N_i為進入網(wǎng)格i的粒子數(shù)量；

粒子輸運時間則通過粒子當前由于權(quán)窗而分裂的次數(shù)進行近似，粒子輸運時間權(quán)重F_time定義為：

公式中T_s為粒子分裂的次數(shù)；

從而單個粒子進入網(wǎng)格i對粒子密度的均勻性的貢獻ΔC_i定義為：

ΔC_i＝F_density×F_time，所述步驟(31)中的C_i為所有粒子產(chǎn)生的貢獻之和：

C_i＝∑_{all particles}ΔC_i。

3.根據(jù)權(quán)利要求1所述的基于粒子密度不均勻性的迭代蒙卡全局權(quán)窗參數(shù)生成方法，其特征在于：所述粒子進入一個網(wǎng)格i后對粒子密度均勻性的貢獻的計算方法，其特征在于：折扣因子f取值范圍為0～1。