逼近函数

yaoyuan4102

function f = Legendre(y,k,x0)
syms t;
  R! _4 v5 N7 v' ^" ]P(1:k+1) = t;
P(1) = 1;
P(2) = t;
/ S$ d4 d, [+ z. N# c' Rc(1:k+1) = 0.0;
  T' n/ x& N; L/ H3 ^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;
9 i+ V4 T) k; G6 N6 o6 o2 uf = 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);
    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
6 M0 P# T: i8 R2 v    f = f + c(i)*P(i);
# n* C# c8 D" y# N4 l, K5 Z4 ^   
    if(i==k+1)
+ e2 ?* V+ u+ B) c$ Y        if(nargin == 3)
            f = subs(f,'t',x0);
        else
7 T% J! G7 C9 \8 h) t            f = vpa(f,6);
        end
* ^6 Y8 B- F; v* ~- P    end
end

Kadyniost

？？？？？

coolistone

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

hjl19890710

lihan厉害厉害厉害厉害厉害

hjl19890710

好 dwadedcsargtg

昝斌

这是什么哦？

zqyzixin

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