逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
5 ~" G: J/ F" Q6 OP(1:k+1) = t;
8 K' u# ^; c* l' d$ pP(1) = 1;
P(2) = t;
c(1:k+1) = 0.0;
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
* @1 k& N* n: K) Y7 v# Uc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
' [1 @2 A# b; o, nf = c(1)+c(2)*t;
for i=3:k+1
" B4 v9 p* T$ y; I- ~0 m2 B- H0 t    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;
% v, J0 l/ v! H4 H; S5 L$ y: S: R    f = f + c(i)*P(i);
   
! ^! m* b; u8 @0 j( E    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
" e& e/ u! R' }( m1 l: [# v        else
            f = vpa(f,6);
5 O) D; o) M+ {' w& c; y8 I5 I4 J3 f        end
    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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