逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
7 R6 t, |1 d* F- xsyms t;
P(1:k+1) = t;
; @3 o8 h  O7 S- @6 p# {9 RP(1) = 1;
) E+ S% Y* @4 mP(2) = t;
4 \6 \! V  M4 m$ ^c(1:k+1) = 0.0;
6 u) ]$ O0 v4 m' Z' jc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
+ p4 f% W8 B# f7 hc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
1 A, e  X9 `0 ?f = c(1)+c(2)*t;
8 k* h1 c7 H4 O6 a$ d8 s0 efor i=3:k+1
2 I0 }6 }$ W! h$ g5 _; E, ^# N    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
6 a7 l8 n+ y5 c' C    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
    f = f + c(i)*P(i);
9 ?! l+ c; X) l. w, i# J* K* h   
+ _) Z, `3 o' v; z; X, g    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
        else
  a$ l; t! S, N. X7 G1 i: H            f = vpa(f,6);
        end
    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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