[½»Á÷] ¡¾ÖÂHÐÖ¡¿ÊýÁÐÌâµÄÍòÄܶàÏîʽÄâºÏÇó·¨ ÒÑÓÐ2È˲ÎÓë

function av = myPolyFit(y)
%A * av = y
    n = length(y);
    y = y(:); %make `y` a column vector
   
    %construct `A`
    A = (1:n)' * ones(1,n);
    for j = n-2 :-1: 1
        A(:,j) = A(:,j) .* A(:,j+1);
    end
    A(:,n) = ones(n,1);
   
    av = A\y;

y=[22 26 20 23 24 21 20 18];
av=myPolyFit(y);
n=1:0.001:9;
val=polyval(av, n);
poly2str(av, 'n')
plot(n, val)

   0.0047619 n^7 - 0.18889 n^6 + 3.025 n^5 - 25.1806 n^4 + 116.3583 n^3
   - 293.6306 n^2 + 365.6119 n - 144

1Â¥: Originally posted by sudo at 2011-09-09 16:11:12:
ÖÐÇï¿ìµ½ÁËÎÞÐŤ×÷£¬¾ÍÓéÀÖһϰÉ~

²»ÖªµÀÊÇÔõô»ØʵijæÓÑÇë¿´Õâ¸öÌû×Ó£º
http://muchong.com/bbs/viewthread.php?tid=3573558&fpage=0&view=&highlight=&page=1

·ÖÎöÒ²ÔÚÄÇ ...

2Â¥: Originally posted by huycwork at 2011-09-09 17:51:44:
Ù¯Õ⶯¾²ÄÖµÃÕâ°ãÕ𺳣¬·ç¿ÚÀ˼â¶ùµÄ
Ò»°ãÇé¿öÏÂÈç¹ûµÈ²îÊýÁУ¬±ÈÈç1 2 3 4 5 6¡£
°³±íʾµ­¶¨£¬Ëµ£¬Õâ¸öÓÃÒ»¸öÏî¾Í¿ÉÒԵóö´ð°¸ÁË£¡
Matlab±íʾµ­¶¨£¬Ï룬sudoͬѧ¸øÁË°³Ò»¸ö¶Î¶ùÄØ£¡
¼ÆËã»ú±íʾµ­¶¨£¬ºÏ ...

   1.8575e-018 n^5 - 2.8578e-017 n^4 + 1.6489e-016 n^3 - 4.2998e-016 n^2
   + 1 n - 2.0817e-016

3Â¥: Originally posted by wangww2011 at 2011-09-09 18:03:00:
ÖÐÇï¿ìÀÖ£¡

ÄãÕâ¸ö´¿¶àÏîʽÄâºÏ²»¿¿Æ×£¬ÒòΪÕâÊǾֲ¿µÄÄâºÏ£¬¶øʵ¼ÊÉÏÄǸöÊýÁпÉÒÔÎÞÏÞ³¤µÄ°¡

4Â¥: Originally posted by huycwork at 2011-09-09 18:39:45:
ÒÀÎÒ¿´£¬´ý¶¨ÏµÊý·¨¾Í²î²»ÀëÁË£¬Äܽâ¾ö90%ÒÔÉÏ´ËÀàÌ⣬·ÖÁ½À࣬һÀàÊÇÖ¸Êýʽ£º
f2n = ab^n
ÁíÍâÒ»ÀàÊǶàÏîʽ£¬Õâ¸ö×î¸ßÏîÓÉÊý¾ÝµÄ¶àÊý¾ö¶¨£¬ÉÏÃæµÄʽ×ÓռȥÁ½¸ö²ÎÊý£¬Ê£Ïµľ͹éËü£º
f1n = an^3+bn^2+cn+d
½« ...

7Â¥: Originally posted by sudo at 2011-09-09 19:47:19:
ÎÒŪÕâ¸öÖ÷ÒªÊÇΪÁ˱íʾ¶ÔÕâÖÖÌâÄ¿µÄÎÞÓ

Èç¹ûÃæÊÔµÄʱºò±»¿¼¹ÙÎʵ½ÀàËƵÄÌâÄ¿£¬¿ÉÒÔÏÈÏêϸµØÎÊÎÊËû¡°ÓÐʲô¾ßÌåµÄÏÞ¶¨Ìõ¼þûÓС±

Èç¹ûËû²»¼ÓÉÏÖîÈç¡°ÄǸöͨÏʽµÄ¶¨ÒåÓòÈç¹û»»³ÉʵÊýÓò£¬ÄÇôËüÔÚʵ ...

8Â¥: Originally posted by sudo at 2011-09-09 20:06:32:
Ŷ£¬²îµãÍüÁË£¬×£¸÷λÖÐÇï¿ìÀÖ~~~