File:Brusselator space.gif
From formulasearchengine
Jump to navigation
Jump to search
Brusselator_space.gif (125 × 125 pixels, file size: 1.44 MB, MIME type: image/gif, looped, 200 frames, 20 s)
This file is from Wikimedia Commons and may be used by other projects. The description on its file description page there is shown below.
Summary
DescriptionBrusselator space.gif |
English: some frames from a simulation of the brusselator model in 2D |
Date | |
Source | Own work |
Author | Jkrieger |
GIF development InfoField | ![]() This diagram was created with MATLAB. |
Source code InfoField | MATLAB codeclear;
% timestep
DeltaT=0.1;
% diffusion coefficients
DX=0.2;
DY=0.02;
% reaction constants
k1=1;
k2=1;
k3=1;
k4=1;
% initial concentrations
X0=1;
Y0=1;
% constant concentrations
A=1;
B=3;
% sigma of the random initial distortion
initRandomSigma=2;
% width/height of the simulation grid
width=100;
% number of iterations before display
iterationsPre=3000;
% number of visible iterations after pre
iterations=300;
iterationsDT=6;
% number of visible iterations before pre
iterationsStart=300;
% number of visible iterations before pre
iterationsInit=300;
% number of invisible iterations before iterationsStart
iterationsStartPre=600;
% number of visible iterations after pre
iterationsPre1=6000;
% number of visible iterations after pre
iterations1=300;
% laplace operator
laplace=0.25*[0 1 0; 1 -4 1; 0 1 0];
% initialize reaction variables
X(1:width,1:width)=X0+(rand(width,width)-0.5)*2*initRandomSigma;
Y(1:width,1:width)=Y0+(rand(width,width)-0.5)*2*initRandomSigma;
X(X<0)=0;
Y(Y<0)=0;
X(X>4.5)=4.5;
Y(Y>4.5)=4.5;
Xi=X;
Yi=Y;
% colormap
rbmap(1:100,1:3)=0;
rbmap(:,1)=((length(rbmap)-1)-(0:(length(rbmap)-1)))/(length(rbmap)-1);
rbmap(:,3)=(0:(length(rbmap)-1))/(length(rbmap)-1);
figure(1)
s=get(gcf, 'Position');
s(3)=250;
s(4)=200;
set(gcf, 'Position', s);
mf=1;
for it=1:(iterationsInit+iterationsStartPre+iterationsStart+iterationsPre+iterations+iterationsPre1+iterations1)
lapX=conv2(X, laplace, 'same');
lapY=conv2(Y, laplace, 'same');
vX=k1*A+k2*X.^2.*Y-k3*B*X-k4*X+DX.*lapX;
vY=k3*B*X-k2*X.^2.*Y+DY.*lapY;
% leave border as distrotion centers
vX(1:width,[1 width])=0;
vX([1 width],1:width)=0;
vY(1:width,[1 width])=0;
vY([1 width],1:width)=0;
if ((it<=iterationsInit) ...
|
Licensing
![w:en:Creative Commons](https://upload.wikimedia.org/wikipedia/commons/thumb/7/79/CC_some_rights_reserved.svg/90px-CC_some_rights_reserved.svg.png)
![attribution](https://upload.wikimedia.org/wikipedia/commons/thumb/1/11/Cc-by_new_white.svg/24px-Cc-by_new_white.svg.png)
This file is licensed under the Creative Commons Attribution 3.0 Unported license.
Attribution: Jan Krieger
- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
Captions
Add a one-line explanation of what this file represents
Items portrayed in this file
depicts
3 August 2012
image/gif
d4e6c595dfd0671268202f6c17def14540cfb977
1,506,406 byte
20.000000000000014 second
125 pixel
125 pixel
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 11:23, 8 August 2012 | ![]() | 125 × 125 (1.44 MB) | wikimediacommons>Jkrieger | more regions |
File usage
There are no pages that use this file.
Retrieved from "https://en.formulasearchengine.com/wiki/File:Brusselator_space.gif"