[发明专利]一种基于改进自适应重要抽样的性能可靠性仿真方法

摘要

一种基于改进自适应重要抽样的性能可靠性仿真方法，该方法依照下述四个阶段进行：第一阶段：抽取初始失效样本，以参数原始均值为采样中心，原始方差的1-3倍为采样协方差，进行随机搜索，并通过迭代找到初始失效样本；第二阶段：抽取一批失效样本，以第一阶段搜索到的失效样本作为第二阶段的初始采样中心进行随机抽样，并在下一次的循环中，调整采样中心，循环至抽取到指定数目的失效样本；第三阶段：自适应重要抽样，根据第二阶段抽取到的失效样本，计算第三阶段的初始采样中心和初始采样协方差，进行循环抽样仿真；第四阶段：统计计算，完成一次抽样后，重新计算采样中心和采样协方差，当失效概率的仿真结果趋于稳定时，停止仿真，完成自适应重要抽样的性能可靠性仿真。

主权项

一种基于改进自适应重要抽样的性能可靠性仿真方法，其特征在于：该方法依照下述四个阶段进行：第一阶段：抽取初始失效样本初始失效样本的选择对于算法的收敛有关键作用，应使其尽可能地靠近失效面上的设计点；以原始概率密度函数的均值μ为采样中心，原方差σ的1‑3倍作为采样协方差，进行随机搜索得随机参数向量样本Xi＝[x1ix2i...xLi]，仿真，得性能结果θi；给定误差限ζ，判断θi和失效阈值θf的距离是否在给定误差限内，即|θi‑θf|＜ζ 若是则判断系统是否失效，是则进入第二个阶段，否则继续循环；在下一次的循环中，对采样中心μ进行调整，μi＝argmin(|θ(μi‑1)‑θf|，|θ(Xi‑1)‑θf|)，即将性能结果|θ(μi‑1)‑θf|和|θ(Xi‑1)‑θf|进行比较；若|θ(μi‑1)‑θf|＞|θ(Xi‑1)‑θf|，说明抽样值Xi‑1比采样中心μi‑1更靠近失效面，则将Xi‑1作为下次抽样的中心；否则仍然以μi‑1作为采样中心；第二阶段：抽取一批失效样本抽取完初始失效样本后，进入第二个阶段；以第一阶段搜索到的失效样本作为第二阶段的初始采样中心，仿真并计算其概率密度f(Xj)；判断系统是否失效，是则令失效样本个数M＝M+1；判断M是否达到指定失效样本个数Nr，是则进入第三阶段，否则继续循环；在下一次的循环中，对采样中心μ进行调整，μi+1＝argmax(f(Xi‑1)，f(μi‑1))，选择概率密度值大者作为新的采样中心，循环至抽取到Nr个失效样本，进入第三阶段；第三阶段：自适应重要抽样根据第二阶段抽取到的Nr个失效样本，分别利用下列式(1)和式(2)计算第三阶段的初始采样中心μf(1)和初始采样协方差σf2(1)，即重要抽样函数h(.)的初始参数；令循环次数L＝1，每次仿真Nf次，失效概率在指定误差内的循环次数S＝0； μ ^ f = 1 N s Σ i = 1 N r x i f - - - ( 1 ) σ ^ f 2 = 1 N r diag { σ ^ 1 , f 2 , σ ^ 2 , f 2 , . . . σ ^ n , f 2 } - - - ( 2 ) 重要抽样函数h(.)为：h(X)＝I(X)f(X)/Pf    (3)其中Pf为系统性能可靠度，I(X)为示性函数，即对于失效阈值θf，若性能结果θ大于(或小于)θf时，系统为失效，此时I(X)为1，否则为0，即 I ( X ) = 1 , θ ≥ θ f 0 , θ < θ f - - - ( 4 ) 采用蒙特卡洛数值求解系统性能可靠度，通常是利用概率密度函数f(X)抽取足够多的样本值，然后用示性函数I(X)的均值来近似估算失效概率Pf，即： P f ≈ P ^ f = 1 N Σ i = 1 N I ( X i ) - - - ( 5 ) 第四阶段：统计计算完成一次抽样后，根据下列式(6)‑(7)计算采样中心μf(L)和采样协方差σf2(L)，并作为下次抽样的抽样中心和协方差；根据下列式(8)计算失效概率Pf(L)，并预先给定容许误差η，计算第L次仿真所得失效率Pf(L)和第L‑S‑1次仿真所得失效率Pf(L‑S‑1)的相对误差δ＝|Pf(L)‑Pf(L‑S‑1)|/Pf(L‑S‑1)；如果δ＞η，说明失效概率Pf(L)和第L‑s‑1次循环所得失效概率Pf(L‑s‑1)的相对误差不在给定容许误差内，则令失效概率在指定误差内的循环次数s＝0，并进入下一次循环；否则，令S＝S+1，进入下一次循环；如果S＝10，即连续10次循环所得失效概率都在容许误差内，说明仿真结果已经趋于稳定，停止仿真，完成自适应重要抽样的性能可靠性仿真；所述公式(6)、(7)、(8)列如下： μ f ( l ) = 1 l ( Σ i l - 1 μ f ( i ) + μ f ( l ) ) = l - 1 l μ f ( l - 1 ) + 1 l μ f ( l ) - - - ( 6 ) σ f 2 ( l ) = 1 l ( Σ i = 1 l - 1 σ f 2 ( i ) + σ f 2 ( l ) ) = l - 1 l σ f 2 ( l - 1 ) + 1 l σ f 2 ( l ) - - - ( 7 ) P f ( l ) = 1 l ( Σ i = 1 l - 1 P f ( i ) + P f ( l ) ) = l - 1 l P f ( l - 1 ) + 1 l P f ( l ) - - - ( 8 ) .