²é¿´: 1660  |  »Ø¸´: 11

chj787

½ð³æ (СÓÐÃûÆø)

[ÇóÖú] matlabÄâºÏÇóÖú

ÓÃmatlabÄâºÏÒ»¸ö·½³Ì£º
    y=(a+bx+cx^2+dx^3)exp(-bx)
x=[0.8,0.85,0.9,0.94,1,1.1,1.14,1.144,1.146,1.148,1.149,1.15,1.151,1.152,1.153,1.157,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2,2.05,2.1,2.15,2.2,2.25,2.3,2.35,2.4];
y=[129698.6607,81517.14907,48462.38869,30160.54842,12625.39888,830.7904906,0,2.343989057,8.695584894,18.42928479,24.54823752,31.4946096,39.26181681,47.84108015,57.22581538,102.7141275,687.9234793,3872.493796,7589.20011,11578.99871,15654.88623,19681.73873,23751.25051,27388.0944,30756.13673,33816.1486,36542.19252,38921.37873,40955.22819,42660.53072,44067.23581,45213.10408,46138.07991,46880.36067,47474.18002,47948.97387,48329.25318,48634.92208,48881.88591,49082.65691];

±äÁ¿£ºx,y
ÄâºÏ³£Êý£º a,b,c,d
Çó¸ßÊÖ°ïÖú£¡¸ø³ömatlabÄâºÏ³ÌÐò¡£

[ Last edited by chj787 on 2012-5-20 at 16:29 ]
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
»ØÌûÖ§³Ö ( ÏÔʾ֧³Ö¶È×î¸ßµÄÇ° 50 Ãû )

cactus586

½ð³æ (ÖøÃûдÊÖ)

Á÷À˹ýµÄÈË×ÜÊÇÏëÔÙ´Î̤ÉÏÂÃ;

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xzhdty: ½ð±Ò+1, »¶Ó­³£À´³ÌÐòÓïÑÔ¿´¿´ 2012-05-21 05:57:43
chj787: ½ð±Ò+30, ¡ïÓаïÖú, 1 2012-05-24 21:55:56
×òÌì¿´µ½Ò»¸öÓÃMatLabÑÝ×࿨ٯµÄ ¿ÉÒÔÕ¾¶ÌÎÊÎÊËû
http://muchong.com/bbs/viewthread.php?tid=4519096
¾²Ë®Á÷Éî
2Â¥2012-05-20 16:40:48
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
ÆÕͨ»ØÌû

chj787

½ð³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by cactus586 at 2012-05-20 16:40:48:
×òÌì¿´µ½Ò»¸öÓÃMatLabÑÝ×࿨ٯµÄ ¿ÉÒÔÕ¾¶ÌÎÊÎÊËû
http://muchong.com/bbs/viewthread.php?tid=4519096

лл£¡
3Â¥2012-05-20 17:59:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lijie169

Í­³æ (ÖøÃûдÊÖ)

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

¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xzhdty: ½ð±Ò+2, »¶Ó­³£À´³ÌÐòÓïÑÔ 2012-05-21 05:58:00
chj787: ½ð±Ò+3, ¡ï¡ï¡ïºÜÓаïÖú 2012-05-29 15:30:03
CODE:
function f=fit(x,y)
f=lsqnonlin(@(a)fun(a,x,y),[1,1,1,1]');
function f1=fun(a,x,y)
f1=y-(a(1)+a(2)*x+a(3)*x.^2+a(4)*x.^3).*exp(-a(2)*x);
end
end

x,y¾ÍÊÇÂ¥Ö÷ÊäÈëµÄ,aÊÇÒ»¸ö¾ØÕ󣬶ÔÓ¦×Åabcd
4Â¥2012-05-20 18:08:21
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lijie169

Í­³æ (ÖøÃûдÊÖ)

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

¶î£¬ÊäÈëfit(x,y)¾Í¿ÉÒÔÁË
5Â¥2012-05-20 18:09:06
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

coolxh

ľ³æ (СÓÐÃûÆø)

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

¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xzhdty: ½ð±Ò+1, »¶Ó­³£À´³ÌÐòÓïÑÔ¿´¿´ 2012-05-21 05:58:22
chj787: ½ð±Ò+2, ¡ï¡ï¡ïºÜÓаïÖú 2012-05-29 15:29:52
Èô²»Ïë×Ô¼ºÐ´´úÂ룬Òà¿ÉÓÃÇúÏßÄâºÏ¹¤¾ß£¬
cftool(x,y)

È»ºó²ÉÓÃ×Ô¶¨Ò幫ʽ£¬°´ÄãÏëÒªµÄ¹«Ê½¶¨Òå±í´ïʽ£¬½øÐÐÄâºÏ

È»ºóÔÙ×Ô¶¯Éú³ÉÄâºÏ´úÂë
6Â¥2012-05-20 21:43:41
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lgycjpcqu

½ð³æ (ÕýʽдÊÖ)

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xzhdty: ½ð±Ò+1, ллÁË 2012-05-22 06:01:48
chj787: ½ð±Ò+15, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, ·Ç³£¸Ðл£¡ÆÀ·ÖÍíÁË£¬Íò·Ö±§Ç¸~£¡ 2012-05-29 15:29:20
x=[0.8,0.85,0.9,0.94,1,1.1,1.14,1.144,1.146,1.148,1.149,1.15,1.151,1.152,1.153,1.157,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2,2.05,2.1,2.15,2.2,2.25,2.3,2.35,2.4];
y=[129698.6607,81517.14907,48462.38869,30160.54842,12625.39888,830.7904906,0,2.343989057,8.695584894,18.42928479,24.54823752,31.4946096,39.26181681,47.84108015,57.22581538,102.7141275,687.9234793,3872.493796,7589.20011,11578.99871,15654.88623,19681.73873,23751.25051,27388.0944,30756.13673,33816.1486,36542.19252,38921.37873,40955.22819,42660.53072,44067.23581,45213.10408,46138.07991,46880.36067,47474.18002,47948.97387,48329.25318,48634.92208,48881.88591,49082.65691];
ybeta=inline('(beta(1)+beta(2)*x+beta(3)*x.^2+beta(4)*x.^3).*exp(-beta(2)*x)','beta','x');
beta=nlinfit(x,y,ybeta,[1e+004,0,-1e+004,1e+004]);
a=beta(1);
b=beta(2);
c=beta(3);
d=beta(4);
yy=(a+b*x+c*x.^2+d*x.^3).*exp(-b*x);
plot(x,yy,x,y,'*')
µ«ÊÇ»áÏÔʾ
Warning: The Jacobian at the solution is ill-conditioned, and some
model parameters may not be estimated well (they are not identifiable).
Use caution in making predictions.ËùÒÔ¿ÉÄÜ»¹ÐèÒª½øÒ»²½ÊÖ¶¯ÓÅ»¯£¬ÎÒ×öÄãºÎʱÓÃnlinfitÒ²¾­³£³öÏÖÕâ¸ö¾¯¸æ£¬ÆÚ´ýÓиßÊÖ½â´ð¡£
ÄâºÏ½á¹ûÊÇ
3866678.20479525       
2.55862906421524       
-8281766.51676900       
4642226.52606192
Ï£ÍûÄÜ°ïµ½Äã
7Â¥2012-05-21 10:31:55
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lijie169

Í­³æ (ÖøÃûдÊÖ)

ÒýÓûØÌû:
4Â¥: Originally posted by lijie169 at 2012-05-20 18:08:21
function f=fit(x,y)
f=lsqnonlin(@(a)fun(a,x,y),');
function f1=fun(a,x,y)
f1=y-(a(1)+a(2)*x+a(3)*x.^2+a(4)*x.^3).*exp(-a(2)*x);
end
end

x,y¾ÍÊÇÂ¥Ö÷ÊäÈëµÄ,aÊÇÒ»¸ö¾ØÕ󣬶ÔÓ¦×Åabcd...

ºÜÓаïÖú²Å¸ø3¸ö½ð±Ò  

» ±¾ÌûÒÑ»ñµÃµÄºì»¨£¨×îÐÂ10¶ä£©

8Â¥2012-05-29 17:14:23
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

chj787

½ð³æ (СÓÐÃûÆø)

ËÍÏÊ»¨Ò»¶ä
ÒýÓûØÌû:
8Â¥: Originally posted by lijie169 at 2012-05-29 17:14:23
ºÜÓаïÖú²Å¸ø3¸ö½ð±Ò  ...

9Â¥2012-05-30 22:10:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lijie169

Í­³æ (ÖøÃûдÊÖ)

ÒýÓûØÌû:
9Â¥: Originally posted by chj787 at 2012-05-30 22:10:20
...

10Â¥2012-05-31 07:04:05
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌøת ÎÒÒª¶©ÔÄÂ¥Ö÷ chj787 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍƼö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[ÂÛÎÄͶ¸å] ·¢±í¿ìµÄËÄÇøsci +3 yangqingxy 2025-02-17 6/300 2025-02-19 02:09 by bobvan
[½Ìʦ֮¼Ò] ¸¸Ä¸ÈÕÒæ±äÀÏ£¬¶àÅã°é¸¸Ä¸°É£¬±ÈÕù¶áԺϵÁìµ¼Ö°Îñ¡¢·¢±íÂÛÎĵÈÓÐÒâÒå¶àÁË +11 ˲ϢÓîÖæ 2025-02-17 12/600 2025-02-18 20:38 by ½£¿Í
[»ù½ðÉêÇë] ¸ø½ÜÇàÓÅÇà¸ÄÃû×ÖÕâÖÖÃ÷ÏÔµÄÑÚ¶úµÁÁåµÄʶùΪɶ»á·¢Éú? +4 babu2015 2025-02-17 7/350 2025-02-18 20:33 by brblmd
[˶²©¼ÒÔ°] ¸Ð¾õÕâ±²×ÓÒѾ­Íêµ°ÁË +27 Ϻ»¬ºÃ³ÔÂð 2025-02-14 41/2050 2025-02-18 17:27 by crusher1994
[ÂÛÎÄͶ¸å] Dalton Transactions 10+3 ÕÅæ¼æ¥zjy 2025-02-16 4/200 2025-02-18 15:32 by maolC
[½Ìʦ֮¼Ò] ÓÐûÓа²»ÕºÍºÓÄϵØÇø¸ßУ²©Ê¿ÇóÖ°½»Á÷Ⱥ +4 Äϼ«10¡æ 2025-02-16 5/250 2025-02-18 15:13 by СÁúϺ2008
[»ù½ðÉêÇë] AIʱ´ú£¬×ÔÈ»»ù½ðÉêÇëÊéÔÙ¸ã°Ë¹ÉÎÄ£¬ÄýÁ¶ÎÊÌâ¾ÍÊǸãЦÁË +5 babu2015 2025-02-13 6/300 2025-02-18 14:35 by adist
[ÂÛÎÄͶ¸å] SCIÒ»Çø×ÜÊDZ»¾Ü¾ø?£¬·³ËÀÁË£¡ +5 ÑâʦÐÖ 2025-02-17 7/350 2025-02-18 14:27 by ÑâʦÐÖ
[¿¼ÑÐ] ѧ¿ÆÓ¢Óïµ÷¼Á +4 Àõ×Ó¼ÓÓÍѼ 2025-02-16 7/350 2025-02-18 06:47 by Ò»ºø×ǾÆ6686
[¿¼ÑÐ] ¿¼ÑÐÕÐÉú +11 ³Ì³¤£¨zhang£© 2025-02-14 21/1050 2025-02-17 22:00 by רҵ¿Î²ÄÁÏÁ¦Ñ§£
[˶²©¼ÒÔ°] ´ó¼Ò¸øÎҵ㽨Òé +16 ɽÓêÒàÆæ 2025-02-16 18/900 2025-02-17 20:16 by huixiong0627
[ÂÛÎÄͶ¸å] »·¾³×¨ÒµÆÚ¿¯Í¶¸å 15+3 »Ê×Ó?±£ 2025-02-14 5/250 2025-02-17 17:01 by Wanggy
[»ù½ðÉêÇë] B06ÈϿɵÄÆÚ¿¯ÓÐÄÄЩ£¿ +3 jqq123456 2025-02-14 6/300 2025-02-17 15:34 by jqq123456
[˶²©¼ÒÔ°] ¸Ã²»¸ÃÌṩÊý¾Ý +17 lupo140155 2025-02-12 33/1650 2025-02-16 20:55 by ÁÁÁÁ¼ÓÓÍ
[˶²©¼ÒÔ°] É격-²©µ¼Í¨¹ýÓʼþѯÎÊÎÊÌâ +4 CG-KY 2025-02-16 4/200 2025-02-16 19:27 by white7wine
[¿¼²©] ²©Ê¿ÉêÇë +4 yuhui51480 2025-02-15 6/300 2025-02-16 17:05 by dmj0610
[¿¼ÑÐ] Çóµ÷¼Á +3 Àõ×Ó¼ÓÓÍѼ 2025-02-16 3/150 2025-02-16 13:33 by 5d94@0eq
[¿¼²©] ¸£ÖÝ´óѧÐÂÄÜÔ´²ÄÁÏÓ빤³ÌÑо¿ÔºÕÐÊÕ¹¤³Ì²©Ê¿Éú +4 °ëÀå 2025-02-13 7/350 2025-02-16 11:02 by Ricoch4t
[¿¼ÑÐ] 320Çóµ÷¼Á +3 ÓÀ²»ÑÔÆúµÄС͸à 2025-02-12 3/150 2025-02-15 10:47 by д´úÂëµÄcoder
[΢Ã׺ÍÄÉÃ×] ¹¤³Ì²©Ê¿ 5+3 541258ldj 2025-02-12 5/250 2025-02-14 12:54 by EGXYRZ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û