宸�鍒嗘柟绋嬪強寰�鍒嗘柟绋嬫暟鍊艰В - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

 

 

1 

涓滃崡澶у�︺�婃暟瀛﹀疄楠屻�嬫姤鍛?/p>

 

 

瀹為獙鍐呭�癸細宸�鍒嗘柟绋嬪強寰�鍒嗘柟绋嬫暟鍊艰В

 

涓�

 

瀹為獙鐩�鐨�

 

 

 

 

 

 

鐔熸倝杩�浠ｆ硶鍙婂井鍒嗘柟绋嬫暟鍊兼柟娉?/p>

 

浜?/p>

 

棰勫�囩煡璇�

 

 

 

锛?/p>

1

锛変簡瑙ｅ樊鍒嗘柟绋嬬ǔ瀹氭�с�佸懆鏈熷垎瑙ｃ�佹贩娌岀瓑鐩稿叧鐭ヨ瘑

 

 

 

锛?/p>

2

锛変簡瑙ｆ�ф媺鏂规硶銆侀緳鏍?/p>

-

搴撶壒鏂规硶銆?/p>

 

涓?/p>

 

瀹為獙鍐呭�逛笌瑕佹�?/p>

 

 

 

 

锛堜竴锛?/p>

V

olterra

鏂圭▼鏁板�艰В

 

鏂圭▼

0

,

,

,

,

,

飥?/p>

飪?/p>

飪?/p>

飪?/p>

飪?/p>

飪?/p>

飪?

飪?/p>

飥?/p>

飥?/p>

飥?/p>

飥?/p>

飥?/p>

d

c

b

a

dxy

cx

dt

dy

bxy

ax

dt

dx

 

鍏朵腑

a=1,b=0.1,c=0.5,d=0.04 

.

2

)

0

(

,

25

)

0

(

飥?/p>

飥?/p>

y

x

 

鍛戒护涓庣粨鏋?/p>

 

鍦ㄥ嚱鏁扮紪杈戝櫒涓�杈撳叆锛�

 

function

 dxdt = euler( t,x ) 

dxdt=[ x(1)*(1-0.1*x(2))    

    x(2)*(-0.5+0.02*x(1))]; 

end

 

 

鍥涢樁榫欐牸

-

搴撳�斿叕寮忥�?/p>

 

鍦ㄥ懡浠ょ獥鍙ｄ腑杈撳叆锛?/p>

 

tspan=[0 15]; 

x0=[25;2]; 

[t,x]=ode45(@euler,tspan,x0); 

plot(t,x(:,1),'r-','LineWidth',0.5); 

hold on; 

plot(t,x(:,2),'g-','LineWidth',0.5); 

hold on; 

 

锛?/p>

1

锛?/p>

 

鐩?/p>

鍥?/p>

 

鍥涢樁榫欐牸

-

搴撳�斿叕寮忥�?/p>