【摘要】:对于校正矩阵E可由具体的公式来计算,不同的公式对应不同的变尺度法,将在下面迸行讨论。但不论哪种变尺度法E必须满足拟牛顿条件M(k+1)y=z即y=z或Ey=z-My满足上式的E有无穷多个,因此上述变尺度法构成一族算法。图4-10 变尺度法的计算程序框图
1)选定初始点X0和收敛精度ε;
2)计算▽f(X(0)),选取初始对称正定矩阵M(0),如M(0)=I;
3)计算搜索方向,S(k)=-M(k)▽f(X(k));
4)沿S(k)方向迸行一维搜索,计算▽f(X(k+1)),S(k)=X(k+1)-X(k),y(k)=▽f(X(k+1))-▽f(X(k));
5)判断是否满足迭代终止准则,若满足X*=X(k+1),则结束,否则转6);
6)当迭代n次后还没找到极小点时,重置M(k)为单位矩阵I,并以当前设计点X(0)←X(k+1)为初始点,返回到2)迸行下一轮迭代,否则转到7);
7)计算矩阵M(k+1)=M(k)+E(k),置k←k+1返回到3)。
对于校正矩阵E(k)可由具体的公式来计算,不同的公式对应不同的变尺度法,将在下面迸行讨论。但不论哪种变尺度法E(k)必须满足拟牛顿条件(www.xing528.com)
M(k+1)y(k)=z(k)
即(M(k)+E(k))y(k)=z(k)
或E(k)y(k)=z(k)-M(k)y(k)
满足上式的E(k)有无穷多个,因此上述变尺度法(属于拟牛顿法)构成一族算法。
变尺度法的计算程序框图如图4-10所示。
图4-10 变尺度法的计算程序框图
免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。