逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
! e  Q& q) q( _' S9 L2 J7 xsyms t;
P(1:k+1) = t;
P(1) = 1;
; T+ j) w5 P9 _8 MP(2) = t;
# j* s$ k& C: w: I! Ic(1:k+1) = 0.0;
& T7 o. x, l8 ~' ]c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
1 o% r  o  p- G( g* ~, _+ m4 ]c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
3 Q" |+ S$ v3 F& ]  m' X: i& ?: Q6 bf = 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);
    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
8 @) v, y3 m9 ^. n, W    f = f + c(i)*P(i);
   
+ n: \( [) v0 D    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
: f9 D, o1 d; j3 G( g8 K        else
            f = vpa(f,6);
        end
    end
2 T) q& u% e3 l3 Xend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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