逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
- u% A. I+ D  FP(1:k+1) = t;
% Y! w5 S/ }) h, v. wP(1) = 1;
: e, K" |% B! E  m6 DP(2) = t;
' _: k  M+ n& Dc(1:k+1) = 0.0;
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
7 R% I; ^# |9 I' Sc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
! f2 m$ N: U' ?% b# c! g: s1 ~7 r( Rf = 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);
/ u! I3 R0 B/ r. [    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
& ^4 ~+ X- D/ m0 v; ]$ @    f = f + c(i)*P(i);
   
3 K; j: a3 o/ j- i* V    if(i==k+1)
        if(nargin == 3)
            f = subs(f,'t',x0);
        else
            f = vpa(f,6);
        end
8 c* t8 x2 g! x" ~    end
* Z* E& n" n' d1 i1 N- Mend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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