[交流] 【致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

6楼: Originally posted by sudo at 2011-09-09 19:42:55:
这个拟合出来的式子也可以无限长啊

哦，莫非你觉得不是整数不爽？简单改改就行：

把n从0开始计数，设给定的数列长度为N-1，同理使用上面的方法得到一个多项式F(n)之后，再令

f(n) = F(n mod N)

不就 ...

-1,1,-1,1,-1,1,-1,1,...

FindSequenceFunction[{-1, 1, -1, 1, -1, 1, -1, 1}, n]

(-1)^n

InterpolatingPolynomial[{-1, 1, -1, 1, -1, 1, -1, 1}, x] // Expand

-255 + (64864 x)/105 - (2856 x^2)/5 + (12056 x^3)/45 - 70 x^4 + (464 x^5)/45 - (4 x^6)/5 + (8 x^7)/315

y(x)=-255 + (64864 x)/105 - (2856 x^2)/5 + (12056 x^3)/45 - 70 x^4 + (464 x^5)/45 - (4 x^6)/5 + (8 x^7)/315

1楼: Originally posted by sudo at 2011-09-09 16:11:12:
中秋快到了无心工作，就娱乐一下吧~

不知道是怎么回事的虫友请看这个帖子：
http://muchong.com/bbs/viewthread.php?tid=3573558&fpage=0&view=&highlight=&page=1

分析也在那 ...