We'll use these feed and kill rates f=.055 Īnd we'll need a couple of other variables. I'm going to do it using circshift with the following function: function out = my_laplacian(in) There are a couple of ways to implement this in MATLAB. This means that cells on the right side were treated as if they were adjacent to cells on the left and cells on the top were treated as if they were adjacent to cells on the bottom. The one which Karl used originally did a "toroidal wrap". MATLAB's del2 operator clamps at the edges of the grid. Where Aold and Bold are the concentrations of the two chemicals at the previous time step and Anew and Bnew are the concentrations at the next time step. You can write this pretty directly in MATLAB like this: Anew = Aold + (da*del2(Aold) - Aold.*Bold.^2 + f*(1-Aold))*dt īnew = Bold + (db*del2(Bold) + Aold.*Bold.^2 - (k+f)*Bold)*dt Where A and B are the concentrations of our two chemicals, and k and f are the "kill" and "feed" rates. We'll start from a wonderful version which Karl Sims wrote for the Connection Machine.
Let's try writing a simulation of the Gray-Scott model in MATLAB. There are a couple of different versions of this now, but one of the most interesting is the Gray-Scott model. This system came to be known as a two-component reaction diffusion system. It turns out that fiddling the reaction rates and the diffusion rates can yield all sorts of different patterns which look very similar to patterns you see in nature. They react with each other, and they each diffuse from areas with high concentrations to areas with low concentrations.
#Circshift matlab skin#
You have two chemicals which are distributed nonuniformly through the skin of the animal. It describes an idea he had for how all of the different patterns on the skins of various animals are formed. It amazed me that he could find time to write papers on subjects like computing the function with all of the other work he was doing during the war.īut I think that my favorite Turing paper is one he wrote after the war called The Chemical Basis of Morphogenesis. But I was also impressed by the collection of papers which Alan Turing wrote. There are lots of interesting things to see there, like old Enigma machines and a working replica of the bombe. A couple of years ago, I visited Bletchley Park.
0 kommentar(er)
