Chaos Synchronization
lorenz_Synchronous.m
function lorenz_Synchronous
% Lorenz 系统同步
%
% Author:LDYU
% Author's email: ustb03-07@yahoo.com.cn
%
[T,Y]=ode45(@lorenzf1,[0,3],[1;-10;1;10;2]);
subplot(2,1,1)
plot(T,Y(:,3),'k')
hold on
plot(T,Y(:,5),'b')
[T,Y]=ode45(@lorenzf2,[0,3],[10;10;10;-20;2]);
subplot(2,1,2)
plot(T,Y(:,3),'k')
hold on
plot(T,Y(:,5),'b')
%-------------------------------------
function dx=lorenzf1(t,x)
% dx= a*(y-x)
% dy= -x*z+r*x-y
% dz= x*y-b*z
% dy1= -x*z1+r*x-y1
% dz1= x*y1-b*z1
a=16;b=4;r=45.92;
dx=[ a*(x(2)-x(1))
-x(1)*x(3)+r*x(1)-x(2)
x(1)*x(2)-b*x(3)
-x(1)*x(5)+r*x(1)-x(4)
x(1)*x(4)-b*x(5)];
%-------------------------------------
function dx=lorenzf2(t,x)
% dx= a*(y-x)
% dy=-x*z+r*x-y
% dz= x*y-b*z
% dx1=a*(y-x1)
% dz1=x1*y-b*z1
a=16;b=4;r=45.92;
dx=[ a*(x(2)-x(1))
-x(1)*x(3)+r*x(1)-x(2)
x(1)*x(2)-b*x(3)
a*(x(2)-x(4))
x(4)*x(2)-b*x(5)];
9 comments:
顶~
ӏ'm really inspired with your writing skills and also with the format for your weblog. Is that this a paid subject matter or did you modify it your self? Anyway stay up the excellent high quality writing, it is rare to peer a nice blog like this one nowadays..
My web page vitrification parquet
I enjοу, lеad to I found just whаt I wаs tаking a look for.
You've ended my 4 day lengthy hunt! God Bless you man. Have a great day. Bye
Take a look at my site www.fenetrepvcsurmesure.com
It's appropriate time to make a few plans for the future and it is time to be happy. I have learn this submit and if I could I wish to suggest you few fascinating issues or suggestions. Maybe you can write subsequent articles regarding this article. I want to read even more things approximately it!
Feel free to visit my homepage: phenix option
Ηmm is аnyone else having problemѕ with thе imаgеs on this blog loaԁing?
I'm trying to determine if its a problem on my end or if it'ѕ the blοg.
Any suggestions ωould be greаtly appгeciated.
Hаve a lοok at my homepage porte de garage sectionnelle hormann
There's certainly a lot to know about this topic. I really like all of the points you've made.
Stop by my page: devis en ligne fenetre pvc
Hello there! Do you know if they make any plugins to help with Search Engine Optimization?
I'm trying to get my blog to rank for some targeted keywords but I'm
not seeing very good success. If you know of any
please share. Many thanks!
Here is my website - motorisation volet roulant
Gooԁ dаy! I ϳust ωould like to give you а hugе thumbs up for yοur grеаt іnfo yоu've got right here on this post. I'll be returning to youг web
ѕite fοr more soon.
Feel frее to vіsіt my ωeb-site
:: prodog
Great post! We are linking to this particularly great article on our website.
Keep up the good writing.
My homepage ... gagner de l'argent
Post a Comment