逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
- t6 {! f  F1 ~) l+ L  Nsyms t;
/ A+ o( k% o- O: k/ ZP(1:k+1) = t;
' l8 b3 T  j( n9 G  k! p' _P(1) = 1;
1 h2 i) }, T  k% |4 h/ n* zP(2) = t;
8 A: |0 Z6 P4 d6 ?5 x- N8 wc(1:k+1) = 0.0;
: }# u. r. @- i" A8 j2 Ic(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
0 O7 w! V& ?# nc(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);
    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
    f = f + c(i)*P(i);
1 E. m' A2 h0 x3 R   
    if(i==k+1)
0 t: _0 |0 D; O) q" @8 y9 m        if(nargin == 3)
            f = subs(f,'t',x0);
        else
            f = vpa(f,6);
        end
    end
, t3 k' Q, O! S2 j' \  H3 a4 Q3 i/ bend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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