逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
' J( R$ N4 }% A6 dsyms t;
! L1 r. b/ t7 i4 Z! C1 TP(1:k+1) = t;
5 z/ i- Z5 d0 YP(1) = 1;
P(2) = t;
c(1:k+1) = 0.0;
# h* P! \- e* N# A2 T3 Y2 gc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
0 j1 X4 O- r6 k8 E5 ~c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
4 u  c4 ~+ s' {% [8 Hf = c(1)+c(2)*t;
" B& f( ^; x: [! s1 u5 \/ mfor i=3:k+1
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
5 q9 G1 U1 _& C8 S/ G4 `7 ~    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
. Y& F) a8 |& B, s! W    f = f + c(i)*P(i);
   
    if(i==k+1)
" I& {& D  w6 k( w        if(nargin == 3)
5 C/ j4 k1 \3 ]+ E- |            f = subs(f,'t',x0);
        else
* \& q3 G3 H; E, v+ f$ J            f = vpa(f,6);
5 E' R) S$ _8 V  n7 M/ k& z        end
$ a" c- a/ [2 l9 `. Q/ B- w* N    end
2 z3 i/ }; Q) |( L: r( m4 R# N1 Jend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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