【正文】 m=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 16 G=m1**m3。 K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 17 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 T=n31*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 18 m12=m1+K0*n1*R。m22=m2+K0*S*n2。m32=m31+K0*T*n31。 VE=(m12m1)/n1。VI=(m22m2)/n2。VP=(m32m31)/n31。 第一次迭代結(jié)束 19 第二次迭代 首先將 P的 I型變量轉(zhuǎn)為正態(tài)變量 m1=m12。n1=*10^7。m2=m22。n2=*10^4。m3=m32。n3=1。 c=。d=。 Fm=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 然后求 KI（可靠度） G=m1**m3。 Q=Gm2*n1*VEm1*n2*VI+*n31*VP。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 20 最后計(jì)算新設(shè)計(jì)驗(yàn)算點(diǎn)的均值 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 T=n3*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 m12=2*10^7+K1**10^7*R。m22=10^4+K1*S**10^4。m32=m31+K1*T*n31。 VE=(m12m1)/n1。VI=(m22m2)/n2。VP=(m32m31)/n31。 第 2次迭代結(jié)果為 21 最后根據(jù)可靠度驗(yàn)算條件編程 while abs(K1K0)*K0 m1=m12。n1=*10^7。m2=m22。n2=*10^4。m3=m32。n3=1。 c=。d=。 K0=K1。 Fm=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 G=m1**m3。 Q=Gm2*n1*VEm1*n2*VI+*n31*VP。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 T=n3*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 m12=2*10^7+K1**10^7*R。m22=10^4+K1*S**10^4。m32=m31+K1*T*n31。 VE=(m12m1)/n1。VI=(m22m2)/n2。VP=(m32m31)/n31。 end 22 經(jīng)過循環(huán)迭代結(jié)果為 23 E正態(tài)、 PI I型分布 這種組合先把 PII型分布轉(zhuǎn)為正太變量，用計(jì)算后的平均值、標(biāo)準(zhǔn)差求可靠度，因?yàn)閮?nèi)容太多 就只寫程序碼和運(yùn)算結(jié)果，不詳細(xì)寫了。 m1=2*10^7。n1=*10^7。m2=10^4。n2=*10^4。m3=4。n3=1。 c=。d=。 Fm=exp(exp(c*(m2d)))。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 m21=m2n21*norminv(Fm,0,1)。 Fm=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 G=m1**m3。 K0=G/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 T=n31*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 24 第二次迭加 m12=m1+K0*n1*R。m22=m21+K0*S*n21。m32=m31+K0*T*n31。 VE=(m12m1)/n1。VI=(m22m21)/n21。VP=(m32m31)/n31。 m1=m12。n1=*10^7。m2=m22。n2=*10^4。m3=m32。n3=1。 c=。d=。 Fm=exp(exp(c*(m2d)))。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 m21=m2n21*norminv(Fm,0,1)。 Fm=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 G=m1**m3。 Q=Gm21*n1*VEm1*n21*VI+*n31*VP。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 T=n3*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 m12=2*10^7+K1**10^7*R。m22= m21+K1*S*n21。m32=m31+K1*T*n31。 VE=(m12m1)/n1。VI=(m22m21)/n21。VP=(m32m31)/n31。 25 循環(huán)迭加 while abs(K1K0)*K0 m1=m12。n1=*10^7。m2=m22。n2=*10^4。m3=m32。n3=1。 c=。d=。 K0=K1。 Fm=exp(exp(c*(m2d)))。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 m21=m2n21*norminv(Fm,0,1)。 Fm=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 G=m1**m3。 Q=Gm21*n1*VEm1*n21*VI+*n31*VP。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 T=n3*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 m12=2*10^7+K1**10^7*R。m22= m21+K1*S*n21。m32=m31+K1*T*n31。 VE=(m12m1)/n1。VI=(m22m21)/n21。VP=(m32m31)/n31。 end 26 最終運(yùn)行結(jié)果為 : 27 EIP都為 I型分布 這種組合應(yīng)該先把 EIP都轉(zhuǎn)化為正態(tài)分布，在利用均值、標(biāo)準(zhǔn)差求可靠度 第一次迭加： m1=2*10^7。n1=*10^7。m2=10^4。n2=*10^4。m3=4。n3=1。 c=。d=。 Fm=exp(exp(c*(m1d)))。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 m11=m1n11*norminv(Fm,0,1)。 Fm=exp(exp(c*(m2d)))。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 m21=m2n21*norminv(Fm,0,1)。 Fm=exp(exp(c*(m3d)))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m31=m3n31*norminv(Fm,0,1)。 G=m1**m3。 K0=G/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*