Znn3bq.jpeg
²é¿´: 1655  |  »Ø¸´: 3

wtisunday

ľ³æ (ÕýʽдÊÖ)

[ÇóÖú] ÈçºÎÔÚmatlabÖÐʵÏÖÏòÁ¿Ïòº¯ÊýµÄ´«µÝ

for i=1:8
        a=k(i);b=k(i+8);
        [x(i),fval]=fsolve(@pp,x0);
end

function f=pp(x)
pp1=x(1)+a;
pp2=x(2)-b;
f=[pp1,pp2]';

ÏëÎÊһϣ¬ÈçºÎʵÏÖÏòÁ¿·ÖÁ¿Ïòº¯ÊýµÄ´«µÝ°¡£¿Ð»Ð»
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

dbb627

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

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
¾Ù¸öÀý×Ó°É
CODE:
function aa
k=1:16;
for i=1:8
    x0=[2 1];
        a=k(i);b=k(i+8);
        [x(i,:),fval]=fsolve(@(x)pp(x,a,b),x0);
        end
x
function f=pp(x,a,b)
pp1=x(1)+a;
pp2=x(2)-b;
f=[pp1,pp2]';

½á¹û

x =

   -1.0000    9.0000
   -2.0000   10.0000
   -3.0000   11.0000
   -4.0000   12.0000
   -5.0000   13.0000
   -6.0000   14.0000
   -7.0000   15.0000
   -8.0000   16.0000
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.
2Â¥2012-02-16 09:50:37
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wtisunday

ľ³æ (ÕýʽдÊÖ)

ÒýÓûØÌû:
Â¥: Originally posted by dbb627 at 2012-02-16 09:50:37:
¾Ù¸öÀý×Ó°É
[code]
function aa
k=1:16;
for i=1:8
    x0=[2 1];
        a=k(i);b=k(i+8);
        [x(i,,fval]=fsolve(@(x)pp(x,a,b),x0);
        end
x
function f=pp(x,a,b)
pp1=x(1)+a;
pp2 ...

Èç¹ûkÊÇÒ»¸öÈ«¾Ö±äÁ¿£¬±ÈÈçk=[a1,a2,a3...,a16]£¬ÓÖÈçºÎʵÏÖÔÚiÑ­»·ÏÂʵÏÖ¶ÔkµÄÑ¡Ôñ°¡
3Â¥2012-02-16 12:48:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wtisunday

ľ³æ (ÕýʽдÊÖ)

ÒýÓûØÌû:
Â¥: Originally posted by dbb627 at 2012-02-16 09:50:37:
¾Ù¸öÀý×Ó°É
[code]
function aa
k=1:16;
for i=1:8
    x0=[2 1];
        a=k(i);b=k(i+8);
        [x(i,,fval]=fsolve(@(x)pp(x,a,b),x0);
        end
x
function f=pp(x,a,b)
pp1=x(1)+a;
pp2 ...

function f=ObjFunc1(x0,k,xexp,yexp)
global c
c=k;
G=[336370,64246.4,44849.7,10649.7,5151.99,1224.43,63.4524,2.67362]';
t=[0.00423,0.01621,0.06216,0.23841,0.91437,3.50692,14.6856,101.343]';
xexp=[0.03162 0.1 0.31623 1 3.16228 10.5635 20.395]';
yexp=[84834.3 84834.3 79215.4 62322.2 52509.5 32501.6 25570.4]';
for j=1:7
    Xexp=xexp(j);Yexp=yexp(j);
    for i=1:8
        x0=[1,0]';
        T=t(i);g=G(i);a=c(i);b=c(i+8);
        [x(i,,fval]=fsolve(@(x)pp(x,a,b,Xexp,T,g),x0);
        fff(i)=(x(:,1)-x(:,2))/Xexp;
    end
    ff(j)=sum(fff')-Yexp;
end
f=sum((ff')^2);
ÒÔÉÏÊÇÎÒµÄÒ»¸ö×Ó³ÌÐò£¬k=[a1,...a16]£¬ÎªÈ«¾Ö±äÁ¿
Äú°ïæ¸ø¿´¿´Ôõô»ØÊ°ɣ¬ÎªÊ²Ã´·µ»ØÏÔʾcÖ»ÓÐ1¡Á7´óÐ¡ÄØ
4Â¥2012-02-16 13:13:43
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ wtisunday µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 070300»¯Ñ§ Çóµ÷¼Á +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤×¨Ë¶329 Çóµ÷¼Á +24 ¶îcc 2026-04-06 25/1250 2026-04-09 16:01 by wp06
[¿¼ÑÐ] 265Çóµ÷¼Á +4 ·ç˵ËýÔçÍüÁË 2026-04-07 4/200 2026-04-09 13:59 by onlyÖÜ
[¿¼ÑÐ] 0860004 Çóµ÷¼Á 309·Ö +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by °¡Àî999
[¿¼ÑÐ] 328Çóµ÷¼Á +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[¿¼ÑÐ] 296Çóµ÷¼Á +3 Íô£¡£¿£¡ 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£©×ö¹ý·Ö×ÓʵÑé +6 ÏàÐűػá¹ââÍòÕ 2026-04-07 7/350 2026-04-08 16:49 by tjzhao
[¿¼ÑÐ] 347²ÄÁÏר˶Çóµ÷¼Á +18 zj8215216 2026-04-06 18/900 2026-04-08 16:27 by luoyongfeng
[¿¼ÑÐ] Ò»Ö¾Ô¸¼ª´ó»¯Ñ§327Çóµ÷¼Á +12 ÍõÍõ°×ʯ 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[¿¼ÑÐ] 388Çóµ÷¼Á +6 ËÄ´¨ÍõÌÎ 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[¿¼ÑÐ] 313Çóµ÷¼Á +3 Ê®Áùʰ½ 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[¿¼ÑÐ] 285Çóµ÷¼Á +15 ŶßϺôo 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[¿¼ÑÐ] ²ÄÁϹ¤³Ì310ר˶µ÷¼Á +14 ÀÌÀÌÎÒ¡­. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[¿¼ÑÐ] Ò»Ö¾Ô¸ºÓ±±¹¤Òµ´óѧ²ÄÁϹ¤³Ì£¬³õÊÔ344Çóר˶µ÷¼Á +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 331Çóµ÷¼Á +8 ÓÚÕ÷yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[¿¼ÑÐ] 22408 ×Ü·Ö320£¬Ò»ÆªÂÛÎĶþ×÷£¬Á½¸ö¹úÈý£¬Çóµ÷¼Á +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[¿¼ÑÐ] 348Çóµ÷¼Á +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by Öí»á·É
[¿¼ÑÐ] 085602µ÷¼Á ³õÊÔ×Ü·Ö335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[¿¼ÑÐ] Äܶ¯µ÷¼Á326ר˶ +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[¿¼ÑÐ] 338Çóµ÷¼Á +7 êɹ¦? 2026-04-03 7/350 2026-04-04 20:37 by À¶ÔÆË¼Óê
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û