df_mthod.m

这是国外用的研究分岔的完整的M程序

function method=df_mthod(kind)

% function method=df_mthod(kind)
% INPUT:
%	kind the kind of default method wanted
% OUTPUT:
%	method default method

% (c) DDE-BIFTOOL v. 2.00, 30/11/2001

% during continuation:
method.continuation.steplength_condition=1; % use steplength condition
method.continuation.plot=1; % plot new points
method.continuation.prediction=1; % use secant prediction
method.continuation.steplength_growth_factor=1.2; % grow steplength with factor 1.2
method.continuation.plot_progress=1; % plot progress gradually
method.continuation.plot_measure=[]; % use default plot measures
method.continuation.halt_before_reject=0; % rejection of points allowed

if kind=='stst'
  % stst
  method.point.newton_max_iterations=5;
  method.point.newton_nmon_iterations=1;
  method.point.halting_accuracy=1e-10;
  method.point.minimal_accuracy=1e-8;
  method.point.extra_condition=0;
  method.point.print_residual_info=0;
  % root
  [alpha,beta]=time_lms('bdf',4);
  method.stability.lms_parameter_alpha=alpha;
  method.stability.lms_parameter_beta=beta;
  method.stability.lms_parameter_rho=time_saf(alpha,beta,0.01,0.01);
  method.stability.interpolation_order=4;
  method.stability.max_number_of_eigenvalues=100;
  method.stability.minimal_time_step=0.01;
  method.stability.maximal_time_step=0.1;
  method.stability.minimal_real_part=[];
  method.stability.max_newton_iterations=6;
  method.stability.root_accuracy=1e-6;
  method.stability.remove_unconverged_roots=1;
  if(nargin('sys_tau'))~=0
    method.stability.delay_accuracy=-1e-8;
  end;
elseif kind=='fold'
  % fold
  method.point.newton_max_iterations=5;
  method.point.newton_nmon_iterations=1;
  method.point.halting_accuracy=1e-9;
  method.point.minimal_accuracy=1e-7;
  method.point.extra_condition=0;
  method.point.print_residual_info=0;
  % root
  [alpha,beta]=time_lms('bdf',4);
  method.stability.lms_parameter_alpha=alpha;
  method.stability.lms_parameter_beta=beta;
  method.stability.lms_parameter_rho=time_saf(alpha,beta,0.01,0.01);
  method.stability.interpolation_order=4;
  method.stability.max_number_of_eigenvalues=100;
  method.stability.minimal_time_step=0.01;
  method.stability.maximal_time_step=0.1;
  method.stability.minimal_real_part=[];
  method.stability.max_newton_iterations=6;
  method.stability.root_accuracy=1e-6;
  method.stability.remove_unconverged_roots=1;
  if(nargin('sys_tau'))~=0
    method.stability.delay_accuracy=-1e-8;
  end;
elseif kind=='hopf'
  % hopf
  method.point.newton_max_iterations=5;
  method.point.newton_nmon_iterations=1;
  method.point.halting_accuracy=1e-9;
  method.point.minimal_accuracy=1e-7;
  method.point.extra_condition=0;
  method.point.print_residual_info=0;
  % root
  [alpha,beta]=time_lms('bdf',4);
  method.stability.lms_parameter_alpha=alpha;
  method.stability.lms_parameter_beta=beta;
  method.stability.lms_parameter_rho=time_saf(alpha,beta,0.01,0.01);
  method.stability.interpolation_order=4;
  method.stability.max_number_of_eigenvalues=100;
  method.stability.minimal_time_step=0.01;
  method.stability.maximal_time_step=0.1;
  method.stability.minimal_real_part=[];
  method.stability.max_newton_iterations=6;
  method.stability.root_accuracy=1e-6;
  method.stability.remove_unconverged_roots=1;
  if(nargin('sys_tau'))~=0
    method.stability.delay_accuracy=-1e-8;
  end;
elseif kind=='psol'
  % point
  method.point.newton_max_iterations=5;
  method.point.newton_nmon_iterations=1;
  method.point.halting_accuracy=1e-8;
  method.point.minimal_accuracy=1e-6;
  method.point.extra_condition=0;
  method.point.print_residual_info=0;
  method.point.phase_condition=1;
  method.point.collocation_parameters=[];
  method.point.adapt_mesh_before_correct=0;
  method.point.adapt_mesh_after_correct=3;
  % mult
  method.stability.max_number_of_eigenvalues=100;
  method.stability.minimal_modulus=0.01;
  method.stability.collocation_parameters=[];
  if(nargin('sys_tau'))~=0
    method.stability.delay_accuracy=-1e-8;
  end;
elseif kind=='hcli'
  % point
  method.point.newton_max_iterations=10;
  method.point.newton_nmon_iterations=1;
  method.point.halting_accuracy=1e-8;
  method.point.minimal_accuracy=1e-6;
  method.point.extra_condition=0;
  method.point.print_residual_info=0;
  method.point.phase_condition=1;
  method.point.collocation_parameters=[];
  method.point.adapt_mesh_before_correct=0;
  method.point.adapt_mesh_after_correct=3;
end;

return;