parrm.m

来自「基于matlab的约束非线性规划算法库」· M 代码 · 共 16 行

function alpha = parrm(v1,v2,v3,alphaIN)
%Call: alpha = parrm(v1,v2,v3,alphaIN)
%
v3n=dot(v3,v3);
coff(1)=1;
coff(2)=(1.5*dot(v2,v3))/v3n;
coff(3)=(dot(v1,v3)+0.5*dot(v2,v2))/v3n;
coff(4)=0.5*dot(v1,v2)/v3n;
%
rr=roots(coff);
rr=rr(find(imag(rr) == 0));
alpha=min(rr(find(rr > 0)));
if alpha > alphaIN
    alpha=alphaIN/2;
end
%pause