转跳到内容

[数学算法挑战]趁有空,悄咪咪出个数列相关的数学问题


推荐贴

:mx040:防止没get到方法,补充了两个小问123263144_LF6DOK3X6PG78Q4FX87P.png.ec37673cb6f312cfe9d1663f9c185206.png

已知对于任意正整数n,有a(n+1)=a(n)^2-2
(I)a(1)=1,求a(n)
(II)a(1)=2,求a(n)
(III)a(1)=3,求a(n)

虽然说这是个非线性递推公式,不过构造技巧足够的话,也不需要什么特别高端的工具呢w

,由Tokur修改
注释
骚男 骚男 40.00节操 不愧是……社长www
链接到点评
  • Tokur将标题更改为[数学算法挑战]趁有空,悄咪咪出个数列相关的数学问题

代换a_n=2cosh(x_n)  (cosh(x)=(e^x+e^-x)/2)  得到一个式子(具体啥可以自己写下,没有latex太麻烦)

利用性质cosh(2x)=2*cosh(x)^2-1 得到cosh(x_{n+1})=cosh(2*x_n)  取同一支就有x_{n+1}=2*x_n

a1=3即e^x1+e^-x1=3,解得两根p,q

a_n=2cosh(2^n-1 *x1)=(e^x1)^(2^n-1)+(e^-x1)^(2^n-1)=p^(2^n-1)+q^(2^n-1)

事实上 对x_{n+1}=ax_n^2+bx_n+c,判别式为2b 或2b+8的都有显的结果 一种配方做 另一种三角换元

其他情况基本没有 可以看看Mandelbrot set

注释
骚男 骚男 10.00节操 对高中题来说,是降维打击了ww
链接到点评

:NEKOMIMI_PARADISE_7:……构造时意识到了一些不对劲(咱已经是多年没碰数学的废物辣w

不过倒是能构造出来

1/2 a_(n+1) = 1/2 a_n ^ 2 -1

然后就是楼上用到的三角换元 

不过完整来看,是要判断1/2 an 的范围,来决定用cosx 还是 coshx ……

:wn002:依稀记得这个-2系数是比较特殊的来着?

链接到点评
16 小时前,huancai77说道:

代换a_n=2cosh(x_n)  (cosh(x)=(e^x+e^-x)/2)  得到一个式子(具体啥可以自己写下,没有latex太麻烦)

利用性质cosh(2x)=2*cosh(x)^2-1 得到cosh(x_{n+1})=cosh(2*x_n)  取同一支就有x_{n+1}=2*x_n

a1=3即e^x1+e^-x1=3,解得两根p,q

a_n=2cosh(2^n-1 *x1)=(e^x1)^(2^n-1)+(e^-x1)^(2^n-1)=p^(2^n-1)+q^(2^n-1)

:mx040:好喔w 答案基本是对的,除了懒没解pq外233(x)p,q解出来是(3±\sqrt{5})/2,:b6:不过说起来咱想问问,ch(x_{n+1})=ch(2x_n)可以直接取x_{n+1}=2*x_n呢x 如果取x_{n+1}=-2*x_n似乎也是可以的 除了使用面向答案证明的数学归纳法,似乎无论是使用三角换元用复数表示还是双曲换元,都躲不开这个问题 @huancai77

16 小时前,huancai77说道:

事实上 对x_{n+1}=ax_n^2+bx_n+c,判别式为2b 或2b+8的都有显的结果 一种配方做 另一种三角换元

其他情况基本没有 可以看看Mandelbrot set

:wn007:是的w 不过有些情况是可以得出类常数列的结果w

,由Tokur修改

Tokur收和谐资源时被小萝莉围观良心发现失去-2节操

链接到点评
9 小时前,骚男说道:

:NEKOMIMI_PARADISE_7:……构造时意识到了一些不对劲(咱已经是多年没碰数学的废物辣w

不过倒是能构造出来

1/2 a_(n+1) = 1/2 a_n ^ 2 -1

然后就是楼上用到的三角换元 

不过完整来看,是要判断1/2 an 的范围,来决定用cosx 还是 coshx ……

:wn002:依稀记得这个-2系数是比较特殊的来着?

:NEKOMIMI_PARADISE_38:揉揉男酱w 其实这个例子是咱在学不动点构造的时候想到的

这个数列能解出不动点2和-1,不过也仅限于此,得到a(n+1)-2=(a(n)-2)(a(n)+2)这样的形式,后面就无法继续推进了

其实如果不限制使用复数表示的话,用cos和cosh都是可以的说,在前者会得到一个复数三角函数表示实数的一个实例ww

链接到点评
2 小时前,Tokur说道:

:NEKOMIMI_PARADISE_38:揉揉男酱w 其实这个例子是咱在学不动点构造的时候想到的

这个数列能解出不动点2和-1,不过也仅限于此,得到a(n+1)-2=(a(n)-2)(a(n)+2)这样的形式,后面就无法继续推进了

其实如果不限制使用复数表示的话,用cos和cosh都是可以的说,在前者会得到一个复数三角函数表示实数的一个实例ww

:NEKOMIMI_PARADISE_1: 龙娘,斯巴拉西!

:wn002: 问了咱搞数学的同学,如果不是-2的话,推不出正常的通项(大概这个意思w

:mx051: 那么龙娘会去参加竞赛这类的嘛  (有端联想,今天是高中数竞的复赛w

链接到点评
4 小时前,Tokur说道:

:mx040:好喔w 答案基本是对的,除了懒没解pq外233(x)p,q解出来是(3±\sqrt{5})/2,:b6:不过说起来咱想问问,ch(x_{n+1})=ch(2x_n)可以直接取x_{n+1}=2*x_n呢x 如果取x_{n+1}=-2*x_n似乎也是可以的 除了使用面向答案证明的数学归纳法,似乎无论是使用三角换元用复数表示还是双曲换元,都躲不开这个问题 @huancai77

:wn007:是的w 不过有些情况是可以得出类常数列的结果w

这是因为把sin 和 sinh分开本来就是一种委曲求全的作法 如果我们在复数域上看就有一种统一的思路

复数域上 sinz=(e^iz-e^-iz)/2i  cosz=(e^iz+e^-iz)/2

明显有z属于R时和原来的三角函数一样 z=-iR时有cosh(R)=cos(-iR) sinh(R)=isin(-iR)    (也就是说ch和原来的cos均是复的cos的一个截面)

此时,和ch的偶函数性质不同 只有cos(z)=cos(z+2kpi)的周期性,该周期性实际上是由e^z带来的(实际上那个负号就隐藏在周期性中)

那么,我们实际上可以令一种新的“数”z_0={所有z_0+2kpi}(这形成了一种等价关系),那实质上就把复数域用2pi长的带子割开了(这实际上是一个求取商集的过程)

那我们就可以做对应:域C 对应 新数 对应 代表元(0<=Im<2pi),这最后解释了为啥我们可以直接取0<=Im<2pi这样一个定义域(cos的一支)在实的就对应是取了0到正无穷那一段

 

注释
Tokur Tokur 66.00节操 太强惹ww
链接到点评
3 小时前,骚男说道:

:NEKOMIMI_PARADISE_1: 龙娘,斯巴拉西!

:wn002: 问了咱搞数学的同学,如果不是-2的话,推不出正常的通项(大概这个意思w

:mx051: 那么龙娘会去参加竞赛这类的嘛  (有端联想,今天是高中数竞的复赛w

:NEKOMIMI_PARADISE_38:咱在网上有查到相应的判别式什么的w

:NEKOMIMI_PARADISE_38:猜猜www

链接到点评
1 小时前,huancai77说道:

这是因为把sin 和 sinh分开本来就是一种委曲求全的作法 如果我们在复数域上看就有一种统一的思路

复数域上 sinz=(e^iz-e^-iz)/2i  cosz=(e^iz+e^-iz)/2

明显有z属于R时和原来的三角函数一样 z=-iR时有cosh(R)=cos(-iR) sinh(R)=isin(-iR)    (也就是说ch和原来的cos均是复的cos的一个截面)

此时,和ch的偶函数性质不同 只有cos(z)=cos(z+2kpi)的周期性,该周期性实际上是由e^z带来的(实际上那个负号就隐藏在周期性中)

那么,我们实际上可以令一种新的“数”z_0={所有z_0+2kpi}(这形成了一种等价关系),那实质上就把复数域用2pi长的带子割开了(这实际上是一个求取商集的过程)

那我们就可以做对应:域C 对应 新数 对应 代表元(0<=Im<2pi),这最后解释了为啥我们可以直接取0<=Im<2pi这样一个定义域(cos的一支)在实的就对应是取了0到正无穷那一段

 

:b10:唔姆w 前面的部分咱可以理解 查了一下,后面的部分是集合论的东西嘛ww 虽然说咱还没学到x,但是感觉很奇妙的说ww 商集什么的还没了解过x 倒是代表元在学初等数论的时候有了解到

,由Tokur修改

Tokur收和谐资源时被小萝莉围观良心发现失去-3节操

链接到点评
20 小时前,Tokur说道:

:NEKOMIMI_PARADISE_38:咱在网上有查到相应的判别式什么的w

:NEKOMIMI_PARADISE_38:猜猜www

 

20 小时前,Tokur说道:

:b10:唔姆w 前面的部分咱可以理解 查了一下,后面的部分是集合论的东西嘛ww 虽然说咱还没学到x,但是感觉很奇妙的说ww 商集什么的还没了解过x 倒是代表元在学初等数论的时候有了解到

:wn002: 快去学离散数学,可有意思了(

:mx051: 玩玩反正无妨ww

链接到点评

这个其实看形式就很容易想到,令an=x+1/x,那么an^2就=(x+1/x)^2=x^2+1/x^2+2。而an+1=an^2-2就恰好=x^2+1/x^2。所以{an}就是可以写成x^2+1/x^2的形式的东西。

那么当a1=a时,x+1/x = a

即x^2-a*x+1 = 0,用求根公式就能求出x^2是多少,再开个跟就是x。

同样当类似的形式,当an+1=an^2-2b,a1=a时,这个公式就是x^2-a*x+b = 0

所以题目中的三问的答案分别是:

(1).  {an} = ((1+sqrt(-3))/2)^(2^(n-1))+((1-sqrt(-3))/2)^(2^(n-1))由于2^n肯定不能整除3,所以当n>=2时这里的x和1/x的幅角肯定是3n+1或3n+2,再加上x的模长为1,所以这两个共轭负数的和必为-1,所以{an}={1, -1, -1...} 

(2).{an} = (1^(2^(n-1))+(1^(2^(n-1))=2

(3).{an} = ((3+sqrt(5))/2)^^(2^(n-1))+((3-sqrt(5))/2)^(2^(n-1))

 

同理,如果题目给的是an+1=an^2-1,那么就设an=sqrt(1/2)*x^(2^(n-1))+sqrt(1/2)/x^(2^(n-1))

如果题目给的是an+1=an^2-200,那么就设an=10*x^(2^(n-1))+10/x^(2^(n-1))

用这个方式也可以算:a1=114514,an+1=an^2-1919810时的通向公式,有想法的话也可以用计算器验证哦~

 

 

以下是修改增添内容(再给你们扩展一下这道题):

顺道一提,如果数列递推公式的形式是an+1=k*an^2+k*b的话,那么可以简单地把这个系数a给提取出来,通向公式就是:an=k^n*(x^(2^(n-1))+(-b/2x)^((2)^(n-1)),其中x为x^2-a1*x+b = 0的解。

例:an+1=3*an^2+8=3*an^2+3*8/3, a1=4,求{an}

解:令an=x-8/3x,则an+1=3*an^2+8=3*(x^2-(8/3x)^2)

故当a1=4时,

令a1=x-8/3x

x-8/3x=4

x^2-4x-8/3=0

x1x2=-8/3(韦达定理)

∴我们要求的x和8/3x其实就是方程的两根x1和x2

解得x1=2+sqrt(20/3)

        x2=2-sqrt(20/3)

所以an=3^n*((2+sqrt(20/3))^(2^n-1)-(2-sqrt(20/3))^(2^(n-1)))

,由makuwa修改
注释
骚男 骚男 30.00节操 虽然有些小瑕疵,有些小臭,但是糖w
链接到点评
  • 骚男锁定了本主题
游客
此主题已关闭。
×
×
  • 新建...

重要消息

为使您更好地使用该站点,请仔细阅读以下内容: 使用条款