±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1790  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 299Çóµ÷¼Á +8 15188958825 2026-03-25 8/400 2026-03-29 01:36 by fmesaito
[¿¼ÑÐ] Ò»Ö¾Ô¸¹þ¶û±õ¹¤Òµ´óѧ²ÄÁÏÓ뻯¹¤·½Ïò336·Ö +9 ³½ãå5211314 2026-03-26 9/450 2026-03-29 01:12 by ÎÒÊÇС¿µ
[¿¼ÑÐ] Ò»Ö¾Ô¸ÏÃÃÅ´óѧ»¯Ñ§Ñ§Ë¶307Çóµ÷¼Á +10 y7czhao 2026-03-26 10/500 2026-03-28 14:23 by ÌÆãå¶ù
[¿¼ÑÐ] 311Çóµ÷¼Á +9 lin0039 2026-03-26 9/450 2026-03-28 13:05 by ÌÆãå¶ù
[¿¼ÑÐ] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by ±£»¤µØÇòÄãÎÒ×öÆ
[¿¼ÑÐ] 291Çóµ÷¼Á +7 ‹üÈA 2026-03-22 7/350 2026-03-28 04:02 by fmesaito
[Óлú½»Á÷] ¸ßθßѹ·´Ó¦ÇóÖú 10+4 chibby 2026-03-25 4/200 2026-03-27 21:08 by BT20230424
[¿¼ÑÐ] Ò»Ö¾Ô¸211£¬335·Ö£¬0856£¬Çóµ÷¼ÁԺУºÍµ¼Ê¦ +4 Çã____Ïô 2026-03-27 5/250 2026-03-27 11:52 by zhshch
[¿¼ÑÐ] ²ÄÁÏѧ˶£¬Çóµ÷¼Á 6+5 ÌǺù«888ll 2026-03-22 10/500 2026-03-27 08:18 by hypershenger
[¿¼ÑÐ] Ò»Ö¾Ô¸±±»¯Çóµ÷¼Á +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv·çÀ×
[¿¼ÑÐ] 321Çóµ÷¼Á +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
[¿¼ÑÐ] 281Çóµ÷¼Á +3 ÑÇ¿ËÎ÷good 2026-03-26 5/250 2026-03-26 19:48 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] ²ÄÁÏ¿¼ÑÐÇóµ÷¼Á +3 Dendel 2026-03-23 6/300 2026-03-26 17:51 by fmesaito
[¿¼ÑÐ] 352Çóµ÷¼Á +4 ´óÃ×·¹£¡ 2026-03-22 4/200 2026-03-26 16:40 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ´óѧ»¯Ñ§Ñ§Ë¶0703×Ü·Ö337Çóµ÷¼Á +7 niko- 2026-03-22 7/350 2026-03-25 20:14 by qingfeng258
[¿¼ÑÐ] 282Çóµ÷¼Á +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[¿¼ÑÐ] 318Çóµ÷¼Á +3 plumÀî×Ó 2026-03-23 3/150 2026-03-25 09:42 by ÎíÉ¢ºóÏàÓölc
[¿¼ÑÐ] 289²ÄÁÏÓ뻯¹¤£¨085600£©BÇøÇóµ÷¼Á +4 ÕâôÃû×ÖÕ¦Ñù 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[¿¼ÑÐ] 085404µç×ÓÐÅÏ¢284·ÖÇóµ÷¼Á +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³£¬×Ü·Ö308Çóµ÷¼Á +7 īīĮ 2026-03-23 8/400 2026-03-23 20:36 by Creta
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û