QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4697|回复: 6
打印 上一主题 下一主题

逼近函数

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

2

主题

0

听众

14

积分

升级  9.47%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2009-4-17 16:02 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
function f = Legendre(y,k,x0)  [8 D1 }, H7 M/ O0 b  O3 f
syms t;  y7 F; i% }( H6 V' ~% x$ q
P(1:k+1) = t;
; ?3 b1 v0 ~2 T0 R* y6 j7 l' iP(1) = 1;% Z/ j. x1 j3 a+ D9 f/ H+ }
P(2) = t;
2 d/ G$ p- I" L$ Q# U) ?c(1:k+1) = 0.0;9 X3 X: n; j& L7 Q! _" u* w$ v
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*P(1),t,-1,1)/2;& P, T6 R& S+ q! o3 _8 ?+ |
c(2)=int(subs(y,findsym(sym(y)),sym('t'))*P(2),t,-1,1)/2;2 n# @$ P# t* _: j+ Z1 k* d! p
f = c(1)+c(2)*t;
, `: G/ g3 G0 J* d9 v2 \. C! e( v9 {for i=3:k+1& O/ Y4 {8 H+ x* M: i1 H
    P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1);2 G# x' m( J7 w
    c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2;5 _, n1 x7 Z, G, w6 x& w
    f = f + c(i)*P(i);- [/ s& X4 @, S: p
   
, q9 E6 q: @; A6 e% J' j: W$ U    if(i==k+1)! Y2 S0 [6 B3 t5 z) K+ u0 m* L1 i6 J
        if(nargin == 3)
! N3 @/ F/ Y5 n$ h& C4 W! p  I, l3 s            f = subs(f,'t',x0);! {' l, z+ ~: E( }; b
        else2 n: @1 Z' ~" h: n6 _% m9 q
            f = vpa(f,6);
) [* B$ a; z6 `( P        end$ }+ c  u. U: ~4 o$ F9 J
    end' \6 B/ i+ }3 V4 I7 {  d% B
end
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
Kadyniost        

0

主题

4

听众

187

积分

升级  43.5%

该用户从未签到

回复

使用道具 举报

1

主题

3

听众

50

积分

升级  47.37%

  • TA的每日心情
    慵懒
    2013-4-11 14:40
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    回复

    使用道具 举报

    0

    主题

    4

    听众

    98

    积分

    升级  97.89%

  • TA的每日心情
    无聊
    2017-1-2 20:55
  • 签到天数: 49 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    0

    主题

    4

    听众

    98

    积分

    升级  97.89%

  • TA的每日心情
    无聊
    2017-1-2 20:55
  • 签到天数: 49 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    昝斌        

    15

    主题

    8

    听众

    1239

    积分

    升级  23.9%

  • TA的每日心情
    开心
    2014-8-4 18:33
  • 签到天数: 359 天

    [LV.8]以坛为家I

    自我介绍
    一个求学中的人

    社区QQ达人

    群组Matlab讨论组

    群组学术交流A

    群组数学建模

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

    群组中国科学院软件研究所

    回复

    使用道具 举报

    zqyzixin 实名认证       

    1

    主题

    5

    听众

    1818

    积分

    升级  81.8%

  • TA的每日心情
    难过
    2013-10-14 10:21
  • 签到天数: 78 天

    [LV.6]常住居民II

    社区QQ达人

    群组小草的客厅

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-3 07:48 , Processed in 0.432302 second(s), 87 queries .

    回顶部