±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1498  |  »Ø¸´: 1

jwj216

½ð³æ (ÕýʽдÊÖ)

[ÇóÖú] ³£Î¢·Ö·½³ÌµÄ²ÎÊýÄâºÏ ÒÑÓÐ1È˲ÎÓë

²ÎÊýÄâºÏÒ»Ö±³ö´í£¬µ«ÊÇÕÒ²»³ö´íÎó¡£ÄÄλÄܰïæ¸ø½â¾öһϣ¿ÏÈллÁË£¡



function Par_SI
clear
clc
ID=[1 3 25 72 222 282 256 233 189 123 70 25 11 4];
X0=[762,1];%³õʼÌõ¼þ
lb=[0,0];%Ô¼ÊøÌõ¼þ£¬²ÎÊýÖµµÄϽç
ub=[0.1,1];%Ô¼ÊøÌõ¼þ£¬²ÎÊýÖµµÄϽç
par1guess=[0.01 0.1];%²ÎÊýÖµµÄ³õʼ
options=optimset('Display','final','MaxIter',2000,'MaxFunEvals',2000);%Çó×îÓŽâµÄÓÅ»¯Ìõ¼þ
[par1,fval]=fmincon(LSmin,par1guess,[],[],[],[],lb,ub,[],options,ID,X0);

%ÀûÓÃÓÅ»¯Ëã·¨µÃµ½µÄ²ÎÊýÄâºÏÊý¾Ý
[T1,X1]=ode45(SImodel,[0 14],X0,[],par1);%ÀûÓùÀ¼Æ²ÎÊýÖØÐÂÔËÐÐÄ£ÐÍ
figure(1)
plot(T1,X1(:,1),'k-');
hold on
plot(T1,X1(:,2),'k-',0:14,ID,'ko')

%¶¨Òå×îС¶þ³ËµÄÄ¿±êº¯Êý
function d=LSmin(par1,ID,X0)
[T,x]=ode45(SImodel,[0:1:14],X0,[],par1)
d=norm(x(:,2)'-ID)^2;

%¶¨Òå΢·Ö·½³Ì
function dy=SImodel(t,y,par1)
dy=zeros(2,1);
beta=par1(1);
gamma=par1(2);
dy=[-beta*y(1)*y(2);beta*y(1)*y(2)-gamma*y(2)];
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

Çà´º¾ÍÊÇÓÃÀ´ÕÛÌڵģ¡
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ÄÚÈÝÒÑɾ³ý
MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
2Â¥2014-08-16 16:59:14
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ jwj216 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] һ־Ը̫ԭÀí¹¤°²È«¹¤³Ì300·Ö£¬Çóµ÷¼Á +3 0857Çóµ÷¼Á. 2026-03-24 3/150 2026-03-27 01:18 by ÇóÖªÈô¿Êlz
[¿¼ÑÐ] Çóµ÷¼Á +5 «lty 2026-03-25 6/300 2026-03-26 20:49 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 085601Çóµ÷¼Á×Ü·Ö293Ó¢Ò»Êý¶þ +4 ¸ÖÌú´óÅÚ 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
[¿¼ÑÐ] 334·Ö Ò»Ö¾Ô¸ÎäÀí ²ÄÁÏÇóµ÷¼Á +4 ÀîÀî²»·þÊä 2026-03-26 4/200 2026-03-26 16:00 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 291 Çóµ÷¼Á +7 »¯¹¤2026½ì±ÏÒµÉ 2026-03-21 8/400 2026-03-26 11:25 by AlenQIN.
[¿¼ÑÐ] 085600 ²ÄÁÏÓ뻯¹¤ 329·ÖÇóµ÷¼Á +9 Mr. Z 2026-03-25 9/450 2026-03-26 10:36 by baoball
[¿¼ÑÐ] 299Çóµ÷¼Á +4 15188958825 2026-03-25 4/200 2026-03-25 22:56 by 418490947
[¿¼ÑÐ] 332Çóµ÷¼Á +6 032500 2026-03-25 6/300 2026-03-25 22:45 by 418490947
[¿¼ÑÐ] 321Çóµ÷¼Á +3 è±Óñ~~ 2026-03-25 3/150 2026-03-25 19:07 by Zhanglab-TJU
[¿¼ÑÐ] ¡¾2026¿¼Ñе÷¼Á¡¿ÖÆÒ©¹¤³Ì 284·Ö ÇóÏà¹Ø×¨Òµµ÷¼ÁÃû¶î +4 Ô¬Û¼Û¼ 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[¿¼ÑÐ] 0854È˹¤ÖÇÄÜ·½ÏòÕÐÊÕµ÷¼Á +4 ÕÂСÓã567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[¿¼ÑÐ] ʳƷר˶ һ־Ը˫һÁ÷ 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[¿¼ÑÐ] Ò»Ö¾Ô¸±±»¯315 Çóµ÷¼Á +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by ÁËÁËÁËÁË¡£¡£
[¿¼ÑÐ] ²ÄÁÏ¿¼Ñе÷¼ÁÉú +3 »ÆÁ»Ò»ÃÎǧÄê 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[¿¼ÑÐ] 291Çóµ÷¼Á +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[¿¼ÑÐ] 276Çóµ÷¼Á +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by ×íÔÚ·çÀï
[¿¼ÑÐ] 260Çóµ÷¼Á +3 ÖìÜÆÁÕ 2026-03-20 4/200 2026-03-22 15:12 by ÖìÜÆÁÕ
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄÏ´ó£¬0703»¯Ñ§£¬·ÖÊý336£¬Çóµ÷¼Á +3 ÊÕµ½VS 2026-03-21 3/150 2026-03-21 18:42 by ѧԱ8dgXkO
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +4 ÄÝÄÝninicgb 2026-03-21 4/200 2026-03-21 18:39 by ѧԱ8dgXkO
[¿¼ÑÐ] Çóµ÷¼Á +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û