逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
3 D2 T  l& B* b0 o* t4 ^" q2 Gsyms t;
P(1:k+1) = t;
$ G3 n  @" p* l+ D- uP(1) = 1;
7 b' y6 ^( a# pP(2) = t;
c(1:k+1) = 0.0;
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
6 @; u* l* i  nf = c(1)+c(2)*t;
! X1 H/ \" W  R( o. \( Vfor i=3:k+1
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
    f = f + c(i)*P(i);
2 o  |+ i" t8 \0 E  C   
    if(i==k+1)
        if(nargin == 3)
" S" L$ ]# `% H1 k            f = subs(f,'t',x0);
. }# A  o. `/ c3 c  @( S; {        else
5 Y) i: }4 n5 [" D1 B' Q            f = vpa(f,6);
        end
    end
end

zqyzixin

呵呵，非常感谢，很好很强大

昝斌

这是什么哦？

hjl19890710

好 dwadedcsargtg

hjl19890710

lihan厉害厉害厉害厉害厉害

coolistone

顶顶顶顶顶顶顶顶顶顶顶顶的

Kadyniost

？？？？？