复变函数知识梳理

历史寻根复变函数的前世今生:

历史上第一个遇到‘虚数’即复数的人是印度的数学家Bhaskara Achary(约114-1185)他在解方程时候认为x^2=-1​是没有意义的。1484年法国数学家N.ChuQuest(约1445-1500)在解方程x^2-3x+4=0​时候得到的根是x=\frac 12 \pm\sqrt{\frac 49-4} ​,他被这个"怪数"弄的不知所措。1545年,意大利数学家G.Cardano(1501-1576)在解方程x(10-x)=40​时候,把这个方程的两个根5\pm\sqrt{(-15)}​,从而引进了复数。与他同时期的另一位意大利数学家Rafael Bombelli(1526-1572) 在其《代数》一书中从已知的实数运算法则类推出了复数的四则运算。
1629年,荷兰科学家A.Girard(1595-1632)在其《代数新发明》一书中引入了符号\sqrt {(-1)}​来表示虚数单位。稍后,法国科学家R.Descartes(1596-1632)用i​来记\sqrt{(-1)}​,并且第一次使用“复数”,“虚数”这些概念。1843年,瑞士科学家L.Euler(1707-1783)发现了著名的欧拉公式e^{i\varTheta}=cos\vartheta+isin\vartheta​。1797年,丹麦数学家C.Wessel在坐标平面上引进了实轴与虚轴,使得复数a+bi​与平面上的点一一对应,从而使得“复数”有了“立足之地”。
此后,爱尔兰数学家W.R.Hamilton(1805-1965)发展了复数的一个代数解释:每个复数都可以用一个实数对(a,b)​表示。18世纪以后,以欧拉为首的数学家们发展起来了一门新的数学分支—复变函数论。19世纪以后,法国数学家柯西,德国数学家黎曼,威尔士特拉斯等人使复变函数论得到巨大的发展,并且广泛地应用到空气动力学,流体力学,电学,热学等等方面


写到这里我要说点题外话了,情不知所起,科学家真的很伟大我真的很佩服他们,无数的科学家呕心沥血前赴后继才有了科学技术才有了人类文明的繁荣昌盛,才让我们人类拥有了改造大自然,利用大自然的能力,我想在这里向伟大的数学家,伟大的物理学家,伟大的生物学家,一切为人类社会做贡献的人致敬!

(笔者大学专业是信息与计算机科学,学了一半数学也学了一半计算机,个人认为我还是对计算机更加的感兴趣,当然数学也感兴趣但是总体而言我对计算机的热情要高过数学,因此笔者学数学主要在吸收数学思想知道如何使用数学工具,不会去像其他人一样把每一个证明过程都牢牢的记在心里,我只注重实用,会忽略到许多细枝末节文章有浅薄之处还请各位高抬贵手放小弟一马!)


复变函数知识树如下:

复变函数与积分变换.png

复数与复变函数

浅谈复数:

复数的一般形式为z=x+iy 具有性质 i^2=-1显然复数由实部(x)和虚部(y)构成实部和虚部(x,y)很自然的对应着直角坐标系上的一点因此我们就可以把复数与向量联系起来,如此复数的加法和减法都可以类比向量的运算,复平面上的一个向量的模长就是其对应复数的模设若负向量与x轴的夹角为\varTheta则有\begin{cases} x=rcos\vartheta\\y=rsin\vartheta\\\end{cases}于是复数又可以表示成z = x + iy=r(cos\vartheta+isin\vartheta)又由Euler公式e^{i\vartheta}=cos\vartheta+isin\vartheta可得z=re^{i\vartheta}这种形式称为复数的指数表达式

何为复变函数?

所谓复变函数形如\omega=f(z),z=x+iy复数的函数\omega=u+vy,复变函数的英文拼写为(Complex Variables),所谓复变函数即为变量为复数的函数。实际上\omega的实部uv往往可以表示为一个含x,y的二元函数即 z=u(x,y)+iv(x,y),这样我们研究复变函数就只要分别研究复变函数所对应的实部函数u(x,y)和虚部函数v(x,y)

研究复变函数的tools

在数学分析或者高等数学中我们研究实变函数式通过研究函数的可导性,连续性,和极限来刻画实变函数的性质的,于是我们伟大的数学家们希望这些成熟的工具拿过来刻画复变函数,实践证明这是可行的,实际上我认为数学家们研究复变函数并不仅仅是因为数学游戏好玩兴趣使然才驱使数学家绞尽脑汁的去研究这么个东西,而是因为19世纪-20世纪物理学的蓬勃发展因为当时的科学技术需要才催生出的这么一门学科,如果解决了复数问题,那么生活中的许多问题都可以迎刃而解了。数学本就是和艺术一样源于生活又高于生活但是却与生活息息相关的一门科学,而不是什么披着宗教外衣的伪科学。

解析函数的概念与解析的充要条件

毫无疑问要想把复变函数像实变函数那样研究并且希望找到实变函数与复变函数桥梁把二者联系起来,那么复变函数的连续性与可导性,复变函数的积分性质以及复变函数在级数上的表实就必须弄清

复变函数的连续性和可导性质

上文已经提到过了要研究复变函数的性质只需要研究复变函数的虚部函数u(x,y)和虚部函数v(x,y)的性质。
复变函数在某一点可导的充要条件:
虚部函数v(x,y)与实部函数u(x,y)在某一点可微且满足Cauchy-Rieman方程:\frac {\partial u}{\partial x}=\frac{\partial v}{\partial y} and \frac{\partial u}{\partial y}=-\frac {\partial v}{\partial x}
复变函数解析的概念:
如果复变函数在一点可导且在这点的一个领域内处处可导,则称复变函数在这一点解析(注意复变函数在一点可导未必解析即可导是解析的必要不充分条件),如果复变函数在区域D内处处可导则称复变函数在区域D内解析。
复变函数在区域D解析的充要条件:
如果复变函数的实部函数与虚部函数在D内处处可微且满足柯西—黎曼方程那么称复变函数f(z)=u(x,y)+iv(x,y)为区域D内的解析函数。
Laplace方程的导出:
因为解析函数满足Cauchy-Rieman条件,由于解析函数可以求N阶导数,把柯西-黎曼方程:
\begin{cases} \frac {\partial u}{\partial x} =\frac {\partial v}{\partial y}\\ \frac {\partial u}{\partial y}=-\frac{\partial v}{\partial x} \end{cases}
两式分别关于x,y求导后可以导出Laplace:\frac {\partial^2 u}{\partial x^2}+\frac {\partial^2 u}{\partial y^2}=0同理\frac {\partial^2 v}{\partial x^2}+\frac {\partial^2 v}{\partial y^2}=0,我们把满足Laplace方程的函数称为调和函数,显然解析函数的实部函数和虚部函数是调和函数,于是数学家们把满足C-R条件的一对调和函数称为共轭调和函数,事实上一对共轭调和函数也可以构成一个解析函数,于是我们又得到复变函数在区域D内解析的一个充要条件即:复变函数的实部函数和虚部函数是一对共轭调和函数。于是我们只要知道一个调和函数我们就可以根据这个调和函数由柯西-黎曼条件我们就可以求出与已知的调和函数共轭的另一个调和函数,并且可以把这一对共轭调和函数组合成一个解析函数。
复初等函数:
因为实变函数与复变函数的主要差别就在与复变函数的变量为复数事变函数的为实数,总所周知在实变函数中许多的函数都是由初等函数复合而成,由此我们不难想象许多的复变函数也是由复初等函数复合而成的,因此认识清楚复变函数的初等函数也是由必要的,下面我只列举出复初等函数的形式,因为比较简单,在这里不做过多的叙述。
\begin{cases} 指数复函数:e^z\\ 对数复函数:z=e^{\omega},\omega=lnz\\ 幂复函数:\omega=z^\alpha,\omega=e^{\alpha lnz}\\ 三角复函数:e^{iy}=cosy+isiny,cosy=\frac 12(e^{iy}+e^{-iy}),siny=\frac {1}{2i}(e^{iy}+e^{-iy}) \end{cases}
上面研究的复变函数的可导性与解析性让我们对复变函数有了一个初步的了解,本质上就是把实变函数的性质在复变函数上推广,数学家们研究复变函数也是为了解决实际问题的,许多在实变函数上面问题用现有的实变函数上面的性质解决的话要么就是过于复杂计算困难,要么就是用实变函数根本无法啊解决,因此数学家门希望从复变函数上面找到突破口,成为解决实际问题的一种新的方法,实际上实变函数于复变函数有着很自然的联系,这正好从哲学上应验了万事万物或多或少都存在着联系
上面的知识都只是把实变函数的某些性质实在复变函数上面做简单的推广只是简单的从函数的角度来刻画复变函数,下面将展示如何从数的角度来刻画复变函数

复变函数的积分性质

积分公式的推导:
容易f(z)=u(x,y)+i(x,y)是C上的连续函数,且复积分\int_c f(z)dz存在,则:\int_cf(z)dz=\int_cu(x,y)dx-v(x,y)dy+i\int_cv(x,y)dx+u(x,y)dy
设曲线C的参数方程是x=x(t),y=y(t),(a\leqslant t \leqslant b)
将参数方程带入导出可得\int_Cf(z)dz=\int^b_af[z(t)]z^\prime(t)dt
单连通区域的柯西积分公式:
如果函数f(z)在单连通区域D解析,则f(z)在D内沿任一简单曲线C的积分:\oint_Cf(z)dz=0
由格林公式:\oint_Cf(z)dz=\oint_Cudx-vdy+i\oint vdx+udy=-\underset{D}\iint(\frac {\partial v}{\partial x}+\frac {\partial u}{\partial y})dxdy+i\underset{D}\iint(\frac {\partial u}{\partial x}-\frac {\partial v}{\partial y})
又由C-R条件则得:\oint_Cf(z)dz=0
柯西积分公式是在单连通区域内才满足,如果是多连通区域又怎么计算简单曲线的积分?如果不是单连通区域我们只需要做辅助线把多连通区域变成单连通区域即可,由此就得到在多连通区域上的复合闭路定理:\oint_Cf(z)dz=\sum\limits_{k=1}^{n}{\oint_{C_k}f(z)dz},特别的如果D是由内外两条闭路C,C_1围成的环形区域,而f(z)在D内及其边界上是解析的,则有\oint_Cf(z)dz=\oint_{C_1}f(z)dz
柯西积分公式:
设函数f(z)在简单闭曲线C上及其D内部是解析的,z_0是D内的任意一点则:f(z_0)=\frac {1}{2\pi i}\oint_C \frac{f(z)}{z-z_0}dz (1)
这个公式是由复合闭路定理推导而来,其思想是以z_0为中心做一个小圆周K其半径为r,由复合闭路定理有:\oint_C \frac{f(z)}{z-z_0}dz=\oint_K \frac{f(z)}{z-z_0}dz,当小圆周K的半径r \rightarrow 0时候\oint_K \frac{f(z)}{z-z_0}dz \rightarrow \oint_K \frac{f(z_0)}{z-z_0}dz\oint_K \frac{1}{z-z_0}dz=2\pi i 即求得一式,当然这只是思想,证明还需要严格的推理,这里不做叙述。
高阶导数公式
f^{(n)}(z_0)=\frac {n!}{2\pi i}\oint_C \frac{f(z)}{(z-z_0)^{n+1}},证略
复数幂级数:
解析函数的泰勒展开定理:
函数f(z)在区域D内解析,z_0是区域内一点R为z_0到D边界的最短距离则当|z-z_0|\lt R时候有
f(z)=\sum\limits_{n=0}^{\infty}a_n(z-z_0)^n,a_n=\frac{1}{2\pi i}\oint_C \frac{f(z)}{(z-z_0)^{n+1}}dz=\frac {f^n(z_0)}{n!}
复幂级数 \sum\limits_{n=0}^{\infty}a_nz^n如果收敛,一定在某一个圆域内收敛于一个解析函数,解析函数在解析域内也能够展开成一个复幂级数,由此可见复幂级数于解析函数之间有着天然的联系
洛朗级数:
解析函数可以展开成幂级数这是一个不争的事实,但是在现实生活中的问题往往是函数在某个点不解析但是在这个点附近的圆盘区域内解析,此时f(z)不能用含有z-z_0的幂级数展开。
我们称形如 \sum\limits_{n=-\infty}^{\infty}a_n(z-z_0)^n的级数称之为洛朗级数。数学家们证明了洛朗级数的收敛区域为圆环而且f(z)在圆环区域D:R_1 \lt |z-z_0| \lt R_2内解析则f(z)=\sum\limits_{n=-\infty}^{\infty}a_n(z-z_0)^n,其中a_n=\frac{1}{2\pi i}\oint_C \frac{f(z)}{(z-z_0)^{n+1}}dz,(n=0,\pm 1,\pm 2,\cdots)这里的C为圆环区域内的任意的圆周:|z-z_0|=R,R_1 \lt R \lt R_2
留数:
解析函数可以在圆环域内展开为幂级数,可以在圆环域内展开为洛朗级数。圆环的一种退化形式是一点的去心领域,当函数在一点的去心领域内解析而在这点不解析的时候这一点就是复变函数的一个孤立奇点,所以洛朗级数就成为研究复变函数孤立奇点的一个有力工具,而解析函数在孤立奇点处的留数是解析函数论中的重要概念之一,且留数在计算上有着巧妙的运用,复变函数在闭曲线上的积分问题可以转化求其孤立奇点的留数问题。
孤立奇点的分类:

如果一个复变函数的在其孤立奇点处的洛朗展开式中不包含z-z_0的负幂项,那么就称这个奇点为孤立奇点,如果负幂项次数绝对值的最大值为m我们就称这个奇点为m级级点,如果有无穷多个负幂项那么就称这个奇点为本性奇点。
留数的定义及计算:
函数f(z)z_0的去心领域内解析但是在z_0不解析,则函数在z_0点可以洛朗展开,f(z)=\sum\limits_{n=-\infty}^{\infty}C_n(z-z_0)^n其中C_n=\frac{1}{2\pi i}\oint_C \frac{f(z)}{(z-z_0)^{n+1}}dz ,C_{-1}=\frac {1}{2\pi i} \oint_C f(z)dz,相当于知道了C_{-1}就等于知道了积分\oint_C f(z)dz的值由此就有了留数的概念,若z_0f(z)的孤立奇点记f(z)z_0处的留数为Res[f(z),z_0]=\frac {1}{2\pi i}\oint_Cf(z)dz=C_{-1}
留数定理:
C是一条正向的简单闭曲线,若f(z)在C上及其C的内部D除去有限个孤立奇点z_1,z_2,\cdots,z_n外处处解析,那么\oint_Cf(z)dz=2\pi i \sum\limits_{k=1}^{n}Res[f(z),z_k]定理证明由复合闭路定理。
极点处的留数:
如果z_0为函数f(z)的m级级点则:
Res[f(z),z_0]=\frac{1}{(m-1)!}\lim_{z=z_0}\frac {d^{m-1}}{dz^{m-1}}[(z-z_0)^mf(z)]
如果z_0为函数f(z)的一级级点则:
Res[f(z),z_0]=\lim_{z=z_0}(z-z_0)f(z)
上面的知识点分别从函数性质以及数的角度刻画了复变函数,如果想要更加形象更加直观的刻画复变函数毫无疑问必须从集合的角度出发,我们可以从变换和映射的角度去考虑复变函数

复变函数的几何特性

在研究许多实际问题中,往往会遇到区域的复杂性,给问题的研究带来困难,那么我们该怎么办嘞?在空间解析几何中我们学习过极坐标变换,椭球面变换,球面变换,这都可以简化问题的难度从而更好的解决问题,在复变函数中我们可以利用解析函数所构成的变换——共形映射来把复杂的问题简单化
共形映射的概念:
共形映射顾名思义,就是经过映射后能够保留之前图形的特性的映射,那么我们如何来描叙原像经过映射后像保留了原像的特性嘞?伟大的数学家们想出了用保持角度不变和伸缩率的不变性(即像与原像的比例)来刻画“共形”这一精妙绝伦的概念,所谓保角性就是原像中有交与一点A的两条曲线,这两条曲线经过映射后任然交于一点,且两条像曲线过交点的夹角与原像曲线中的夹角一致(注意是一致不是相等,一致包括方向大小一样)那么就称这种映射在A是保角的,所谓伸缩率的不变性就是原像经过映射之后像在某一点的长度与原像在某一点的长度的极限值为一个定值与原像无关。综合以上两种特性就有了共形映射的概念:
如果\omega=f(z)z_0的领域内是一一的,在z_0处具有保角性和伸缩率的不变性那么称映射在z_0是共形的如果映射在D的没一点都共形即原像与像的相对位置保持不变,图形的比列保持不变则称\omega=f(z)在区域D内是共形的。
解析函数与共形映射:
如果函数\omega=f(z)z_0解析且f^\prime (z_0)\neq0那么映射\omega=f(z)z_0处是共形状的,如果解析函数的导数处处不为零那么映射在D内是共形映射。
几个初等函数所构成的映射:
幂函数:\omega=z^n 此映射在将角形域映射成一个角形域,且在z=0处不保角。
指数函数:\omega=e^z吃映射将X=常数直线映射成圆周,将Y=常数映射成射线,将水平带形域映射成角形域
研究复变函数的原本目的在于解决实变函数所解决不了的问题或者是简化实变函数的问题,上面的知识都只是把复变函数的基础体系建立起来,而Fourier变换和Laplace变换才是把实变函数与复变函数联系起来的接口

傅里叶变换与拉普拉斯变换

傅里叶变换

我们曾经在数学分析和高等代数中学过Fourier级数,一个以L为周期的函数f_L(t),如果在区间[-\frac{L}{2},\frac{L}{2}]上面连续那么在[-\frac{L}{2},\frac{L}{2}]上可以展开成傅里叶级数:

f_L(t)=\frac{a_0}{2}+\sum\limits_{n=1}^{\infty}(a_ncosn\omega t+b_nsinn\omega t)

\omega=\frac{L}{2\pi},a_0=\frac L2\int_{-\frac L2}^{\frac L2}f_L(t)dt, a_n=\frac L2\int_{-\frac L2}^{\frac L2}f_L(t)cosn\omega tdt,n=1,2,3,\cdots

b_n=\frac L2\int_{-\frac L2}^{\frac L2}f_L(t)sinn\omega tdt,n=1,2,3,\cdots

由欧拉公式:cost=\frac{e^{it}+e^{-it}}{2},sint=\frac{e^{it}-e^{-it}}{2i} 将两个式子带入f_L(t)可以导出:

f_L(t)=\sum\limits_{n=-\infty}^{\infty}c_ne^{in\omega t},c_n=\frac1L\int^{\frac2L}_{-\frac2L}f_L(t)e^{-in\omega t}dt,n=0,\pm1,\pm2,\cdots

上面研究的是周期函数事实上对于任何一个非周期函数f(t) 都可以看成是一个由某个周期函数L的函数f_L(t)当周期L\rightarrow + \infty时候转化而来的,当L\rightarrow +\infty时候可以导出一个等式:

f(t)=\frac{1}{2\pi}\int_{-\infty}^{+\infty}\mathcal{F}[f(t)](\omega)e^{i\omega t}d\omega,\mathcal{F}[f(t)](\omega)=\int_{-\infty}^{+\infty}f(t)e^{-i\omega t}dt

\mathcal F[f(t)](\omega)称为函数f(t)的傅里叶变换,f(t)=\frac{1}{2\pi}\int_{-\infty}^{+\infty}\mathcal{F}[f(t)](\omega)e^{i\omega t}d\omega 称为傅里叶逆变换

傅里叶变换的的四个性质:线性性质,位移性质,微分性质,积分性质。

卷积的傅里叶变换:

卷积的定义:设f,g定义在(-\infty,\infty)上,若任意的x\in(-\infty,+\infty),积分

\int^{+\infty}_{-\infty}f(y)g(x-y)dy收敛。则称改积分为函数f,g的卷积记为f*g ,卷积满足交换律与乘法对加法的分配律

卷积傅里叶定理:

\mathcal F[f*g](\omega)=\mathcal F[f](\omega)\bullet\mathcal F[g](\omega),\mathcal F^{-1}[\mathcal F[f](\omega)\bullet\mathcal F[g](\omega)](t)=f*g(t)

拉普拉斯变换:

如果想要求一个函数的傅里叶变换还需要函数在区间上满足狄氏条件,以及函数在无限区间上绝对可积,但是事实上这个条件很强很多的函数都不满足,因此傅里叶变换的应用范围受到的限制很大,为了然积分变换更具一般性,数学家们就引入了一个新的概念Laplace变换,实际上Laplace变换是一种特殊的傅里叶变换,他只不过是把一个不满足傅里叶条件的函数通过乘上一个衰减因子(使得函数收敛)以及乘上一个跃阶函数(控制收敛区间)来使得函数满足傅里叶变换的条件,然后把经过改造后的函数在实行傅里叶变换于是就得到了Laplace变换的概念:

函数f(t)t\geqslant 0时候由定义而且积分\int_{0}^{+\infty}f(t)e^{-st}dt在复数s的某一个区域内收敛则此积分所确定的函数记为F(s)=\mathcal L[f(t)](s)=\int^{+\infty}_{0}f(t)e^{-st}dt称为函数f(t)的拉普拉斯变换,f(t)=\mathcal L^{-1}[F(s)](t)为拉普拉斯逆变换

利用留数求拉普拉斯的逆变换:

其思想大致为将乘以跃阶函数与衰减函数之后的改造函数进行傅里叶变换逆变换之后就导出式子:

f(t)=\frac{1}{2\pi i}\int^{\sigma+i\infty}_{\sigma-i\infty}F(s)e^{st}ds\cdots(1)

f(t)的积分区域变成闭曲线积分区域然后将f(t)在新的闭曲线积分区域上积分,此时可以将积分分为两部分,一部分为原有的积分区间积分,另一部分为新增的使得积分区间变成闭曲线的积分曲线的积分,然后在证明第二部分当积分区间趋于无穷时候极限为零,再由函数f(t)在简单闭曲线上的积分等于其奇点的留数值由此我们可以导出式子:

f(t)=\frac{1}{2\pi i}\int^{\sigma+i\infty}_{\sigma-i\infty}F(s)e^{st}ds=\sum\limits_{k=1}^{n}Res[F(s)e^{st},s_k],s_k为f(t)在闭曲线内部的全部奇点

Laplace变换的性质:

F(s)=\mathcal L[f(t)](s),G(s)=\mathcal L[g(t)](s)

线性性质:\mathcal L[\alpha f(t)+\beta f(t)](s)=\alpha F(s)+\beta G(s)

\mathcal L^{-1}[\alpha F(s)+\beta F(s)](t)=\alpha f(t)+\beta g(t)

延迟性质:对t_0 \gt0\mathcal L[f(t-t_0)](s)=e^{-st_0}F(s)

\mathcal L^{-1}[e^{-st_0}F(s)](t)=f(t-t_0)

位移性质:\mathcal L [e^{s_0t}f(t)](s)=F(s-s_0)

微分性质:\mathcal L[f^{\prime}(t)](s)=sF(s)-f(0^+)

F^{(n)}(s)=\mathcal L[(-t)^nf(t)](s)

\mathcal L[f^{(n)}](s)=s^nF(S)-s^{n-1}f(0^+)-s^{n-2}f^{\prime}(0^+)-\cdots-f^{n-1}(0^+)

积分性质:\mathcal L[\int^t_0f(t)dt](s)=\frac 1sF(s)且若积分\int^{\infty}_sF(s)ds收敛则\frac{f(t)}{t}的拉普拉斯变换存在且\mathcal L[\frac{f(t)}{t}](s)=\int^{\infty}_sF(s)ds

初值定理:若\lim\limits_{s\rightarrow+\infty}sF(s)存在则有f(0^+)=\lim\limits_{s\rightarrow+\infty}sF(s)

卷积定理和傅里叶变换一致。


利用Laplace求变系数线性常微分方程的解

用Laplace变换求变系数微分方程

ty^{\prime\prime}-(1+t)y^{\prime}+2y=t-1,y(0)=0,y^{\prime}(0)=0的解

解:对微分方程两边求拉式变换可得

\mathcal L[ty^{\prime\prime}]-\mathcal L[(1+t)y^{\prime}]+\mathcal L[2y]=\mathcal L[t-1]\mathcal L[y(t)]=Y(s) ,并且利用微分性质可得

-[s^2Y(s)-sy(0)-y^{\prime}(0)]^{\prime}-sY(s)+y(0)+[sY(s)-y(0)]^{\prime}+2Y(s)=\dfrac{1}{s^2}-\dfrac 1s

将条件带入可得:

-[s^2Y(s)]^{\prime}-sY(s)+[sY(s)]^{\prime}+2Y(s)=\dfrac {1}{s^2}-\dfrac 1s

方程整理可得:

Y^{\prime}(s)(s-s^2)+3(1-s)Y(s)=\dfrac {1-s}{s^2}

可解得:Y^{\prime}(s)+\dfrac 3sY(s)=\dfrac {1}{s^3},利用一阶线性非齐次微分方程公式,得

Y(s)=e^{-\int\frac{3}{s}ds}[\int\frac {1}{s^3}e^{\frac 3s}ds+c]=\frac{1}{s^2}+\frac{c}{s^3}

显然Y(s)只有一个三级极点s=0

对方程求拉普拉斯逆变换可得:

y(t)=\sum\limits_{k=1}^{n}Res[Y(s)e^{st},s_k]=Res[Y(s)e^{st},s=0]=\dfrac {1}{2!}\lim\limits_{s\rightarrow0}\dfrac{d^2}{ds^2}[\dfrac {1}{s^2}+\dfrac {c}{s^3}]e^{st}(s-0)^3=t+\dfrac c2t^2

求微分方程组:

\begin{cases} x^{\prime} +y+z^{\prime}=1\\x+y^{\prime}+z=0\\y+4z^{\prime}=0 \end{cases}

满足x(0)=0,y(0)=0,z(0)=0的解

解: 对方程组两边分别去拉普拉斯变换令\mathcal L[x(t)]=X(s),\mathcal L[y(t)]=Y(s),\mathcal L[z(t)]=Z(s)并且考虑满足的初始条件,由微分性质可得像函数满足的方程组为:

\begin{cases} sX(s)+Y(s)+sZ(s)=\dfrac 1s\\X(s)+sY(s)+Z(s)=0\\Y(s)+4sZ(s)=0\end{cases}

由此解方程得:

X(s)=\dfrac{4s^2-1}{4s^2(s^2-1)}

Y(s)=\dfrac {-1}{s(s^2-1)}

Z(s)=\dfrac {1}{4s^2(s^2-1)}

对每一像函数取拉普拉斯逆变换可得:

x(t)=\mathcal L^{-1}[X(s)]=\dfrac 14\mathcal L^{-1}[\dfrac {3}{s^2-1}+\dfrac{1}{s^2}]=\dfrac 14(3sinht+t)

y(t)=\mathcal L^{-1}[Y(s)]=\mathcal L^{-1}[\dfrac 1s-\dfrac {s}{s^2-1}]=1-cosht

z(t)=\mathcal L^{-1}[Z(s)]=\dfrac 14\mathcal L^{-1}[\dfrac {1}{s^2-1}-\dfrac {1}{s^2}]=\dfrac 14(sinht-t)

注:双曲正余弦函数

sinht=\dfrac {e^t-e^{-t}}{2}

cosht=\dfrac {e^t+e^{-t}}{2}


文章下载地址:https://github.com/luozekun1230/Math-Inspiration.git

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 205,386评论 6 479
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,939评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,851评论 0 341
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,953评论 1 278
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,971评论 5 369
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,784评论 1 283
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,126评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,765评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 43,148评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,744评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,858评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,479评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,080评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,053评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,278评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,245评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,590评论 2 343

推荐阅读更多精彩内容