Znn3bq.jpeg
²é¿´: 1318  |  »Ø¸´: 7
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

xwndf250

Òø³æ (СÓÐÃûÆø)

[ÇóÖú] matlab´¦Àí³£Î¢·Ö·½³Ì×÷ͼÎÊÌâ

Á÷Ðв¡Ä£ÐÍ sir
dI/dt=a*S*I-b*I
dS/dt=-a*S*I
dR/dt=b*I
S+I+R=1ÇÒ0 ÏëÒª×öÒ»¸öºá×ø±êt=[0,50]£¬ ×Ý×ø±êdI/dt¡¢dS/dt¡¢dR/dtµÄ3ÌõÇúÏߣ¬Çó¾ßÌå³ÌÐò£¬±¾ÈËдµÄ³ÌÐò
mÎļþ£ºfunction y=SIR(t,x)
a=0.2;b=0.1;
y=[a*x(1)*x(2)-b*x(1);
-a*x(1)*x(2);
b*x(3)];
end
ÃüÁî:Ó¦¸ÃÔõÑùд£¬²ÅÊÇtÓëdI/dt¡¢dS/dt¡¢dR/dtµÄͼÏñ£¨¾ÍÊÇS¡¢I¡¢RµÄ±ä»¯ÂÊËæÊ±¼äµÄ±ä»¯£©£¬×¢Ò⣬²»ÊÇtÓëS¡¢I¡¢RµÄͼÏñ£¡×îºÃÓÃmatlabÓïÑÔ¡£
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

xwndf250

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
5Â¥: Originally posted by csgt0 at 2013-02-27 10:37:52
¾ÍÊÇt=0ÊǵÄISR°¡£¬¾ßÌå¶àÉٵÿ´ÄãµÄʵ¼ÊÇé¿ö£¬ÎÒÖ»ÊÇËæÒâдµÄ¸öÊý¡£...

ºÃµÄ ллÄã
6Â¥2013-02-27 11:09:25
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 8 ¸ö»Ø´ð

csgt0

ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)

²ÊÉ«¹Òͼ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
dbb627: ½ð±Ò+2, ¸ÐлӦÖú 2013-02-26 14:07:58
CODE:
function xwn
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
intvalue=[0.1 0.4 0.5];       %t=0ʱµÄ³õÖµ
[T,Y] = ode45(@rigid,[0 50],intvalue,options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
title('T-Yͼ')
for i=1:length(T)     
vdy(:,i)=rigid(T(i),Y(i,:));
end
figure
plot(T,vdy(1,:),'-',T,vdy(2,:),'-.',T,vdy(3,:),'.')
title('T-dYͼ')
end

function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
a=0.2;
b=0.1;
dy(1) = a*y(2)*y(1)-b * y(1);
dy(2) = -a*y(2) * y(1);
dy(3) = b * y(1);
end

showmethemoney
2Â¥2013-02-26 11:06:11
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

xwndf250

Òø³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by csgt0 at 2013-02-26 11:06:11
function xwn
options = odeset('RelTol',1e-4,'AbsTol',);
intvalue=;       %t=0ʱµÄ³õÖµ
= ode45(@rigid,,intvalue,options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
title('T-Yͼ')
for i=1: ...

function dy = rigid(t,y)ºóÃæµÄÎÒ¿´¶®ÁË£¬µ«ÊÇÇ°Ãæfunction xwnûÓп´¶®£¬ÎҼǵÃÇóµ¼Ã²ËÆÓÃdiffº¯ÊýµÄ£¬Äܲ»ÄÜÏêϸ½²½âÏÂodeset£¨£©º¯Êý£¬intvalue×÷Ó㬡°for i=1:length(T) ºÍvdy(:,i)=rigid(T(i),Y(i,);¡±µÄÒâ˼¡£ÎÒÊÇС°×£¬ÇóÖ¸½Ì¡£ºÃµÄ»°ÎÒ¶à¸ø½ð±Ò¡£
3Â¥2013-02-26 21:06:43
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

csgt0

ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)

²ÊÉ«¹Òͼ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
xwndf250: ½ð±Ò+15, ¡ï¡ï¡ïºÜÓаïÖú, ·Ç³£¸ÐлÄúµÄ°ïÖú£¬¸øÄú½ð±Ò£¬Ë³±ãÎÊÏ£¬t0³õÖµÔõôȷ¶¨µÄ£¿Ò²¾ÍÊÇ0.1£¬0.4£¬0.5ÊÇÔõôȷ¶¨µÄ£¿ 2013-02-27 10:00:25
fegg7502: ½ð±Ò+1, Ó¦ÖúÖ¸Êý+1, ¹ÄÀø½»Á÷ 2013-04-02 09:30:14
odeset£¨£©º¯ÊýÉèÖÃÊýÖµ½âµÄ¼ÆË㾫¶È£¬Æäʵһ°ã²»ÓÃÒ²¿ÉÒÔ
intvalueÊdzõÖµ°¡£¬t=0ʱµÄ³õÖµ
¡°for i=1:length(T) ºÍvdy(:,i)=rigid(T(i),Y(i,);¡±µÄΪÁËÑ­»·¼ÆËãÿ¸ötϵÄ3¸öµ¼ÊýÖµ
showmethemoney
4Â¥2013-02-27 09:54:46
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ʳƷÓëÓªÑø£¨0955£©271Çóµ÷¼Á +7 Éý¸ñ°¢´ï 2026-04-12 7/350 2026-04-12 22:56 by yuyin1233
[¿¼ÑÐ] ҩѧÇóµ÷¼Á +3 RussHu 2026-04-12 4/200 2026-04-12 17:49 by ³ÂƤƤ
[¿¼ÑÐ] 26×ÔÈ»µØÀíѧ303·ÖÇóµ÷¼Á +6 Ò»Õ½³É˶°¡°¡°¡° 2026-04-06 11/550 2026-04-11 21:27 by Åáºêΰ
[¿¼ÑÐ] µ÷¼Á +10 Ö»ÐðÀë±ð´Ç 2026-04-09 12/600 2026-04-11 20:57 by ÄæË®³Ë·ç
[¿¼ÑÐ] 283Çóµ÷¼Á 086004¿¼Ó¢¶þÊý¶þ +17 ÄǸöàà×Ó 2026-04-10 18/900 2026-04-11 16:27 by Ã÷Ô´ËʱÓÐ
[¿¼ÑÐ] Ò»Ö¾Ô¸211£¬»¯Ñ§310·Ö£¬±¾¿ÆÖصãË«·Ç£¬Çóµ÷¼Á +23 ŬÁ¦·Ü¶·112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[¿¼ÑÐ] 263ÄÜÔ´¶¯Á¦×¨Ë¶Çóµ÷¼Á +4 ¼Ó´óºÅ·¹ºÐ´ü 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[¿¼ÑÐ] 265Çóµ÷¼Á +12 ·ç˵ËýÔçÍüÁË 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[¿¼ÑÐ] 071000ÉúÎïѧµ÷¼ÁÇóÖú +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by Çóµ÷¼Ázz
[¿¼ÑÐ] 293µ÷¼Á +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by ²ñС°×
[¿¼ÑÐ] Ò»Ö¾Ô¸»ª¶«Ê¦·¶ÉúÎïѧ326·Ö£¬Çóµ÷¼Á +8 Áõīī 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[¿¼ÑÐ] Òѵ÷¼Á +18 ²ñ¿¤Ã¨_ 2026-04-09 19/950 2026-04-09 22:10 by ²ñ¿¤Ã¨_
[¿¼ÑÐ] Çó»úеר˶297µÚ¶þÅúµ÷¼Á +5 ʰÆâ12¡£ 2026-04-08 5/250 2026-04-09 16:43 by Ôʵ±ÊʶÈ
[¿¼ÑÐ] ²ÄÁÏ307·ÖÇó´óÀÐ×éÊÕÁô +17 Hllºú 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[¿¼ÑÐ] 283µç×ÓÐÅÏ¢Çóµ÷¼Á +4 ÈýʯWL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤×¨Ë¶306·ÖÕÒºÏÊʵ÷¼Á +27 ²×º£ÇáÖÛe 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á 348·Ö +14 °¦ÎÒ³¬ÕæÃ»ÕÐÁË 2026-04-06 15/750 2026-04-08 19:16 by ÎÒ¼õ·Ê1
[¿¼ÑÐ] 304Çóµ÷¼Á +10 ËØÄê¼ÀÓï 2026-04-06 17/850 2026-04-08 09:05 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 344Çóµ÷¼Á +11 κ×Óper 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[¿¼ÑÐ] 085602µ÷¼Á ³õÊÔ×Ü·Ö335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û