快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

逼近函数

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

2

主题

0

听众

14

积分

升级  9.47%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-4-17 16:02 |只看该作者 |正序浏览
|招呼Ta 关注Ta
function f = Legendre(y,k,x0)
4 z" ?" H$ f6 p0 C: Ysyms t;- A! p6 }! m" N* C
P(1:k+1) = t;9 g% E  ^$ ~5 A. [
P(1) = 1;
+ \8 h6 r* e: m4 o( @) t4 i5 JP(2) = t;
, o: Z6 s8 K9 o7 J" lc(1:k+1) = 0.0;3 y& c7 d  `0 |2 f
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;
: S. O( x) T! k: V& pc(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;
/ \# z$ c, |2 J. U% f/ Qf = c(1)+c(2)*t;
' R3 D4 ~' U9 ^9 Ffor i=3:k+1
3 j5 |( L. O  G' x7 C    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);
5 y  }% a! W9 s: ?% M    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;
) k" S  u/ m/ ?- m0 L& s  W    f = f + c(i)*P(i);; p* H9 \) |0 R$ I5 x
   
" C- A& a6 E3 q- A. n4 N    if(i==k+1)" E3 p- V3 E1 E" M( Z1 p; g, t4 P
        if(nargin == 3)
& i, r6 v* D$ Z" z4 v4 b            f = subs(f,'t',x0);3 ]; x7 e' P% h9 b
        else
: b8 ~; J2 R- N" t            f = vpa(f,6);' D" q3 B) \% ~8 Q; [0 c7 }
        end' e+ e; u, S* y$ K2 m
    end( O4 @3 l+ B4 r% \! `
end
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
回复

使用道具 举报

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
    呵呵，非常感谢，很好很强大
    回复

    使用道具 举报

    昝斌        

    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
    这是什么哦？
    回复

    使用道具 举报

    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
    回复

    使用道具 举报

    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厉害厉害厉害厉害厉害
    回复

    使用道具 举报

    coolistone        

    1

    主题

    3

    听众

    50

    积分

    升级  47.37%

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

    • 签到天数: 3 天

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

    使用道具 举报

    Kadyniost        

    0

    主题

    4

    听众

    187

    积分

    升级  43.5%

    该用户从未签到
    2#
    发表于 2009-8-9 23:57 |只看该作者 |招呼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-4 02:44 , Processed in 0.737490 second(s), 87 queries .

    回顶部