[求助]大家帮忙看看这段matlab程序

dudu84

以下是一段解偏微分方程组的程序，谁能告诉我用的是什么方法解的，我怎么查了好多资料都看不出来啊，谢谢大家了！

( j! \/ |0 H3 V3 ^ 5 C! a- T7 Z! `& o

%
% Programme pour leopard
%

global d K F N

N = 64

. @8 ?: y8 C/ j+ f; W) G* s0 ^

dx = 1/N;

d=10
K=0.01
F=0.1

5 h1 D0 E" n+ y

dt=0.001;

. V' t" ?* F2 W+ o: l2 M! R

U=rand(N);
V=rand(N);

: \% y _$ Y/ @+ @9 R3 b7 i: O7 ]

FourU=fft2(U);
FourV=fft2(V);

for i=1:N,
for j=1:N,
i0=i-1;
j0=j-1;
if i0 > N/2
i0=1+N-i;
end
if j0 > N/2
j0=1+N-j;
end

FourU(i,j)=FourU(i,j)/(1+(i0*i0+j0*j0)/8);
FourV(i,j)=FourV(i,j)/(1+(i0*i0+j0*j0)/8);
end
end

' X0 O) K1 S" e9 e6 c

U=real(ifft2(FourU));
V=real(ifft2(FourV));

/ I0 D1 ?! g1 q

compteur=0;

for t=0:dt:1,

8 h: I. I& b H4 n5 V }2 [; d

format long
max(max(U))
max(max(V))
min(min(V))

7 W5 }+ g7 E2 W9 [" j) ?* |& f" [' [

FourU=fft2(U);
FourV=fft2(V);

for i=1:N,
for j=1:N,
i0=i-1;
j0=j-1;
if i0 > N/2
i0=1+N-i;
end
if j0 > N/2
j0=1+N-j;
end

FourU(i,j)=FourU(i,j)*exp(-dt*(i0*i0+j0*j0)/10);
FourV(i,j)=FourV(i,j)*exp(-d*dt*(i0*i0+j0*j0)/10);
end
end

: o9 u, e+ F5 L! s) j7 h# S

U=real(ifft2(FourU));
V=real(ifft2(FourV));

, \; Y% `2 {8 }4 m

for i=1:N,
for j=1:N,
P(i,j)=-U(i,j)*V(i,j)^2+F*(1-U(i,j));
G(i,j)=U(i,j)*V(i,j)^2-(F+K)*V(i,j);
end
end

9 C: K: l9 _' j* u0 A

U = U + dt*P;
V = V + dt*G;

for i=1:N,
for j=1:N,
W(i,j)=V(i,j);
M(i,j)=U(i,j);
end
end

) g V1 H+ d" g0 z2 t7 H

if rem(compteur,10) == 0
clf;
if W>0.6
image(100*W-60);
else
image(100*W)
end
if M>0.6
image(100*M-60);
else
image(100*M);
end
drawnow;
end

9 U! ]8 y; P F5 y$ u4 G) X/ f/ l

compteur = compteur +1;

end
%------------------------------------------

争气的败家子

MATLAB 吧 里面有rand随机数生成，有fft傅立叶算法，得等，可能不会错吧。。