步骤

操作

第一步

定义初始值 $c,d,h,s,{\lambda }_1,{\lambda }_2,{\mu }_1,{\mu }_2,{\mu }_{2v},\theta ,{\theta }_1,\eta \left(\eta ={10}^{-10}\right)$ ，且设 $R=0$

第二步

输入矩阵 $A$ ， $B$ ， $C$

第三步

令 $R_n=R$ ，可求出 $R_{n+1}=-\left(R_n^{2}B+C\right)A^{-1}$

第四步

若 $\Vert R_{n+1}-R_n\Vert >\eta$ ，则令 $n=n+1$ 返回第三步继续迭代，

若 $\Vert R_{n+1}-R_n\Vert <\eta$ ，停止迭代，进入第五步。

第五步

得到率阵 $R$ 的近似解： $R=R_{n+1}$