逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
+ j) U. e) N9 `) j+ b2 @" dP(1:k+1) = t;
$ @5 d* C5 G) Q8 k: q; Q: C4 XP(1) = 1;
; [1 }% p( i6 IP(2) = t;
c(1:k+1) = 0.0;
- d! z0 @0 Z/ pc(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;
2 B* z' M* N5 a+ P9 sf = 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;
    f = f + c(i)*P(i);
" H0 |9 y& l* r) j, s   
& I0 b& z/ P7 U* L: [    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
        else
' j0 ~8 G5 z. ^# v2 p. _% x            f = vpa(f,6);
        end
    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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