逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
% L" `/ H7 B1 C0 Vsyms t;
( `8 s" }8 j# T( sP(1:k+1) = t;
P(1) = 1;
P(2) = t;
' r$ o/ E8 c/ E6 J; Jc(1:k+1) = 0.0;
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
: M2 Q8 t7 P9 w6 _- t: g% b+ Mc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
3 J" l6 L- m0 r& n* g% ?! O9 O8 `f = c(1)+c(2)*t;
for i=3:k+1
. }$ \* Q! t4 P! c! e( t7 w- v    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
6 r# u4 O+ f6 C, p- X    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
1 p4 V" U! F  U$ {& a6 S    f = f + c(i)*P(i);
   
( F& A3 y, N4 g( T0 ?& c; O    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
7 T. v$ G  x2 o0 [' E        else
            f = vpa(f,6);
        end
. w. [6 m& s& K6 ]    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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