Znn3bq.jpeg
ÉÇÍ·´óѧº£Ñó¿ÆÑ§½ÓÊܵ÷¼Á
²é¿´: 2752  |  »Ø¸´: 7

lihongru

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

[ÇóÖú] matlab·ÇÏßÐÔ²ÎÊýÄâºÏÎÊÌâ

ÏÖÒÑÖªxºÍyµÄһϵÁÐʵÑéÖµ£¬ÇÒÒÑÖªxºÍy¾ßÓÐÈçϹØÏµ£º
y=kln(x-ky)+a,ÆÚÖÐkºÍaÊÇ´ýÄâºÏµÄ²ÎÊý£¬ÎÒÒ»¿ªÊ¼²ÉÓÃÏßÐÔÄâºÍµÄ·½·¨£¬ÏȸøÀ¨ºÅÖеÄk¸³³õÖµk0£¬¶ÔyºÍln(x-k0)ÏßÐԻعé³öеÄkºÍa£¬È»ºó¶Ôk½øÐеü´úºÍÏßÐÔÄâºÏ£¬Ö±µ½kÊÕÁ²£¬ÏÖÔÚµÄÎÊÌâÊÇÓеÄÇé¿öÏÂÄܹ»ÊÕÁ²³ök£¬ÓÐʱ²»ÊÕÁ²£¬ÇëÎÊÓÐûÓиüºÃµÄ·ÇÏßÐԻع鷽·¨£¬ÇóµÃkºÍa²ÎÊýÄØ£¿
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

dingd

Ìú¸Ëľ³æ (Ö°Òµ×÷¼Ò)

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

lihongru(½ð±Ò+10): лл£¬ÎÒ¶Ô1stOpt²»ÊÇÌØ±ðÊìϤ£¬»ØÍ·ÊÔһϡ£ 2011-10-24 20:58:45
ÄãµÄÎÊÌâÒѱä³ÉÒþº¯ÊýÄâºÏÁË£¬¿ÉÒÔÊÔÏÂ1stOpt£¬´¦ÀíµÄºÜºÃ¡£
2Â¥2011-10-24 09:36:44
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yqx1985

ľ³æ (ÖøÃûдÊÖ)

ÔÆÖÐÏÉ

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

¡ï ¡ï
΢³¾¡¢ÃÎÏë(½ð±Ò+2): лл²ÎÓëÓ¦Öú 2011-10-24 17:46:23
lihongru(½ð±Ò+10): лл£¬ÎÒÊÔ×ÅÔËÐÐһϡ£ 2011-10-24 21:00:32
ÄãµÄ±í´ïʽ£¬¾­¹ýת»¯²»¾ÍÊÇÏÂÃæµÄÐÎʽÂð£¿
x=e^[(y-a)/k]+ky;
Äã¿ÉÒÔÔÚmatlabÖн¨Á¢Ò»¸öº¯Êý
function y=fun(x,p)
a=p(1);k=p(2);
y=e^((x-a)./k)+k*x;
end
½ø¶øÄã¿ÉÒÔÓÃnlinfitº¯Êý½øÐзDzÎÊýÄâºÏ
options = statset('Robust','on');
[beta,r,J,Sigma,mse] = nlinfit(x,y,@fun,rand(1,2),options);
¼´¿ÉµÃµ½²ÎÊýµÄ¹À¼ÆÖµ£¬µ±È»£¬ÄãÒ²¿ÉÒÔÖØÐ¶¨Ò庯Êý£¬²ÉÓÃ×îС¶þ³ËµÄ·½·¨½øÐвÎÊý¹À¼Æ¡£
ÈËÃæ²»ÖªºÎʱȥ£¬ÌÒ»¨ÒÀ¾ÉЦ´º·ç
3Â¥2011-10-24 11:24:04
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lihongru

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

ÎÒÊÔ×ÅÓÃÁíÒ»ÖÖ·½·¨½âÁËһϣ¬¾ÍÊÇÎÒÖªµÀkµÄ´ó¸Å·¶Î§£¬È»ºó³ÌÐòÈçÏ£º
k=6:0.1:15%ÁíÒ»ÖÖÇókÖµµÄ·½·¨
g=length(k);
for i=1:g
d0=d-k(i)*s;
lnd0=log(d0);
X0=[lnd0,DT,cons];
p=X0\lns;
m=p(1)-k(i);
f=abs(m);
if f>=delt
    return
else
end
end
p
¿ÉÊÇÔËÐв»³ö½á¹û£¬Ò²²»Öª³ÌÐò´íÔÚʲôµØ·½¡£
4Â¥2011-10-24 21:04:03
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lihongru

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

ԭʼµÄÒªÄâºÏµÄ·½³ÌÊÇlns=kln(d-k*s)+a*DT+c£¬s¡¢dºÍDTÊÇÒÑÖªµÄÈý¸öÁÐÏòÁ¿£¬ÒªÄâºÏk¡¢aºÍc£¬ÈôDT¹Ì¶¨ºó¾Í³ÉÁËÒ»¸öÒþº¯ÊýÁË¡£²»Öª²ÉÓÃÉϱߵķ½·¨¿ÉÐв»¿ÉÐС£
5Â¥2011-10-24 21:13:27
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lihongru

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

ÉÏÃæ³ÌÐòÔËÐв»³öÀ´£¬ÓÃÈçϳÌÐò½â¾öÁË
delt=1*10^-4;%É趨µü´úÖÕÖ¹Öµ
k=6:0.1:15;%ÁíÒ»ÖÖÇókÖµµÄ·½·¨
f=1;
i=0;
if f>=delt
    i=i+1;
    d0=d-k(i)*s;
    lnd0=log(d0);
    X0=[lnd0,DT,cons];
    p=X0\lns;
    m=p(1)-k(i);
    f=abs(m);
end
p
²»ÖªÊDz»ÊÇÕýÈ·
6Â¥2011-10-24 22:04:24
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yqx1985

ľ³æ (ÖøÃûдÊÖ)

ÔÆÖÐÏÉ

¡ï
jjdg(½ð±Ò+1): ¸Ðл֧³Ö 2011-10-25 00:16:11
ÒýÓûØÌû:
3Â¥: Originally posted by yqx1985 at 2011-10-24 11:24:04:
ÄãµÄ±í´ïʽ£¬¾­¹ýת»¯²»¾ÍÊÇÏÂÃæµÄÐÎʽÂð£¿
x=e^[(y-a)/k]+ky;
Äã¿ÉÒÔÔÚmatlabÖн¨Á¢Ò»¸öº¯Êý
function y=fun(x,p)
a=p(1);k=p(2);
y=e^((x-a)./k)+k*x;
end
½ø¶øÄã¿ÉÒÔÓÃnlinfitº¯Êý½øÐзDzÎÊýÄâºÏ
opt ...

have a try
ÈËÃæ²»ÖªºÎʱȥ£¬ÌÒ»¨ÒÀ¾ÉЦ´º·ç
7Â¥2011-10-24 23:54:50
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lihongru

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

×îºóÓÃÈçϳÌÐò½â¾öÁË
delt=1*10^-4;%É趨µü´úÖÕÖ¹Öµ
k=6:0.00001:15;%ÁíÒ»ÖÖÇókÖµµÄ·½·¨
g=length(k);
f=1;
i=0;
while f>=delt
  i=i+1;
  d0=d-k(i)*s;
  lnd0=log(d0);
  X0=[lnd0,DT,cons];
  p=X0\lns;
  m=p(1)-k(i);
  f=abs(m);
end
p
i
¾ÍÊÇÕûÌå¸Ð¾õ±È½Ï±¿ÖØ£¬Èç¹ûÓ÷ÇÏßÐÔÄâºÏµÄ»°£¬¿ÉÄÜ»á¸üºÃ£¬¸Ðл¥ÉÏÁ½Î»µÄ°ïÖú£¡
8Â¥2011-10-26 08:05:48
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ lihongru µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 291 Çóµ÷¼Á +37 »¯¹¤2026½ì±ÏÒµÉ 2026-04-09 38/1900 2026-04-15 00:27 by lihaoda1994
[»ù½ðÉêÇë] ɽ¶«Ê¡»ù½ð2026 +7 jerry681 2026-04-08 9/450 2026-04-14 22:01 by °¢ÓµÄÀºÇò
[¿¼ÑÐ] 272·Ö²ÄÁÏ×ÓÇóµ÷¼Á +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[¿¼ÑÐ] »¯¹¤Ñ§Ë¶294·Ö£¬Çóµ¼Ê¦ÊÕÁô +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[¿¼ÑÐ] ¿ÒÇëÓÐѧУÊÕÁô +3 ¿Â»´È» 2026-04-12 3/150 2026-04-14 16:25 by ÄæË®³Ë·ç
[¿¼ÑÐ] ²ÄÁÏרҵ344Çóµ÷¼Á +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[¿¼ÑÐ] ¸÷λÀÏʦºÃ£¬Çóµ÷¼Á£¬±¾¿Æ211£¬Ò»Ö¾Ô¸Ìì½ò´óѧÉúÎïÓëҽҩѧ˶£¬²îÁ½Ãû¼ȡ¡£ +11 ·ÁùÁùjjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[¿¼ÑÐ] ²ÄÁÏ085601µ÷¼Á +32 ºÎÈó²É123 2026-04-10 34/1700 2026-04-14 08:47 by ľľmumu¡«
[¿¼ÑÐ] ¿¼ÑжþÂÖµ÷¼Á +10 ¹ÊÈË?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[¿¼ÑÐ] Ò»Ö¾Ô¸Õã´óÉúÎï325·ÖÇóµ÷¼Á +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[¿¼ÑÐ] 0860004 Çóµ÷¼Á 309·Ö +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[¿¼ÑÐ] ҩѧר˶µ÷¼Á +8 ? һ·Éú?»¨? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[¿¼ÑÐ] 293Çóµ÷¼Á +6 ÓÂÔ¶¿â°®314 2026-04-08 6/300 2026-04-11 10:08 by zhq0425
[¿¼ÑÐ] »¯Ñ§¹¤³ÌÓë¼¼Êõ324µ÷¼Á +23 Ëï³£»ª 2026-04-09 25/1250 2026-04-11 00:07 by ÆïÅ£¶Éº®½­
[¿¼ÑÐ] 22408 366·Ö£¬±¾¿Æ211£¬Ò»Ö¾Ô¸Î÷¹¤´ó +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[¿¼ÑÐ] Ò»Ö¾Ô¸»¦9£¬326ÉúÎïѧÇóÏà¹Ø×¨Òµµ÷¼Á +4 Áõīī 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[¿¼ÑÐ] 296Çóµ÷¼Á +6 Íô£¡£¿£¡ 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[¿¼ÑÐ] ²ÄÁÏר˶(0856) 339·ÖÇóµ÷¼Á +9 ¹þ¹þ¹þ¶ì¹þ¹þ¹þ 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[¿¼ÑÐ] µ÷¼Á +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[¿¼ÑÐ] 083200 ³õÊÔ305·Ö Çóµ÷¼Á Ôݲ»¿¼ÂÇ¿çרҵ +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û