数学建模社区-数学中国
标题:
逼近函数
[打印本页]
作者:
yaoyuan4102
时间:
2009-4-17 16:02
标题:
逼近函数
function f = Legendre(y,k,x0)
4 w3 ^( L) x! V
syms t;
- Z/ E8 S: O; `
P(1:k+1) = t;
- d: S+ Z% P- U+ x: i% h# {
P(1) = 1;
5 m- N$ y2 r% h, X* e# h( ]5 M
P(2) = t;
$ r) O6 e2 ?6 S) z) t" D* E" q
c(1:k+1) = 0.0;
# @+ @( \- s+ _% I
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
0 L) g) e6 q+ t9 U2 I
c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
, L! H, w: F Q- Q" t
f = c(1)+c(2)*t;
" r, H4 b- t. Q* S+ t$ Q
for i=3:k+1
2 b& w( ~4 U# \1 G, \( K9 M* l6 ?
P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
5 @( ?" s5 u9 E" a2 [0 e/ ~% [) Q
c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
* S+ I$ k/ J# L6 `: z" y$ U# C
f = f + c(i)*P(i);
" r& g+ B8 M# K: N, b/ p ~
# T: {, T+ q$ r6 G/ D3 S/ v0 n2 o
if(i==k+1)
- b% `- T0 }$ d) V
if(nargin == 3)
! ]: v3 \# ~) K1 D
f = subs(f,'t',x0);
- u+ r: e. u, X+ t) M* \" X" {) Q
else
# D( d" R1 C4 F$ i; R3 X# Z8 m
f = vpa(f,6);
" m% p5 o" i. N# m0 f
end
3 S# j, v$ g! S. m) r% i1 \; F
end
: }2 N8 c1 s4 k! X" u
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