数列a(n+1)=sqrt(S(n)^2-S(n)+2)，其中S(n)是数列前n项和，求a(n)的表达式。
数列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)比较，稳合。
這是高中的題吧？
页:
[1]