逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
% C) g2 F- P/ h$ G$ DP(1:k+1) = t;
P(1) = 1;
* b* c# o' D1 l5 `/ m2 NP(2) = t;
1 [1 f7 b) L) c6 ^c(1:k+1) = 0.0;
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
: A5 e7 h6 ^1 Y3 Y  qc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
& w1 J- f2 ^$ t" s. ]0 |f = c(1)+c(2)*t;
for i=3:k+1
, e  J( e5 b! {& h, L8 y' s$ M    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
  Y; D, [/ C$ ]$ ^0 ~( N    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
/ a+ w# D' C; m  n) o    f = f + c(i)*P(i);
   
    if(i==k+1)
9 _! O1 g5 W/ ^( C" P        if(nargin == 3)
6 D' u* c* k+ ?$ E, H2 z& @            f = subs(f,'t',x0);
        else
& {: w! Q* d4 M+ T" Y            f = vpa(f,6);
        end
3 q# P% S" L; x$ A1 ]    end
8 U. m7 M( z: xend

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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