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

ssmile37

Ìú³æ (СÓÐÃûÆø)

[ÇóÖú] matlab ÔËÐгö´í£¬Çó°ïæÀ²~~¡¤

ÓÃ'Levenberg-Marquardt'¼ÆË㣬ÔËÐкóÌáʾOptimization terminated: no further progress can be made.Cannot generate point that reduces the sum of squares.
Problem may be ill-conditioned.
²»ÖªÊǺÎÔ­Òò£¬ÏְѳÌÐò·¢ÉÏÀ´£¬Çó´ó¼Ò°ïÖú£¡£¡£¡£¡



function sh
options = optimset('Algorithm','MaxFunEvals',5000);
[x,RESNORM,RESIDUAL,EXITFLAG]=lsqnonlin(@sh1,[1,1,1],[],[],options);
function f=sh1(beta)
data=[286 0.0646 0.999354;297 0.1079 0.998921;303 0.1330 0.998670;323 0.2508 0.997492;353 0.8079 0.991921;369 1.5109 0.984891];
T=data(:,1)';
x1=data(:,2)'*0.01;
x2=data(:,3)';
new1=84.6;
new2=18.0;
date1=24.77;
date2=48.15;
for i=1:6
r1(i)=1/x1(i)*exp(103.3*1000/8.314*(1/459.6-1/T(i)));
l12(i)=beta(1)+beta(2)*T(i)+beta(3)*T(i)^2;
fai2(i)=new2*x2(i)/(new1*x1(i)+new2*x2(i));
r1(i)=new1*((date1-date2)^2+2*l12(i)*date1*date2)*fai2(i)^2;
x1c(i)=1/r1(i)*103.3*1000/8.314*(1/459.6-1/T(i));
end
f=sum((x1-x1c).^2);
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

Ö»ÒªÏàÐÅ£¬×Ü»áÓÐÆæ¼£
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ssmile37

Ìú³æ (СÓÐÃûÆø)

ÎÊʲôûÓÐÈË»ØÓ¦
Ö»ÒªÏàÐÅ£¬×Ü»áÓÐÆæ¼£
2Â¥2012-06-20 13:50:09
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

csgt0

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

²ÊÉ«¹Òͼ

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

¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ºØÒÇ: ½ð±Ò+1, лл½¨Ò飡 2012-06-22 22:32:24
ssmile37: ½ð±Ò+3, ¡ïÓаïÖú 2012-06-26 21:41:57
ÄãÏÈ¿´¿´ÄãµÄ³ÌÐòûÎÊÌâÂð£¿r1¼ÆËãÁ½´Î£¬³£Á¿ÏÈËã³öÀ´£¬¿´ÆðÀ´¼òµ¥µã£¬±È½ÏºÃ·ÖÎö¡£
´ËÍâʹÓÃf=x1-x1c¾Í¿ÉÒÔÁË£¬lsqnonlin×Ô¼º¿ÉÒÔÇ󯽷½ºÍ£¬²»È»Ì«´ó£¬Ð¡µÄbeta¸Ä±äÒýÆð´óµÄf¸Ä±ä¾Í»áÔì³É²¡Ì¬ÎÊÌâ¡£
ÔÙÓУ¬ÄãµÄ³õÖµÒ²¿ÉÄܾÍÊÇÒ»¸ö¼«Ð¡Öµµã
showmethemoney
3Â¥2012-06-21 09:52:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ssmile37

Ìú³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
3Â¥: Originally posted by csgt0 at 2012-06-21 09:52:51
ÄãÏÈ¿´¿´ÄãµÄ³ÌÐòûÎÊÌâÂð£¿r1¼ÆËãÁ½´Î£¬³£Á¿ÏÈËã³öÀ´£¬¿´ÆðÀ´¼òµ¥µã£¬±È½ÏºÃ·ÖÎö¡£
´ËÍâʹÓÃf=x1-x1c¾Í¿ÉÒÔÁË£¬lsqnonlin×Ô¼º¿ÉÒÔÇ󯽷½ºÍ£¬²»È»Ì«´ó£¬Ð¡µÄbeta¸Ä±äÒýÆð´óµÄf¸Ä±ä¾Í»áÔì³É²¡Ì¬ÎÊÌâ¡£
ÔÙÓУ¬ÄãµÄ³õ ...

лл~ÎÒÒѾ­Ð޸ĺÃÁË£¬ÔËÐнá¹û³öÀ´ºóÊÇOptimization terminated: first-order optimality less than OPTIONS.TolFun,
and no negative/zero curvature detected in trust region model.²»ÖªµÀÕâ¸öÊDz»ÊÇÓÐÎÊÌâµÄ
Ö»ÒªÏàÐÅ£¬×Ü»áÓÐÆæ¼£
4Â¥2012-06-23 10:38:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ssmile37 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] »¯Ñ§357·Ö£¬¿¼Ñе÷¼Á +9 .Starry. 2026-04-04 10/500 2026-04-05 09:46 by MinkyXMQ
[¿¼ÑÐ] Çóµ÷¼Á +9 111623 2026-04-04 9/450 2026-04-05 09:34 by ÈçÃÎÎÞºÛ1
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +9 Ò¼·¡·¡ÒÚ 2026-04-04 9/450 2026-04-05 08:00 by qlm5820
[¿¼ÑÐ] 306·Ö²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +6 Àè°ÉÀ²À²ÄãºÜÓÐà 2026-04-03 6/300 2026-04-04 23:59 by ¹û¶³´óÍõ
[¿¼ÑÐ] ²ÄÁÏר˶322·Ö +11 ¹þ¹þ¹þºðºðºð¹þ 2026-04-02 11/550 2026-04-04 23:37 by ÓÀ×ÖºÅ
[¿¼ÑÐ] 298·Ö 070300Çóµ÷¼Á +13 zwen03 2026-04-02 13/650 2026-04-04 23:13 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 278Çóµ÷¼Á +14 ·¶æÃÄÈ 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[¿¼ÑÐ] 338Çóµ÷¼Á +7 êɹ¦? 2026-04-03 7/350 2026-04-04 20:37 by À¶ÔÆË¼Óê
[¿¼ÑÐ] ÇóÉúÎïѧµ÷¼Á +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[¿¼ÑÐ] Ò»Ö¾Ô¸°²»Õij211 0703»¯Ñ§×Ü·Ö339Çóµ÷¼Á +6 Íí·ç²»Íí 2026-04-04 6/300 2026-04-04 20:11 by dongzh2009
[¿¼ÑÐ] 332Çóµ÷¼Á +10 Lyy930824@ 2026-03-29 10/500 2026-04-04 18:41 by macy2011
[¿¼ÑÐ] »¯¹¤Çóµ÷¼Á +11 ÀóÏã֥ʿҬÄÌ 2026-04-03 11/550 2026-04-03 22:06 by à£à£à£0119
[¿¼ÑÐ] »¯¹¤µ÷¼Á303·Ö£¬¹ýËļ¶ +28 ÆÜÎà´ý·ç 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[¿¼ÑÐ] ÊÞÒ½µ÷¼Á +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[¿¼ÑÐ] ũѧ¿¼ÑÐÇóµ÷¼Á +3 dkdkxm 2026-04-01 3/150 2026-04-02 16:04 by wangjagri
[¿¼ÑÐ] 337Çóµ÷¼Á +11 ¡¶Ê÷¡· 2026-03-29 11/550 2026-04-02 10:20 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 298Çóµ÷¼Á +4 ʲôÊÇÅÖÍ·Óã 2026-03-30 6/300 2026-04-01 22:06 by ¿Í¶ûÃÀµÂ
[¿¼ÑÐ] °²È«¹¤³Ì 285 Çóµ÷¼Á +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by ¾²¾²¾²¾²¾²¾²¾²¾
[¿¼ÑÐ] Çóµ÷¼Á£¬Ò»Ö¾Ô¸±±ÁÖʳƷÓëÓªÑø095500£¬301·Ö£¬ÒѹýÁù¼¶£¬ÓпÆÑо­Àú +4 ¿ìÀÖ´¢Ðî¹Þ 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á£ºÒ»Ö¾Ô¸£ºÄϾ©´óѧ רҵ£º0705 ×Ü·Ö320 £¬±¾¿Æ985£¬ËÄÁù¼¶Òѹý +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û