幾類特殊遞推數(shù)列的矩陣算法(前篇)
前言
擱置了許久,終于回來更新以兌現(xiàn)前周立下的flag了。
在學(xué)習(xí)數(shù)列時,我們會遇到這種遞推式,其中一種算法是用特征根,而特征根之前做過一期專欄寫了個人的奇妙理解方法,感興趣的可以參考:
而這期專欄,我們可以換種視角,從線性代數(shù)的角度來解決。
另外,在做一些數(shù)列的題中難免會遇到像,
這種惡心的數(shù)列。答案是通過"暴力"反復(fù)迭代得出的周期,不禁讓人摸不著頭腦!
"這誰能想到???"[惱]
那么命題人的心機(jī)何在呢?今兒就來一探究竟~
另外,此篇專欄會涉及到不少矩陣運(yùn)算的知識,萌新可以先參考3b1b的視頻講解:
【官方雙語/合集】線性代數(shù)的本質(zhì) - 系列合集
此篇文章不針對應(yīng)試,而是給數(shù)學(xué)愛好者們拓展下思路

整式的線性遞推
先來講講其概念,所謂線性遞推,也即(~表系數(shù))形式,這里的
次數(shù)均為1
比如為線性遞推,而像
這種有次數(shù)不為1的就不是線性遞推了
再來引入矩陣表示式以便良好地銜接
對于線性方程組,如:
我們想更直觀地表示與
的映射關(guān)系,可將其化為矩陣形式:
幾何意義上,即向量由向量
作線性變換
得到,其中矩陣
就類似于函數(shù)(function),只不過函數(shù)
的功能是將一個x映射為一個y,而矩陣的功能則是將一個向量映射為另一個向量。由于向量起點(diǎn)均默認(rèn)于原點(diǎn),故也可視為將一個坐標(biāo)(向量終點(diǎn))映射為另一個坐標(biāo)

先來看一道例題:已知,求通項(xiàng)公式?
由于筆者能力暫有限,想不到更好地銜接鋪墊了,所以提前劇透(
有了上述的鋪墊,如果我們能將遞推式化為該多好!
為什么好呢?因?yàn)檫@樣就有:,
我們便驚奇地發(fā)現(xiàn),這運(yùn)算不就和等比數(shù)列很相近么?
于是得到
因此,要求,只需對初值構(gòu)成的向量
作n次矩陣A的變換即可
而我們要求,那么就需作n-1次矩陣A的變換

回到題目,先選取已知遞推式,再選取另一條式子
寫成矩陣形式,即:
于是可形成遞推:
則
下面計算
對角化后,得:
代入化簡得:
即得:

如果是3階遞推又怎樣呢?
比如①
方法是一樣的,則考慮構(gòu)建3階矩陣:
則還需選?。?/p>
②
③
由①②③,寫成矩陣形式,即:
遞推即得:
后面就也是對角化處理了
不過求特征值時,特征多項(xiàng)式的次數(shù)跟矩陣的階數(shù)是對應(yīng)的,因此上述矩陣特征多項(xiàng)式就是3次方程。而上面的系數(shù)是隨便選取,因此解析解理論可求,只是可能不平凡。這也是出題時最多只出到二階線性遞推的主要原因。

推廣到任意階線性遞推
其中表系數(shù)
考慮構(gòu)建n階矩陣:
則還需選取:
寫成矩陣形式,即:
觀察系數(shù)矩陣可知,其有如下特點(diǎn):
第一行的數(shù)字即遞推式中對應(yīng)的由到
的系數(shù)
從第二行開始,從左往右依次寫1,其余寫0。也即從第1列第2行起,沿"左上--右下"走向處數(shù)字為1,其余為0
舉個例子當(dāng)練習(xí)
對應(yīng)矩陣形式遞推,即:

有了以上背景,我們便可以從更高觀點(diǎn)來證明這個奇葩式子了
對應(yīng)矩陣形式遞推,即:
遞推得:
下面求解
求得矩陣特征根為一對共軛復(fù)數(shù):
ps:特征值為復(fù)數(shù)則對應(yīng)伸縮+旋轉(zhuǎn)(模長對應(yīng)伸縮,輻角對應(yīng)旋轉(zhuǎn))矩陣運(yùn)算在復(fù)數(shù)域成立已有嚴(yán)格證明,這里就先不作拓展了
那么對角化后,則有:
其中
于是有
說明對向量作6次J的變換后回到原位,因此周期為6
到此悟性好的讀者腦海中已經(jīng)掀起了滔天巨浪!
這正是特殊的輻角具有旋轉(zhuǎn)周期性造成的!??!

原來這令人頭大的結(jié)論背后有如此絕妙的背景!
這便是數(shù)學(xué)!純真美妙的數(shù)學(xué)!不被應(yīng)試名利銅臭玷污的數(shù)學(xué)!

總結(jié)
這期專欄主要講解了用矩陣方法求解整式線性遞推數(shù)列,順帶發(fā)掘了一些奇葩結(jié)論背后美妙的命題背景。
其中矩陣次方的處理采用了矩陣對角化的方法。而這種方法對于二階矩陣而言需要滿足特征根為2不等實(shí)根或一對共軛復(fù)根時使用。換而言之,還存在特征根為重根的情況,這種情況矩陣不可對角化。由于個人尚未了解充分,因此先把尚未解決的這種特殊情況先遺留與此。后續(xù)掌握后再補(bǔ)充討論。
整式線性遞推寫完也花了不少篇幅了,因此一階分式線性遞推就放后續(xù)再講了~
另外,感覺缺少例題的講解,有空翻翻好像還沒丟的一輪書找找例題附上[滑稽]