线性代数及其应用 第四章

第 4 章 向量空间

4.1 向量空间和子空间

定义 一个向量空间是由一些被称为向量的对象构成的非空集合 V,在这个集合上定义了两种运算,称为加法和标量乘法(标量取实数),服从以下公理(或法则),这些公理必须对 V 中所有向量 u,v,w 及所有标量(或数)cd 均成立。

1. u,v 之和(表示为 u+v)属于 V

2. u+v=v+u

3. (u+v)+w=u+(v+w)

4. V 中存在一个向量 0,使得 u+0=u

5.V 中每个向量 u,存在 V 中一个向量 u,使得 u+(u)=0

6. u 与标量 c 的标量乘法(记为 cu)属于 V

7. c(u+v)=cu+cv

8. (c+d)u=cu+du

9. c(du)=(cd)u

10. 1u=u

利用这些公理可以证明零向量唯一,且 u负向量 u 也是唯一的。有如下简单事实:

V 中每个向量 u 和任意标量 c,有

0u=0

c0=0

u=(1)u

注:更专业地,V 是一个实数向量空间,本章所有理论对负向量空间(标量取复数)也成立。

关于向量空间可以给出许多例子。例如对 n0,次数最高为 n 的多项式集合 Pn 也是向量空间。

子空间

定义 向量空间 V 的一个子空间V 的一个满足以下三个性质的子集 H

a. V 中的零向量在 H 中。

b. H 对向量加法封闭,即对 H 中任意向量 u,v,和 u+v 仍在 H 中。

c. H 对标量乘法封闭,即对 H 中任意向量 u 和任意标量 c,向量 cu 仍在 H 中。

这些性质(公理 1,4,6)保证 H 本身对 V 中定义的向量空间运算而言是一个向量空间。这样,每个子空间都是一个向量空间。反之,每个向量空间都是一个子空间(针对本身或其他更大的空间而言)。

有一些简单的例子:

  • 向量空间 V 中仅由零向量组成的集合是 V 的一个子空间,称为零子空间,写成 {0}

  • P 为全体实系数多项式的集合,则 P 是定义在 R 上的全体实值函数的空间的一个子空间。再者,对每个 n0PnP 的子空间。

  • 注意向量空间 R2 不是 R3 的子空间,因为 R2 不是 R3 的子集,但集合 H={[st0]:st}R3 的子空间。

由一个集合生成的子空间

定理 1

v1,v2,,vp 在向量空间 V 中,则 Span{v1,v2,,vp}V 的一个子空间。

Span{v1,v2,,vp} 是由 {v1,v2,,vp} 生成(或张成的子空间。给定 V 的任一子空间 HH生成(或张成是集合 {v1,v2,,vp}H,满足 H=Span{v1,v2,,vp}

现有

H={(a3b,ba,a,b):a,bR}

证明 HR4 的子空间。


H 中向量写完列向量,则

[a3bbaab]=a[1110]+b[3101]

令右边两个向量为 v1,v2,这说明 H=Span{v1,v2}。由定理 1 知 HR4 的子空间。

4.2 零空间、列空间、行空间和线性变换

矩阵的零空间

我们在第二章中已接触过。

定义 矩阵 A 的列空间写成 NulA,是齐次方程 Ax=0 的全体解的集合,用集合符号表示,即

NulA={x:xRn,Ax=0}

定理 2

m×n 矩阵 A 的零空间是 Rn 的一个子空间。等价地,m 个方程、n 个未知数的齐次线性方程组 Ax=0 的全体解的集合是 Rn 的一个子空间。

首先 NulARn 的子集,且 0NulA。然后取 u,vNulA,验证 NulA 对向量加法和标量乘法封闭即可。

NulA 的一个显式刻画

NulA 中向量与 A 的元素之间无明显关系,称 NulA隐式地定义。解出 Ax=0 可以得到 NulA显式刻画。用自由变量作权,假设已解得 x=xiu+xjv+xkw,则 u,v,w 的每一个线性组合都是 NulA 中的一个元素,反之亦然。从而 {u,v,w}NulA 的一个生成集。

有以下简单事实:

  • 上述的 {u,v,w} 必然线性无关。仅当 xi=xj=xk=0xiu+xjv+xkw 为零。

  • NulA 包含非零向量时,它的生成集中向量的个数等于方程 Ax=0 中自由变量的个数。

矩阵的列空间

定理 3

m×n 矩阵 A 的列空间是 Rm 的一个子空间。

A 的列生成 Rm 当且仅当方程 Ax=b 对任一个 b 有解,重述可得:

m×n 矩阵 A 的列空间等于 Rm 当且仅当方程 Ax=bRm 中每个 b 有一个解。

m×n 矩阵 ANulAColA 有许多显然的区别,需要注意的是 NulARn 的一个子空间,而 ColARm 的一个子空间。

An×n 矩阵,且 ColA=NulA,证明 NulA2=Rn


设任意 xRn,显然 AxColA。由题知 A(Ax)=A2x=0,即 xRnxNulA2。得证。

行空间

m×n 矩阵 A 的行向量的所有线性组合的集合称为 A 的行空间,记为 RowA,它是 Rn 的一个子空间。也可以用 ColAT 代替 RowA

线性变换的核与值域

我们推广 1.8 节中线性变换的定义。

定义 由向量空间 V 映射到向量空间 W线性变换 T 是一个规则,它将 V 中每个向量 x 映射成 W 中唯一向量 T(x),且满足:

(i) T(u+v)=T(u+v),对 V 中所有 u,v 均成立。

(ii) T(cu)=cT(u),对 V 中所有 u 和所有数 c 均成立。

线性变换 T(或零空间)是 V 中所有满足 T(u)=0 的向量 u 的集合(0W 中的零向量)。T值域W 中所有具有形式 T(x)(任意 xV)的向量的集合。若 T(x)=AxT 的核与值域恰好是前面定义的 A 的零空间和列空间。

T 的核是 V 的子空间,值域是 W 的子空间。

应用中,一个子空间往往由一个适当的线性变换的核或值域产生。例如一个齐次线性微分方程的全部解的集合最终是一个线性变换的核。典型地,这样一个线性变换用关于一次函数的一阶或高阶导数描述。以下有两个例子。

V 为定义在区间 [a,b] 上的所有连续可导的实函数 f 构成的向量空间,令 W[a,b] 上所有连续函数构成的向量空间 C[a,b],且令 D:VW 是将 Vf 变为其导数 f 的变换。由微积分中两个简单的微分法则,

D(f+g)=D(f)+D(g),D(cf)=cD(f)

于是 D 是一个线性变换,可以证明 D 的核是 [a,b] 上的常函数的集合,D 的值域是 [a,b] 上所有连续函数的集合 W

有微分方程

y+ω2y=0

ω 是常数,常常用来描述物理系统中的一个变化过程。该式的解集是将函数 y=f(t) 映成函数 f(t)+ω2f(t) 的线性变换的核。寻找这个向量空间的一个显式刻画是微分方程中的一个问题。

4.3 线性无关集和基

定义在向量空间 V 中的线性无关、线性相关等定义与 Rn 中一致。下列定理与 2.8 中已经讨论过的大致相同。

定理 4

两个或多个向量组成的向量集合 {v1,v2,,vp}(如果 v10)是线性相关的,当且仅当某 vj(j>1) 是其前面向量 v1,v2,,vj1 的线性组合。

一般向量空间中的线性相关与 Rn 中的线性相关与 Rn 的不同之处在于,当向量不是 n 元组时,齐次方程

c1v1+c2v2++cpvp=0

通常不能写为一个 n 个方程的线性方程组。即为研究方程 Ax=0,向量不能从一个矩阵 A 中得到,反而要依靠线性相关的定义和定理 4。

例如 p1(t)=1,p2(t)=t,p3(t)=4tp 是函数),有 p3=4p1p2,从而 {p1,p2,p3} 是线性相关的;集合 {sint,cost}C[0,1] 中是线性无关的,因为不存在数 c 使得 cost=csint 对任意 t[0,1] 成立,而 {sintcost,sin2t} 是线性相关的。

定义 H 是向量空间 V 的一个子空间。V 中向量的指标集 B 称为 H 的一个,如果

(i) B 是一线性无关集。

(ii)B 生成的子空间与 H 相同,即 H=SpanB

这显然允许 H=VV 的一个基是生成 V 的一个线性无关集。当 HV 时,条件 (ii) 蕴涵 B 中每个向量都属于 H

S={1,t,t2,,tn},证明 SPn 的一个基,称为 Pn 的标准基。


显然 S 生成 Pn,为证 S 线性无关,假设 c0,c1,,cn 满足

c01+c1t++cntn=0(t)

该式对任意 t 成立,于是左式为零多项式,c0=c1==cn=0,故 S 线性无关且是 Pn 的一个基。

定理 5(生成集定理)

S={v1,v2,,vp} 是向量空间 V 中的向量集,H=Span{v1,v2,,vp}

a.S 中某一个向量 vkS 中其余向量的线性组合,则 S 中去掉 vk 后形成的集合仍然可以生成 H

b.H{0},则 S 的某一子集是 H 的一个基。

定理 6

矩阵 A 的主元列构成 ColA 的一个基。

行变换会改变矩阵的列空间,而行空间不受影响。

定理 7

如果矩阵 A 行等价于矩阵 B,那么它们的行空间是相同的。如果 B 是阶梯形,那么 B 的非零行构成了 A 的行空间的基,也是 B 的行空间的基。

4.4 坐标系

在 2.9 节中已接触过。

定理 8(唯一表示定理)

B={b1,b2,,bn} 是向量空间 V 的一个基,则对 V 中每个向量 x,存在唯一的一组数 c1,c2,,cn 使得

x=c1b1+c2b2++cnbn

定义 假设 B={b1,b2,,bn} 是向量空间 V 的一个基,xV 中,x 相对于基 B 的坐标(或 xB坐标)是使得 x=c1b1+c2b2++cnbn 的权 c1,c2,,cn

c1,c2,,cnxB坐标,则 Rn 中的向量

[x]B=[c1c2cn]

x(相对于 B)的坐标向量xB 坐标向量,映射 x[x]B 称为(由 B 确定的)坐标映射

{e1,e2}{e1,e2,e3} 分别作为 R2R3 的基是一个十分简单的具有几何意义的例子。

对于 Rn,设其一基为 B,则求向量 xB 坐标向量只需解方程 x=PB[x]B,其中

PB=[b1b2bn]

PB 为从 BRn 中标准基的坐标变换矩阵。它是可逆的,于是

PB1x=[x]B

PB1 产生的映射 x[x]B 即前面的坐标映射,此坐标映射是一个由 RnRn 的一对一线性变换。

坐标映射

定理 9

B={b1,b2,,bn} 是向量空间 V 的一个基,则坐标映射 x[x]B 是一个由 V 映射到 Rn 上的一对一的线性变换。

如 1.8 节,将坐标映射的线性性质推广到线性组合:即若 u1,u2,,up 属于 Vc1,c2,,cp 是数,则

[c1u1+c2u2++cpup]B=c1[u1]B+c2[u2]B++cp[up]B

定理 9 的坐标映射是一个由 VRn同构的重要例子。一般而言,从一个向量空间 V 映射到另一个向量空间 W 的一对一线性变换称为从 VW 上的同构(isomorphism(希腊语中 iso 表相同,morph 表形状或结构)。每一个在 V 中的向量空间的计算可以完全相同地出现在 W 中,反之亦然

B={1,t,t2,t3}P3 的标准基。P3 中的一个典型元素 p 具有形式

p(t)=a0+a1t+a2t2+a3t3

断定

[p]B=[a0a1a2a3]

于是坐标映射 p[p]B 是一个 P3R4 的同构。如果考虑将 P3R4 分别在两个计算机展现出来,两个显示屏上的向量看起来不同,但作为向量的“作用”是完全相同的,因为若两个显示屏由坐标变换相联系,P3 的每一个向量空间的运算被正确地复制到 R4 的一个对应的向量运算。

4.5 向量空间的维数

定理 9 表明,向量空间 V 的基 B 若含有 n 个向量,则 VRn 同构。数 nV 的一个内在性质(称为维数),它不依赖基的选择。

定理 10

若向量空间 V 具有一组基 B={b1,b2,,bn},则 V 中任意包含多于 n 个向量的集合一定线性相关。

V 中有向量集 {u1,u2,,up}p>n,则坐标向量 [u1]B,[u2]B,,[up]B 线性相关,存在不全为 0 的数 c1,c2,,cp 使得

c1[u1]B+c2[u2]B++cp[up]B=0

[c1u1+c2u2++cpup]B=0c1u1+c2u2++cpup=PB0=0

因为 ci 不全为零,{u1,u2,,up} 线性相关。得证。

定理 10 可应用到 V 的无限集中。一个无线集称为线性相关的,如果其中某一有限集是线性相关的,否则这个无限集是线性无关的。

定理 11

若向量空间 V 有一组基含有 n 个向量,则 V 的每一组基一定恰好含有 n 个向量。

设有两个不同的基 B1,B2B1n 个向量。因为 B1 是一个基,B2 中的向量个数不超过 n;同理因 B2 是一个基,B2 中的向量个数不少于 n。得证。

若非零向量空间 V 由有限集 S 生成,由生成集定理,S 的一个子集是 V 的一个基。由此定理 11 保证下列定义有意义。

定义 若向量空间 V 由一个有限集生成,则 V 称为有限维的V 的维数写成 dimV,是 V 的基中向量的个数。零向量空间 {0} 的维数定义为零。如果 V 不能由一有限集生成,则 V 称为无穷维的

显然地,dimRn=ndimPn=n+1

有限维空间的子空间

下一定理是生成集定理的一个自然配对。

定理 12

H 是有限维向量空间 V 的子空间。若有必要的话,H 中任一个线性无关集均可以扩充成为 H 的一个基。H 也是有限维的并且

dimHdimV

定理 13(基定理)

V 是一个 p 维向量空间,p1V 中任意含有 p 个元素的线性无关集必然是 V 的一个基。任意含有 p 个元素且生成 V 的集合自然是 V 的一个基。

NulA,ColA,RowA 的维数

定义 m×n 矩阵 A 矩阵 A 的秩是列空间的维数,A 的零维是零空间的维数。

A 的秩就是主元列个数。RowA 的一个基可从 A 的行最简阶梯形的主元行中找到,所以 RowA 的维数也等于 A 的秩。A 的零维是自由变量的数量。

该定理与 2.9 节定理 14 相同。

定理 14(秩定理)

一个 m×n 矩阵 A 的列空间的维数和 A 的零空间的维数满足下面的方程:

rankA+dimNulA=A

应用到方程组

秩定理是处理线性方程组信息的一个有力工具。

已对一个 40 个方程 42 个变量的齐次方程组求出了两个解。它们不是倍数关系,且其他所有解均能表示为这两个解的适当倍数之和。能否确定一个相应的非齐次方程组有解?


A 为方程组的 40×42 系数矩阵。条件蕴涵两个解线性无关且能生成 NulA,故 dimNulA=2。由秩定理,rankA=422=40。由于 R40R40 唯一的维数是 40 的子空间,一定有 ColA=R40,这表明每个非齐次方程组 Ax=b 有一个解。

4.6 基的变换

对一个向量空间 V,考虑两个基 B={b1,b2}C={c1,c2},满足 b1=4c1+c2b2=6c1+c2。假设 x=3b1+b2,即 [x]B=[31],求 [x]C

x=3b1+b2[x]C=3[b1]C+[b2]C,将它写为矩阵方程:

[x]C=[[b1]C[b2]C][31]

确定矩阵的列即可得到 [x]C,由题有 [b1]C=[41][b2]C=[61],于是

[x]C=[4611][31]=[64]


将前述矩阵方程推广可产生如下结果。

定理 15

B={b1,b2,,bn}C={c1,c2,,cn} 是向量空间 V 的基,则存在一个 n×n 矩阵 PCB 使得

[x]C=PCB[x]B

PCB 的列是 B 中向量的 C 坐标向量,即

PCB=[[b1]C[b2]C[bn]C]

矩阵 PCB 称为BC 的坐标变换矩阵,乘以 PCB 的运算将 B 坐标变为 C 坐标。

PCB 的列是线性无关的,因为它们是线性无关集 B 的坐标向量。于是 PCB 可逆,有

(PCB)1[x]C=[x]B

(PCB)1 是将 C 坐标变为 B 的矩阵,即 (PCB)1=PBC

Rn 中基的变换

BRn 的一个基,E 为标准基 {e1,e2,,en}(这个字母是花体的 E,不是 ε),则 [bi]E=bi,故 PEB 与坐标变换矩阵 PB 相同,即

PEB=PB=[b1b2bn]

定理 15 为在 Rn 中两个非标准基间变换坐标提供帮助。

B,CR2 中的基,b1=[91],b2=[51],c1=[14],c2=[35],求 BC 的坐标变换矩阵。


求矩阵 PCB 涉及 b1b2C 坐标向量。设 [b1]C=[x1x2],[b2]C=[y1y2],即

[c1c2][x1x2]=b1,[c1c2][y1y2]=b2

同步解出这两个方程组,将 b1,b2 增加到系数矩阵并行化简:

[c1c2b1b2]=[13954511][10640153]

于是 [b1]C=[65],[b2]C=[43],所求坐标变换矩阵为

PCB=[[b1]C[b2]C]=[6453]

PCB 已经出现在行化简后的矩阵内,因为 PCB 的第 1 列是将 [c1c2b1] 行化简为 [I[b1]C],第 2 列同理,于是

[c1c2b1b2][IPCB]

显然地能够通过 PBC=(PCB)1 求出由 CB 的坐标变换矩阵。

另外地,对 Rn 中的 x,有

PB[x]B=x

PC[x]C=x[x]C=PC1x

于是

[x]C=PC1x=PC1PB[x]B

显然 PCB=PC1PB,直接对 [PCPB] 行化简可以较快地计算它。

4.7 数字信号处理(Digital Signal Processing,DSP)

此为应用性内容,探讨的是离散时间信号空间 S 及其子空间的一些性质,以及如何用线性变换来处理、过滤和合成信号中包含的数据。

离散时间信号

离散时间信号的向量空间 S 中的一个信号是一个无限的序列 {yk}k 取遍所有整数。以下为若干例子。(不清楚下面的记号字母是不是对的)

  • δ(记号 δ),描述为 {dk}dk={1k=00k0 向量形式 (,0,0,0,1,0,0,0,)

  • 单位阶梯(记号 ν),描述为 {uk}uk={1k00k<0 向量形式 (,0,0,0,1,1,1,1,)

  • 常数(记号 χ),描述为 {ck}ck=1 向量形式 (,1,1,1,1,1,1,1,)

  • 交替(记号 α),描述为 {ak}ak=(1)k 向量形式 (,1,1,1,1,1,1,1,)

  • 斐波那契(记号 F),描述为 {fk}fk={0k=01k=1fk1+fk2k>1fk+2fk+1k<0 向量形式 (,2,1,1,0,1,1,2,)

  • 指数(记号 εc),描述为 {ek}ek=ck 向量形式 (,c2,c1,c0,c1,c2,)

周期信号也是常用信号,其中特别的是信号 {pk},其存在一个正整数 q 使得 kZ,pk=pk+q。正弦信号由周期函数 σf,θ={cos(fkπ+θπ)} 刻画,f,θ 为固定有理数。

线性时不变变换

采用线性时不变(Linear Time Invariant,LTI)变换处理信号,一种方法是在需要时创建信号而非用空间储存信号。

如下例,对仅有的一个信号重复应用移位 LTI 变换可以创建与 δ 类似的信号。

S 是将信号中的每个元素根据 S({xk})={yk}yk=xk1 向右移动的变换,记为 S({xk})={xk1}。将一个信号向左移动,考虑 S1({xk})={xk+1}。显然有 S1S({xk})=xkS1S=SS1=I 是一个恒等变换,即 S 为可逆变换。对 δ 应用移位信号得到的结果是简单的。

容易验证 S 满足线性变换的性质。对任何一个整数 qS({xk+q})={xk1+q},可以把这个性质看作时不变性质。

定义 (线性时不变(LTI)变换)

如果一个变换 T:SS 满足以下条件,称为线性时不变

(i) T({xk+yk})=T({xk})+T({yk}) 对所有信号 {xk},{yk} 都成立。

(ii) T(c{xk})=cT({xk}) 对任意常数 c 和信号 {xk} 成立。

(iii) 如果 T({xk})={yk},则 T({xk+q})={yk+q} 对任意整数 q 和信号 {xk} 成立。

定理 16(LTI 变换为线性变换)

在信号空间 S 上的线性时不变变换是一种特殊类型的线性变换。

数字信号处理

对任意 mN+,周期 m 的移动平均 LTI 变换为

Mm({xk})={yk},yk=1mj=km+1kxk

易证它是一个线性变换。

{pk}S 表示在很长一段时期内每天记录的股票价格的集合(可以对不在研究时间段内的 k 假设 pk=0)。创建一个两天的移动平均转换 M2,它的核 {pk} 满足对于任意 kpk+pk12=0,于是 pk=pk1。从 k=0 开始,核中任意信号都可写成 pk=(1)kp0,一个由 (1)k 描述的交替信号。由于两天移动平均函数的核由所有交替序列的倍数组成,它平滑了每天的波动,而没有趋平总体趋势。

随着 m 增加,应用 Mm 可使信号更平滑。

另一种类型的 DSP 与平滑或过滤相反,它通过组合信号来增加它们的复杂性。音频化是一种用于娱乐行业的处理过程,可使虚拟生成的声音具有更高的音质。

生成 S 的子空间的基

长为 n 的信号集 Sn 被定义为,当 k<0k>n 时,所有满足 yk=0 的信号集合 {yk}

定理 17

同构于 Rn+1 的集合 SnS 的子空间,且信号集合 Bn={δ,S(δ),,Sn(δ)} 形成了 Sn 的一个基。

由于 Sn 有一个有限基,其中任何向量都可以表示为 Rn+1 中的向量。

S2 的基 B2={δ,S(δ),S2(δ)} 表示信号 {yk},其中

yk={0,k<0k>32,k=03,k=11,k=2


{yk}=2δ+3S(δ)+(1)S2(δ)

[{yk}]B2=[231]

有限支持信号的集合 Sf 是信号 {yk} 的集合,其中只有有限多的项是非零的。SfS 的子空间,通过记录股票的每日价格增长而产生的信号猛然是有限的支持,这些信号属于 Sf 而不属于任何特定的 Sn

定理 18

集合 Bf={Sj(δ),jZ} 是无限维向量空间 Sf 的基。

4.8 在差分方程中的应用

信号空间 S 中的线性无关性

考虑一仅包含三个信号 {uk},{vk},{wk} 的集合 S,当方程

c1uk+c2vk+c3wk=0k

蕴含 c1=c2=c3=0 时,{uk},{vk},{wk} 线性无关。

则其同时蕴含

c1uk+1+c2vk+1+c3wk+1=0k

c1uk+2+c2vk+2+c3wk+2=0k

从而

[ukvkwkuk+1vk+1wk+1uk+2vk+2wk+2][c1c2c3]=[000]k

该方程组的系数矩阵称为信号的 Casorati 矩阵,其行列式称为 {uk},{vk},{wk}Casorati 行列式。若对至少一个 kCasorati 矩阵可逆,则上述方程蕴含 c1=c2=c3=0,则三个信号线性无关。一般代入 k=0 较为方便。

Casorati 矩阵不可逆,则相应信号是线性相关或线性无关的。但是若这些信号是同一个齐次差分方程的所有解,则 Casorati 矩阵对所有 k 是可逆的且这些信号线性无关,否则其对所有 k 都不可逆且这些信号线性相关。

线性差分方程

给定数 a0,a1,,an 和信号 {zk}a0an 非零,方程

a0yk+n+a1yk+n1++an1yk+1+anyk=zkk

称为一个 n 阶线性差分方程(或线性递归关系)。为简化,a0 通常取 1。若 {zk} 是零序列,方程是齐次的,否则是非齐次的

在数字信号处理(DSP)中,类似上述的差分方程用来描述一个线性时不变LTI滤波器a0,a1,,an 称为滤波器系数,这里用它们来描述与线性差分方程相关的 LTI 滤波器。

定义

T=a0Sn+a1Sn+1++an1S1+anS0

k,{zk}=T{yk},上述方程描述了两个信号中各项之间的关系。

对滤波器系数为 24,12,24 的滤波器输入两个信号:

{yk}={,cos(0),cos(π4)cos(2π4),cos(3π4),}

它由连续信号 y=cos(πt4) 在整数值 t 抽样而生成,而 {wk}y=cos(3πt4) 以同种方式生成。

会发现 {yk} 能通过而 {wk} 不行,因为输入 {wk} 的结果是零序列。若一个滤波器使 {yk} 能通过而截断高频的 {wk},称该滤波器为低通滤波器

齐次差分方程的解通常具有形式 {yk}={rk} 对某 r 成立。

求解方程

yk+32yk+25yk+1+6yk=0k


rk 代替 {yk} 并因式分解,可得

rk(r1)(r+2)(r3)=0

{1k},{(2)k},{3k} 都是方程的解。

一般而言,一个非零信号 {rk} 满足齐次差分方程

yk+n+a1yk+n1++an1yk+1+anyk=0k

当且仅当 r辅助方程

rn+a1rn1++an1r+an=0

的一个根,不考虑 r 是重根的情形。当这个辅助方程有复根时,差分方程具有形如 {skcoskω}{sksinkω} 的解,s,ω 为常数。

线性差分方程的解集

对于 LTI 转换 T:SS

T=a0Sn+a1Sn+1++an1S1+anS0

将信号 {yk} 变换为信号 {wk},故齐次方程

yk+n+a1yk+n1++an1yk+1+anyk=0k

的解集是 T 的核,且描述了被滤掉或转换为零信号的信号。一个齐次方程的解集也是 S 的子空间,任何解的线性组合仍然是解。

定理 19

an0{zk} 给定,只要 y0,y1,,yn1 给定,方程

yk+n+a1yk+n1++an1yk+1+anyk=zkk

有唯一解。

定理 20

n 阶齐次差分线性方程

yk+n+a1yk+n1++an1yk+1+anyk=0k

的解集 H 是一个 n 维向量空间。

由于 H 是一个线性变换的核,它是 S 的子空间。对 H{yk},设 F{yk}Rn 中向量 (y0,y1,,yn1)。易证 F:HRn 是一个线性变换。对 Rn 中任意向量 (y0,y1,,yn1),由定理 19 知存在 H 中唯一信号 {yk} 使得 F{yk}=(y0,y1,,yn1) 成立。这说明 F 是由 HRn 上的一对一线性变换,即 F 是一个同构,dimH=dimRn=n

回到差分方程

yk+32yk+25yk+1+6yk=0k

若要求出解集的一个基,由于解空间恰好是三维的,{1k},{(2)k},{3k} 构成解空间的一个基。

描述这种差分方程的“通解”的标准方法是对所有解构成的子空间给出它的一个基,这样的基称为它的基础解系

非齐次方程

非齐次差分方程

yk+32yk+25yk+1+6yk=zkk

的通解能写成它的一个特解加上对应的齐次差分方程的一个基础解系的任意线性组合,就像 Ax=bAx=0 的关系一样。它们有相同的意义,xAx{yk}{zk} 的映射都是线性的。

化简成一阶方程组

研究 n 阶齐次线性差分方程的现代方法是用等价的一阶差分方程组代替它,其中一阶差分方程如下:

xk+1=Axkk

其中 xkRnAn×n 矩阵。、

将下列差分方程写成一个一阶方程组:

yk+32yk+25yk+1+6yk=0k


对每个 k

xk=[ykyk+1yk+2]

yk+3=6yk+5yk+1+2yk+2

xk+1=[yk+1yk+2yk+3]=[0+yk+1+00+0+yk+26yk+5yk+1+2yk+2]=[010001652][ykyk+1yk+2]

xk+1=Axk 对所有 k 成立,其中 A=[010001652]

一般而言,方程

yk+n+a1yk+n1++an1yk+1+anyk=0k

可重写成 xk+1=Axk,对所有 k 成立,其中

xk=[ykyk+1yk+n1],A=[010000100001anan1an2a1]

本文作者:SError

本文链接:https://www.cnblogs.com/SError0819/p/18229781

版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。

posted @   SError  阅读(48)  评论(0编辑  收藏  举报
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
收起