在Von-Mises屈服准则及正交流动准则的前提下，建立了循环载荷下叠加型A-F(Armstrong-Frederick)非线性随动强化模型的迭代算法，并根据塑性应变增量的收敛控制实现内部的平衡迭代。为验证本文数值方法的正确性，以Chaboche和Ohno-Abde-Karim随动硬化模型为例，将本文方法的计算结果与通用有限元软件ANSYS的分析结果及试验数据进行了比较，均吻合良好，验证了本文算法的可靠性。

近，我们这两周进行光电综合实验，要求使用Matlab的GUI编程。有一道题目是用Fox_Li数值迭代法对平行平面腔的自再现模求解近似结果，绘制最终自再现模的振幅和相位分布曲线。我在网上搜到了等于是该题答案的源代码，很开心。代码的可读性还算勉强，但遗憾的在表达第1次、第299次、第300次迭代时竟然使用的是几乎同样的代码段。为什么不在Matlab的函数里再设置一个子函数来实现迭代的过程呢？这样可以明显提高代码的简洁性和可读性嘛！而且该程序在变量命名方面也显得稍微有些混乱。于是，针对这两点，我在自己的源程序里给出了优化，并且经过多次调试，理论上自至臻完美。另外，由于Fox_Li数值迭代法涉及循环和积分，耗时非常长。我在程序里又增加了一个计时器，迭代300次，大约需要150s的时间，还在我们的tolerance范围内。当然，我的程序是属于Windows Application的范畴，有个漂亮的、友好的界面。

平行平面镜--自在现模形成过程--Fox-Li数值迭代法 % 参数初始化 clc;clear; lambda=600e-9; %波长 L=100*lambda; %腔长 a=25*lambda; %腔镜线宽 k=2*pi/lambda; %波矢 x1=linspace(-a,a,1000); %取1000个点积分；linspace用于产生-a,a之间的N点行矢量，相邻数据跨度相同。

非稳腔的FOX-LI数值迭代解法：利用FOX—LI数值迭代解法，对非稳腔的模式进行数值求解。用Matlab软件程序模拟出腔内光场的振荡情况，给出腔内光场振幅和相位分布，计算出每次迭代能量的损耗。

数值迭代算法基本原理，及其Matlab具体实现，内附一具体实例，绝对货真价实。