复数的矩阵表示

复数

复数z定义:a+bi

其中a,bRi2=1i又称为圆复数(虚数)单位,a为实部Re(z)=ab为虚部Im(z)=b,复数域记作C

0.复数三角形式和指数形式

  • z=a+bi=r(cosθ+isinθ)=reiθ

1.复数域是实数域的代数(加法和乘法)闭包,定义复数域内的加法和乘法

  • (a+bi)+(c+di)=(a+c)+(b+d)i
  • (a+bi)×(c+di)=(acbd)+(bc+ad)i
  • 加法单位元为0,乘法单位元为1
  • 加法逆元唯一,除零乘法逆元唯一
  • 对加法和乘法满足交换律,结合律,分配律

2.复数的模

  • z=a+bi的模:|z|=a2+b2

3.复数的辐角

  • 指数形式下,z的模等于rθ称为z的辐角,记作Arg(z)
  • θ[π,π],则称为z的主辐角,记作arg(z)

4.共轭复数

  • z=a+bi的共轭复数:|z^|=abi
  • z与其共轭复数的模相等,二者乘积为模长的平方

复数的矩阵映射

除了圆复数外,对偶数和双曲复数都具有同构矩阵

对于复数z=a+bia,bR,复数到特定形式二阶矩阵的同构映射

a+bi=(abba)

显然,该类矩阵对加法和乘法满足交换律,结合律,分配律;纯虚数对应的矩阵为反对称矩阵(AT=A,主对角线元素必为0)

(a+bi)+(c+di)=(a+c(b+d)b+da+c)=(c+di)+(a+bi)

(a+bi)×(c+di)=(abba)(cddc)=(acbd(bc+ad)bc+adacbd)=(acbd)+(bc+ad)i

1.虚数i,0,1

i=(0110)0=(0000)1=(1001)

2.共轭复数=矩阵的转置

z^=abi=(abba)=(abba)T

对于负单位模长(假设存在)x2+y2=1,取P=(xyyx)=PT,那么:
P1=(xx2+y2yx2+y2yx2+y2xx2+y2)=(xyyx)=P

PZP1=Z^,满足群元素的共轭,那么P对应的复数(又一种复数?)又是怎样的?

因为x2+y2=1,令x=icosθ,y=isinθθ为实数,则

P=(xyyx)=(icosθisinθisinθicosθ)=PT

T=( cosθ sinθ sinθ cosθ),发现T=T1=TTT2=E

T的转置等于其逆,且T为正交矩阵,矩阵T有何特殊之处

共轭复数与群元素共轭是两个不同的概念

  • 复数域满足交换律,按照群元素共轭概念来来,只有实数同时具有元素共轭且满足共轭复数的定义
  • 若群也满足交换律,那么群内元素的共轭为该元素本身

3.加法逆元

abi=(abba)

4.模长的平方=行列式

  • 模=行列式开方

|z|2=|zz^|=|abba|=a2+b2

5.乘法逆元=转置矩阵/行列式的平方

z1=(abba)/(a2+b2)2

6.三角形式

a+bi=r(cosθsinθsinθcosθ)

所以复数(z=a+bia,bR)对加法和乘法构成域,同构映射到特定形式的矩阵((abba))对矩阵加法和矩阵乘法构成矩阵域

复数只包含两个参数,2×2矩阵包含四个参数,所有必可以找到一类矩阵与复数同构

7.矩阵的对角矩阵或相似矩阵

a+bi(abba)=(ii11)(aib00a+ib)(i212i212)=P(aib00a+ib)P1

(abba)有两个特征向量(特征值):v1=(i1), 特征值 λ1=aibv2=(i1), 特征值 λ2=a+ib

其中,P=(ii11)P1=(i212i212)

a+bi=abi(abba)=(ii11)(aib00a+ib)(i212i212)=Q(aib00a+ib)Q1
(abba)有两个特征向量(特征值):v1=(i1), 特征值 λ3=aibv2=(i1), 特征值 λ4=a+ib

其中,Q=(ii11)=2(P1)TQ1=(i212i212)=12PT

8.幂

(a+bi)2=(abba)2=(a2b22ab2aba2b2)

(a+bi)3=(abba)3=(a33ab2b33a2bb3+3a2ba33ab2)

(a+bi)4=(abba)4=(a4+b46a2b24ab34a3b4ab3+4a3ba4+b46a2b2)

(a+bi)n=(abba)n=(ii11)(aib00a+ib)n(i212i212)

9.秩

rank(abba)={2;a0,b01;a0,b=02;a=0,b00;a=0,b=0

10.矩阵指数

ea+bi=exp(abba)=(eacos(b)easin(b)easin(b)eacos(b))

11.矩阵乘法性质

  • 乘法结合律: (AB)C=A(BC)
  • 乘法左分配律:(A+B)C=AC+BC
  • 乘法右分配律:C(A+B)=CA+CB
  • 对数乘的结合性:k(AB)=(kA)B=A(kB)
  • 转置: (AB)T=BTAT

四元数的矩阵映射

哈密顿放弃乘法交换律,i扩充成三个虚单位元i,j,k,且有i2=j2=k2=1,ijk=1,ij=ji,四元数定义为q=x+yi+zj+wkx,y,z,wR

1.四元数的模长

|q|=x²+y²+z²+w²

模长为1的四元数为单位四元数

2.四元数的共轭

q^=xyizjwk

qq^=|q|2

3.没有良好的三角形式,只有方向余弦

q=r(cosα+icosβ+jcosγ+kcosθ),其中αβγθq与四个坐标系的夹角,rq的模长;所以使用四元数表示旋转时,没有旋转矩阵的角度的意义显然

4.两种矩阵表示

4.1 表示方法一,2\*2复矩阵
  • q=x+yi+zj+wk(xwiy+ziy+zix+wi)

  • 那么四元数的加法与矩阵的加法相容,四元数的乘法也与矩阵的乘法相容,而四元数的模长为矩阵行列式的平方根,四元数的共轭就是矩阵的共轭转置。

  • z=w=0时,四元数q=x+yi+zj+wk退化为复数x+yi

  • i(0110)j(0ii0)k(i00i)1(1001)

    4.1 表示方法二,44实矩阵

  • q=x+yi+zj+wk(xywzyxzwwzxyzwyx)

  • 四元数的模长是矩阵行列式的平方根,四元数的共轭是矩阵的转置

  • k=ij,可以得到q=x+yi+zj+wij=(x+yi)+(z+wi)j,记a=x+yib=z+wi,即q=a+bj,其中a,bC

5.四元数的乘法

  • 不满足交换律
  • 结合叉乘和点乘

平移、缩放、旋转变换矩阵

平移矩阵(Translation)

(x,y,z)表示当前位置, (x,y,z)表示新的位置, (dx,dy,dz)平移的量。

2D平移

(xy)=(xy)+(dxdy)

(xy1)=(10dx01dy001)(xy1)=(1(x)+0(y)+dx(1)0(x)+1(y)+dy(1)0(x)+0(y)+1(1))=Txy(xy1)

3D平移

(xyz)=(xyz)+(dxdydz)

(xyz1)=(100dx010dy001dz0001)(xyz1)=(1(x)+0(y)+0(z)+dx(1)0(x)+1(y)+0(z)+dy(1)0(x)+0(y)+1(z)+dz(1)0(x)+0(y)+0(z)+1(1))=Txyz(xyz1)


缩放矩阵(Scaling)

(x,y,z)表示未缩放前的原始位置,(x,y,z)表示经过缩放后的新位置,SxSySz分别表示在x轴、y轴和z轴方向上的缩放因子。

注意:计算多个点的缩放,需要将每个点位置分别代入公式计算。

2D缩放

(xy1)=(Sx000Sy0001)(xy1)=(Sx(x)+0(y)+0(1)0(x)+Sy(y)+0(1)0(x)+0(y)+1(1))=Sxy(xy1)

3D缩放

(xyz1)=(Sx0000Sy0000Sz00001)(xyz1)=(Sx(x)+0(y)+0(z)+0(1)0(x)+Sy(y)+0(z)+0(1)0(x)+0(y)+Sz(z)+0(1)0(x)+0(y)+0(z)+1(1))=Sxyz(xyz1)


旋转矩阵(rotation)

(x,y,z)表示未旋转前的原始位置,(x,y,z)表示一个点经过旋转后的新位置, $$为在右手坐标系中,逆时针旋转角度(编程中以弧度为单位) ,

注意:计算多个点的旋转,需要将每个点位置分别代入公式计算。

2D旋转

(xy1)=(cosθsinθ0sinθcosθ0001)(xy1)=R2d(xy1)

注:θ为平面内绕原点逆时针旋转角度

绕Z轴3D旋转,Z轴指向自己,X正半轴旋转90°到Y正半轴为逆时针

(xyz1)=(cosθsinθ00sinθcosθ0000100001)(xyz1)=Rz(xyz1)

绕X轴3D旋转,X轴指向自己,Y正半轴旋转90°到Z正半轴为逆时针

(xyz1)=(10000cosθsinθ00sinθcosθ00001)(xyz1)=Rx(xyz1)

绕Y轴3D旋转,Y轴指向自己,Z正半轴旋转90°到X正半轴为逆时针

(xyz1)=(cosθ0sinθ00100sinθ0cosθ00001)(xyz1)=Ry(xyz1)

坐标轴具有轮换对称性,Rz通过行line或列column交换(偶次数)可变成Rx,具体操作如下,两次行交换+两次列交换(保持坐标系的手性进行,对空间向量进行线性变换等价于对坐标系进行逆变换):

  • L3L2
  • L1L2
  • C3C2
  • C1C2

同理,Rz通过行line或列column交换(偶次数)可变成Ry,具体操作如下:

  • L3L2
  • C3C2
  • L1L3
  • C1C3

此外,还有倾斜或者剪切矩阵,就是将向量某一投影方向进行缩放,即对其分量(点乘基向量)进行缩放

复变函数

复变函数:以自变量为复数,因变量亦为复数的函数

复数z=x+yiw=f(z)=u(x,y)+iv(x,y),其中x,yRu(x,y),v(x,y)是实变函数,那么w=f(z)就是复变函数

一个复变函数的值域取决于以下几个因素:

  1. 定义域:定义域是复变函数自变量取值的所有可能点构成的集合。一个函数的值域首先受限于其定义域的选取。
  2. 连续性:一个复变函数如果在定义域内连续,那么其值域至少是连续的。连续性保证了函数在某个区域内不会出现跳跃,从而有助于确定值域。
  3. 解析性:复变函数如果在其定义域内是解析的,即存在导数,那么其值域通常更加广泛。解析函数具有很好的局部性质,其值域可以包含复平面上的任意点。

详细地,对于具体的复变函数,可以通过以下方法探究其值域:

  • 对于简单的复变函数,如多项式函数(f(z)=ki=0aizi ),其值域是整个复平面,因为多项式函数可以取到复平面上任意点。
  • 对于有理函数,其值域则可能是一个去掉了一些点的复平面。这是因为有理函数在某些点可能无定义或者趋于无穷。
  • 对于指数函数、对数函数和三角函数等初等函数,它们的值域通常是复平面上的某个区域。

复数可视化

复变函数w=f(z)是将复数映射到复数,即输入z包含实部和虚部的二维信息,输出w包含实部和虚部的二维信息,使用三维空间是不足以描述其图像的

1. 3D plots

不同颜色的三维曲面

在复平面上增加一个轴R,复平面表示输入,轴R表示输出的模长,点采用颜色表示输出的相位(辐角),色相环正好对应[0,2π)

当然,轴R表示输出的实部,点采用颜色表示输出的虚部——没有上述巧妙

可以描述连续变化

2. Vector fields

不同颜色和长度的二维向量

在复平面上标记有限个点(一般选取整数格点),格点表示输入,使用带颜色的不同模长的矢量表示输出

3. Domain colouring

不同颜色和亮度的平面

用平面上点的位置代表自变量的实部和虚部,用亮度和色相分别代表函数值的模长和辐角。

4. z-w planes

弯曲坐标系,把自变量的点弯曲到函数值

images/复数的矩阵表示-20240612161419168.webp

images/复数的矩阵表示-20240612161334639.webp

儒可夫斯基函数w=12(z+1z)

images/复数的矩阵表示-20240612161554167.webp

5. Riemann sphere

平面作为自变量,球面作为因变量,把平面放在球面上

images/复数的矩阵表示-20240612161641673.webp

6.自动控制原理中的根轨迹

传递函数一般是分式线性映射,讨论根的分布

复变函数的矩阵映射

1.非线性多项式复变函数的矩阵映射

f(z)=z2+z+2=(x+yi)2+(x+yi)+2=(x2y2+x+2)+(2xy+y)i

那么

f(A)=A2+A+2E=(xyyx)2+(xyyx)+2(1001)=(x2y2+x+2(2xy+y)2xy+yx2y2+x+2)

例如,

1+2if6i

(1221)f(0660)

能否将f变换成矩阵乘法,假设可以,那么这里等价于f(A)=PAQ,由于矩阵乘法是线性变换,所以

f(A)=PAQ=(m11x+n11ym12x+n12ym21x+n21ym22x+n22y),其中m,n,bR

显然,(m11x+n11ym12x+n12ym21x+n21ym22x+n22y)(x2y2+x+2(2xy+y)2xy+yx2y2+x+2)

所以非线性变换无法等价于矩阵乘法;P的元素pij均为常数,A的元素aij均为变量,则PA的元素必为线性pkjajk(不是一次,没有截距)

2.一次复变函数的矩阵映射

f(z)=kz+a+bi=k(x+yi)+by(x+yi)+abxy ,其中x,y,k,a,bR,y0

y0f(z)=kz+b形式,故以此为通用的线性复变函数形式

那么

f(A)=kA+bE=(kx+bkykykx+b)

同上,二阶矩阵A通过乘法,PAQ结果中的元素不会出现常数项b

显然,对平移变量b,需要进行升维才能处理

z=x+yi(xy0yx0001)

f(z)=kz+b=kx+b+kyi(kx+bky0kykx+b0001),其中k,bR

P(xy0yx0001)Q=(kx+bky0kykx+b0001)

同样不能找到PQ,因为f(z)不是线性函数(有疑问,线型函数是描述直线方程的一个比较恰当的同义词),仿射变换(其常数项不为0),其集合空间中的元素为矩阵——一次函数不一定是线性的

  • f(z1+z2)=k(z1+z2)+bf(z1)+f(z2) f(A1+A2)=k(A1+A2)+bEf(A1)+f(A2)
  • f(mz1)=k(mz1)+bmf(z1)f(A1)=k(A1)+bEkf(A1)

矩阵空间

线性空间V中的一个变换A,要验证它是否为一个线性变换

只要看对于V中任意的元素αβ和数域P中任意k,是否都有

  • A(α+β)=A(α)+A(β)
  • A(kα)=kA(α)

其中A(α)称为α在变换A下的像,可以表示为 ()=,其中A是一个矩阵,x是一个向量

线性空间是由向量集合V、域P、加法运算+和标量乘法(数乘)组成的模类代数结构

则称代数系统(V,+,,P)V关于+,构成P上的一个线性空间P为线性空间的基域V中元素称为 向量P中元素称为标量。当域P为实数域时,称为实线性空间。当域P为复数域时,称为复线性空间。

矩阵空间是由矩阵集合T、域P、加法运算+和标量乘法(数乘)组成的模类代数结构

只要看对于T中任意的元素AB和数域P中任意k,是否都有(P为矩阵作用)

  • P(A+B)=P(A)+P(B):乘法分配律
  • P(kA)=kP(A)

显然,矩阵空间T是线性空间,所以函数作用f(A)=PA=nA,结果必为线性

而且复数只能表示二维平面上的点,不如直接矩阵来得直观;且复变函数包含平移、缩放和旋转等操作,用矩阵来研究复变函数多有不便

直接映射的想法是行不通的(哪里出差错?,一次函数不一定是线性变换,所以若要使得函数作用在矩阵上等价于矩阵乘法,该函数必须为线性函数)

增广矩阵,解非齐次方程,仿射变换,将平移量放在第三个基上[0,0,1]T,用于保存b的信息

缩放矩阵和平移矩阵来描述复变函数

z=x+yi(xy1)

f(z)=kz+b+ci=kx+b+(ky+c)i(kx+bky+c1),其中k,bR

f(z)=kz+b+ci表示对z同方向缩放k倍后,再向实轴正方向平移b,向虚轴正方向平移c

(kx+bky+c1)=TxySxy(xy1)=(10b01c001)(k000k0001)(xy1)=(k0b0kc001)(xy1)

所以只能把复数当作一维向量看待,才能完成同构,此时

f(k0b0kc001)

3.复变函数的保角映射(共形映射)

映射w=f(z)是区域D上的保角映射的充分必要条件f(z)D内解析, 并且f(z)0

images/复数的矩阵表示-20240612095032588.webp

3.1 分式线性映射

w=az+bcz+d,其中a,b,c,d为复常数,且adbc0

显然,倒数映射(又叫反演映射)w=1z是保角映射

平移变换(w=z+b,b为复常数),等比例缩放(w=kzk0)和旋转变换(w=eiθz)都是保角映射

分式线性映射是平移、旋转、相似、反演映射 的复合映射

性质

  • 一一映射,双方单值映射
  • 保角性
  • 保圆性,圆周映射后还是圆周(直线看作半径为无穷大的圆周)
  • 保对称性

3.2 指数函数的保角映射

w=ez在复平面内处处解析,且w0

3.3 儒可夫斯基函数,用于解决机翼截面的绕流

w=12(z+1z)

参考资料与注释

作者:invo

出处:https://www.cnblogs.com/invo/p/18243532

版权:本作品采用「署名-非商业性使用-相同方式共享 4.0 国际」许可协议进行许可。

posted @   Invo1  阅读(782)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· DeepSeek 开源周回顾「GitHub 热点速览」
· 记一次.NET内存居高不下排查解决与启示
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· .NET10 - 预览版1新功能体验(一)
more_horiz
keyboard_arrow_up light_mode palette
选择主题
点击右上角即可分享
微信分享提示