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

[11.其他] 如何求出一无穷三角数列之和的收敛值?

[复制链接]
发表于 2011-1-3 14:29:43 | 显示全部楼层 |阅读模式 来自 湖南湘潭
本帖最后由 lin2009 于 2011-1-3 14:58 编辑

1) 下面无穷数列之和能否收敛,若能,收敛的确切值是多少?


上式是
  1. fabnv := proc (q, beta, nu) options operator, arrow; sin((1/2)*nu*beta*Pi)^2*sin((1/6)*nu*Pi)^2/(q^2*sin((1/6)*nu*Pi/q)^2*nu^2) end proc
复制代码
在beta=1,q=3时的特例。

maple不能得出该无穷数列的收敛值。但从下面的有限项之和可以看出
该无穷数列应该是收敛的,其值应在0.06879左右。

> evalf(subs({q = 3, beta = 1}, sum((sin((1/2)*nv*beta*Pi)*sin((1/6)*nv*Pi)/(q*sin((1/6)*nv*Pi/q)*nv))^2, nv = 2 .. 10000)));
print(`output redirected...`); # input placeholder
                         0.06878334966
> evalf(subs({q = 3, beta = 1}, sum((sin((1/2)*nv*beta*Pi)*sin((1/6)*nv*Pi)/(q*sin((1/6)*nv*Pi/q)*nv))^2, nv = 2 .. 20000)));
print(`output redirected...`); # input placeholder
                         0.06879168540
> evalf(subs({q = 3, beta = 1}, sum((sin((1/2)*nv*beta*Pi)*sin((1/6)*nv*Pi)/(q*sin((1/6)*nv*Pi/q)*nv))^2, nv = 2 .. 30000)));
print(`output redirected...`); # input placeholder
                         0.06879446199
> evalf(subs({q = 3, beta = 1}, sum((sin((1/2)*nv*beta*Pi)*sin((1/6)*nv*Pi)/(q*sin((1/6)*nv*Pi/q)*nv))^2, nv = 2 .. 40000)));
print(`output redirected...`); # input placeholder
                         0.06879585084

> evalf(subs({q = 3, beta = 1}, sum((sin((1/2)*nv*beta*Pi)*sin((1/6)*nv*Pi)/(q*sin((1/6)*nv*Pi/q)*nv))^2, nv = 2 .. infinity)));
                          原式输出

2)如何求出该无穷数列之和关于q和beta的一般表达式。

  1. f :=  (q, beta, nu) -> sin((1/2)*nu*beta*Pi)^2*sin((1/6)*nu*Pi)^2/(q^2*sin((1/6)*nu*Pi/q)^2*nu^2)
复制代码
已知q为正整数,2/3 <= beta <= 1
发表于 2011-1-4 12:39:56 | 显示全部楼层 来自 新加坡
Simdroid开发平台
v=18的时候,式子是不是无意义呀?
回复 不支持

使用道具 举报

 楼主| 发表于 2011-1-4 17:08:05 | 显示全部楼层 来自 湖南湘潭
2# zsq-w
v=18及18倍数时,f(q,beta,v)应为0;
单项计算时,会提示出错信息:
f(3, 1, 18)
Error, (in f) numeric exception: division by zero;

但在求和时Maple应该是按照极限值来计算的:
.
回复 不支持

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-4-28 08:36 , Processed in 0.034881 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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