hom_main.m

来自「同伦方程组求非线性方程组所有解」· M 代码 · 共 33 行

%Homotopy_Optimization_Method
clear;
global x00 lambda;
m=20;deta=1/m;
format long;
%%%%%%%%%%%%%%%%%%
%%%%%%%%%%(1)
n=4;
for i=1:n
    x1(i,1)=rand(1);
end
lambda0(1)=0;x00=x1;

xtemp1=x1;
for i=2:m+1
    lambda0(i)=lambda0(i-1)+deta;lambda=lambda0(i);
    
    x1=mulNewton_H(x1);BB='牛顿法(基本方法)';
    %x1=mulDNewton_H(x00);BB='牛顿下山法';
    % x1=mulGSND(x0);BB='高斯牛顿法';
   %  x1=mulDamp(x0);BB='阻尼最小二乘法';
    %x1=mulConj(x0);BB='共轭梯度法???????????????????';
    %x1=mulFastDown(x0);BB='最速下降法??????????';
    %x1=mulNumYT(x0);BB='数值延拓法!???/????????';
   %x1=newtoniteration6(x0);BB='牛顿迭代法改进《机床与液压》';
    x1=x1;f1=myfun(x1);
  end
%%%%%%%%%%
format short;
x1',f1'