用于研究非线性系统matlab模块

用于研究非线性系统的混沌性质。如洛伦兹系统

function varargout = matds（varargin）
% MATDS Application M-file for matds.fig
%    FIG = MATDS launch matds GUI.
%    MATDS（‘callback_name‘ ...） invoke the named callback.

% Last Modified by GUIDE v2.0 27-Jan-2003 12:54:54

% DS-system
global DS;

% Calculation constants
global session_time;
global aboutpress;
global WasChanged;
global calculation_progress first_call;


% Windows
global nbrwin;
global session_windows;
global driver_window;

% Parameters
global P;

% Path
global matdspath;

% Trajectory bufer
global TRJ_bufer Time_bufer;
global bufer_i;

% Initialization
  if nargin == 0
     calculation_progress=0;
     WasChanged = 0;
     matdspath.main = pwd;    
     tmpi=findstr（matdspath.main（ （length（pwd）-5）:length（pwd））‘matds‘）;
     if isempty（tmpi）==1
        fprintf（‘\n The MATDS-directory must be a current!‘）;
        fprintf（‘\n The work is impossible sorry...‘）;
        return
     end;

% Path update

     addpath（matdspath.main）;
     matdspath.systems=[matdspath.main filesep ‘systems‘];
     addpath（matdspath.systems）;
     matdspath.gui=[matdspath.main filesep ‘gui‘];
     addpath（matdspath.gui）;
     matdspath.temp=[matdspath.main filesep ‘temp‘];
     addpath（matdspath.temp）;
     matdspath.maths=[matdspath.main filesep ‘maths‘];
     addpath（matdspath.maths）;

     session_time = clock;
     fprintf（‘         Welcome to dynamical system investigation‘）;

     DS（1）.currfile = [matdspath.systems filesep ‘lorenz‘];

     session_windows.sysdef = [];
     session_windows.editsys = [];
     session_windows.vars = [];
     session_windows.params = [];
     session_windows.intdata = [];


    P（1）=30;
    P（2）=8/3;
    P（3）=10.0;

     DS（1）.Val_param=P;

%  Integration data
     DS（1）.time_start = 0.0;
     DS（1）.time_end   = 100.0;
     DS（1）.max_step = 0.1;
     DS（1）.abs_error = 0.000001;
     DS（1）.rel_error = 0.00001;
     DS（1）.method_int = 7;

%  DS data
     DS（1）.name = ‘lorenz‘;
     DS（1）.vars = {‘x‘ ‘y‘ ‘z‘};
     DS（1）.param ={‘r‘ ‘b‘ ‘sigma‘};
     DS（1）.periodic={0 0 0};
     DS（1）.periodic_value=[0 0 0];
     DS（1）.equations = {‘-sigma*（x-y）‘ ...
                   ‘-x*z+r*x-y‘ ...
                   ‘x*y-b*z‘ };


% Define current trajectory as initial point

     DS（1）.currenttrajectory = 0;

% Define text output regime. 0 - no text output to main MATLAB window; 1 -
% trajectory points output to main MATLAB window.
% Always equilibria cycles and other values print to main MATLAB window;
     DS（1）.trj_text_out = 0;

% Regime of time output to main window
     DS（1）.time_indication = 0;

% Define equation of poincare section {1} - lhs depending on variables {2}
% rhs const. In symbolic format.
     DS（1）.poincare_map = [];
% Switch variable for poincare map calculation:   0 - no section; 1 - section on time; 
%                               2 - section by plane in both direction
%                               3 - section by plane in positive direction
%                               4 - section by plane in

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2004-05-19 19:11  matds\gui\
     文件         934  2003-11-14 16:55  matds\gui\Check_And_Delete.m
     文件        3787  2004-05-13 13:31  matds\gui\Check_And_Load.m
     文件        1211  2003-11-14 16:55  matds\gui\Check_And_Write.m
     文件         150  2003-11-14 16:55  matds\gui\ColorPlot2Dtraj.m
     文件         220  2003-11-14 16:55  matds\gui\ColorPlot3Dtraj.m
     文件         217  2003-11-14 16:55  matds\gui\ColorPlotTraj.m
     文件        1036  2003-11-14 16:55  matds\gui\Delete2D.m
     文件        1025  2003-11-14 16:55  matds\gui\Deletewin.m
     文件         396  2003-11-14 16:55  matds\gui\EquilibriaPlot2D.m
     文件         396  2003-11-14 16:55  matds\gui\EquilibriaPlot3D.m
     文件        1435  2003-11-14 16:55  matds\gui\Load_point.m
     文件        2556  2003-11-14 16:55  matds\gui\Open2DWin.m
     文件        2334  2003-11-14 16:55  matds\gui\Open3DWin.m
     文件         392  2003-11-14 16:55  matds\gui\PoincarePlot2D.m
     文件         354  2003-11-14 16:55  matds\gui\PoincarePlot3D.m
     文件         549  2003-11-14 16:55  matds\gui\PointSizeSet.m
     文件         124  2003-11-14 16:55  matds\gui\RegionPlot2D.m
     文件         124  2003-11-14 16:55  matds\gui\RegionPlot3D.m
     文件        1243  2003-11-14 16:55  matds\gui\Save_As.m
     文件         543  2003-11-14 16:55  matds\gui\Save_point.m
     文件         362  2003-11-14 16:55  matds\gui\TimeIndication.m
     文件         867  2003-11-14 16:55  matds\gui\TrjPlot2D.m
     文件         898  2003-11-14 16:55  matds\gui\TrjPlot3D.m
     文件       17528  2003-11-14 16:55  matds\gui\about.fig
     文件         803  2003-11-14 16:55  matds\gui\about.m
     文件       11664  2003-11-14 16:55  matds\gui\calc_driver.fig
     文件        3660  2004-05-17 18:42  matds\gui\calc_driver.m
     文件         347  2003-11-14 16:55  matds\gui\checksymb.m
     文件         174  2003-11-14 16:55  matds\gui\clear2Dtraj.m
     文件         174  2003-11-14 16:55  matds\gui\clear3Dtraj.m
............此处省略119个文件信息