逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
P(1:k+1) = t;
( D9 e6 A! q0 ]4 p% pP(1) = 1;
. g# {( ^+ H) P  J5 V/ T/ M3 sP(2) = t;
3 o" D! l  n% M/ j, [3 a" Hc(1:k+1) = 0.0;
3 e8 u4 @$ y' s$ e' Oc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
! g) E0 n, d* ^9 c. u* L/ xc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
  Y" P4 ~0 m; [2 `f = c(1)+c(2)*t;
for i=3:k+1
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
% x* J; [' a/ c  w    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
    f = f + c(i)*P(i);
/ Z$ O3 V7 U1 n; Y   
( s4 J- T. K; t. w# w9 y    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
/ h4 E7 r6 w, u: p: a1 V        else
) ?( @% C6 W9 {8 [            f = vpa(f,6);
        end
    end
end

zqyzixin

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

昝斌

这是什么哦？

hjl19890710

好 dwadedcsargtg

hjl19890710

lihan厉害厉害厉害厉害厉害

coolistone

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

Kadyniost

？？？？？