特征方程法求解

合集下载

高阶常微分方程的解法

高阶常微分方程的解法

高阶常微分方程的解法高阶常微分方程是数学中重要的研究对象,它描述了自变量以及其各阶导数之间的关系。

解决这类方程对于理解自然科学以及工程应用都具有重要意义。

本文将介绍一些常见的高阶常微分方程的解法。

一、分离变量法分离变量法是解决一阶常微分方程常用的方法,对于高阶常微分方程也同样适用。

首先,我们需要将方程重写为关于各阶导数的方程。

例如,考虑一个形如 $y''(x) + p(x)y'(x) + q(x)y(x) = 0$ 的二阶常微分方程。

我们可以将其表示为 $D^2y + pDy + qy = 0$,其中 $D =\frac{d}{dx}$ 是微分算子。

接下来,我们可以假设解为 $y(x) =u(x)v(x)$ 的形式,通过代入原方程并合并同类项,可以得到一阶关于$u$ 和 $v$ 的方程组。

然后我们利用这个方程组进行求解,最终得到$y$ 的解。

二、特征方程法特征方程法常用于线性常系数齐次高阶常微分方程的求解。

这类方程的特点是方程中仅包含自变量 $x$ 及其各阶导数,没有出现 $y$ 和它的导数。

我们以二阶方程 $y''(x) + py'(x) + qy(x) = 0$ 为例,首先使用代数技巧将方程转化为特征方程 $r^2 + pr + q = 0$。

通过求解特征方程得到两个根 $r_1$ 和 $r_2$,从而得到方程的通解形式 $y(x) =C_1e^{r_1x} + C_2e^{r_2x}$,其中 $C_1$ 和 $C_2$ 是常数。

三、常数变易法常数变易法是一种适用于非齐次高阶常微分方程的解法。

对于方程$D^ny = f(x)$,我们可以先求解对应的齐次方程 $D^ny = 0$ 的通解$y_c(x)$,接着我们假设非齐次方程的解为 $y_p(x)$,通过代入原方程并解得非齐次方程的特解 $y_p(x)$,最终得到原方程的通解 $y(x) =y_c(x) + y_p(x)$。

常系数线性微分方程

常系数线性微分方程

常系数线性微分方程常系数线性微分方程是微分方程中一类重要的特殊形式,其特点是方程中的系数是常数。

本文将介绍常系数线性微分方程的定义、求解方法以及相关性质。

一、常系数线性微分方程的定义常系数线性微分方程又称为齐次线性微分方程,其一般形式为:\[a_ny^{(n)}+a_{n-1}y^{(n-1)}+...+a_1y'+a_0y=0\]其中,n为方程的阶数,\(y^{(n)}\)表示y的n阶导数。

二、常系数线性微分方程的求解方法1. 特征方程法通过设定方程的解为\(y=e^{mx}\),将其代入原方程中,得到特征方程:\[a_nm^n+a_{n-1}m^{n-1}+...+a_1m+a_0=0\]解特征方程,可得到n个不同的解,分别是\(m_1, m_2,..., m_n\)。

则原方程的通解为:\[y=c_1e^{m_1x}+c_2e^{m_2x}+...+c_ne^{m_nx}\]其中,\(c_1, c_2,..., c_n\)为常数。

2. 变量分离法对于一些特殊的常系数线性微分方程,可以通过变量转换将其化为可分离变量的形式,从而简化求解过程。

三、常系数线性微分方程的性质1. 零解的存在唯一性对于常系数线性微分方程,其零解必然存在且唯一。

2. 齐次性质如果y1(x)是常系数线性微分方程的一个解,那么ky1(x)(k为常数)也是该微分方程的解。

3. 叠加性质如果y1(x)和y2(x)分别是常系数线性微分方程的解,那么y(x)=y1(x)+y2(x)也是该微分方程的解。

4. 线性性质设y1(x)和y2(x)分别是齐次常系数线性微分方程的两个解,c1和c2为常数,则c1y1(x)+c2y2(x)也是该微分方程的解。

总结:常系数线性微分方程作为微分方程中的重要形式,在工程、物理学以及其他科学领域中具有广泛的应用。

求解常系数线性微分方程的方法多种多样,特征方程法和变量分离法是常用的求解方法。

同时,常系数线性微分方程满足一系列重要性质,这些性质使得我们可以更加灵活地利用微分方程进行问题的建模和求解。

微分方程中的常系数齐次线性方程求解

微分方程中的常系数齐次线性方程求解

微分方程中的常系数齐次线性方程求解在微积分学中,常系数齐次线性方程是一类常见的微分方程。

它们的解可以通过一定的方法得到。

在本文中,我们将介绍如何求解常系数齐次线性方程。

一、什么是常系数齐次线性方程常系数齐次线性方程是指形如y″+ay′+by=0的微分方程,其中a和b为常数。

它们的特点是方程中的未知函数及其导数的系数都是常数。

二、求解常系数齐次线性方程的方法1. 特征方程法特征方程法是求解常系数齐次线性方程的一种常用方法。

具体步骤如下:(1)写出微分方程的特征方程,特征方程就是对应的代数方程。

对于y″+ay′+by=0,其特征方程为r²+ar+b=0。

(2)解特征方程,求得特征根。

设特征根为r₁和r₂,则特征方程的解为r₁和r₂。

根的个数和重根的情况会影响方程的解形式。

(3)根据特征根求解原方程的解。

当r₁和r₂为不同的实根时,原方程的通解可以表示为y=C₁e^(r₁x)+C₂e^(r₂x),其中C₁和C₂为常数。

当r₁和r₂为不同的复数根时,通解可以表示为y=e^(αx)(C₁cos(βx)+C₂sin(βx)),其中α为实部,β为虚部。

2. 代入法代入法也是一种常用的求解常系数齐次线性方程的方法。

具体步骤如下:(1)设定未知函数的形式。

根据方程的阶数,设定未知函数的形式,如y=e^(mx)。

(2)将未知函数及其导数带入微分方程,消去常数,得到相应的代数方程。

(3)解代数方程,得到未知函数的表达式。

根据代数方程的解,确定未知函数的形式。

(4)确定未知函数的常数。

根据给定的初始条件,确定未知函数中的常数值。

3. 傅里叶级数法对于特定的边界条件,常系数齐次线性方程还可以通过傅里叶级数法进行求解。

该方法主要适用于周期性边界条件的问题。

三、实例分析为了更好地理解求解常系数齐次线性方程的方法,我们来看一个具体的实例。

例题:求解方程y″+3y′+2y=0.解法:首先写出特征方程r²+3r+2=0,解得特征根r₁=-1,r₂=-2.特征根不相等,所以方程的通解为y=C₁e^(-x)+C₂e^(-2x)。

特征方程特征根法求解数列通项公式

特征方程特征根法求解数列通项公式

特征方程特征根法求解数列通项公式
1、将数列的前两项给出,在此基础上推雅可比数列,得到数列的递
推公式;
2、将递推公式化为特征方程,且特征方程只包含未知数x;
3、求解特征方程的特征根,得到特征根为{r1,r2,…,rm};
4、使用特征根构造数列的通项公式:利用特征根构造出原数列的通
项公式,即an = A1*r1^(n-1) + A2*r2^(n-1) + … + Am*rm^(n-1)(此
处n>=1);
5、求解参数A1,A2,…,Am,即将特征根对应的数列项代入原数列,解方程组求出所有参数;
6、给出最终的数列通项公式:将前面求出的所有参数代入数列通项
公式中,得到最终的数列通项公式。

二、实例演示
下面以解决下列特征方程求数列的通项公式为例,详细介绍特征方程
特征根法的求解:
原特征方程:x^2-x-6=0;特征根:r1=3,r2=2;推出数列:a1=4,
a2=10;
求数列通项公式:
1、根据特征方程求出特征根:
原特征方程:x^2-x-6=0;
解之,得:x=3,2;
即特征根为r1=3,r2=2;。

计算方法之计算矩阵的特征值和特征量

计算方法之计算矩阵的特征值和特征量

计算方法之计算矩阵的特征值和特征量计算矩阵的特征值和特征向量是线性代数中的一个重要问题,它在科学研究和工程应用中有着广泛的应用。

本文将介绍计算矩阵特征值和特征向量的方法,包括特征方程法、幂法、反幂法和QR方法。

一、特征值和特征向量的定义给定一个n阶方阵A,如果存在一个非零向量x和一个标量λ,满足以下方程:Ax=λx其中,x被称为A的特征向量,λ被称为A的特征值。

二、特征方程法特征方程法是计算矩阵特征值和特征向量的一种常用方法,其基本思想是通过求解矩阵的特征方程来求得特征值。

对于一个n阶方阵A,其特征方程为:A-λI,=0其中,I是n阶单位矩阵,A-λI,表示A-λI的行列式。

解特征方程可以得到n个特征值λ₁,λ₂,...,λₙ。

然后,将这些特征值带入原方程组(A-λI)x=0,求解线性方程组得到n个特征向量x₁,x₂,...,xₙ。

三、幂法幂法是一种通过迭代来计算矩阵最大特征值和对应的特征向量的方法。

首先,随机选择一个非零向量b₀,并进行归一化,得到单位向量x₀=b₀/,b₀。

然后,通过迭代的方式,计算xₙ₊₁=Axₙ,其中xₙ为第k次迭代得到的向量。

在迭代过程中,向量xₙ的模长会逐渐趋近于最大特征值对应的特征向量。

当迭代收敛后,xₙ就是矩阵A的最大特征值对应的特征向量。

四、反幂法反幂法是一种通过迭代来计算矩阵最小特征值和对应的特征向量的方法。

首先,随机选择一个非零向量b₀,并进行归一化,得到单位向量x₀=b₀/,b₀。

然后,通过迭代的方式,计算xₙ₊₁=(A-σI)⁻¹xₙ,其中σ为待求的特征值。

在迭代过程中,向量xₙ的模长会逐渐趋近于特征值σ对应的特征向量。

当迭代收敛后,xₙ就是矩阵A的特征值为σ的特征向量。

五、QR方法QR方法是一种通过迭代来计算矩阵特征值和特征向量的方法。

首先,将矩阵A进行QR分解,得到矩阵A=QR,其中Q是正交矩阵,R是上三角矩阵。

然后,计算矩阵B=RQ,重复以上步骤,直到矩阵B收敛。

特征值与特征向量的求解

特征值与特征向量的求解

特征值与特征向量的求解特征值和特征向量是线性代数中一对重要的概念,广泛应用于物理学、工程学和计算机科学等领域。

在本篇文章中,我们将深入探讨特征值和特征向量的定义、性质以及求解方法。

一、特征值与特征向量的定义在介绍特征值与特征向量的求解方法之前,我们先来了解它们的定义。

在一个n维向量空间V中,若存在一个n阶方阵A和一个非零向量X,使得下式成立:AX = λX其中,λ为标量,称为矩阵A的特征值;X为矩阵A的特征向量。

特征值与特征向量的求解方法有多种,下面将介绍其中两种常用的方法。

二、特征值与特征向量的求解方法1. 特征方程法特征方程法是求解特征值和特征向量的一种常用方法。

假设A是一个n阶方阵,我们的目标是求解它的特征值和特征向量。

首先,我们将上述特征方程AX = λX两边同时左乘一个单位矩阵I,得到:(A-λI)X = 0其中,I为n阶单位矩阵,0为n维零向量。

由于X为非零向量,所以矩阵(A-λI)必须是奇异矩阵,即其行列式为0:|A-λI| = 0这就是特征方程。

接下来,我们需要求解特征方程|A-λI| = 0的根λ,即矩阵A的特征值。

求解得到的特征值λ可以有重根。

然后,将每个特征值λ带入原特征方程(A-λI)X = 0,解得对应的特征向量X。

注意,对于每个不同的特征值,我们都可以对应多个特征向量。

通过特征方程法,我们可以求解出矩阵A的所有特征值和对应的特征向量。

2. 幂法幂法是求解矩阵特征值和特征向量的一种迭代方法,适用于大规模稀疏矩阵。

幂法的基本思想是:通过迭代将初始向量不断与矩阵A进行乘法运算,使得向量的模不断增大,趋向于对应最大特征值的特征向量。

具体做法是:1) 先选择一个非零向量X0作为初始向量。

2) 迭代计算X(k+1) = AX(k),其中k表示迭代次数。

3) 归一化向量X(k+1),即X(k+1) = X(k+1) / ||X(k+1)||,其中||X(k+1)||表示向量X(k+1)的模。

特征方程法求解递推关系中的数列通项

特征方程法求解递推关系中的数列通项

特征方程法求解递推关系中的数列通项一、(一阶线性递推式)设已知数列{a n }的项满足a j = b,a n 4 = ca n • d ,其中c = 0, c = 1,求这个数列的通项公式。

采用数学归纳法可以求解这一问题,然而这样做太过繁琐,而且在猜想通项公式中容易出错,本文提出一种易于被学生掌握的解法一一特征方程 法:针对问题中的递推关系式作出一个方程 x =cx • d,称之为特征方程;借助这个特征方程的根快速求解通项公式•下面以定理形式进行阐述.定理1:设上述递推关系式的特征方程的根为 x 0,则当x 0 = a 4时,a n为常数列,即a n 二a i ;当X o 二a i 时,a^ b n ' x o ,其中{b n }是以c 为公比 的等比数列,即 b n = b 4c n J,b 4 =a 4-x 0.pl证明:因为c = 0,1,由特征方程得x 0——.作换元b n = a n - x 0,贝U 1 -c n 1当X 。

=a 1时,b 1 =0 ,数列{b n }是以c 为公比的等比数列, 故b n =b1C _; 当 x ° 二a 1 时,d =0 , {b n }为 0 数列,故 a * =a 1,n • N.(证毕) 下面列举两例,说明定理 1的应用.1例1•已知数列{a n }满足:a n^^a -2,- N,a—,求a n.13 解:作方程x x -2,则x 0. 3 2b"a n「x0 © d—注乂a .cd1 -c二 c(a n -X °) = cb n . 11一2 -3 一2 +X — a-fl等的比公为11 1 n4丁 3) ,a n-3b n —3叫-」)n‘, n N. 2 2 2 3b n列是例2.已知数列{a n}满足递推关系:a n ^(2a n - 3)i, n,N,其中i为虚数3单位。

当a i 取何值时,数列{a .}是常数数列?a^ :-,a 2二:给出的数列:a n 爲方程x 2- px -q =0,叫做数列 :a n / 的特征方程。

10第十讲:多自由度系统特征值问题求解方法(I)

10第十讲:多自由度系统特征值问题求解方法(I)

k 23 k 2 2k1
k3=245MN/m
k 33 k 2 k 3 2k1 2.5k1 4.5k1
2
2015/12/7
第十讲:多自由度系统特征值问题求解方法(I) 第十讲:多自由度系统特征值问题求解方法 (I)
一、特征方程法(Characteristic Equation Method)
φ (3) 1 .000 , 1 . 247 ,
0.802
0 . 555
T
T
节点
k3=245MN/m
频率 方程 为
k11 ω2 j m1 k21 k31
k12 k22 ω2 j m2 k32
k13 k23 0 k33 ω2 m j 3
第十讲:多自由度系统特征值问题求解方法(I) 第十讲:多自由度系统特征值问题求解方法 (I)
频率方程:
K 2 M k1 m1 2 1 k 1 1 0 1 m1 2 3 1. 5 k 1 2 0 2 m 4 . 5 1 .5 1 2 k1 0
x1 : x2 : : xn A1 : A2 : : An
(K i M ) A ( i ) 0
i 1, 2, , n
由于上述方程组的系数行列式为零,该方程是降秩的;对于没有 重根的情况,只有n-1个方程是独立的。
并不随时间而变化,也就是说在任何时刻结构的振动都保持同一形状, 整个结构就像一个单自由度结构一样在振动。这种多自由度结构按任一 自振频率 k 进行的简谐振动称为主振动,与其相应的特定振动形式称为 主振型(振型) 无阻尼自由振动的解为
(2)系统的运动方程由刚度矩阵表示 设系统的广义坐标列向量为x,则运动微分方程为 为二阶齐次方程 按微分方程理论,设齐次方程的特解为 x j A j sin(ωt θ ) j 1,2, ,n

第八讲 特征值问题求解方法

第八讲 特征值问题求解方法

同济大学土木工程防灾国家重点实验室、桥梁工程系
小结
特征方程法—精确方法、求全部特征值、适用于低阶矩阵
相似变换法—精确方法、求全部特征值、适用于高阶矩阵
同济大学土木工程学院桥梁工程系
高等结构动力学
主讲教师:葛耀君 博士 教授 曹曙阳 博士 教授
同济大学 土木工程防灾国家重点实验室 桥梁工程系
同济大学土木工程学院桥梁工程系
第七章 特征值问题求解
第一节 特征方程法 第二节 相似变化法 第三节 矢量迭代法 第四节 Sturm 排序法
同济大学土木工程防灾国家重点实验室、桥梁工程系
经过多次旋 转正交使矩 阵对角化
1954年以前唯一的相似变换法 严格按照相似变换关系计算出矩阵[N ]的元素 通过相似变换使矩阵[A ] 为对角矩阵
(2) Givens 三对角化方法(Givens Tridiagonalization)
坐标旋转变 换,使矩阵 三对角化
通过相似变换使矩阵[A ] 为三对角元素矩阵
{v }=ω [D ]{v }
(1) 1 2 1 (0) 1
忽略频率影响后的一次近似值
{v }= [D]{v }
(1) 1 ( 0) 1
同济大学土木工程防灾国家重点实验室、桥梁工程系
3.1 一阶振型迭代 (续)
一次近似值标准化
{ v } {v }= ({ }) ref v
(1) 1 (1) 1 (1) 1
同济大学土木工程防灾国家重点实验室、桥梁工程系
2.2 相似变换法分类(续)
(5) LR变换法(LR Transformation Method) 通过相似变换使矩阵[A ] 为上三角元素矩阵 相似变换公式:[A ] = [L ][A][R ] 改进特点:不必求逆运算 (6) QR或QL变换法(QR or QL Transformation Method) 每次相似变换使矩阵[A ] 右下角4元素分离,即

三阶常系数线性非齐次微分方程特解的两种解法

三阶常系数线性非齐次微分方程特解的两种解法

三阶常系数线性非齐次微分方程特解的两种解法一种解法是求解特征方程,另一种解法是采用逐步求解法。

1、求解特征方程法:设三阶常系数线性非齐次微分方程为:y'''+a2y''+a1y'+a0y=f(x)其中a2,a1,a0为常数,f(x)为右端函数。

(1)求解特征方程:设特征根为λ1,λ2,λ3,则特征方程为:λ3+a2λ2+a1λ+a0=0求解特征方程,得到特征根:λ1,λ2,λ3(2)求解特解:令特解为y=C1eλ1x+C2eλ2x+C3eλ3x代入方程,得:C1eλ1x+C2eλ2x+C3eλ3x+a2(C1eλ1x+C2eλ2x+C3eλ3x)+a1(C1eλ1x+C2eλ2x+C3eλ3x)+a0(C1eλ1x+C2eλ2x+C3eλ3x)=f(x)即:(C1λ3+C2λ2+C3λ1)eλ1x+(C1λ3+C2λ2+C3λ1)eλ2x+(C1λ3+C2λ2 +C3λ1)eλ3x=f(x)化简得:C1λ3+C2λ2+C3λ1=f(x)解得:C1=f(x)λ3/(λ3-λ2)(λ3-λ1)C2=f(x)λ2/(λ2-λ1)(λ2-λ3)C3=f(x)λ1/(λ1-λ2)(λ1-λ3)故特解为:y=f(x)λ3/(λ3-λ2)(λ3-λ1)eλ1x+f(x)λ2/(λ2-λ1)(λ2-λ3)eλ2x+f(x)λ1/(λ1-λ2)(λ1-λ3)eλ3x2、逐步求解法:设三阶常系数线性非齐次微分方程为:y'''+a2y''+a1y'+a0y=f(x)(1)求解一阶线性微分方程:设y1(x)为一阶线性微分方程的解,则有:y1'+a2y1=0解得:y1=C1e-a2x(2)求解二阶线性微分方程:设y2(x)为二阶线性微分方程。

特别解析特征方程法求解递推关系中的数列通项

特别解析特征方程法求解递推关系中的数列通项

特别解析:特征方程法求解递推关系中的数列通项一、一阶线性递推式设已知数列}{n a 的项满足d ca a b a n n +==+11,,其中,1,0≠≠c c 求这个数列的通项公式;定理1:设上述递推关系式的特征方程的根为0x ,则当10a x =时,n a 为常数列,即0101,;x b a a x a a n n n +===时当,其中}{n b 是以c 为公比的等比数列,即01111,x a b c b b n n -==-.证明:因为,1,0≠c 由特征方程得.10cdx -=作换元,0x a b n n -=则.)(110011n n n n n n cb x a c ccdca c d d ca x a b =-=--=--+=-=--当10a x ≠时,01≠b ,数列}{n b 是以c 为公比的等比数列,故;11-=n n c b b当10a x =时,01=b ,}{n b 为0数列,故.N ,1∈=n a a n 证毕例1.已知数列}{n a 满足:,4,N ,23111=∈--=+a n a a n n 求.n a解:作方程.23,2310-=--=x x x 则 当41=a 时,.21123,1101=+=≠a b x a 数列}{n b 是以31-为公比的等比数列. 于是:.N ,)31(2112323,)31(211)31(1111∈-+-=+-=-=-=---n b a b b n n n n n n例2.已知数列}{n a 满足递推关系:,N ,)32(1∈+=+n i a a n n 其中i 为虚数单位;当1a 取何值时,数列}{n a 是常数数列 解:作方程,)32(i x x +=则.5360i x +-=要使n a 为常数,即则必须.53601ix a +-== 二、二阶线性递推式定理2:对于由递推公式n n n qa pa a +=++12,βα==21,a a 给出的数列{}n a ,方程02=--q px x ,叫做数列{}n a 的特征方程;若21,x x 是特征方程的两个根,当21x x ≠时,数列{}n a 的通项为1211--+=n n n Bx Ax a ,其中A,B 由βα==21,a a 决定即把2121,,,x x a a 和2,1=n ,代入1211--+=n n n Bx Ax a ,得到关于A 、B 的方程组;当21x x =时,数列{}n a 的通项为11)(-+=n n x B A a ,其中A,B 由βα==21,a a 决定即把2121,,,x x a a 和2,1=n ,代入11)(-+=n n x Bn A a ,得到关于A 、B 的方程组;例3:已知数列{}n a 满足),0(0253,,1221N n n a a a b a a a n n n ∈≥=+-==++,求数列{}n a 的通项公式;解法一待定系数、迭加法由025312=+-++n n n a a a ,得)(32112n n n n a a a a -=-+++, 且a b a a -=-12;则数列{}n n a a -+1是以a b -为首项,32为公比的等比数列, 于是:11)32)((-+-=-n n n a b a a ;把n n ,,3,2,1⋅⋅⋅=代入,得:a b a a -=-12, )32()(23⋅-=-a b a a , ••• ,21)32)((---=-n n n a b a a ;把以上各式相加,得:])32()32(321)[(21-+⋅⋅⋅+++-=-n n a b a a )(321)32(11a b n ---=-; a b b a a a b a n n n 23)32)((3)]()32(33[11-+-=+--=∴--;解法二特征根法:数列{}n a :),0(025312N n n a a a n n n ∈≥=+-++, b a a a ==21,的特征方程是:02532=+-x x ;32,121==x x , ∴1211--+=n n n Bx Ax a 1)32(-⋅+=n B A ; 又由b a a a ==21,,于是:⎩⎨⎧-=-=⇒⎪⎩⎪⎨⎧+=+=)(32332b a B a b A B A b BA a 故1)32)((323--+-=n n b a a b a三、分式递推式定理3:如果数列}{n a 满足下列条件:已知1a 的值且对于N ∈n ,都有hra qpa a n n n ++=+1其中p 、q 、r 、h 均为常数,且r h a r qr ph -≠≠≠1,0,,那么,可作特征方程hrx q px x ++=. 1当特征方程有两个相同的根λ称作特征根时,若,1λ=a 则;N ,∈=n a n λ若λ≠1a ,则,N ,1∈+=n b a n n λ其中.N ,)1(11∈--+-=n r p rn a b n λλ特别地,当存在,N 0∈n 使00=n b 时,无穷数列}{n a 不存在;2当特征方程有两个相异的根1λ、2λ时,则112--=n n n c c a λλ,,N ∈n 其中).(,N ,)(211212111λλλλλ≠∈----=-a n rp r p a a c n n 其中例3、已知数列}{n a 满足性质:对于,324,N 1++=∈-n n n a a a n 且,31=a 求}{n a 的通项公式.解:依定理作特征方程,324++=x x x 变形得,04222=-+x x 其根为.2,121-==λλ故特征方程有两个相异的根,使用定理2的第2部分,则有:∴.N ,)51(521∈-=-n c n n ∴.N ,1)51(521)51(52211112∈----⋅-=--=--n c c a n n n nn λλ 即.N ,)5(24)5(∈-+--=n a nn n 例5.已知数列}{n a 满足:对于,N ∈n 都有.325131+-=+n n n a a a1若,51=a 求;n a 2若,31=a 求;n a 3若,61=a 求;n a 4当1a 取哪些值时,无穷数列}{n a 不存在解:作特征方程.32513+-=x x x 变形得,025102=+-x x特征方程有两个相同的特征根.5=λ依定理2的第1部分解答.1∵∴=∴=.,511λa a 对于,N ∈n 都有;5==λn a 2∵.,311λ≠∴=a a ∴λλr p rn a b n --+-=)1(11令0=n b ,得5=n .故数列}{n a 从第5项开始都不存在,当n ≤4,N ∈n 时,51751--=+=n n b a n n λ. 3∵,5,61==λa ∴.1λ≠a ∴.,811)1(11N n n r p r n a b n ∈-+=--+-=λλ令,0=n b 则.7n n ∉-=∴对于.0b N,n ≠∈n∴.N ,7435581111∈++=+-+=+=n n n n b a nn λ 4、显然当31-=a 时,数列从第2项开始便不存在.由本题的第1小题的解答过程知,51=a 时,数列}{n a 是存在的,当51=≠λa 时,则有.N ,8151)1(111∈-+-=--+-=n n a r p r n a b n λλ令,0=n b 则得N ,11351∈--=n n n a 且n ≥2.∴当11351--=n n a 其中N ∈n 且N ≥2时,数列}{n a 从第n 项开始便不存在. 于是知:当1a 在集合3{-或,:1135N n n n ∈--且n ≥2}上取值时,无穷数列}{n a 都不存在.定理3证明:分式递推问题:如果数列}{n a 满足下列条件:已知1a 的值且对于N ∈n ,都有hra q pa a n n n ++=+1其中p 、q 、r 、h 均为常数,且r ha r qr ph -≠≠≠1,0,,那么,可作特征方程hrx qpx x ++=.1当特征方程有两个相同的根λ称作特征根时,若,1λ=a 则;N ,∈=n a n λ若λ≠1a ,则,N ,1∈+=n b a n n λ其中.N ,)1(11∈--+-=n r p r n a b n λλ特别地,当存在,N 0∈n 使00=n b 时,无穷数列}{n a 不存在.2当特征方程有两个相异的根1λ、2λ称作特征根时,则112--=n n n c c a λλ,,N ∈n 其中).(,N ,)(211212111λλλλλ≠∈----=-a n rp r p a a c n n 其中证明:先证明定理的第1部分. 作交换N ,∈-=n a d n n λ, 则λλ-++=-=++h ra q pa a d n n n n 11hra hq r p a n n +-+-=λλ)( h d r h q r p d n n ++-+-+=)())((λλλλλλλλr h rd q p h r r p d n n -+--+--=])([)(2 ①∵λ是特征方程的根,∴λ.0)(2=--+⇒++=q p h r hr qp λλλλ将该式代入①式得.N ,)(1∈-+-=+n rh rd r p d d n n n λλ ②将rpx =代入特征方程可整理得,qr ph =这与已知条件qr ph ≠矛盾.故特征方程的根λ,rp≠于是.0≠-r p λ ③ 当01=d ,即λ+=11d a =λ时,由②式得,N ,0∈=n b n 故.N ,∈=+=n d a n n λλ 当01≠d 即λ≠1a 时,由②、③两式可得.N ,0∈≠n d n 此时可对②式作如下变化:.1)(11rp rd r p r h r p d r h rd d n n n n λλλλλ-+⋅-+=--+=+ ④由λ是方程h rx q px x ++=的两个相同的根可以求得.2r hp -=λ ∴,122=++=---+=-+h p p h rrh p p rr h p h r p r h λλ将此式代入④式得.N ,111∈-+=+n rp rd d n n λ 令.N ,1∈=n d b n n 则.N ,1∈-+=+n rp rb b n n λ故数列}{n b 是以r p r λ-为公差的等差数列.∴.N ,)1(1∈-⋅-+=n rp rn b b n λ其中.11111λ-==a db 当0,N ≠∈n b n 时,.N ,1∈+=+=n b d a nn n λλ当存在,N 0∈n 使00=n b 时,λλ+=+=0001n n n b d a 无意义.故此时,无穷数列}{n a 是不存在的. 再证明定理的第2部分如下:∵特征方程有两个相异的根1λ、2λ,∴其中必有一个特征根不等于1a ,不妨令.12a ≠λ于是可作变换.N ,21∈--=n a a c n n n λλ故21111λλ--=+++n n n a a c ,将hra qpa a n n n ++=+1代入再整理得N ,)()(22111∈-+--+-=+n hq r p a hq r p a c n n n λλλλ ⑤由第1部分的证明过程知r p x =不是特征方程的根,故.,21rp r p ≠≠λλ 故.0,021≠-≠-r p r p λλ所以由⑤式可得:N ,2211211∈--+--+⋅--=+n rp h q a r p hq a rp r p c n n n λλλλλλ ⑥∵特征方程hrx q px x ++=有两个相异根1λ、2λ⇒方程0)(2=--+q p h x rx 有两个相异根1λ、2λ,而方程xrp xh q x --=-与方程0)(2=---q p h x rx 又是同解方程.∴222111,λλλλλλ-=---=--rp hq r p h q将上两式代入⑥式得当,01=c 即11λ≠a 时,数列}{n c 是等比数列,公比为rp rp 21λλ--.此时对于N ∈n 都有当01=c 即11λ=a 时,上式也成立. 由21λλ--=n n n a a c 且21λλ≠可知.N ,1∈=n c n所以.N ,112∈--=n c c a n n n λλ证毕注:当qr ph =时,h ra q pa n n ++会退化为常数;当0=r 时,hra qpa a n n n ++=+1可化归为较易解的递推关系,在此不再赘述.求数列通项公式的方法很多,利用特征方程的特征根的方法是求一类数列通项公式的一种有效途径.1.已知数列{}n a 满足1n n n a a b a c a d+⋅+=⋅+......① 其中*0,,c ad bc n N ≠≠∈.定义1:方程ax bx cx d+=+为①的特征方程,该方程的根称为数列{}n a 的特征根,记为,αβ. 定理1:若1,a αβ≠且αβ≠,则11n n n n a a a c a a c a αααβββ++---=⋅---.定理2: 若1a αβ=≠且0a d +≠,则1121n n c a a d a αα+=+-+-.例109·江西·理·22各项均为正数的数列{}n a ,12,a a a b ==,且对满足m n p q +=+的正数,,,m n p q 都有(1)(1)(1)(1)p q m nm n p q a a a a a a a a ++=++++. 1当14,25a b ==时,求通项n a ;2略. 例2 已知数列{}n a 满足*1112,2,n n a a n N a -==-∈,求通项n a . 例 3 已知数列{}n a 满足11122,(2)21n n n a a a n a --+==≥+,求数列{}n a 的通项n a例4已知数列{}n a 满足*11212,()46n n n a a a n N a +-==∈+,求数列{}n a 的通项n a2.已知数列{}n a 满足2112n n n a c a c a ++=+② 其中12,c c 为常数,且*20,c n N ≠∈. 定义2:方程212x c x c =+为②的特征方程,该方程的根称为数列{}n a 的特征根,记为12,λλ.定理3:若12λλ≠,则1122n nn a b b λλ=+,其中12,b b 常数,且满足111222221122a b b a b b λλλλ=+⎧⎨=+⎩. 定理4: 若12λλλ==,则12()nn a b b n λ=+,其中12,b b 常数,且满足1122212()(2)a b b a b b λλ=+⎧⎨=+⎩. 例5已知数列{}n a 满足*12212,3,32()n n n a a a a a n N ++===-∈,求数列{}n a 的通项n a 例6已知数列{}n a 满足*12211,2,44()n n n a a a a a n N ++===-∈,求数列{}n a 的通项n a例7:已知数列{}n a 满足12212,8,44n n n a a a a a ++===-,求通项n a .。

特征方程法求数列地通项公式(1)

特征方程法求数列地通项公式(1)

特征方程法求数列的通项公式求数列通项公式的方法很多,利用特征方程的特征根的方法是求一类数列通项公式的一种有效途径.{}n a 满足1n n n a a ba c a d+⋅+=⋅+......① 其中*0,,c ad bc n N ≠≠∈.定义1:方程ax bx cx d+=+为①的特征方程,该方程的根称为数列{}n a 的特征根,记为,αβ. 定理1:假如1,a αβ≠且αβ≠,如此11n n n n a a a c a a c a αααβββ++---=⋅---. 证明:2()0,ax b a d bx cx d a x b cx d c cαβαβ+-=⇒+--=⇒+==-+(),d a c b c αβαβ∴=-+=-11()()()()()()()()n n n n n n n n n n n n aa ba ca d aab ca d ac a bd aa b a aa b ca d a c a b d ca dαααααβββββ+++--++-+-+-∴===+-+-+-+--+ ()[()]()()()[()]()()n n n n a c a c a c c a c a a c a c a c a c c a c a a c ααβαβααααβαβαβββββ-+-------==-+------- n n a a c a c a ααββ--=⋅-- 证毕 定理2: 假如1a αβ=≠且0a d +≠,如此1121n n c a a d a αα+=+-+-. 证明:22,d a c b c αα=-=-111()()()n n n n n n n n ca d ca d aa b a aa b ca d a c a b dca dααααα+++∴===+-+-+-+--+22222()(2)()()()2n n n n n nca a c ca a c ca a ca d a c a c a c a c a a αααααααααα+-+-+-===+--+---- 2242(2)2()()()()()()()()n n n n n n ca a c ca a c d c a a d a d a a d a a d a αααααα+-+-+-++===+-+-+- 21n c a d a α=++- 证毕 例〔09·某某·理·22〕各项均为正数的数列{}n a ,12,a a a b ==,且对满足m n p q +=+的正数,,,m n p q 都有(1)(1)(1)(1)p q m nm n p q a a a a a a a a ++=++++. (1)当14,25a b ==时,求通项n a ;(2)略. 解:由(1)(1)(1)(1)p q m nm n p q a a a a a a a a ++=++++得121121(1)(1)(1)(1)n n n n a a a a a a a a --++=++++ 将14,25a b ==代入上式化简得11212n n n a a a --+=+ 考虑特征方程212x x x +=+得特征根1x =± 所以11111121112112113112n n n n n n n n a a a a a a a a ------+--+-==⋅+++++ 所以数列11n n a a ⎧⎫-⎨⎬+⎩⎭是以111113a a -=-+为首项,公比为13的等比数列 故11111()()1333n nn n a a --=-⋅=-+ 即3131n n na -=+ 例 数列{}n a 满足*1112,2,n n a a n N a -==-∈,求通项n a . 解: 考虑特征方程12x x=-得特征根1x = 111111111111111(2)11n n n n n n a a a a a a -----====+------ 所以数列11n a ⎧⎫⎨⎬-⎩⎭是以1111a =-为首项,公差为1的等差数列 故11n n a =- 即1n n a n += 例 数列{}n a 满足11122,(2)21n n n a a a n a --+==≥+,求数列{}n a 的通项n a解:其特征方程为221x x x +=+,化简得2220x -=,解得121,1x x ==-,令111111n nn n a a c a a ++--=⋅++ 由12,a =得245a =,可得13c =-,∴数列11n n a a ⎧⎫-⎨⎬+⎩⎭是以111113a a -=+为首项,以13-为公比的等比数列,1111133n n n a a --⎛⎫∴=⋅- ⎪+⎝⎭,3(1)3(1)n n n n na --∴=+- 例数列{}n a 满足*11212,()46n n n a a a n N a +-==∈+,求数列{}n a 的通项n a解:其特征方程为2146x x x -=+,即24410x x ++=,解得1212x x ==-,令1111122n n c a a +=+++由12,a =得2314a =,求得1c =, ∴数列112n a ⎧⎫⎪⎪⎨⎬⎪⎪+⎩⎭是以112152a =+为首项,以1为公差的等差数列,123(1)11552n n n a ∴=+-⋅=-+, 135106n n a n -∴=-2.数列{}n a 满足2112n n n a c a c a ++=+② 其中12,c c 为常数,且*20,c n N ≠∈.定义2:方程212x c x c =+为②的特征方程,该方程的根称为数列{}n a 的特征根,记为12,λλ. 定理3:假如12λλ≠,如此1122n n n a b b λλ=+,其中12,b b 常数,且满足111222221122a b b a b b λλλλ=+⎧⎨=+⎩. 定理4: 假如12λλλ==,如此12()nn a b b n λ=+,其中12,b b 常数,且满足1122212()(2)a b b a b b λλ=+⎧⎨=+⎩. 设)(11-+-=-n n n n ta a s ta a ,如此11)(-+-+=n n n sta a t s a , 令⎩⎨⎧-==+q st pt s 〔*〕(1) 假如方程组〔*〕有两组不同的解),(),,(2211t s t s , 如此)(11111-+-=-n n n n a t a s a t a ,)(12221-+-=-n n n n a t a s a t a ,由等比数列性质可得1111211)(-+-=-n n n s a t a a t a ,1212221)(1-+-=-n n n s a t a a t a ,,21t t ≠ 由上两式消去1+n a 可得()()()n n n s t t s a t a s t t s a t a a 21221221121112..-----=.(2) 假如方程组〔*〕有两组相等的解⎩⎨⎧==2121t t s s ,易证此时11t s =,如此()()112112112111111)(a t a s a t a s a t a s a t a n n n n n n n -==-=-=-----+ ,211121111s a t a s a s a nn n n -=-∴++,即⎭⎬⎫⎩⎨⎧n n s a 1是等差数列, 由等差数列性质可知()21112111.1s a t a n s a s a nn --+=, 所以n n s n s a t a s a t a s a a 1211122111211.⎥⎥⎦⎤⎢⎢⎣⎡-+⎪⎪⎭⎫⎝⎛--=.例数列{}n a 满足*12212,3,32()n n n a a a a a n N ++===-∈,求数列{}n a 的通项n a 解:其特征方程为232x x =-,解得121,2x x ==,令1212n n n a c c =⋅+⋅,由1122122243a c c a c c =+=⎧⎨=+=⎩,得12112c c =⎧⎪⎨=⎪⎩, 112n n a -∴=+例数列{}n a 满足*12211,2,44()n n n a a a a a n N ++===-∈,求数列{}n a 的通项n a解:其特征方程为2441x x =-,解得1212x x ==,令()1212nn a c nc ⎛⎫=+ ⎪⎝⎭,由1122121()121(2)24a c c a c c ⎧=+⨯=⎪⎪⎨⎪=+⨯=⎪⎩,得1246c c =-⎧⎨=⎩, 1322n n n a --∴=例:数列{}n a 满足12212,8,44n n n a a a a a ++===-,求通项n a .解: 考虑特征方程244x x =-得特征根2λ=如此12()2n n a b b n =+ 其中1211222()2024(2)81n n b b b a n b b b +==⎧⎧⇒⇒=⎨⎨+==⎩⎩常见递推数列通项的求解方法高考中的递推数列求通项问题,情境新颖别致,有广度,创新度和深度,是高考的热点之一。

求矩阵特征值的工具

求矩阵特征值的工具

求矩阵特征值的工具求矩阵特征值是线性代数中一个重要的问题,用于解决许多实际问题,如物理系统模型的求解、网络分析、图像处理、数据压缩等。

在过去的几十年里,研究人员开发了许多用于计算矩阵特征值的工具和算法。

本文将介绍一些常用的求矩阵特征值的工具,并对其原理和特点进行详细的阐述。

1.解特征方程法:这是最直接的一种方法,它是通过将矩阵特征值问题转换为求解特征方程根的问题。

对于一个n阶矩阵A,它的特征方程为,A-λI,=0,其中λ是一个标量,I是单位矩阵。

解特征方程可以得到所有的特征值。

这种方法适用于已知特征方程的情况,但对于高阶矩阵来说,求解特征方程是一个复杂和耗时的过程。

2.特征向量法:特征向量法是一种直接计算特征值和特征向量的方法。

它利用特征向量的线性相关性,将特征向量排列成一个阵,通过求解线性方程组来计算特征值。

特征向量法的优点是可以同时求解所有的特征值和特征向量,但需要对矩阵求逆,当矩阵非常大时,计算复杂度较高。

3.幂法:幂法是一种迭代方法,用于求解矩阵的最大特征值和对应的特征向量。

它利用矩阵的最大特征值对特征向量进行放大,从而逼近真正的特征向量。

幂法收敛于最大特征值对应的特征向量,但只能求解最大特征值。

为了提高收敛速度和准确性,可以使用位移策略和正交化技术。

4.反幂法:反幂法是幂法的变形,用于求解矩阵的最小特征值和对应的特征向量。

它利用矩阵的最小特征值对特征向量进行放大,从而逼近真正的特征向量。

反幂法收敛于最小特征值对应的特征向量,但只能求解最小特征值。

同样,为了提高收敛速度和准确性,可以使用位移策略和正交化技术。

5.QR方法:QR方法是一种迭代方法,可以同时求解所有的特征值和特征向量。

它将矩阵分解为正交矩阵和上三角矩阵的乘积,通过反复迭代使得上三角矩阵收敛为对角矩阵,对角元素即为矩阵的特征值。

QR方法的收敛速度较快,但在一些情况下可能会出现不收敛或者收敛到错误的结果的问题。

以上是常用的求矩阵特征值的工具和算法,每种方法都有其适用范围和特点。

特征方程法求行列式

特征方程法求行列式

特征方程法求行列式
特征方程法是一种求解行列式的技巧,对于学生来说,是一项十分重要的高等
数学学习内容,需要认真掌握。

特征方程法指的是将行列式拆分成n项乘积,其中一项是相应阶数的特征方程,即行列式的指数是n的方程,其他项是与特征方程次数相乘的行列式的元素。

换言之,特征方程法就是将复杂的求解过程简化,不仅可以快速准确地求出行列式,而且它还有助于深入了解该行列式的结构特征。

特征方程法的基本思想是“分而治之”,即将难以解决的问题分解为若干个容
易解决的问题,然后逐步解决各个问题,最终达到解题目标。

这个过程类似于分治策略。

在这种分解过程中,特征方程是最重要的。

它为求解行列式提供了便利,使得整个求解过程得以高效进行。

在实际应用中,特征方程法可以解决不同阶数的行列式求解问题,无论是一阶
行列式还是n阶行列式,特征方程法都可以快速解决。

它的优势在于对计算量的节省,可以有效地减小计算工作量,使用特征方程法可以节省很长的时间。

通过掌握特征方程法,学生可以准确快速地求解行列式,为学习和工作中遇到
的类似问题提供解决方案,也为深入理解数学领域的学习提供技能支撑。

掌握此方法是学生做好高等数学学习必备的知识和技能,在今后的学习中也可以灵活运用,帮助提升学习效率和实践能力。

求特征方程

求特征方程

求特征方程
要求特征方程,首先需要明确特征方程的定义和求法。

特征方程是根据线性常微分方程组来定义的。

给定一个n阶线性常微分方程组
dydx=A(x)y(x)\frac{dy}{dx} = A(x)y(x)dxdy=A(x)y(x)其中A(x)是一个n ×n矩阵,我们定义特征方程为
|λE−A(x)=0\lambda E - A(x) = 0λE−A(x)=0|其中λ\lambdaλ是我们要找的数,E是单位矩阵。

具体求解特征方程的过程涉及到对矩阵和常微分方程组的深入理解,包括如何求解矩阵的逆,如何求解矩阵的行列式等。

如果你可以提供更多关于你想求解特征方程的线性常微分方程组的信息,我会很乐意帮助你进一步求解。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

特征方程法求解递推关系中的数列通项考虑一个简单的线性递推问题.设已知数列}{n a 的项满足其中,1,0≠≠c c 求这个数列的通项公式.采用数学归纳法可以求解这一问题,然而这样做太过繁琐,而且在猜想通项公式中容易出错,本文提出一种易于被学生掌握的解法——特征方程法:针对问题中的递推关系式作出一个方程,d cx x +=称之为特征方程;借助这个特征方程的根快速求解通项公式.下面以定理形式进行阐述.定理 1.设上述递推关系式的特征方程的根为0x ,则当10a x =时,n a 为常数列,即0101,;x b a a x a a n n n +===时当,其中}{n b 是以c 为公比的等比数列,即01111,x a b c b b n n -==-.证明:因为,1,0≠c 由特征方程得.10c d x -=作换元,0x a b n n -= 则.)(110011n n n n n n cb x a c ccd ca c d d ca x a b =-=--=--+=-=-- 当10a x ≠时,01≠b ,数列}{n b 是以c 为公比的等比数列,故;11-=n n c b b当10a x =时,01=b ,}{n b 为0数列,故.N ,1∈=n a a n (证毕)下面列举两例,说明定理1的应用.例1.已知数列}{n a 满足:,4,N ,23111=∈--=+a n a a n n 求.n a解:作方程.23,2310-=--=x x x 则 当41=a 时,.21123,1101=+=≠a b x a 数列}{n b 是以31-为公比的等比数列.于是.N ,)31(2112323,)31(211)31(1111∈-+-=+-=-=-=---n b a b b n n n n n n 例2.已知数列}{n a 满足递推关系:,N ,)32(1∈+=+n i a a n n 其中i 为虚数单位.当1a 取何值时,数列}{n a 是常数数列?解:作方程,)32(i x x +=则.5360i x +-=要使n a 为常数,即则必须.53601i x a +-== 现在考虑一个分式递推问题(*).例3.已知数列}{n a 满足性质:对于,324,N 1++=∈-n n n a a a n 且,31=a 求}{n a 的通项公式. 将这问题一般化,应用特征方程法求解,有下述结果.定理2.如果数列}{n a 满足下列条件:已知1a 的值且对于N ∈n ,都有hra q pa a n n n ++=+1(其中p 、q 、r 、h 均a 1=ba n+1=ca n +d为常数,且r h a r qr ph -≠≠≠1,0,),那么,可作特征方程hrx q px x ++=. (1)当特征方程有两个相同的根λ(称作特征根)时,若,1λ=a 则;N ,∈=n a n λ若λ≠1a ,则,N ,1∈+=n b a n n λ其中.N ,)1(11∈--+-=n r p r n a b n λλ特别地,当存在,N 0∈n 使00=n b 时,无穷数列}{n a 不存在.(2)当特征方程有两个相异的根1λ、2λ(称作特征根)时,则112--=n n n c c a λλ,,N ∈n 其中).(,N ,)(211212111λλλλλ≠∈----=-a n rp r p a a c n n 其中 证明:先证明定理的第(1)部分.作交换N ,∈-=n a d n n λ 则λλ-++=-=++hra q pa a d n n n n 11 hra h q r p a n n +-+-=λλ)( h d r h q r p d n n ++-+-+=)())((λλλλ λλλλr h rd q p h r r p d n n -+--+--=])([)(2 ① ∵λ是特征方程的根,∴λ.0)(2=--+⇒++=q p h r h r q p λλλλ 将该式代入①式得.N ,)(1∈-+-=+n r h rd r p d d n n n λλ ② 将r p x =代入特征方程可整理得,qr ph =这与已知条件qr ph ≠矛盾.故特征方程的根λ,r p ≠于是.0≠-r p λ ③当01=d ,即λ+=11d a =λ时,由②式得,N ,0∈=n b n 故.N ,∈=+=n d a n n λλ当01≠d 即λ≠1a 时,由②、③两式可得.N ,0∈≠n d n 此时可对②式作如下变化:.1)(11rp r d r p r h r p d r h rd d n n n n λλλλλ-+⋅-+=--+=+ ④由λ是方程h rx q px x ++=的两个相同的根可以求得.2rh p -=λ ∴,122=++=---+=-+h p p h r r h p p r r h p h r p r h λλ 将此式代入④式得.N ,111∈-+=+n rp r d d n n λ 令.N ,1∈=n d b n n 则.N ,1∈-+=+n rp r b b n n λ故数列}{n b 是以r p r λ-为公差的等差数列. ∴.N ,)1(1∈-⋅-+=n rp r n b b n λ 其中.11111λ-==a d b 当0,N ≠∈n b n 时,.N ,1∈+=+=n b d a n n n λλ 当存在,N 0∈n 使00=n b 时,λλ+=+=0001n n n b d a 无意义.故此时,无穷数列}{n a 是不存在的. 再证明定理的第(2)部分如下: ∵特征方程有两个相异的根1λ、2λ,∴其中必有一个特征根不等于1a ,不妨令.12a ≠λ于是可作变换.N ,21∈--=n a a c n n n λλ 故21111λλ--=+++n n n a a c ,将h ra q pa a n n n ++=+1代入再整理得 N ,)()(22111∈-+--+-=+n hq r p a h q r p a c n n n λλλλ ⑤ 由第(1)部分的证明过程知r p x =不是特征方程的根,故.,21r p r p ≠≠λλ 故.0,021≠-≠-r p r p λλ所以由⑤式可得:N ,2211211∈--+--+⋅--=+n r p h q a r p h q a r p r p c n n n λλλλλλ ⑥ ∵特征方程h rx q px x ++=有两个相异根1λ、2λ⇒方程0)(2=--+q p h x rx 有两个相异根1λ、2λ,而方程xrp xh q x --=-与方程0)(2=---q p h x rx 又是同解方程. ∴222111,λλλλλλ-=---=--r p h q r p h q 将上两式代入⑥式得N ,2121211∈--=--⋅--=-n c rp r p a a r p r p c n n n n λλλλλλ 当,01=c 即11λ≠a 时,数列}{n c 是等比数列,公比为r p r p 21λλ--.此时对于N ∈n 都有 .))(()(12121111211------=--=n n n rp r p a a r p r p c c λλλλλλ 当01=c 即11λ=a 时,上式也成立. 由21λλ--=n n n a a c 且21λλ≠可知.N ,1∈=n c n 所以.N ,112∈--=n c c a n n n λλ(证毕)注:当qr ph =时,h ra q pa n n ++会退化为常数;当0=r 时,h ra q pa a n n n ++=+1可化归为较易解的递推关系,在此不再赘述.现在求解前述例3的分类递推问题)(*. 解:依定理作特征方程,324++=x x x 变形得,04222=-+x x 其根为.2,121-==λλ故特征方程有两个相异的根,使用定理2的第(2)部分,则有.N ,)221211(2313)(11212111∈⋅-⋅-⋅+-⋅--⋅--=--n r p r p a a c n n n λλλλ ∴.N ,)51(521∈-=-n c n n ∴.N ,1)51(521)51(52211112∈----⋅-=--=--n c c a n n n nn λλ 即.N ,)5(24)5(∈-+--=n a nn n 例4.已知数列}{n a 满足:对于,N ∈n 都有.325131+-=+n n n a a a(1)若,51=a 求;n a (2)若,31=a 求;n a (3)若,61=a 求;n a (4)当1a 取哪些值时,无穷数列}{n a 不存在? 解:作特征方程.32513+-=x x x 变形得,025102=+-x x特征方程有两个相同的特征根.5=λ依定理2的第(1)部分解答.(1)∵∴=∴=.,511λa a 对于,N ∈n 都有;5==λn a(2)∵.,311λ≠∴=a a ∴λλr p r n a b n --+-=)1(11 51131)1(531⋅-⋅-+-=n ,8121-+-=n 令0=n b ,得5=n .故数列}{n a 从第5项开始都不存在,当n ≤4,N ∈n 时,51751--=+=n n b a n n λ. (3)∵,5,61==λa ∴.1λ≠a ∴.,811)1(11N n n r p r n a b n ∈-+=--+-=λλ 令,0=n b 则.7n n ∉-=∴对于.0b N,n ≠∈n ∴.N ,7435581111∈++=+-+=+=n n n n b a n n λ (4)显然当31-=a 时,数列从第2项开始便不存在.由本题的第(1)小题的解答过程知,51=a 时,数列}{n a 是存在的,当51=≠λa 时,则有.N ,8151)1(111∈-+-=--+-=n n a r p r n a b n λλ令,0=n b 则得N ,11351∈--=n n n a 且n ≥2. ∴当11351--=n n a (其中N ∈n 且N ≥2)时,数列}{n a 从第n 项开始便不存在. 于是知:当1a 在集合3{-或,:1135N n n n ∈--且n ≥2}上取值时,无穷数列}{n a 都不存在.。

相关文档
最新文档