发明名称 基于伯努利分布的贝叶斯模型混合预测电路成品率方法
摘要 本方法属于集成电路领域,涉及一种基于伯努利分布的贝叶斯模型混合预测电路成品率的方法。该方法通过结合在集成电路设计的不同阶段的信息,加快对只具有“通过—不通过”两种状态的电路的成品率估计过程。该方法为“通过—不通过”的输出结果建立一个伯努利模型,将先验成品率设定为beta分布,并利用最大似然法确定beta分布中的超参数。再使用该超参数,结合比较少量的后验信息,估算出集成电路的成品率。该方法相比传统的蒙特-卡洛方法估计成品率,在达到同一精度的情况下,需要的后验信息少了很多,能明显节省进行后仿真或者进行新一次测试的时间。
申请公布号 CN104978448A 申请公布日期 2015.10.14
申请号 CN201410146481.7 申请日期 2014.04.14
申请人 复旦大学 发明人 曾璇;李昕;杨帆;方晨蕾
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 上海元一成知识产权代理事务所(普通合伙) 31268 代理人 吴桂琴
主权项 一种基于伯努利分布的贝叶斯模型混合预测电路成品率的方法,其特征在于,其包括步骤:步骤201:读取前阶段与后阶段的数据,该些数据经预先处理编码为“0—1”的格式,其中,0代表未通过测试,1代表测试通过;步骤202:将得到的数据分为两组,其中:对于流片前仿真的应用情况,对应的是布局布线前仿真(前阶段)和布局布线后仿真(后阶段);对于流片后测试的应用,对应的是较早一批的测试结果(前阶段)和最新一批的测试结果(后阶段);步骤203:结合后阶段得到的数据,利用最大似然法得到超参数a的值:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&beta;</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>|</mo><mi>&beta;</mi><mo>)</mo></mrow><mi>p</mi><mrow><mo>(</mo><mi>&beta;</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><msup><mi>&beta;</mi><mrow><mi>M</mi><mo>+</mo><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&beta;</mi><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub></mrow></msup></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mn>1</mn><mrow><mo>1</mo><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000490395690000011.GIF" wi="1326" he="288" /></maths>上式中,x代表取得的后阶段数据点,M表征其中通过测试(结果为“1”)的数据点个数,N为总的数据点个数,将上式对β从0到1积分:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mn>1</mn></munderover><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&beta;</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>d&beta;</mi></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mi>a</mi><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000490395690000012.GIF" wi="1436" he="349" /></maths>为了使得上式最大(也即取得观测结果x的可能最大),使用线性搜索的方法获得一个最优的超参数a满足下式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mi>max</mi><mi>a</mi></munder><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mi>a</mi><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000490395690000013.GIF" wi="1413" he="198" /></maths>其中,超参数的大小表征了对于先验信息的确信程度;步骤204:使用最大后验法得到后验成品率根据贝叶斯公式,参数β的后验分布正比于先验分布和似然函数的乘积:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>&beta;</mi><mo>|</mo><mi>x</mi><mo>)</mo></mrow><mo>&Proportional;</mo><msup><mi>&beta;</mi><mrow><mi>M</mi><mo>+</mo><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&beta;</mi><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mi>&beta;E</mi></mrow></msup></mrow>]]></math><img file="FDA0000490395690000021.GIF" wi="1105" he="115" /></maths>经过归一化得到:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>&beta;</mi><mo>|</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mrow><mi>N</mi><mo>+</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mo>/</mo><msup><mi>&beta;</mi><mrow><mi>M</mi><mo>+</mo><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&beta;</mi><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub></mrow></msup></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mi>a</mi><mo>)</mo></mrow><mo></mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000490395690000022.GIF" wi="1459" he="203" /></maths>上式为参数β的后验分布表达式,可以看到,该式也是一个beta分布,其值取得最大时的参数β的值即为BMF方法估计的成品率:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mi>MAP</mi></msub><mo>=</mo><mfrac><mrow><mi>M</mi><mo>+</mo><mi>a</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>+</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>&beta;</mi><mi>E</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000490395690000023.GIF" wi="589" he="172" /></maths>步骤205:输出估计得到的成品率数值β<sub>MAP</sub>。
地址 200433 上海市杨浦区邯郸路220号