发明名称 一种基于近似概率转换的电路性能可靠性的估计方法
摘要 本发明涉及一种基于近似概率转换的电路性能可靠性的估计方法,属于电子线路设计技术领域。对于给定的电路性能函数,将其输入变量表示为区间证据;将该区间证据通过性能函数映射到输出,获得输出量的区间证据;利用可传递信度模型中的近似概率转换,将输出变量的区间证据转化为近似累积概率分布;将该近似分布作为输出变量真实累积概率分布的近似,用其估计电路性能可靠度是否达到要求。本发明方法在区间采样次数确定的情况下,只需实施一次仿真过程即可得到置信水平为100%的确定性估计误差,在同样的估计误差下,本方法所需计算量远远小于已有的蒙特卡罗方法。
申请公布号 CN101980220B 申请公布日期 2012.06.20
申请号 CN201010515429.6 申请日期 2010.10.15
申请人 清华大学 发明人 周东华;徐晓滨;吉吟东;孙新亚
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 1.一种基于近似概率转换的电路性能可靠性的估计方法,其特征在于该方法包括以下各步骤:(1)设电路系统的系统性能函数为:y=g(u),y=g(u)需要满足以下条件:y=g(u)是连续的并且其偏导数也是连续的,y=g(u)关于每个变量u<sub>1</sub>,…,u<sub>i</sub>,…u<sub>n</sub>是严格单调的,其中u=(u<sub>1</sub>,…,u<sub>i</sub>,…u<sub>n</sub>),i=1,…n,u<sub>i</sub>表示电路系统中第i个元件的参数,设定u<sub>i</sub>的标称值u<sub>i_0</sub>和公差Δu<sub>i</sub>,则u<sub>i</sub>的容差区间为I<sub>i</sub>=[u<sub>i_0</sub>-Δu<sub>i</sub>,u<sub>i_0</sub>+Δu<sub>i</sub>],取u<sub>i</sub>为在I<sub>i</sub>上满足截断型正态分布的随机变量N(u<sub>i_0</sub>,σ<sub>i</sub>),其中u<sub>i_0</sub>为u<sub>i</sub>的标称值,σ<sub>i</sub>为u<sub>i</sub>的标准差,且各随机变量之间相互独立,则电路系统的性能函数y为一随机变量;给定性能函数y的容差区间D=[y<sub>A</sub>,y<sub>B</sub>],则y∈D的概率P<sub>y</sub>(D)为:P<sub>y</sub>(D)=F<sub>y</sub>(y<sub>B</sub>)-F<sub>y</sub>(y<sub>A</sub>)其中,F<sub>y</sub>(y)为y的累积概率分布函数,F<sub>y</sub>(y<sub>A</sub>)和F<sub>y</sub>(y<sub>B</sub>)分别为y在y<sub>A</sub>和y<sub>B</sub>点的累积概率分布函数值;设定电路系统的性能可靠度阈值为Υ,则电路性能可靠性评估准则为:若P<sub>y</sub>(D)>Υ,则判断电路系统性能为可靠;若P<sub>y</sub>(D)<Υ,则判断电路系统性能为不可靠;若P<sub>y</sub>(D)=Υ,则判断电路系统性能为临界状态;(2)将上述步骤(1)中性能函数y的输入u表示为区间证据<img file="FSB00000714726600011.GIF" wi="155" he="47" />其过程为:将步骤(1)中各输入变量的容差区间I<sub>i</sub>表示为笛卡尔积U=I<sub>1</sub>×…×I<sub>n</sub>,则u=(u<sub>1</sub>,…,u<sub>i</sub>,…,u<sub>n</sub>)∈U,将I<sub>i</sub>划分为d<sub>i</sub>个互不相交的子区间,则区间证据<img file="FSB00000714726600012.GIF" wi="126" he="47" />中的<img file="FSB00000714726600013.GIF" wi="43" he="40" />是U的一个划分,表示为<img file="FSB00000714726600014.GIF" wi="454" he="68" />其中<maths num="0001"><![CDATA[<math><mrow><msub><mi>A</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>C</mi><mn>1</mn><mi>k</mi></msubsup><mo>&times;</mo><msubsup><mi>C</mi><mn>2</mn><mi>k</mi></msubsup><mo>&times;</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>&times;</mo><msubsup><mi>C</mi><mi>n</mi><mi>k</mi></msubsup><mo>,</mo><msubsup><mi>C</mi><mi>i</mi><mi>k</mi></msubsup><mo>&SubsetEqual;</mo><msub><mi>I</mi><mi>i</mi></msub><mo>,</mo></mrow></math>]]></maths><img file="FSB00000714726600016.GIF" wi="50" he="57" />表示对I<sub>i</sub>进行划分后得到的一个子区间,<img file="FSB00000714726600017.GIF" wi="59" he="66" />表示<img file="FSB00000714726600018.GIF" wi="43" he="40" />的势,即<img file="FSB00000714726600019.GIF" wi="44" he="40" />中元素A<sub>k</sub>的个数,<img file="FSB000007147266000110.GIF" wi="318" he="65" /><img file="FSB000007147266000111.GIF" wi="255" he="65" />区间证据<img file="FSB000007147266000112.GIF" wi="131" he="47" />中的m表示A<sub>k</sub>的概率赋值函数,当输入u中的u<sub>1</sub>,…,u<sub>n</sub>是相互独立的随机变量时,m由下式给出:<maths num="0002"><![CDATA[<math><mrow><mi>m</mi><mo>=</mo><mi>m</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mo>&Integral;</mo><msub><mi>A</mi><mi>k</mi></msub></msub><mi>f</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mi>du</mi><mo>=</mo><msub><mo>&Integral;</mo><mrow><msubsup><mi>C</mi><mn>1</mn><mi>k</mi></msubsup><mo>&times;</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>&times;</mo><msubsup><mi>C</mi><mi>i</mi><mi>k</mi></msubsup><mo>.</mo><mo>.</mo><mo>.</mo><mo>&times;</mo><msubsup><mi>C</mi><mi>n</mi><mi>k</mi></msubsup></mrow></msub><mi>f</mi><mrow><mo>(</mo><msub><mi>u</mi><mn>1</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>u</mi><mi>n</mi></msub><mo>)</mo></mrow><msub><mi>du</mi><mn>1</mn></msub><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>du</mi><mi>i</mi></msub><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>du</mi><mi>n</mi></msub></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mo>=</mo><mrow><mo>(</mo><msub><mo>&Integral;</mo><msubsup><mi>C</mi><mn>1</mn><mi>k</mi></msubsup></msub><mi>f</mi><mrow><mo>(</mo><msub><mi>u</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>du</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><msub><mo>&Integral;</mo><msubsup><mi>C</mi><mi>i</mi><mi>k</mi></msubsup></msub><mi>f</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>du</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><msub><mo>&Integral;</mo><msubsup><mi>C</mi><mi>n</mi><mi>k</mi></msubsup></msub><mi>f</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>n</mi></msub><mo>)</mo></mrow><msub><mi>du</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>=</mo><msub><mi>m</mi><mn>1</mn></msub><mrow><mo>(</mo><msubsup><mi>C</mi><mn>1</mn><mi>k</mi></msubsup><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msubsup><mi>C</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>m</mi><mi>n</mi></msub><mrow><mo>(</mo><msubsup><mi>C</mi><mi>n</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>其中<maths num="0005"><![CDATA[<math><mrow><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msubsup><mi>C</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mo>&Integral;</mo><msubsup><mi>C</mi><mi>i</mi><mi>k</mi></msubsup></msub><mi>f</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>du</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>为输入u<sub>i</sub>落入<img file="FSB00000714726600025.GIF" wi="49" he="57" />的概率,且<img file="FSB00000714726600026.GIF" wi="267" he="132" />(3)将上述步骤(2)得到的区间证据<img file="FSB00000714726600027.GIF" wi="130" he="47" />输入到性能函数中,获得输出量y的区间证据<img file="FSB00000714726600028.GIF" wi="155" he="46" />过程如下:首先获得区间证据<img file="FSB00000714726600029.GIF" wi="127" he="46" />中的<img file="FSB000007147266000210.GIF" wi="66" he="40" /><img file="FSB000007147266000211.GIF" wi="533" he="54" />其中g(A<sub>k</sub>)={g(u)|u∈A<sub>k</sub>};获得区间证据<img file="FSB000007147266000212.GIF" wi="128" he="48" />中的ρ:ρ(R<sub>j</sub>)=∑{m(A<sub>k</sub>)|R<sub>j</sub>=g(A<sub>k</sub>)},其中<img file="FSB000007147266000213.GIF" wi="306" he="65" /><img file="FSB000007147266000214.GIF" wi="62" he="66" />为<img file="FSB000007147266000215.GIF" wi="46" he="40" />中元素R<sub>j</sub>的个数,且<img file="FSB000007147266000216.GIF" wi="268" he="136" />(4)将上述步骤(3)中求得的R<sub>j</sub>表示为区间R<sub>j</sub>=[a<sub>j</sub>,b<sub>j</sub>],则R<sub>j</sub>的概率赋值函数表示为ρ(R<sub>j</sub>)=ρ([a<sub>j</sub>,b<sub>j</sub>]),通过可传递信度模型中的近似概率转换,将<img file="FSB000007147266000217.GIF" wi="126" he="47" />转换为输出变量y的近似累积概率分布为<img file="FSB000007147266000218.GIF" wi="524" he="138" />上式中的求和项BetF(y)<sub>j</sub>为:当y≥b<sub>j</sub>,则BetF(y)<sub>j</sub>=ρ([a<sub>j</sub>,b<sub>j</sub>]),当a<sub>j</sub><y<b<sub>j</sub>,则<img file="FSB000007147266000219.GIF" wi="611" he="122" />当y≤a<sub>j</sub>,则BetF(y)<sub>j</sub>=0;(5)根据上述步骤(1)的电路性能可靠性评估准则,用上述BetF(y)代替电路性能y的累积概率分布F<sub>y</sub>(y),得到电路性能可靠度的估计值<img file="FSB000007147266000220.GIF" wi="98" he="68" /><maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>&gamma;</mi><mo>^</mo></mover><mi>Bet</mi></msub><mo>=</mo><mi>BetF</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>BetF</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>A</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>则<img file="FSB000007147266000222.GIF" wi="76" he="65" />即为y∈D的近似累积概率;(6)将上述<img file="FSB000007147266000223.GIF" wi="76" he="65" />与设定的性能可靠度阈值Υ比较,若<img file="FSB000007147266000224.GIF" wi="184" he="65" />则电路系统性能为可靠;若<img file="FSB000007147266000225.GIF" wi="186" he="65" />则电路系统性能为不可靠;若<img file="FSB000007147266000226.GIF" wi="181" he="65" />则电路系统性能为临界状态。
地址 100084 北京市海淀区清华园1号