复变函数可视化以及代数基本定理

22 篇文章 0 订阅

复变函数:

我们生活的世界是三维的,需要三个坐标信息才能清晰描述一个物体的位置,所以,人类很容易就能够理解三维以下的物体形状,二维笛卡尔坐标系和三维立体坐标系就是用来分析此类问题的数学工具。

但对于复变函数,情况有所不同。

以下面的复函数为例

f(z)=z^3 \quad z\in C

就无法像实变函数那样,在二维或者三维笛卡尔坐标系下绘制其图像。因为想要绘制出一对复数(z,f(z)),需要4个维数,两个用来画z=x+iy,两个用来画f(z)=u+iv.

将函数展开,令

\\f(z)=X+Yi \\ z=x+yi,

则可以得到:

\\X+Yi=(x+yi)^3=(x^2-y^2+2xyi)(x+yi)\\=x^3-xy^2+2x^2yi+x^2yi-y^3i-2xy^2=(x^3-3xy^2)+(3x^2y-y^3)i

所以:

 \left\{\begin{matrix} X(x,y)=x^3-3xy^2\\ Y(x,y)=3x^2y-y^3 \end{matrix}\right.

所以很明显,一元复变函数可以分解成两个二元函数,所以,无论在二维还是三维笛卡尔坐标系中,都无法表示四元

x,y,X(x,y),Y(x,y)

坐标的位置。

二维复空间可以解释做实的,但是四维的空间. 空间维度本来就不够用,笛卡尔坐标系用两个实轴来表示,委实有点浪费,实际上一维的实数空间用一个轴就够了。

在matlab中,是通过分解函数的实部和虚部来实现的,实部用三维笛卡尔坐标系的z轴表示,虚部用颜色表示,xOy表示整个复平面,绘制出来是这个样子:

颜色分析并不直观,有没有办法可以将虚轴的值也通过坐标截距表示出来呢?

问题转化为,如果用三维空间笛卡尔坐标系标识四维函数的问题,方法有两个:

1.用两个笛卡尔坐标系, 分别绘制X(x,y),Y(x,y)

2.在一个笛卡尔坐标系内,分开绘制X(x,y),Y(x,y)

这有点像变成领域的,空间时间互换问题,二者不可兼得,这也是大自然的普遍规律。

下面采用第二种,用geogebra 将

\left\{\begin{matrix} X(x,y)=x^3-3xy^2\\ Y(x,y)=3x^2y-y^3 \end{matrix}\right.

分别绘制在一个笛卡尔坐标系内:

图有些乱,但大体上能够看出绿色和品红色的是两个不同函数的图形,绿色是X(x,y),品红色是Y(x,y)左边是整个立体图形在XOY平面上的投影,立体上,两个图形完全一样,只是旋转了角度。

XOY平面上一条条的是图形的等高线,也是图像上具有相同实部和虚部的点在XOY平面上的投影,每个投影有三条分支,围绕中心原点在圆上呈现对称分布,由于两个图形完全相似,所以XOY平面有两组一共六条投影线,

每两条分别来自不同函数的等高线的交点在复平面上的坐标表示方程

f(z)=z_0, \quad z_0\in C

的一个解,也就是图中等边三角形的三个顶点,所以可以看到,对于任何z_0,方程均有三个解,围绕圆周均匀分布。

z_0=0的时候,也就是

\left\{\begin{matrix} x^3-3xy^2=0\\ 3x^2y-y^3=0\end{matrix}\right.=>\left\{\begin{matrix} x(x-\sqrt{3}y)(x+\sqrt{3}y)=0\\ y(\sqrt{3}x-y)(\sqrt{3}x+y)=0\end{matrix}\right.=>\left\{\begin{matrix} x=0,y=\frac{\sqrt{3}}{3}x, y=-\frac{\sqrt{3}}{3}x\\ x=0,y=\sqrt{3}x, y=-\sqrt{3}x \end{matrix}\right.

得出的方程是每个图像的各三条渐近线.

模曲面:

从几何的角度看问题,在每个复数点z上安置立坐标垂线,它的长度等于多项式f(z)在这一点的模,这些立坐标的端点形成某一曲面M,我们把它叫做多项式f(z)的模曲面,这样就可以将复变函数的实部和虚部一起考虑,而不用分别

绘制实和虚部,模曲面实际上将实部和虚部作为一个整体来看,这样就减少了一个维度,四位塌缩成了三维,就可以画在三维坐标系上了。

模曲面有一些性质:

1.在全部复平面内,模曲面不会下降到复平面以下,模不可能是负的。

2.对于复平面上任何一点z,这个曲面上有一点,并且仅有一点位于这个点的上面或者就是这个点本身,说明与复平面接触。

3.当z在复平面上连续运动时,这个曲面的点的立坐标\left \| f(z)\right \|连续变动。

4.由于多项式的根要么是实根,要么是共轭复根,因为假如复数a是多项式的一个根,那么:

  f(a)=0=>f(\bar{a})=\bar{f}(a)=\bar{0}=0

\bar{a}也是一个根,所以M曲面一定关于实轴对称.

z=x+yi

对于一次函数

f(z)=z+a=x+yi+a=(x+a)+yi

模曲面与复平面有一个触点,并且固定在x轴上,坐标为

(-a,0)

对于二次复变函数:

f(z)=z^2+az+b=x^2-y^2+ax+b +(2xy+ay)i

与复平面的触点有两个,要么是共轭复根,要么是两个实根,这也符合二次函数根的分布情况。

下面绘制

f(z)=z^3+az^2+bz+c

的模曲面:

则:

\\f(z)=(x+yi)^3+a(x+yi)^2+b(x+yi)+c=\\(x^3-3xy^2)+(3x^2y-y^3)i+a(x^2-y^2)+a(2xy)i + bx+byi +c =\\ (x^3-3xy^2+ax^2-ay^2+bx+c)+(3x^2y-y^3+2axy+by)i

所以

\\ \left \| f(z) \right \|=\left \| (x^3-3xy^2+ax^2-ay^2+bx+c)+(3x^2y-y^3+2axy+by)i \right \|=\sqrt{(x^3-3xy^2+ax^2-ay^2+bx+c)^2+(3x^2y-y^3+2axy+by)^2}

图形为:

三个实根的情况:

octave求得算数解:

对于四次的情况:

\\f(z)=z^4+az^3+bz^2+cz+d=\\x^4+4x^3yi+ax^3-6x^2y^2+3ax^2yi+bx^2-4xy^3i-3axy^2+2bxyi+cx+y^4-ay^3i-by^2+cyi+d=\\ (x^4+y^4-6x^2y^2+ax^3-3axy^2+bx^2-by^2+cx+d)+(4x^3y-4xy^3+3ax^2y-ay^3+2bxy+cy)i

其对应的模曲面形状如下所示,有四个点和复平面接触,为其零点。

图中可见,方程的次数和根的个数似乎有某种对应关系,事实上,次数等于根树,这是代数基本定理,模曲面方程和复平面的n个接触点为零点,就好像有n根支柱将整个曲面支撑起来,不断变化a,b,c,d的值,这三个支柱相对位置在复平面上会发生变化,但始终和复平面有接触,根据达朗贝尔引理(了解就好,对理解关系不大),曲面上的任何一个极小值点都是0,没有一个极小值是正的非零值。

这表示复平面上始终有三个点,使\left \| f(z)\right \|=0,由于复数等于零,与它的模等于零是同义的,\left \| f(z)\right \|为0的充分必要条件是

\dpi{120} \left \| f(z)\right \|=0<=>f(z)=0

等价于,多项式

f(z)=0

有三个根。

这是由代数基本定理决定的:代数基本定理是说,任何复系数一元n次多项式 方程在复数域上至少有一根(n≥1),由此推出,n次复系数多项式方程在复数域内有且只有n个根(重根按重数计算)。代数基本定理在代数乃至整个数学中起着基础作用。 据说,关于代数学基本定理的证明,现有200多种证法。

 代数基本定理强调的有一个根,而不是n个根,之所以这样,是因为多根可以由1个根推导出来,设:

f(z)=z^n+a_1z^{n-1}+\cdots++a_n

z_0是f(z)的一个根,用

z-z_0

去除以f(z),由于除式是1次的,因此余数就是一个常数R(根据多项式整除规则,如果遗留带z的多项式做余数,则整除规则可以继续下去,直到只剩下常数项),所以:

f(z)=(z-z_0)f_1(z)+R

式中f_1(z)是一个n-1次的多项式,而R是常数,用z_0代入:

f(z_0)=(z_0-z_0)f_1(z)+R=R

由于

f(z_0)=0

所以

R=0

f(z)=(z-z_0)f_1(z)

同理,如果代数基本定理正确的话,必定

f_1(z)=(z-z_1)f_2(z)

..............

所以,如果存在1个根,就一定存在n个根,代数基本定理只需要说明1根的情况即可。

这n个根之间没有要沿实轴在整个复平面上呈现对称分布,除此之外,没有其它形状要求(比如是否一定要共圆?下图五个根在实部为-2的复直线上,显然不可能共圆).

尝试以几何的方式去理解:

对于形如:

f(z)=z^n+c_1z^{n-1}+\cdots+c_{n-1}z+c_n \qquad c_1\quad c_2 \quad \cdots \quad c_n 为复数

将形式改写为:

 f(z)=z^n(1+\frac{c_1}{z}+\frac{c_2}{z^2}+\cdots+\frac{c_{n-1}}{z^{n-1}}+\frac{c_{n}}{z^{n}})     

明确几点:

1.f(z)是从复数域到复数域的一个映射.

 2.当\left \| z \right \|增大时,\left \| f(z) \right \|无限接近\left \| z^n\right \|,这是由于

\\\left \| f(z) \right \|=\left \| z^n(1+\frac{c_1}{z}+\frac{c_2}{z^2}+\cdots+\frac{c_{n-1}}{z^{n-1}}+\frac{c_{n}}{z^{n}}) \right \|<=\\ \left \| z^n \right \|(1+\left \| \frac{c_1}{z} \right \|+\left \| \frac{c_2}{z^2} \right \|+\cdots+\left \| \frac{c_{n-1}}{z^{n-1}} \right \|+\left \| \frac{c_{n}}{z^{n}} \right \|)\approx \left \| z^n \right \|  

接下来,设复平面上一圆C_r,圆心在原点,半径r足够大。

由于r足够大,当z在C_r上运动时,f(z)主要由其高次项决定,也会非常接近

  g(z)=z^n

当z围绕C_r旋转一圈时,根据复数乘法法则

\\z=r(cos(\theta)+isin(\theta))\\ g(z)=r^n(cos(n\theta)+isin(n\theta))

\theta=2\pi时候,g(z)的幅角运动2n\pi,也就是旋转n 圈。

由于r足够大时,f(z)g(z)足够接近,所以也旋转n 圈。

然后,想象一下

r\rightarrow 0

当z的辐径逐渐收缩到复平面原点的时候,

f(z)\rightarrow c_n

由于r足够大的时候,原点和c_n都包含在C_r内部(总是可控的,只要r足够大,一定满足这个条件) ,那么f(z)在收缩过程中,必定扫过原点,也就是说

存在最少一次的
f(z)=0

得证。

结合M曲面图,

\left \| f(0) \right \|=\left \| c_n \right \|

\left \| f(r(cos\theta+isin\theta)) \right \|=足够大
并且上面已经证明了会经过原点。

所以,在收缩过程中,M曲面必定有点与复平面“摩擦"接触!

联想到自动控制原理中介绍跟轨迹问题,也可以用M曲面理论去理解.

例如:

对于开环传递函数为

  H(s)=\frac{K}{s(0.5s+1)}

的控制系统,如下图所示:

其闭环传递函数为:

G(s)=\frac{H(s)}{1+H(s)}=\frac{2K}{s^2+2s+2K}

所以特征方程

s^2+2s+2K=0

的根是:

\\s_1=-1+\sqrt{1-2K} \\s_2=-1-\sqrt{1-2K}  

如果想知道随着开环增益K变化,控制系统的稳定性如何变化,可以绘制特方程的M曲面:

s=x+yi

f(s)=(x+yi)^2+2(x+yi)+2K=(x^2-y^2+2x+2K)+(2xy+2y)i

所以,M曲面方程为:

M(x,y)=\sqrt{(x^2-y^2+2x+2K)^2+(2xy+2y)^2}

K=0,其图像是两个不同的实根:

K=0.5,两个相同的实根

K>0.5的时候,是两个分离的复数根.

所以,复平面上的跟轨迹实际上是:

kaihuanfenzi=[0,2];
kaihuanfenmu=[1,2,0];
sys=tf(kaihuanfenzi,kaihuanfenmu);
rlocus(sys);
axis([-8 2 -2 2]);

所以控制系统的稳定性问题实际上是多项式在复平面上的根分布问题.

数形结合能够反映代数和几何之间的深层联系,geogebra是数形结合分析问题的利器。

结束!

  • 10
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

papaofdoudou

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值