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