逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
- l" I6 g  u6 Q8 F7 `6 {$ x' ^syms t;
" {4 e' y$ Z# ^( UP(1:k+1) = t;
; M: T6 m& c7 L' g& j# `* m& OP(1) = 1;
+ z/ f* j2 ^2 Z! AP(2) = t;
  e( V& s6 L2 q$ _c(1:k+1) = 0.0;
- [6 \9 `% z4 ]% Mc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
! j, }, U" O6 Gc(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
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
) h$ K' U, X% r+ Y. X7 O- r    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
9 O/ F4 R8 L% w, _3 V    f = f + c(i)*P(i);
   
    if(i==k+1)
        if(nargin == 3)
9 D5 O; @' K  R' h7 ~5 @            f = subs(f,'t',x0);
        else
. d9 E4 [& X: N7 Y  q! X            f = vpa(f,6);
        end
& H. J* W3 ?. n" A2 x3 b    end
+ C$ G4 Q8 ]+ p' e+ c  oend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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