快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

留学>>
返回列表 发新帖
查看: 4338|回复: 6
打印 上一主题 下一主题

逼近函数

[复制链接]
字体大小: 正常 放大
yaoyuan4102        

2

主题

0

听众

14

积分

升级  9.47%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-4-17 16:02 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
function f = Legendre(y,k,x0)
6 T* }0 x, `( Osyms t;( `% ~6 J' e2 x
P(1:k+1) = t;
  i* [* @: O( U7 e# ?P(1) = 1;) l- V- \9 X' ~: W% o% W7 c. x
P(2) = t;8 p5 j( Q; A2 l, [6 k5 e) j$ \+ `
c(1:k+1) = 0.0;$ j, x- N; x- t- X1 d9 C  x6 Z
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
5 u* K7 s" a0 p5 ]* q5 Q: e# ?c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;0 Z0 y( K  d8 i4 K* ]: R
f = c(1)+c(2)*t;8 `7 F3 ~' g; ]  Z$ B
for i=3:k+16 `# n4 d  }' C4 ^) G& [/ A7 w
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
: A; C. t* q  G: C; i! T    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;9 i: l' A; N: K4 ], r: p
    f = f + c(i)*P(i);0 j0 i/ E* O$ h' K6 n
    " v% w9 g( b9 E' f; q
    if(i==k+1)( P9 B) p! k9 W4 s# S; c& |4 z
        if(nargin == 3)# a2 D) M+ R6 a$ e; d9 M# s
            f = subs(f,'t',x0);  k  D1 Q3 `7 D7 s2 ^: R
        else
/ O3 t: d) P" l, A$ t- G            f = vpa(f,6);0 |, ~3 y3 @9 }  l: d9 M/ k9 Z1 R
        end
7 B- H, Z3 [( Q0 @* ?7 h    end% J# z: d  i) a5 D( m, f" N- E3 s
end
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
回复

使用道具 举报

Kadyniost        

0

主题

4

听众

187

积分

升级  43.5%

该用户从未签到
2#
发表于 2009-8-9 23:57 |只看该作者 |招呼Ta 关注Ta
？？？？？
回复

使用道具 举报

coolistone        

1

主题

3

听众

50

积分

升级  47.37%

  • TA的每日心情
    慵懒
    2013-4-11 14:40

    • 签到天数: 3 天

    [LV.2]偶尔看看I
    3#
    发表于 2012-1-31 17:52 |只看该作者 |招呼Ta 关注Ta
    顶顶顶顶顶顶顶顶顶顶顶顶的
    回复

    使用道具 举报

    hjl19890710        

    0

    主题

    4

    听众

    98

    积分

    升级  97.89%

  • TA的每日心情
    无聊
    2017-1-2 20:55

    • 签到天数: 49 天

    [LV.5]常住居民I
    4#
    发表于 2012-2-15 20:27 |只看该作者 |招呼Ta 关注Ta
    lihan厉害厉害厉害厉害厉害
    回复

    使用道具 举报

    hjl19890710        

    0

    主题

    4

    听众

    98

    积分

    升级  97.89%

  • TA的每日心情
    无聊
    2017-1-2 20:55

    • 签到天数: 49 天

    [LV.5]常住居民I
    5#
    发表于 2012-2-15 20:29 |只看该作者 |招呼Ta 关注Ta
    好 dwadedcsargtg
    回复

    使用道具 举报

    昝斌        

    15

    主题

    8

    听众

    1239

    积分

    升级  23.9%

  • TA的每日心情
    开心
    2014-8-4 18:33

    • 签到天数: 359 天

    [LV.8]以坛为家I

    自我介绍
    一个求学中的人

    社区QQ达人

    群组: Matlab讨论组

    群组: 学术交流A

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等

    群组: 中国科学院软件研究所
    6#
    发表于 2012-8-27 18:29 |只看该作者 |招呼Ta 关注Ta
    这是什么哦？
    回复

    使用道具 举报

    zqyzixin 实名认证       

    1

    主题

    5

    听众

    1818

    积分

    升级  81.8%

  • TA的每日心情
    难过
    2013-10-14 10:21

    • 签到天数: 78 天

    [LV.6]常住居民II

    社区QQ达人

    群组: 小草的客厅
    7#
    发表于 2012-12-2 22:39 |只看该作者 |招呼Ta 关注Ta
    呵呵，非常感谢，很好很强大
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    B Color Image Link Quote Code Smilies
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩

    • 电话咨询

    • 04714969085

    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问：王兆丰

    GMT+8, 2025-8-2 14:26 , Processed in 0.956755 second(s), 86 queries .

    回顶部