Lozi混沌系统相关分析的MATLAB程序

一个论坛上的程序，自己运行调试后通过。包括了Lozi混沌系统的映射图、分岔图、最大Lyapunov指数计算和混沌同步。

function Lozi_system_analysis
% Author: Thomas Lee
% E-mail: lixf1979@126.com
% Corresponding: School of Mathematics Physics and Software Engineering Lanzhou Jiaotong University Lanzhou 730070 China
% if you want to get more information please refer to one of the published articles of the author :
%李险峰等.Lozi混沌映射的线性反馈控制.河北师范大学学报.200731（4）:479-483.
clear all
close all
clc
%=======================================================Lozi系统映射图
% X=[1;0.1;1.7;.5];
% %前面两个赋初值后面两个给出了系统中的
% %常数项的值！与newx=Lozi（x）中给出的列向量相对应！
% Y=[]; %给出一个空数组，以便存储得出的迭代解
% for i=1:10000 %循环即迭代次数
%     X=feval（@LoziX）;%调用主函数
%     Y（i:）=X（1:21）;
%     %将每一次迭代后的数组X的上面两个变量存储！
%     %因为下面的两个常数项进行迭代时保持原来的数值不便！
% end
% plot（Y（:1）Y（:2）‘.‘‘markersize‘1）;
% %得到Y=[Y（:1）Y（:2）]分别指代每一次迭代后的
% %数组X的上面两个变量然后点画图.
% title（‘Lozi映射图‘）%加上图像标题
% xlabel（‘x‘）ylabel（‘y‘）
% %分别给对应的向量随处的坐标命名
% %另外还可以使用其它的语句对图像进行进一步的调整

% %====================================================分岔图1
% Z=[];%给出一个空数组，以便存储得出的迭代解！
% for p=linspace（01.71700）;
%     %将分岔参数用linspace在区间[01.7]
%     %等分成1700份
%     x=[1;0;p;.5];
%     %前面两个赋初值后面两个给出了系统中的
%     %常数项的值！与newx=Lozi（x）中给出的列‘
%     %向量相对应其中值得注意的是p是连续变化
%     %的向量.
%     for k=1:500 %循环即迭代次数
%         x=Lozi（x）; %调用主函数
%         if k>400
%            %取迭代后的最后一百个点这里认为
%            %前400个点为迭代过程中的瞬态响应！
%            Z=[Zp+x（1）*i];
%            %将P的值以及在这个值所对应最后400次迭代振动资讯R/mrEK$xBFi ms5MR
%            %后的数组X的第一个变量x存储为一个二维列向量组！振动资讯8lr#Iv;OlZ2xM:_ |~
%         end
%     end
% end
% plot（Z‘.‘‘markersize‘1）
% %得到Z=[Z（:1）Z（:2）]分别指代P的值以及在这个
% %值所对应最后400次迭代后的数组X的第一个变量x然后点画图.
% title（‘Lozi映射分岔图‘）%加上图像标题
% xlabel（‘p‘）ylabel（‘x‘）
% %分别给对应的向量随处的坐标命名
% %另外还可以使用其它的语句对图像进行进一步的调整

% %==========================================================分岔图2
% P=[];%给出一个空数组，以便存储P的值！
% Z=[];%给出一个空数组，以便存储得出的迭代解！
% for p=linspace（01.71700）;
%     %将分岔参数用linspace在区间[01.7]
%     %等分成1700份
%     x=[1;0;p;.5];
%     %前面两个赋初值后面两个给出了系统中的
%     %常数项的值！与newx=Lozi（x）中给出的列振动资讯“b&qaS[dc
%     %向量相对应其中值得注意的是p是连续变化的向量.
%     for k=1:500
%         %循环即迭代次数但是遗憾的是这里的迭代次数太少振动资讯+oR-V I7[
%         %因为如果太多的话由于上面的p点太多运行的时间
%         %就越长但是为了得到较为精确的迭代解建议适当的增加迭代次数
%         %估计在10000左右而这时候取最后的1000个点即可！
%         x=Lozi（x）; %调用主函数
%         if k>400
%            %取迭代后的最后一百个点这里认为
%            %前400个点为迭代过程中的瞬态响应！
%            Z=[Z;x（1）];
%            %这个值所对应最后400次迭代
%            %后的数组X的第一个变量x存储为一个二维列向量组！
%            P=[P;p]; %将P的值存储到空变量！
%          end
%     end
% end
% plot（PZ‘.‘‘markersize‘1）%点画图振动资讯&J;S l~7iQ
% title（‘Lozi映射分岔图‘）%加上图像标题
% xlabel（‘p‘）ylabel（‘x‘）
% %分别给对应的向量随处的坐标命名
% %另外还可以使用其它的语句对图像进行进一步的调整

% %==================================================最大Lyapunov指数图
% d0=1e-8;%定义极小的扰动还可以定义更小
% Z=[];%给出一个空数组，以便存储得出的迭代解！
% P=[];%给出一个空数组，以便存储P的值！
% for p=linspace（01.71700）
%     %将分岔参数用linspace在区间[01.7]
%     %等分成1700份振动资讯]“dR4~:v$h‘Lpp h
%     le=0;%定

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6282  2009-05-18 10:19  Lozi混沌系统\Lozi_system_analysis.m

     文件     256440  2009-05-18 09:25  Lozi混沌系统\Lozi混沌映射的线性反馈控制[1].caj

     目录          0  2009-05-25 16:42  Lozi混沌系统

----------- ---------  ---------- -----  ----

               262722                    3