主权项 |
一种基于伯努利分布的贝叶斯模型混合预测电路成品率的方法,其特征在于,其包括步骤:步骤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>β</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>|</mo><mi>β</mi><mo>)</mo></mrow><mi>p</mi><mrow><mo>(</mo><mi>β</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mi>Γ</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>β</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><msup><mi>β</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>β</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>β</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>β</mi><mi>E</mi></msub></mrow></msup></mrow><mrow><mi>Γ</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mn>1</mn><mrow><mo>1</mo><mo>-</mo><msub><mi>β</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>β</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>∫</mo><mn>0</mn><mn>1</mn></munderover><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>β</mi><mo>|</mo><mi>a</mi><mo>)</mo></mrow><mo>·</mo><mi>dβ</mi></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mi>Γ</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msub><mi>β</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mi>a</mi><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>β</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>β</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>Γ</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>β</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>β</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>Γ</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>β</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>Γ</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>β</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mi>a</mi><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>β</mi><mi>E</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>β</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>β</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>Γ</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>Γ</mi><mrow><mo>(</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>β</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>β</mi><mi>E</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>Γ</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>β</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>β</mi><mo>|</mo><mi>x</mi><mo>)</mo></mrow><mo>∝</mo><msup><mi>β</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>β</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>β</mi><mi>E</mi></msub><mo>)</mo></mrow><mo>·</mo><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mi>β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>β</mi><mo>|</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>Γ</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>β</mi><mi>E</mi></msub><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mo>/</mo><msup><mi>β</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>β</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>β</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>β</mi><mi>E</mi></msub></mrow></msup></mrow></mrow><mrow><mi>Γ</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mi>a</mi><mo>)</mo></mrow><mo></mo><mi>Γ</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>M</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>β</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>β</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>β</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>β</mi><mi>E</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000490395690000023.GIF" wi="589" he="172" /></maths>步骤205:输出估计得到的成品率数值β<sub>MAP</sub>。 |