- 积分
- 32
- 注册时间
- 2009-7-29
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 lin2009 于 2011-1-3 14:58 编辑
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的一般表达式。
- 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 |
|