逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
# ?; w4 @; @4 VP(1:k+1) = t;
P(1) = 1;
! [; M( y$ B$ e+ O0 MP(2) = t;
5 ]9 y; F: m1 p6 f( v7 ic(1:k+1) = 0.0;
0 @: M/ b7 R" `: W4 D7 a% {# R) Jc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
7 v" Z) z% {4 D- S% f( Ec(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
; s# T7 d* Q- a0 T2 k    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;
& t( l! n% q$ R7 K1 {5 f    f = f + c(i)*P(i);
   
) F, L8 \3 F* [% ^8 N- p    if(i==k+1)
        if(nargin == 3)
" ^  k8 }  g' \; a5 I$ n2 y* j            f = subs(f,'t',x0);
        else
            f = vpa(f,6);
        end
# ^, A) f! S# ]' z8 P: a0 q% X    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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