²é¿´: 1774  |  »Ø¸´: 9

zhangluyoa

Òø³æ (³õÈëÎÄ̳)

[ÇóÖú] Çó½âÀ©É¢·½³ÌÎÊÌâ

ÓÐÒ»¸ö·½³ÌÈçͼ
Ô´ÂëÊÇ
function AA
D0=3.55e-19;
options=optimset('tolx',1e-9);
[D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,[],[],options)

function F=fei(D)
C0=0.687;
d=0.02;
Ceq=20;
t=60;
C=4.136;
m=1-(C0-C)/(C0-Ceq);
s=[];
k=0:1000;
s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2)));
s
F=sum((m-s').^2);


¶øµÃµ½µÄ½á¹ûÔõôҲ²»¶Ô ÈçÏÂ
µÃµ½µÄDÖµÔõô¶¼ÊǸú¸³µÄD0ÖµÊÇÒ»ÑùµÄ¶î
s =

    0.9998


s =

    0.2158


s =

   Inf


s =

    0.2447


s =

   Inf


s =

    0.2161

Optimization terminated: norm of the current step is less
than OPTIONS.TolX.

D =

  3.5500e-019


resnorm =

    0.0010


residual =

    0.0318


exitflag =

     2


output =

    firstorderopt: 7.1546e+005
       iterations: 2
        funcCount: 6
     cgiterations: 1
        algorithm: 'large-scale: trust-region reflective Newton'
          message: [1x77 char]


lambda =

    lower: 0
    upper: 0


jacobian =

   (1,1)     2.2480e+007
µ½µ×ÔõôŪµÃ£¬Ð¡µÜÊDzËÄñ£¬Ï£ÍûÓиßÊÖ¸úºÃÐÄÈ˰ïÎÒ5555
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

ya634984094

Í­³æ (СÓÐÃûÆø)

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
³ôË®¹µ: Ó¦ÖúÖ¸Êý-1, ·ÇÓ¦ÖúÌùÇ벻ҪѡÔñ¡°Ó¦Öú»ØÌû¡±~ 2012-04-03 14:47:18
ÄãÕâÊÇÓÃuds±àµÄÊä³öÂ𣿣¿
2Â¥2012-04-03 10:54:59
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhangluyoa

Òø³æ (³õÈëÎÄ̳)

ÒýÓûØÌû:
2Â¥: Originally posted by ya634984094 at 2012-04-03 10:54:59:
ÄãÕâÊÇÓÃuds±àµÄÊä³öÂ𣿣¿

ʲôÊÇubs°¡£¬ÓõÄÊÇmatlab£¬ÎÒÊDzËÄñ
3Â¥2012-04-03 11:35:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ya634984094

Í­³æ (СÓÐÃûÆø)

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

¡ï
³ôË®¹µ: ½ð±Ò+1, лл½»Á÷~~ 2012-04-04 11:15:14
ÒýÓûØÌû:
3Â¥: Originally posted by zhangluyoa at 2012-04-03 11:35:39:
ʲôÊÇubs°¡£¬ÓõÄÊÇmatlab£¬ÎÒÊDzËÄñ

ßÀ~~matlabÎÒ²»ÊǺÜÊ죬udsÊÇfluentÓÃÀ´Çó½â·½³ÌµÄ~~
4Â¥2012-04-03 12:34:11
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhangluyoa

Òø³æ (³õÈëÎÄ̳)

ÎÒ¾ÍÊÇÏëÇóDÖµ£¬µ«ÊDz»ÖªµÀΪʲô¸ø¸³D0ÖµÊÇÒ»ÑùµÄ
5Â¥2012-04-03 14:58:46
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

dbb627

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

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
³ôË®¹µ: ½ð±Ò+3, лл½»Á÷~~ 2012-04-04 11:15:22
zhangluyoa: ½ð±Ò+4, ¡ïÓаïÖú 2012-04-11 16:35:41
CODE:
D0=3.55e-6;
options=optimset('tolx',1e-9);
[D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,[],[],options)

function F=fei(D)
C0=0.687;
d=0.02;
Ceq=20;
t=60;
C=4.136;
m=1-(C0-C)/(C0-Ceq);
s=[];
k=0:5000;
s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2[color=red]*t[/color]/(d^2)));
F=sum((m-s').^2);

Local minimum possible.

lsqnonlin stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.




D =

  5.0774e-008


resnorm =

  1.1370e-007


residual =

  3.3720e-004


exitflag =

     3


output =

    firstorderopt: 39.0251
       iterations: 13
        funcCount: 28
     cgiterations: 0
        algorithm: 'large-scale: trust-region reflective Newton'
          message: [1x458 char]


lambda =

    lower: 0
    upper: 0


jacobian =

   (1,1)     1.1573e+005
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.
6Â¥2012-04-03 17:56:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

dbb627

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

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

¡ï ¡ï ¡ï
³ôË®¹µ: ½ð±Ò+1, лл½»Á÷~~ 2012-04-04 11:15:29
zhangluyoa: ½ð±Ò+2, ¡ïÓаïÖú 2012-04-11 16:35:57
s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2))); ¶à³ËÁË60£¿
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.
7Â¥2012-04-03 17:58:17
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zhangluyoa

Òø³æ (³õÈëÎÄ̳)

ÒýÓûØÌû:
7Â¥: Originally posted by dbb627 at 2012-04-03 17:58:17:
s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2))); ¶à³ËÁË60£¿

²»ÊÇ£¬tÊÇ·ÖÖÓ£¬³ËÒÔ60±ä³ÉÃë
8Â¥2012-04-04 19:08:11
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

luyao_zhang

Òø³æ (³õÈëÎÄ̳)

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
zhangluyoa: ½ð±Ò+44, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, ÖÕÓÚ½â¾öÁË 2013-07-16 17:35:01
xiegangmai: ½ð±Ò+2, лл²ÎÓ룡 2013-07-24 21:03:08
¾«¶È̫С£¬ÒѾ­Ð¡ÓÚ¼ÆËã»úÓÅ»¯¾«¶È£¬ËùÒÔÓÅ»¯¼ÆËãÍ£Ö¹£¬°ÑD Ôö´ó10^8±¶£¬·½³ÌÖмõС10^8±¶£¬¿É˳ÀûÇó½â

D0=3.55;

[D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,[],[])

function F=fei(D)
C0=0.687;
d=0.02;
Ceq=20;
t=60;
C=4.136;
m=1-(C0-C)/(C0-Ceq);
s=[];
k=0:1000;
s=sum(8./(pi*(2*k+1)).^2.*exp(-D*0.00000001*(pi*(2*k+1)).^2*t*60/(d^2)));
s
F=sum((m-s').^2);




½á¹û
D =

   0.04933156266209


resnorm =

    1.200260299201833e-008

D¾ÍµÈÓÚ4.93e-10
ÈȰ®Í´¿à
9Â¥2013-07-16 17:34:30
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

1657956529

гæ (³õÈëÎÄ̳)

ÕâÊǽâµÄʲô·½³Ì£¬ÎÒÒ²Ïë¾ÈÖúһϣ¬fickµÚ¶þ¶¨ÂÉÔõôÇó½â·½³Ì
10Â¥2016-11-22 18:58:42
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zhangluyoa µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 070303 ×Ü·Ö349Çóµ÷¼Á +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[¿¼ÑÐ] 277²ÄÁÏ¿ÆÑ§Ó빤³Ì080500Çóµ÷¼Á +3 ×ÔÓɼå±ý¹û×Ó 2026-03-16 3/150 2026-03-16 14:10 by ÔËÆøyunqi
[¿¼ÑÐ] »·¾³¹¤³Ìµ÷¼Á +4 ´ó¿Édigkids 2026-03-16 4/200 2026-03-16 14:00 by laoshidan
[¿¼ÑÐ] 326Çóµ÷¼Á +3 ŵ±´¶û»¯Ñ§½±êéê 2026-03-15 3/150 2026-03-16 13:51 by houyaoxu
[¿¼ÑÐ] Ò»Ö¾Ô¸985£¬±¾¿Æ211£¬0817»¯Ñ§¹¤³ÌÓë¼¼Êõ319Çóµ÷¼Á +4 Liwangman 2026-03-15 4/200 2026-03-16 08:01 by wang_dand
[¿¼ÑÐ] 327Çóµ÷¼Á +6 ʰ¹âÈÎȾ 2026-03-15 11/550 2026-03-15 22:47 by ʰ¹âÈÎȾ
[¿¼²©] »¶Ó­É격ͬѧÁªÏµ +3 ÌìµÀ³êÇÚ2026686 2026-03-10 7/350 2026-03-15 19:03 by ÌìµÀ³êÇÚ2026686
[¿¼ÑÐ] 265Çóµ÷¼Á +4 Íþ»¯±ý07 2026-03-12 4/200 2026-03-14 17:23 by userper
[¿¼ÑÐ] Ò»Ö¾Ô¸¹þ¹¤´ó²ÄÁÏ324·ÖÇóµ÷¼Á +5 ãÆÐñ¶« 2026-03-14 5/250 2026-03-14 14:53 by ľ¹Ï¸à
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤£¨0856£©304ÇóBÇøµ÷¼Á +7 Çñgl 2026-03-10 11/550 2026-03-14 12:18 by Çñgl
[¿¼ÑÐ] µ÷¼Á +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[¿¼ÑÐ] 337Ò»Ö¾Ô¸»ªÄÏÀí¹¤0805²ÄÁÏÇóµ÷¼Á +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[¿¼ÑÐ] 0703£¬333·ÖÇóµ÷¼Á Ò»Ö¾Ô¸Ö£ÖÝ´óѧ-ÎïÀí»¯Ñ§ +3 Àîħٶ·Åñ 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[¿¼ÑÐ] 308Çóµ÷¼Á +5 ÊÇLupa°¡ 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[¿¼ÑÐ] £Û0860£Ý321·ÖÇóµ÷¼Á£¬abÇø½Ô¿É +4 ±¦¹óÈÈ 2026-03-13 4/200 2026-03-13 22:01 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] Çóµ÷¼Á +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[¿¼ÑÐ] 307Çóµ÷¼Á +5 ³¬¼¶ÒÁ°º´óÍõ 2026-03-12 5/250 2026-03-13 15:56 by °ô°ôÇòÊÖ
[¿¼ÑÐ] 274Çóµ÷¼Á +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[¿¼ÑÐ] 328»¯¹¤×¨Ë¶Çóµ÷¼Á +4 ¡££¬¡££¬¡££¬¡£i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[¿¼ÑÐ] 0817»¯Ñ§¹¤³ÌÓë¼¼Êõ¿¼ÑÐ312·Öµ÷¼Á +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û