找回密码
 欢迎注册
查看: 583|回复: 1

[原创] 数列a(n+1)=sqrt(S(n)^2-S(n)+2),其中S(n)是数列前n项和,求a(n)的表达式。

[复制链接]
发表于 2023-9-16 21:43:35 来自手机 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
数列a(n+1)=sqrt(S(n)^2-S(n)+2),其中S(n)是数列前n项和,求a(n)的表达式。

得到:a(n)=a3*sin(x)/sin(x/2^(n-3)) n>=3。
解:
如果a(1)=sqrt(2),可以构造三角形利用余弦定理,有:
a(n)=sqrt(14)/2/sin((pi-asin(sqrt(14)/4))/2^(n-1))。
这题对于任意a(1)都可以有解析表达式的。 可以写出连乘的出来。
先去掉S(n):
a(n+1)^2=S(n)^2-S(n)+2
a(n+2)^2=S(n+1)^2-S(n+1)+2
a(n+2)^2-a(n+1)^2=(S(n+1)+S(n))*a(n+1)-a(n+1)
=(2S(n+1)-a(n+1))*a(n+1)-a(n+1)
所以:
2S(n+1)=(a(n+2)^2+a(n+1))/a(n+1)
2S(n+2)=(a(n+3)^2+a(n+2))/a(n+2)
2a(n+2)=a(n+3)^2/a(n+2)-a(n+2)^2/a(n+1)
a(n+3)^2/a(n+2)^2=a(n+2)/a(n+1)+2
⚠️⚠️⚠️注意这个等式成立条件是n>=1
除了a1,a(n)都大于0。 可以验证一般情况下a(3)^2/a(2)^2 != a(2)/a(1)+2。
a(1)=sqrt(2)/4的时候,这个等式刚好成立。
设a(n+1)/a(n)=b(n),有:
b(n+1)^2=b(n)+2
b(n+1)=sqrt(2+b(n))
设b(1)=2*cosx
x=acos(b1/2)
cosx=2(cos(x/2))^2-1
所以:
b(n)=2cos(x/2^(n-1))
有:
a1=a1
a2=sqrt(a1*a1-a1+2)
a3=sqrt((a1+a2)^2-(a1+a2)+2)
b1=a(3)/a(2)
x=acos(b1/2)
a(n)=a2* 2^(n-2)*∏cos(x/2^(i-1)
(i取1到n-2) n>=3
a(n)=a3*sin(x)/sin(x/2^(n-3)) n>=3
验证Python代码:
from math import *
rn="\r\n"
num=10
a1=sqrt(2)/2
s=a1
a2=sqrt(s*s-s+2)
s+=a2
a3=sqrt(s*s-s+2)
s+=a3
b1=a3/a2
b=acos(b1/2)
print(rn,"a1,a2,a3=",rn,a1,a2,a3,rn)
print(rn,a2/a1,a3/a2,sqrt(a2/a1+2),rn)
w=pi-acos(sqrt(2)/4)
x=acos(b1/2)
a=a3
for i in range(num):
an=sqrt(s*s-s+2)
s=s+an

b=b/2
c=2*cos(b)
print(an/a,c)

a=an
print(rn,a,a3*sin(x)/sin(x/2**(i+1)))
# print(a1*sqrt(14)/4/sin(w/2**(i+3)))


an是按递推公式计算得出来的,和an公式得到的吻合。a(n+1)/a(n)比值和公式b(n)比较,稳合。
C85036AE-868D-49CB-B4A6-D276B9E80F4A.jpeg
5A035729-BD67-4963-96E2-0DD7D462F50E.png
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-9-17 06:37:16 | 显示全部楼层
這是高中的題吧?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-5-3 00:46 , Processed in 0.056342 second(s), 25 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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