逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
. N7 ^2 U* O) Q+ a* usyms t;
P(1:k+1) = t;
P(1) = 1;
P(2) = t;
' a7 i& W4 o( ]6 B' mc(1:k+1) = 0.0;
+ k( N) [- S; P" |' O. Kc(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
% \4 ~0 H# g& O" i" Gc(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
2 t( T  \1 S/ |1 b2 |( z( _    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;
6 T. K0 D$ T  F$ \0 P, E6 x; ~    f = f + c(i)*P(i);
   
5 E# X% L$ E$ P: L    if(i==k+1)
        if(nargin == 3)
+ u1 o& H! F# p+ Z            f = subs(f,'t',x0);
        else
* c6 z, ~+ ^6 c/ T+ o& c, B/ \7 f            f = vpa(f,6);
5 @, t& n8 }  I; M6 D        end
$ r7 h  n; |9 y& _( V    end
end

zqyzixin

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

昝斌

这是什么哦？

hjl19890710

好 dwadedcsargtg

hjl19890710

lihan厉害厉害厉害厉害厉害

coolistone

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

Kadyniost

？？？？？