逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
P(1:k+1) = t;
P(1) = 1;
! N" ?' m. Y; d4 `* }P(2) = t;
6 D2 k; E/ W/ x- g7 ?$ k  u. Uc(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;
. t& {$ [& q3 b. t! ?; |for i=3:k+1
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
0 O7 _5 A; o% X4 N4 @    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
    f = f + c(i)*P(i);
   
4 u) p- y; Q& M& h! S    if(i==k+1)
        if(nargin == 3)
" B; d0 l6 C: Q) |" H            f = subs(f,'t',x0);
        else
  R7 W8 G% c' L            f = vpa(f,6);
        end
/ {) t8 h, H) f, I& E5 ?; A9 \) Y    end
; u, O, h" y+ cend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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