逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
' b' n& F/ }, J+ oP(1:k+1) = t;
P(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;
& B& S( J3 l4 G' E' M6 X, ac(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
f = c(1)+c(2)*t;
, j& m- _. [9 Z' `for 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;
  d  A5 E: W, |/ R2 M( P    f = f + c(i)*P(i);
7 B% Z$ r1 j( `" P   
    if(i==k+1)
6 [( z5 k6 Q7 P4 |$ m+ c+ Q        if(nargin == 3)
            f = subs(f,'t',x0);
) i* R7 p4 F+ w8 P" F( x        else
2 @4 F6 C, y4 P  z2 u            f = vpa(f,6);
  q1 t  O! q& s4 o4 \' ?6 o2 X0 n        end
    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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