²é¿´: 629  |  »Ø¸´: 3

sara_ecust

гæ (СÓÐÃûÆø)

[ÇóÖú] matlabÄâºÏ²ÎÊý ÒÑÓÐ1È˲ÎÓë

ÒÑÖªPA0 T x t
dx/dt=k*(PA0*(1-x))^n
¸ù¾ÝÒÑÖªµÄÊýÄâºÏk ºÍn
³ÌÐòÔÚ¸½¼þÖУ¬¾ÍÒ»Ö±³öÏÖ
¡®too many input arguments¡¯
Çó½â¾ö
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

ÔÂÖ»À¶

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

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
sara_ecust: ½ð±Ò+5, ¡ïÓаïÖú 2014-07-14 16:28:37
sara_ecust: ½ð±Ò+5, ¡ïÓаïÖú 2014-12-25 09:18:49
µ¥¸ö´íÎóÌáʾ¿ÉÄÜ»á¶ÔÓ¦¶à¸öÔ­Òò¡£
Åųý³ÌÐòÎÊÌ⣬»¹Êǵÿ´Ô­³ÌÐò¡£
¸ø³öÔ­³ÌÐòÀ´¿´¿´¡£
MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
2Â¥2014-07-06 18:17:41
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

sara_ecust

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by ÔÂÖ»À¶ at 2014-07-06 18:17:41
µ¥¸ö´íÎóÌáʾ¿ÉÄÜ»á¶ÔÓ¦¶à¸öÔ­Òò¡£
Åųý³ÌÐòÎÊÌ⣬»¹Êǵÿ´Ô­³ÌÐò¡£
¸ø³öÔ­³ÌÐòÀ´¿´¿´¡£

function biaoguan
clear all;clc
format long;
%global T PA0 tspan x x0
T=[613 633 663];
F0=[15 17.5 20 22.5 25];
%t=0.03/0.78*134/0.004./F0;
t=[85.8974 73.6264 64.4231 57.2650 51.5385];

x1=[0.632 0.581 0.540 0.496 0.468];
x2=[0.600 0.551 0.515 0.480 0.451];
x3=[0.551 0.492 0.459 0.427 0.402];
x=[x1;x2;x3];
x0=[0.632 0.600 0.551];
PA0=[0.001865374 0.002289233 0.00239448];

beta0=[1.1e5 100000 1];
lb=[0 0 0];
ub=[inf inf inf];
tspan =t;
[beta,resnorm,resid,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@Opt,beta0,[],[],[],tspan,x0,T,PA0,x)
ci = nlparci(beta,resid,jacobian);
beta

function f = Opt(beta,tspan,x0,T,PA0,x)
[t xc1] = ode45(@Kin,tspan,x0(1),[],beta,T(1),PA0(1),x(1,);
[t xc2] = ode45(@Kin,tspan,x0(2),[],beta,T(2),PA0(2),x(2,);
[t xc3] = ode45(@Kin,tspan,x0(3),[],beta,T(3),PA0(3),x(3,);

f1= xc1' - x(1,;
f2= xc2' - x(2,;
f3= xc3' - x(3,;
f=[f1;f2;f3]

% ------------------------------------------------------------------
function dxdt = Kin(t,beta,T,PA0,x)
dxdt=beta(1).*exp(-beta(2)./(8.3145*T)).*(PA0.*(1-x)).^beta(3);

» ±¾ÌûÒÑ»ñµÃµÄºì»¨£¨×îÐÂ10¶ä£©

3Â¥2014-07-06 20:30:46
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

sara_ecust

гæ (СÓÐÃûÆø)

Ëͺ컨һ¶ä
ÒýÓûØÌû:
3Â¥: Originally posted by sara_ecust at 2014-07-06 20:30:46
function biaoguan
clear all;clc
format long;
%global T PA0 tspan x x0
T=;
F0=;
%t=0.03/0.78*134/0.004./F0;
t=;

x1=;
x2=;
x3=;
x=;
x0=;
PA0=;

beta0=;
lb=;
ub=;
tspan =t;
= .. ...

Õâ¸öЦÁ³ÊÇ¡°£º¡±
4Â¥2014-07-06 20:31:52
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ sara_ecust µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ×Ü·Ö322ÇóÉúÎïѧ/Éú»¯Óë·Ö×Ó/ÉúÎïÐÅϢѧÏà¹Øµ÷¼Á +5 ÐdzÁuu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[¿¼ÑÐ] µ÷¼ÁÍÆ¼ö +4 Ç徯714 2026-03-26 5/250 2026-03-26 17:49 by fmesaito
[¿¼ÑÐ] 085601Çóµ÷¼Á×Ü·Ö293Ó¢Ò»Êý¶þ +4 ¸ÖÌú´óÅÚ 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
[¿¼ÑÐ] 07»¯Ñ§303Çóµ÷¼Á +5 î£08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[¿¼ÑÐ] ²ÄÁÏѧ˶333Çóµ÷¼Á +4 ±±µÀÏï 2026-03-24 4/200 2026-03-25 14:16 by mapenggao
[¿¼ÑÐ] 293Çóµ÷¼Á +7 ¼ÓÒ»Ò»¾Å 2026-03-24 7/350 2026-03-25 12:02 by userper
[¿¼ÑÐ] Çóµ÷¼Á Ò»Ö¾Ô¸ ±¾¿Æ ±±¿Æ´ó »¯Ñ§ 343 +4 13831862839 2026-03-24 5/250 2026-03-25 09:47 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 07»¯Ñ§280·ÖÇóµ÷¼Á +7 722865 2026-03-23 7/350 2026-03-25 09:29 by aa331100
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á£¬Çóµ¼Ê¦ÊÕ +7 ÌìÌìºÃÔËÀ´Éϰ¶° 2026-03-24 7/350 2026-03-24 20:26 by peike
[¿¼ÑÐ] ²ÄÁÏר˶ÕÒµ÷¼Á +5 ¹þ¹þ¹þºðºðºð¹þ 2026-03-23 5/250 2026-03-24 19:07 by ÁËÁËÁËÁË¡£¡£
[¿¼ÑÐ] 307Çóµ÷¼Á +3 ÓàÒâÇä 2026-03-21 6/300 2026-03-24 15:03 by ÓàÒâÇä
[¿¼ÑÐ] Ò»Ö¾Ô¸¼ª´ó»¯Ñ§322Çóµ÷¼Á +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by ´÷Χ²±µÄСÎÃ×Ó
[¿¼ÑÐ] 298Çóµ÷¼Á +8 Éϰ¶6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[¿¼ÑÐ] 291Çóµ÷¼Á +5 ‹üÈA 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ»¯Ñ§£¨0703£©×Ü·Ö337Çóµ÷¼Á +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[¿¼ÑÐ] 285Çóµ÷¼Á +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] ³õÊÔ 317 +7 °ëÀ­Ô±û 2026-03-20 7/350 2026-03-21 22:26 by peike
[¿¼ÑÐ] ²ÄÁÏÇóµ÷¼Á +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[¿¼ÑÐ] ¿¼Ñе÷¼Á +3 ºôºô£¿~+123456 2026-03-21 3/150 2026-03-21 20:04 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 297Çóµ÷¼Á +3 ϲ»¶»¹ÊDz»¸ÊÐÄ 2026-03-20 3/150 2026-03-21 18:33 by ѧԱ8dgXkO
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û