逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
) z6 R. l' V4 e9 J& t9 [syms t;
P(1:k+1) = t;
P(1) = 1;
0 t- o0 Z- |( D7 x+ `/ O5 GP(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;
f = c(1)+c(2)*t;
for i=3:k+1
9 D* [' j0 n( J* s' C8 n    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
/ p) {% K6 `( I2 D" P    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
  x4 K$ r) p- [5 C    f = f + c(i)*P(i);
. d1 T0 ~3 t0 ^9 h; |  [2 O   
    if(i==k+1)
        if(nargin == 3)
. k- G( J$ R. D& k8 B            f = subs(f,'t',x0);
& |* _7 Y; e& J: T; c        else
) Z# d6 }  L. f- V- P4 d; \            f = vpa(f,6);
        end
    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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