資源簡介
一個論壇上的程序,自己運行調(diào)試后通過。包括了Lozi混沌系統(tǒng)的映射圖、分岔圖、最大Lyapunov指數(shù)計算和混沌同步。
代碼片段和文件信息
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混沌映射的線性反饋控制.河北師范大學(xué)學(xué)報.200731(4):479-483.
clear?all
close?all
clc
%=======================================================Lozi系統(tǒng)映射圖
%?X=[1;0.1;1.7;.5];
%?%前面兩個賦初值后面兩個給出了系統(tǒng)中的
%?%常數(shù)項的值!與newx=Lozi(x)中給出的列向量相對應(yīng)!
%?Y=[];?%給出一個空數(shù)組,以便存儲得出的迭代解
%?for?i=1:10000?%循環(huán)即迭代次數(shù)
%?????X=feval(@LoziX);%調(diào)用主函數(shù)
%?????Y(i:)=X(1:21);
%?????%將每一次迭代后的數(shù)組X的上面兩個變量存儲!
%?????%因為下面的兩個常數(shù)項進行迭代時保持原來的數(shù)值不便!
%?end
%?plot(Y(:1)Y(:2)‘.‘‘markersize‘1);
%?%得到Y(jié)=[Y(:1)Y(:2)]分別指代每一次迭代后的
%?%數(shù)組X的上面兩個變量然后點畫圖.
%?title(‘Lozi映射圖‘)%加上圖像標(biāo)題
%?xlabel(‘x‘)ylabel(‘y‘)
%?%分別給對應(yīng)的向量隨處的坐標(biāo)命名
%?%另外還可以使用其它的語句對圖像進行進一步的調(diào)整
%?%====================================================分岔圖1
%?Z=[];%給出一個空數(shù)組,以便存儲得出的迭代解!
%?for?p=linspace(01.71700);
%?????%將分岔參數(shù)用linspace在區(qū)間[01.7]
%?????%等分成1700份
%?????x=[1;0;p;.5];
%?????%前面兩個賦初值后面兩個給出了系統(tǒng)中的
%?????%常數(shù)項的值!與newx=Lozi(x)中給出的列‘
%?????%向量相對應(yīng)其中值得注意的是p是連續(xù)變化
%?????%的向量.
%?????for?k=1:500?%循環(huán)即迭代次數(shù)
%?????????x=Lozi(x);?%調(diào)用主函數(shù)
%?????????if?k>400
%????????????%取迭代后的最后一百個點這里認為
%????????????%前400個點為迭代過程中的瞬態(tài)響應(yīng)!
%????????????Z=[Zp+x(1)*i];
%????????????%將P的值以及在這個值所對應(yīng)最后400次迭代振動資訊R/mrEK$xBFi?ms5MR
%????????????%后的數(shù)組X的第一個變量x存儲為一個二維列向量組!振動資訊8lr#Iv;OlZ2xM:_?|~
%?????????end
%?????end
%?end
%?plot(Z‘.‘‘markersize‘1)
%?%得到Z=[Z(:1)Z(:2)]分別指代P的值以及在這個
%?%值所對應(yīng)最后400次迭代后的數(shù)組X的第一個變量x然后點畫圖.
%?title(‘Lozi映射分岔圖‘)%加上圖像標(biāo)題
%?xlabel(‘p‘)ylabel(‘x‘)
%?%分別給對應(yīng)的向量隨處的坐標(biāo)命名
%?%另外還可以使用其它的語句對圖像進行進一步的調(diào)整
%?%==========================================================分岔圖2
%?P=[];%給出一個空數(shù)組,以便存儲P的值!
%?Z=[];%給出一個空數(shù)組,以便存儲得出的迭代解!
%?for?p=linspace(01.71700);
%?????%將分岔參數(shù)用linspace在區(qū)間[01.7]
%?????%等分成1700份
%?????x=[1;0;p;.5];
%?????%前面兩個賦初值后面兩個給出了系統(tǒng)中的
%?????%常數(shù)項的值!與newx=Lozi(x)中給出的列振動資訊“b&qaS[dc
%?????%向量相對應(yīng)其中值得注意的是p是連續(xù)變化的向量.
%?????for?k=1:500
%?????????%循環(huán)即迭代次數(shù)但是遺憾的是這里的迭代次數(shù)太少振動資訊+oR-V?I7[
%?????????%因為如果太多的話由于上面的p點太多運行的時間
%?????????%就越長但是為了得到較為精確的迭代解建議適當(dāng)?shù)脑黾拥螖?shù)
%?????????%估計在10000左右而這時候取最后的1000個點即可!
%?????????x=Lozi(x);?%調(diào)用主函數(shù)
%?????????if?k>400
%????????????%取迭代后的最后一百個點這里認為
%????????????%前400個點為迭代過程中的瞬態(tài)響應(yīng)!
%????????????Z=[Z;x(1)];
%????????????%這個值所對應(yīng)最后400次迭代
%????????????%后的數(shù)組X的第一個變量x存儲為一個二維列向量組!
%????????????P=[P;p];?%將P的值存儲到空變量!
%??????????end
%?????end
%?end
%?plot(PZ‘.‘‘markersize‘1)%點畫圖振動資訊&J;S?l~7iQ
%?title(‘Lozi映射分岔圖‘)%加上圖像標(biāo)題
%?xlabel(‘p‘)ylabel(‘x‘)
%?%分別給對應(yīng)的向量隨處的坐標(biāo)命名
%?%另外還可以使用其它的語句對圖像進行進一步的調(diào)整
%?%==================================================最大Lyapunov指數(shù)圖
%?d0=1e-8;%定義極小的擾動還可以定義更小
%?Z=[];%給出一個空數(shù)組,以便存儲得出的迭代解!
%?P=[];%給出一個空數(shù)組,以便存儲P的值!
%?for?p=linspace(01.71700)
%?????%將分岔參數(shù)用linspace在區(qū)間[01.7]
%?????%等分成1700份振動資訊]“dR4~:v$h‘Lpp?h
%?????le=0;%定
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????6282??2009-05-18?10:19??Lozi混沌系統(tǒng)\Lozi_system_analysis.m
?????文件?????256440??2009-05-18?09:25??Lozi混沌系統(tǒng)\Lozi混沌映射的線性反饋控制[1].caj
?????目錄??????????0??2009-05-25?16:42??Lozi混沌系統(tǒng)
-----------?---------??----------?-----??----
???????????????262722????????????????????3
評論
共有 條評論