24СʱÈÈÃŰæ¿éÅÅÐаñ    

±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1496  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 325Çóµ÷¼Á +4 Àî¼Îͼ¡¤S¡¤Â· 2026-03-23 4/200 2026-03-26 23:08 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] Çóµ÷¼Á Ò»Ö¾Ô¸ ±¾¿Æ ±±¿Æ´ó »¯Ñ§ 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 081200-11408-276ѧ˶Çóµ÷¼Á +3 ´Þwj 2026-03-26 3/150 2026-03-26 20:33 by cy128
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©»¯¹¤´óѧ²ÄÁÏÓ뻯¹¤£¨085600£©296Çóµ÷¼Á +9 µ¾ÆÞС±à 2026-03-26 9/450 2026-03-26 16:16 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 0856µ÷¼Á +4 ÇóÇóÈÃÎÒÓÐÊé¶Á° 2026-03-26 5/250 2026-03-26 15:54 by dick_runner
[¿¼ÑÐ] 296Çóµ÷¼Á +5 www_q 2026-03-20 5/250 2026-03-26 12:56 by 3Strings
[¿¼ÑÐ] Ò»Ö¾Ô¸ÉϺ£½»´óÉúÎïÓëҽҩר˶324·Ö£¬Çóµ÷¼Á +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[¿¼ÑÐ] 321Çóµ÷¼Á +3 è±Óñ~~ 2026-03-25 3/150 2026-03-25 19:07 by Zhanglab-TJU
[¿¼ÑÐ] 26¿¼ÑÐ-291·Ö-ÏÃÃÅ´óѧ£¨085601£©-ÈáÐÔµç×ÓѧԺ²ÄÁϹ¤³ÌרҵÇóµ÷¼Á +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[¿¼ÑÐ] ¸÷λÀÏʦÄúºÃ£º±¾È˳õÊÔ372·Ö +5 jjÓ¿77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[¿¼ÑÐ] 0854µç×ÓÐÅÏ¢Çóµ÷¼Á 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[¿¼ÑÐ] 086003ʳƷ¹¤³ÌÇóµ÷¼Á +6 íµíµ111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[¿¼ÑÐ] ʳƷר˶ һ־Ը˫һÁ÷ 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[¿¼ÑÐ] ²ÄÁÏר˶ÕÒµ÷¼Á +5 ¹þ¹þ¹þºðºðºð¹þ 2026-03-23 5/250 2026-03-24 19:07 by ÁËÁËÁËÁË¡£¡£
[¿¼ÑÐ] ²ÄÁÏ292µ÷¼Á +8 éÙËÌ˼ÃÀÈË 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[¿¼ÑÐ] 315·Ö£¬³ÏÇóµ÷¼Á£¬²ÄÁÏÓ뻯¹¤085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[¿¼ÑÐ] Çóµ÷¼ÁÒ»Ö¾Ô¸º£´ó£¬0703»¯Ñ§Ñ§Ë¶304·Ö£¬Óдó´´ÏîÄ¿£¬Ëļ¶Òѹý +6 ÐÒÔËÁ¨Á¨ 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[¿¼ÑÐ] 311Çóµ÷¼Á +3 26ÑÐ0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +4 ÄÝÄÝninicgb 2026-03-21 4/200 2026-03-21 18:39 by ѧԱ8dgXkO
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³×Ü·Ö308Çóµ÷¼Á +7 īīĮ 2026-03-20 7/350 2026-03-21 16:36 by barlinike
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û