QQ登录

只需要一步,快速开始

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

傅里叶级数(fourier series)展开的Matlab实现

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

2620

主题

162

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2015-3-12 15:35
  • 签到天数: 207 天

    [LV.7]常住居民III

    社区QQ达人 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组第六届国赛赛前冲刺培

    群组国赛讨论

    群组2014美赛讨论

    群组2014研究生数学建模竞

    群组数学中国试看培训视频

    跳转到指定楼层
    1#
    发表于 2014-9-18 12:25 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    傅里叶级数(fourier series)展开的Matlab实现
    MATLAB和maple语言均未直接提供求解fourier级数的系数的直接函数,好我们自己动手丰衣足食,下面提供了一个代码,大家可以参照下


    file:///http://latex.codecogs.com/svg.latex?12\,\sin \left( x \right) +3/2\,\sin \left( 2\,x \right) +4/9\,\sin \left( 3\,x \right) +3/16\,\sin \left( 4\,x \right) +{\frac {12}{125}}\,\sin \left( 5\,x \right) +1/18\,\sin \left( 6\,x \right) +{\frac {12}{343}}\,\sin \left( 7\,x \right) +{\frac {3}{128}}\,\sin \left( 8\,x \right) +{\frac {4}{243}}\,\sin \left( 9\,x \right) +{\frac {3}{250}}\,\sin \left( 10\,x \right) +{\frac {12}{1331}}\,\sin \left( 11\,x \right) +{\frac {1}{144}}\,\sin \left( 12\,x \right)给个例子说明下:将函数y=x*(x-pi)*(x-2*pi),在(0,2*pi)的范围内傅里叶级数展开

    syms x
    fx=x*(x-pi)*(x-2*pi);
    [an,bn,f]=fseries(fx,x,12,0,2*pi)%前12项展开
    latex(f)%将f转换成latex代码
    复制代码
    得到如下结果

    an =
    [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

    bn =
    [      -12,      3/2,     -4/9,     3/16,  -12/125,     1/18,  -12/343,    3/128,   -4/243,    3/250, -12/1331,    1/144]

    f =
    12*sin(x)+3/2*sin(2*x)+4/9*sin(3*x)+3/16*sin(4*x)+12/125*sin(5*x)+1/18*sin(6*x)+12/343*sin(7*x)+3/128*sin(8*x)+4/243*sin(9*x)+3/250*sin(10*x)+12/1331*sin(11*x)+1/144*sin(12*x)

    ans =
    12\,\sin \left( x \right) +3/2\,\sin \left( 2\,x \right) +4/9\,\sin \left( 3\,x \right) +3/16\,\sin \left( 4\,x \right) +{\frac {12}{125}}\,\sin \left( 5\,x \right) +1/18\,\sin \left( 6\,x \right) +{\frac {12}{343}}\,\sin \left( 7\,x \right) +{\frac {3}{128}}\,\sin \left( 8\,x \right) +{\frac {4}{243}}\,\sin \left( 9\,x \right) +{\frac {3}{250}}\,\sin \left( 10\,x \right) +{\frac {12}{1331}}\,\sin \left( 11\,x \right) +{\frac {1}{144}}\,\sin \left( 12\,x \right)
    复制代码
    使用论坛的LaTex代码显示如下





    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    深V礼        

    246

    主题

    24

    听众

    8718

    积分

  • TA的每日心情
    慵懒
    2014-11-8 19:08
  • 签到天数: 55 天

    [LV.5]常住居民I

    管理员

    超级版主

    版主

    发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    回复

    使用道具 举报

    张银银 实名认证       

    4

    主题

    13

    听众

    1268

    积分

    创造思考

    升级  26.8%

  • TA的每日心情
    开心
    2018-9-20 23:23
  • 签到天数: 394 天

    [LV.9]以坛为家II

    群组华南理工大学

    回复

    使用道具 举报

    0

    主题

    11

    听众

    41

    积分

    升级  37.89%

  • TA的每日心情
    开心
    2015-9-20 08:23
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    回复

    使用道具 举报

    2802

    主题

    160

    听众

    8882

    积分

  • TA的每日心情
    开心
    2017-4-26 10:25
  • 签到天数: 491 天

    [LV.9]以坛为家II

    自我介绍
    即使不开心也不要皱眉,因为你永远不知道有谁会爱上你的微笑!

    社区QQ达人 发帖功臣 新人进步奖 最具活力勋章

    群组数学中国试看培训视频

    群组2017美赛两天强训

    群组2015司守奎matlab培训

    群组2016国赛优秀论文解析

    群组国赛护航思路养成班

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2025-5-22 23:30 , Processed in 0.551668 second(s), 75 queries .

    回顶部