逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
$ i9 Y0 F& Y# nP(1:k+1) = t;
7 _( V2 a$ L& L9 Z4 P2 KP(1) = 1;
* Q. X! S' O$ W* ]6 gP(2) = t;
0 T2 B' m' v+ X2 D9 Wc(1:k+1) = 0.0;
$ V' k* A! {% mc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
9 h3 v' B+ g3 z3 f8 b# _c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
/ |7 Y$ S$ h8 e7 O/ K% I  Tf = c(1)+c(2)*t;
/ n  V9 d9 }& T6 u4 d; Mfor i=3:k+1
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
. `3 f( i- a! ]9 S. |    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
! v! r; ]% P1 K$ a! p    f = f + c(i)*P(i);
1 p( L" X# u' a3 ~  @   
    if(i==k+1)
* t) X9 P7 S2 X; }        if(nargin == 3)
1 ]- ~6 x  v8 U, p            f = subs(f,'t',x0);
        else
7 A. n9 R3 H; ]            f = vpa(f,6);
        end
& S9 P8 b( T  c* R; j    end
- d! J8 U& ?* C7 h8 H% P6 Zend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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