数学建模社区-数学中国
标题:
逼近函数
[打印本页]
作者:
yaoyuan4102
时间:
2009-4-17 16:02
标题:
逼近函数
function f = Legendre(y,k,x0)
. s! O" e; D4 l3 y, y
syms t;
! w! n: e, D; J& t# [8 T9 F7 v- v/ E$ E
P(1:k+1) = t;
. @& ^' y: L# U& K: S9 k
P(1) = 1;
/ P. d' M" Z& n9 W( m8 ]& W
P(2) = t;
; d$ @+ Y+ F1 e; G+ T, Q
c(1:k+1) = 0.0;
4 o1 U) Z2 n1 N% [% T d) m
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
5 r$ x B. q0 T. k" M
c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
7 W/ {7 W9 ~- k3 ^6 \
f = c(1)+c(2)*t;
& {' g" {8 d# |) L
for i=3:k+1
7 o. h* }5 m. S0 t7 F% F
P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
8 @+ U( ?% A7 ]9 G) I/ K
c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
k( f: B+ A8 @( I8 y8 ^2 d
f = f + c(i)*P(i);
; j7 h9 [: [* p9 v) H
& }/ T/ D$ O0 d* Q( ~( f
if(i==k+1)
9 P! n) {; M* d( C' T+ h7 P
if(nargin == 3)
+ O* O Z3 K% o/ P: d
f = subs(f,'t',x0);
5 s4 x( p, y N! n. o8 d
else
( s- M- F. X3 ^0 O
f = vpa(f,6);
/ N6 P/ E& Q7 {' r4 C; I
end
) @2 ~4 `3 G8 y! ~
end
" u- X& k- ] L, X; T
end
作者:
Kadyniost
时间:
2009-8-9 23:57
？？？？？
作者:
coolistone
时间:
2012-1-31 17:52
顶顶顶顶顶顶顶顶顶顶顶顶的
作者:
hjl19890710
时间:
2012-2-15 20:27
lihan厉害厉害厉害厉害厉害
作者:
hjl19890710
时间:
2012-2-15 20:29
好 dwadedcsargtg
作者:
昝斌
时间:
2012-8-27 18:29
这是什么哦？
作者:
zqyzixin
时间:
2012-12-2 22:39
呵呵，非常感谢，很好很强大
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5