QQ登录

只需要一步,快速开始

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

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

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

1

主题

0

听众

17

积分

升级  12.63%

该用户从未签到

新人进步奖

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

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

. v$ \5 }7 l$ O& ~4 [, R3 A ?, o$ D) I1 h6 I) I. F0 v4 [

%
% Programme pour leopard
%

: P1 H) J v3 V8 ~. f2 \& ?

global d K F N

2 ]. W& X# t7 v! o d- f% b9 t

N = 64

* [9 @$ m" B8 a( R0 m* U

dx = 1/N;

& R r; C0 o+ ]8 W0 x" m

d=10
K=0.01
F=0.1

3 [* o& G) [7 z: }3 h8 R


dt=0.001;

5 m1 t, D1 J, g. p: i: H [2 k' G9 T

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

, x0 K2 e+ o2 x6 _; Y9 S" P

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

# m; B; i, m6 p! H

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

$ N# _6 Z& _3 M

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

) u4 |, ?; J n2 m, v$ k

compteur=0;

7 A- {( J" I3 L3 ~6 K- G+ U5 ^

for t=0:dt:1,

, I2 w/ ]/ x" K3 T! t

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

0 Q' _- D! g! k9 W! p- L+ s

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

% F- ?+ u ?4 w& d

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

{* ^7 d! f5 Y

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

+ s( u6 H, ^; P

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

% w3 j/ d- t. T2 G+ V. e

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

' Q$ X- ^# x# `8 E- @4 S) s

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

6 d( ?1 @% g5 I1 x: e; L

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

! X5 I5 C' ~ i4 b+ _0 ]

compteur = compteur +1;

: w. q0 D/ u$ Y+ A" ]. x5 q

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 13:36 , Processed in 0.404572 second(s), 63 queries .

回顶部