逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
( L, ~0 q8 V% H* h0 w& a6 C1 _3 B5 GP(1:k+1) = t;
1 C% o' E) t  ^) F4 d# zP(1) = 1;
; F" M: x. [* v0 d: P" Z' p1 dP(2) = t;
c(1:k+1) = 0.0;
c(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;
f = c(1)+c(2)*t;
, D' P* v6 t" h+ }+ l" m8 xfor i=3:k+1
  J4 N2 q1 P0 a2 x& `    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
2 T, |7 N7 b. {9 L    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
    f = f + c(i)*P(i);
   
    if(i==k+1)
7 X6 @/ x- m8 S$ f7 Q# D7 P        if(nargin == 3)
            f = subs(f,'t',x0);
- _$ D, \! k+ m7 \, a' p" K+ x        else
  O# n* Z! M; Y; L; y            f = vpa(f,6);
8 a6 j) E0 Y3 U. O  O7 C8 i: D        end
    end
& u- o5 x# }8 Vend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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