数学建模社区-数学中国
标题:
逼近函数
[打印本页]
作者:
yaoyuan4102
时间:
2009-4-17 16:02
标题:
逼近函数
function f = Legendre(y,k,x0)
( D8 i% S6 s+ u+ l
syms t;
9 S* U6 w9 b9 d) Y: R
P(1:k+1) = t;
3 _/ e3 D1 i' k' _$ _* _
P(1) = 1;
0 \& G" Q$ r. T1 x5 i! u% U) U) c
P(2) = t;
' z7 q9 \6 Z' ^& U' @
c(1:k+1) = 0.0;
/ a t1 ^: ]7 S: D6 c& g
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
/ p6 c; E& s! f9 m! g/ x2 p
c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
6 r- p- l( R! `; M* f% e- Z7 y8 ?
f = c(1)+c(2)*t;
3 E1 \, F9 [" g1 N& c+ T
for i=3:k+1
( e$ T; N- v7 P1 u M6 P8 U' |
P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
! d) g, ? N, x+ Q: ~
c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
7 f6 n: b d# v% f- y: v' U
f = f + c(i)*P(i);
7 i9 y% a N4 }/ v6 y' B0 r
8 Y: Q4 q: l) G
if(i==k+1)
9 @+ i' A( K# V
if(nargin == 3)
( ]" K% A2 u( ~" o$ Y1 E
f = subs(f,'t',x0);
; A4 h' H( C+ v* g- t- V
else
8 f, j8 i O! R4 T P/ X
f = vpa(f,6);
. S* a0 c3 h% n4 M" E
end
/ W7 f5 \" z$ @6 S) H! |- u
end
- C, q D# x( c" k
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