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

Ê·µÏ×Ðand

гæ (³õÈëÎÄ̳)

[ÇóÖú] ÇúÏßÄâºÏ·½·¨ ÒÑÓÐ1È˲ÎÓë

´ó¼ÒºÃ£¬ÎÒ×î½üÔÚ×öÒ»¸ö·ÇÏßÐÔÄâºÏµÄ¿ÎÌ⣬²ÎÊýÖ»ÒªÒ»±ä¶à¾ÍÄâºÏ²»³öÀ´£¬ÎÒÖªµÀÈ·¶¨µÄ¹«Ê½f=x.^2./(a(1)+sqrt(a(1).^2-(a(2)+1)*x.^2))+a(3)*x.^4+a(4)*x.^6+a(5)*x.^8+a(6)*x.^10+a(7)*x.^12+a(8)*x.^14+a(9)*x.^16+a(10)*x.^18+a(11)*x.^20);ÄâºÏ³öa£¬×Ô¼ºÐ´µÄ³ÌÐò

[x,y]=textread('·ÇÇòÃæÊý¾Ý.txt','%f %f');
plot(x,y,'o')
hold on
s = fitoptions('Method','NonlinearLeastSquares','Algorithm', 'levenberg-marquardt','Lower',[-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1], 'Upper',[Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf]);
f = fittype('x.^2./(a(1)+sqrt(a(1).^2-(a(2)+1)*x.^2))+a(3)*x.^4+a(4)*x.^6+a(5)*x.^8+a(6)*x.^10+a(7)*x.^12+a(8)*x.^14+a(9)*x.^16+a(10)*x.^18+a(11)*x.^20','options',s);
[c,gof] = fit(x,y,f,'StartPoint',[4,0,0,0,0,0,0,0,0,0,0]);

ÕæÐÄÏ£ÍûÄܹ»µÃµ½ÄãÃǵİïÖú£¬×Ô¼º±àдµÄLMËã·¨£¬½á¹ûµ÷ÊÔûÎÊÌ⣬һÔËÐоͱÀÀ£ÁË£¬ÏÖÔÚ´øÓÃËû×Ô¼ºµÄº¯Êý£¬ÓÖ×ö²»³öÀ´£¬²»ÖªµÀÕâÑùµÄ·½³ÌÓÃÔõÑùµÄÄâºÏ·½·¨Ð§¹û×îºÃÄØ£¿matlabÀïÃæÓкöàÖÖÄâºÏËã·¨Âð£¬×Ô¼º¿´ÁË¿´»ù±¾¶¼ÊǾֲ¿×îÓŽ⣬ÓÐûÓÐÈ«¾Ö×îÓÅÄØ£¿Ï£ÍûÓо­ÑéµÄÈËÄܸøÒÔÖ¸µ¼¡£1stoptÈç¹ûÄâºÏµÄ»°ÄÜÉú³ÉexeÎļþÂð£¬ÎÒ×îÖÕÏë×ö³ÉÒ»¸öͨÓõijÌÐò£¬¶ø²»ÊÇÖ»ÒªÒ»´Î½á¹û¡£ÏÈлл¸÷λÁË
»Ø¸´´ËÂ¥

» ±¾Ìû¸½¼þ×ÊÔ´Áбí

  • »¶Ó­¼à¶½ºÍ·´À¡£ºÐ¡Ä¾³æ½öÌṩ½»Á÷ƽ̨£¬²»¶Ô¸ÃÄÚÈݸºÔð¡£
    ±¾ÄÚÈÝÓÉÓû§×ÔÖ÷·¢²¼£¬Èç¹ûÆäÄÚÈÝÉæ¼°µ½ÖªÊ¶²úȨÎÊÌ⣬ÆäÔðÈÎÔÚÓÚÓû§±¾ÈË£¬Èç¶Ô°æÈ¨ÓÐÒìÒ飬ÇëÁªÏµÓÊÏ䣺xiaomuchong@tal.com
  • ¸½¼þ 1 : ·ÇÇòÃæÊý¾Ý.txt
  • 2014-12-12 17:29:24, 13.58 K

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

¶¯Á¦Ñ§ÄâºÏ

» ²ÂÄãϲ»¶

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

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

Ê·µÏ×Ðand

гæ (³õÈëÎÄ̳)

ÒýÓûØÌû:
2Â¥: Originally posted by feixiaolin at 2014-12-12 20:36:34
¿¼ÂÇÕý½»¶àÏîʽ£¬»òÇбÈÑ©·òÄâºÏʲôµÄ£¬ÄãµÄ½×ÊýÌ«¸ßÁË¡£

ºÃ£¬ÎÒ¿´¿´£¬matlabÖÐÔõÑùʹÄâºÏµÄ½á¹ûֻȡʵÊý½âÄØ£¬ÎÒÄâºÏ³öÀ´½á¹ûºÃ¶àÐé²½°¡
3Â¥2014-12-12 22:28:23
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

feixiaolin

ÈÙÓþ°æÖ÷ (ÎÄ̳¾«Ó¢)

ÓÅÐã°æÖ÷

¿¼ÂÇÕý½»¶àÏîʽ£¬»òÇбÈÑ©·òÄâºÏʲôµÄ£¬ÄãµÄ½×ÊýÌ«¸ßÁË¡£
2Â¥2014-12-12 20:36:34
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

dbb627

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

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
È«¾ÖËã·¨
CODE:
function Insphere_fit
[xdata,ydata]=textread('·ÇÇòÃæÊý¾Ý.txt','%f %f');
plot(xdata,ydata,'o')
hold on
myfun=@(x,xdata)xdata.^2./(x(1)+sqrt(x(1).^2-(x(2)+1)*xdata.^2))+x(3)*xdata.^4+x(4)*xdata.^6+x(5)*xdata.^8+x(6)*xdata.^10+x(7)*xdata.^12+x(8)*xdata.^14+x(9)*xdata.^16+x(10)*xdata.^18+x(11)*xdata.^20;
tic
x0=[4 rand(1,10)];
format shortg
ms=MultiStart;
F =@(x)norm(xdata.^2./(x(1)+sqrt(x(1).^2-(x(2)+1)*xdata.^2))+x(3)*xdata.^4+x(4)*xdata.^6+x(5)*xdata.^8+...
    x(6)*xdata.^10+x(7)*xdata.^12+x(8)*xdata.^14+x(9)*xdata.^16+x(10)*xdata.^18+x(11)*xdata.^20-ydata);
opts=optimset('Algorithm', 'interior-point');
problem=createOptimProblem('fmincon','x0',x0,...
    'objective',F,'lb',[-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1],...
    'ub',[Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf],'options',opts);
[xminm,fminm,flagm,outptm,manyminsm]=run(ms,problem,20)
t2=toc
figure;plot(xdata,ydata,'r-*',xdata,myfun(xminm,xdata),'bo-')
legend('real','model')
Text2=['ÀûÓÃÈ«¾ÖËã·¨ multistartºÍfmincon¹À¼ÆµÃµ½µÄ²ÎÊý,ºÄʱ',num2str(t2),'s'];
title(Text2)

MultiStart completed the runs from all start points.

All 20 local solver runs converged with a positive local solver exit flag.

xminm =

  Columns 1 through 5

       132.11       1497.3       12.388     -0.99669     -0.98238

  Columns 6 through 10

     -0.99812      0.68971     -0.16267     0.017723  -0.00074323

  Column 11

-4.2727e-006


fminm =

       788.62


flagm =

     1


outptm =

                funcCount: 15871
         localSolverTotal: 20
       localSolverSuccess: 20
    localSolverIncomplete: 0
    localSolverNoSolution: 0
                  message: [1x128 char]


manyminsm =

  1x20 GlobalOptimSolution

  Properties:
    X
    Fval
    Exitflag
    Output
    X0




t2 =

       59.298
ÇúÏßÄâºÏ·½·¨
untitled.jpg

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
4Â¥2014-12-14 23:52:47
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[Óлú½»Á÷] ¶þ±½¼×ͪËáÀàÑÜÉúÎï 50+3 С°×°®Ö÷ÈË 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[¿¼ÑÐ] 320Çóµ÷¼Á +4 Éî½¼akm 2026-04-17 4/200 2026-04-17 17:32 by Ò»ÇÐOK
[¿¼ÑÐ] 297£¬¹¤¿Æµ÷¼Á? +4 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-14 4/200 2026-04-16 22:52 by wulijun2012
[¿¼ÑÐ] »¯Ñ§070300 Çóµ÷¼Á +28 ¹þ¹þ¹þ^_^ 2026-04-12 28/1400 2026-04-16 21:36 by ´óÁ¦Ë®ÊÖÁ¦´óÎÞÇ
[¿¼ÑÐ] 291Çóµ÷¼Á +11 ¹ØÒä±±. 2026-04-14 11/550 2026-04-16 15:18 by jiahl2024
[¿¼ÑÐ] 26ҩѧר˶105500Çóµ÷¼Á +6 à¶¹þ¼ÓÓÍ 2026-04-13 7/350 2026-04-16 14:31 by zhouxiaoyu
[»ù½ðÉêÇë] RY£ºÖйú²ú³öµÄ¿ÆÑ§À¬»øÂÛÎÄ£¬¾ø¶ÔÊýÁ¿ºÍ±ÈÀý¶¼ÊÀ½çµÚÒ» +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by »¶ÀÖËÌÒ¶Ýè
[¿¼ÑÐ] 322Çóµ÷¼Á +8 123°²¿µ 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[¿¼ÑÐ] 290µ÷¼ÁÉúÎï0860 +38 ÍÛ¹þ¹þ£¬¡£ 2026-04-11 44/2200 2026-04-16 09:52 by cuisz
[¿¼ÑÐ] 279ѧ˶ʳƷרҵÇóµ÷¼ÁԺУ 20+7 ¹Â¶ÀµÄÀǰ®³ÔÑò 2026-04-12 29/1450 2026-04-16 09:00 by screening
[¿¼ÑÐ] 289 ·Ö105500ҩѧר˶Çóµ÷¼Á(ÕÒBÇøÑ§Ð£) +4 °×ÔÆ123456789 2026-04-13 4/200 2026-04-16 00:18 by ·ÛÇßÈô³¾
[¿¼ÑÐ] 310Çóµ÷¼Á +16 666ÕæºÃ 2026-04-11 18/900 2026-04-15 13:28 by ºÚ¿Æ¼¼¿óÒµ
[¿¼ÑÐ] ²ÄÁϹ¤³Ì281»¹Óе÷¼Á»ú»áÂð +43 xaw. 2026-04-11 44/2200 2026-04-15 12:46 by Î÷±±Íû¡ª·çɳ
[¿¼ÑÐ] 211±¾¿Æ²ÄÁÏ»¯¹¤Çóµ÷¼Á +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[¿¼ÑÐ] 305Çóµ÷¼Á +8 Â꿨°Í¿¨boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[¿¼ÑÐ] Çóµ÷¼Á£¬985²ÄÁÏÓ뻯¹¤348·Ö +9 º­ÖñÁõ 2026-04-11 14/700 2026-04-13 22:26 by º­ÖñÁõ
[¿¼ÑÐ] 2026˶ʿµ÷¼Á_Äܶ¯_ºÓÄÏũҵ´óѧ +4 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[¿¼ÑÐ] 302Çóµ÷¼Á +10 Ò×£¡? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[¿¼ÑÐ] 344 ²ÄÁÏרҵ Çóµ÷¼Á211 ÎÞµØÓòÒªÇó +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[¿¼ÑÐ] 359Çóµ÷¼Á +5 θ¾·ÂÎÀÛÁË 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û