找回密码
 注册
Simdroid-非首页
查看: 90|回复: 0

[1.符号] 简单的分段符号函数积分(方波傅立叶级数幅值计算)的速度问题

[复制链接]
发表于 2010-7-23 18:27:37 | 显示全部楼层 |阅读模式 来自 湖南湘潭


如图,设0< alpha[1]<Pi/2。波形关于Pi/2偶对称,关于Pi奇对称。
则分段函数及在其上定义傅立叶级数幅值函数分别为:

  1. restart;
  2. yfun := piecewise(t < alpha[1], 0,
  3.                   t < Pi-alpha[1], 1,
  4.                   t < Pi+alpha[1], 0,
  5.                   t < 2*Pi-alpha[1], -1,
  6.                   t <= 2*Pi, 0);
  7. a :=(n) ->int(yfun*cos(n*t), t = 0 .. 2*Pi))/Pi ;#omega = 1;
复制代码


定义函数时Maple不进行积分运算,而当执行a(n)就进行积分运算。

  1. a(n); # 输出 a(n)
复制代码


可以发现在maple中,这个积分运算很难算出!
但是手工推导却很简单。即a(n)  = 1/n*(4/Pi)*cos(n*alpha[1]), n = 2*k+1, k = 0,1,2,3,...(即为奇数)

这是什么问题呢,有无好的解决方法?
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-27 23:18 , Processed in 0.037480 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表