求解留数,函数表达式如下

导读:复变函数与实变函数在MATLAB中的计算有着相似之处,因为不管自变量是实数还是复数,都是将自变量的值直接代入函数表达式中去计算。而MATLAB对复变函数和实变函数运算时最大的区别在于MATLAB只对复变函数的主值进行计算。
关键词:MATLAB,复变函数

复变函数理论以其完美的理论与精湛的技巧成为数学的一个重要组成部分。它推动了许多学科的发展,在解决某些实际问题中也是强有力的工具,因此《复变函数》课程是数学与应用数学专业、信息与计算科学专业的主要基础课。和其他数学课程一样,它的学习有着较为枯燥的一面。如何把枯燥的内容变得生动有趣,这是每个授课教师必须要正视的问题。

MATLAB是一种具有强大数值计算,分析和图形处理功能的科学计算语言,其应用领域极为广泛,而且使用方便、调试容易,代码少、效率高,有人称为第四代程序设计语言,本文则把《复变函数》的教学过程和MATLAB结合起来,旨在提高学生学习数学的兴趣,减轻教师的负担,优化学习环境,缩短课时,实现低价高效的教学效果。

1.MATLAB在复变函数计算中的应用

复变函数与实变函数在MATLAB中的计算有着相似之处,因为不管自变量是实数还是复数,都是将自变量的值直接代入函数表达式中去计算。而MATLAB对复变函数和实变函数运算时最大的区别在于MATLAB只对复变函数的主值进行计算。我们学习实变函数里的一些初等函数时,总是先用描点法等作出函数的图像,然后根据图像得出函数的相关性质,而复变函数同样可以采取这样的方法,以增加学生的创新思维和学习兴趣,下面举几个例子加以说明。

例一:计算的函数值、函数值的实部、虚部、辐角、模、共轭函数,并作出函数图像,MATLAB程序如下:

从图中可以看出,为单值函数,的绝对值可以大于1,在图形上轴所表示的函数的实部已经几乎达到60.

例二:计算函数在处的留数。

解:因为在扩充复平面有三个极点,分别为1,-1,,MATLAB程序如下:

2.MATLAB绘图功能在复变函数中的应用

随着计算机处理数据和图形的功能越来越强,复变函数和计算机的结合已经成为必然的选择。比如从定理的推导证明到繁杂的运算,单调乏味,十分影响学习的兴趣。和计算机结合起来就不同了,学生可以把复变函数的一些初等函数,利用绘图功能将该函数用图形直观的表达出来,由此可以通过图形来观察出函数图形的一些性质,这使得教学过程直观生动。

但由于复变函数的自变量是复数,函数值也是复数,就需要有四个量来表示,而MATLAB表现四位数据的方法是用三个空间坐标再加上颜色。下面将通过一些具体算例来说明通过复变函数的图形可以帮助我们理解复变函数课程中一些比较难于理解的抽象的函数概念和性质。

例四:用MATLAB作出复变函数和函数的图像。

图2函数的图像图3函数函数的图像

从图2中可以看出,自变量的取值在水平面的单位园内,x轴是实轴,y轴是虚轴。画函数时,是以坐标系的z轴表示函数的实部,其大小变化范围为,上面的每一个横条都有相同的实部值,因为平面上的颜色表示虚部,从颜色轴对应的数值看出变化范围也是,还可以从图中看出,函数是一个单值函数,它所形成的曲面对应三个高峰和三个低谷,对应函数的实部有三个极大值和三个极小值。并且,还可以从图形中看出该函数是单值的。

从图3中可以看出函数为多值函数,0和为该函数的两个支点。

以上是《复变函数》中一些问题的MATLAB求解。从中可以看出,利用MATLAB求解这些问题具有规范、简洁、灵活等特点;大大简化了数学问题的求解过程,便于求解一些实际应用中较为复杂的数学问题;对于理解掌握《复变函数》理论知识也具有一定的辅助作用。


[1] 彭芳麟.数学物理方程的MATLAB解法与可视化[M].北京:清华大学出版社,2004.
[2]张莹.浅谈MATLAB在复变函数论中的应用[J].沈阳教育学院学报。2005,(7).
[3]路见可,钟寿国,刘士强.复变函数.[M].武汉:武汉大学出版社。2001。

}

1. 微分算子法适用于求常系数线性非齐次微分方程的特解

(A)思想是将求导运算看成线性算符。右边非齐次项仍然是函数,就等价于求一个算符的逆的问题,同时在辅助特征值与特征函数理论,可以求解非齐次项是多项式,指数(三角函数通过欧拉公式化为指数)的形式。

(B)计算中通常结合位移定理以及级数,因式分解,短除法等方法使用。

(C)算子法的优点是能快速得到非齐次部分的特解而不需要特殊记忆特解的形式进行待定系数运算。这种方法能快速得到整个微分方程的通解。

注意: 算子法的题目比较简单无脑,所以近年来直接考通解的题目已经很少,逐渐变为在边界条件下求解的问题,这样算子法的优势就不如后面介绍的拉普拉斯变换方法。

2. 拉普拉斯变换适用于求解带初值条件的常系数线性非齐次微分方程

注意:常用的单边拉普拉斯变换求解的微分方程定义域只能在[0,+无穷)上

(A)拉普拉斯变换的思想是通过一个积分变换,将方程左右两边同时进行变换,使得微分方程变成一个代数方程,求解这个代数方程在做一个反变换,就能得到原方程的解。

(B)拉普拉斯变换的左边,形式与算子法很像,但会多出一个系数部分,为了完成非齐次项的拉普拉斯变换,需要记住常用函数的变换公式。最后还需要反变换回去。

(C)拉普拉斯变换法的优点在于能将题目所给的初值条件直接融入到等式中,而不需要最后再来待定系数。这种方法在某些问题上具有奇效。

来看拉普拉斯变换的定义等式

显然,要计算的就是F(0)

注意: 对这类问题拉普拉斯变换有奇效是因为巧妙的回避了反变换问题,仅在s域就解决了问题,如果要进行反变化,并不一定特别快速。但优点是反变换之后得到的就是原初值问题的解,而不需要在解线性方程确定系数。

1. 格林函数方法的思想是非齐次项可以通过叠加原理分解为一连串的冲击响应,只要求解系统的冲击响应,那么整个响应可以通过冲击响应与输入函数卷积得到。格林函数方法可以直接给出方程的形式上的积分解。这点也可以从拉普拉斯变换看出来,两个拉普拉斯变化的乘积,返回到时域就是卷积。

2. 直接求解系统的冲击响应并不是一件容易的事,而且格林函数的定义其实严格说来也依赖于初值条件。所以这个问题可以又转化为拉普拉斯变换。

3. 可以这么理解,格林函数实际给出了方程的通解的形式,这种形式在特定问题下往往很有用。

(1)  y = x -1, 这是可以直接观察出来或者用微分算子法快速求出来.

所以方程的通解是Ae(-x) +x-1

(2)构造一个方程的特解如下(空降)

然后证明这是一个周期函数

注意:实际问题中格林函数的坑(因果性)

}

我要回帖

更多关于 MATLAB表达式 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信