²é¿´: 2052  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸Î人Àí¹¤²ÄÁϹ¤³Ìר˶µ÷¼Á +5 Doleres 2026-03-19 5/250 2026-03-19 20:14 by ÖÆ¶ÈµÄ
[¿¼ÑÐ] Áº³ÉΰÀÏʦ¿ÎÌâ×é»¶Ó­ÄãµÄ¼ÓÈë +9 һѼѼӴ 2026-03-14 11/550 2026-03-19 17:22 by £¡±¾°µÒ»´Î£¡
[¿¼²©] ¶«»ªÀí¹¤´óѧ»¯²Äרҵ26½ì˶ʿ²©Ê¿ÉêÇë +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by ÇáËɲ»ÉÙËæ
[¿¼ÑÐ] ÕÐÊÕµ÷¼Á˶ʿ +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by Óà÷ëÓà
[¿¼ÑÐ] ²ÄÁÏר˶274Ò»Ö¾Ô¸ÉÂÎ÷ʦ·¶´óѧÇóµ÷¼Á +8 Ñ¦ÔÆÅô 2026-03-13 8/400 2026-03-19 15:36 by haoshis
[¿¼ÑÐ] 0703»¯Ñ§ 305Çóµ÷¼Á +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[¿¼ÑÐ] 304Çóµ÷¼Á +6 ˾¿Õ. 2026-03-18 6/300 2026-03-18 23:03 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 330Çóµ÷¼Á +3 С²Ä»¯±¾¿Æ 2026-03-18 3/150 2026-03-18 21:55 by ÎÞи¿É»÷111
[¿¼ÑÐ] 295Çóµ÷¼Á +3 Ò»Ö¾Ô¸¾©Çø211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[¿¼ÑÐ] 311Çóµ÷¼Á +11 ¶¬Ê®Èý 2026-03-15 12/600 2026-03-18 14:36 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 278Çóµ÷¼Á +5 ÑÌ»ðÏÈÓÚ´º 2026-03-17 5/250 2026-03-18 08:43 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] 290Çóµ÷¼Á +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[¿¼ÑÐ] 275Çóµ÷¼Á +4 Ì«Ñô»¨ÌìÌ쿪ÐÄ 2026-03-16 4/200 2026-03-17 10:53 by ¹¦·ò·è¿ñ
[¿¼ÑÐ] 302Çóµ÷¼Á +4 С¼Öͬѧ123 2026-03-15 8/400 2026-03-17 10:33 by С¼Öͬѧ123
[¿¼ÑÐ] 318Çóµ÷¼Á +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[¿¼ÑÐ] 085601²ÄÁϹ¤³Ì315·ÖÇóµ÷¼Á +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[¿¼ÑÐ] 288Çóµ÷¼Á +4 Ææµã0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[¿¼ÑÐ] 266Çóµ÷¼Á +4 ѧԱ97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[¿¼ÑÐ] ÕÐÊÕ0805£¨²ÄÁÏ£©µ÷¼Á +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%¡¢
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û