QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2900|回复: 1
打印 上一主题 下一主题

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

[复制链接]
字体大小: 正常 放大
dudu84        

1

主题

0

听众

17

积分

升级  12.63%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2003-1-15 11:40 |只看该作者 |正序浏览
|招呼Ta 关注Ta

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

. C3 @* q' F2 r# J) a' |6 R+ H) [) R+ w, [# e' d

%
% Programme pour leopard
%

. i. U: q" I' x$ d

global d K F N

* v, x5 w& A$ v. b! [! U [+ O

N = 64

2 z% G( j, Q! v2 ]0 M* {

dx = 1/N;

/ D/ j* O K% Z

d=10
K=0.01
F=0.1

7 ?* D" p8 O2 U' }0 {


dt=0.001;

) P" x* q8 w8 A( _% k& p5 J( e( q B

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

, A/ H) O5 m* ?6 K

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

1 d% Z% F. r$ u3 K% _3 i b7 J) w

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

. k D Y' Y. Y% d

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

/ K! Y' }& q0 X5 h2 j

compteur=0;

; ]8 Z$ N9 G ]$ H# b( W9 b

for t=0:dt:1,

3 W2 z8 |9 r" ~; L

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

1 j( ^# p/ {! q) u( Y

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

6 d: G1 ~; N) |8 s8 H: }8 x

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

! i0 q6 P' Z# j4 d* U" g$ z5 m

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

9 i! y) q7 u0 A( B8 f

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

b) C# d! x; L# T8 g+ t

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

8 n5 q5 i& _: C; {% ^2 }1 S8 C* z% _

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

9 J% J \+ S, ]' F0 q0 A

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

8 i* D0 V$ ?! F# G

compteur = compteur +1;

: r+ d. `6 i! t5 F# R5 u

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

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

0

主题

3

听众

20

积分

升级  15.79%

该用户从未签到

新人进步奖

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-8-6 14:39 , Processed in 0.426226 second(s), 63 queries .

回顶部