MATLAB在复变函数与积分变换的应用

1、用Matlab解决复变函数与积分变换中的问题摘要：复变函数与积分变换理论性较强,又是解决实际问题的强有力的工具.该课程已深入到数学的各个分支,如微分方程、积分方程、概率论和数论等多个学科.然而该课程的很多内容比较抽象,学起来比较枯燥且难学.本文利用MATLAB讨论了复变函数与积分变换中的复数运算、泰勒级数的展开、留数、有理函数展开、Fourier变换、Laplace变换和复变函数图形绘制等几个问题.这样不仅提高和完善复变函数与积分变换方法的实用性,同时可以培养学习者运用MATLAB语言编程的能力,对学习者以后的专业课及工作中使用数学软件进行数据处理有很大帮助.关键词：MATLAB; 复变函数;

2、 积分变换目 录1 引言(12 复常数的运算(12.1 求复数的实部、虚部、模、幅角、共轭复数(12.2 对于两个复常数之间进行乘法、除法运算及复方程求根(23 泰勒级数的展开(34 留数计算及积分计算和有理函数的部分分式展开(44.1 留数计算及积分计算(44.2 有理函数的部分分式展开(55 Fourier变换及其逆变换(66 Laplace换变换及其逆变换(87 复变函数图形绘制(9参考文献（10）1 引言复变函数与积分变换是电力工程、控制领域和通讯等理工科必备的重要课程,同时在解决实际问题中也有十分重要的作用.但是大多数人在学习这门课程时都会感觉内容抽象,学起来感觉枯燥且难学.如何应用

3、现代高科技信息技术,让比较难理解的理论与繁杂枯燥的内容变得生动有趣,激发学习的兴趣,以及可以提高计算能力、实践能力就相当重要.在国际学术界,MATLAB已经被接受为一种准确、可靠的标准计算软件.用户可以直接在Command Window内输入执行命令，或者可以建立一个M文件，输入较大应用程序，编译完成后一起运行.现在常用的MATLAB语言是基于最为流行的C+语言基础之上的，因此语法与C+语言有很大的相识，而且较C+语言更加简单，更符合研究人员对数学表达式的书写格式.使之更便利与非专业人员的使用.并且这种语言可拓展性极强，具有良好的可移植性，这也是在各个领域流行MATLAB的重要原因.本文把复变

4、函数与积分变换的学习过程和MATLAB结合起来，把复杂的计算交于计算机，目的是为了提高学生学习的兴趣与爱好同时也可以减轻学习的负担，缩短学习时间，大大提高了教学效果与质量.2 复常数的运算2.1 求复数的实部、虚部、模、幅角、共轭复数在MATLAB中的求解格式为：real(x %回车x的实部imag(x %回车x的虚部abs(x %回车x的模angle(x %回车x的幅角conj(x %回车x的共轭复数例1 求下列复数的实部、虚部、模、幅角、共轭复数.(1 (2 (3 解：在编辑器中建立M文件001.m如下：format rat X=5/4+7i,3*exp(2i*pi/5,i7+i(3/7+

5、5re=real(Xim=imag(Xab=abs(Xan=angle(Xco=conj(X运行结果如下：Z = 5/4 + 7i 305/329 + 2565/899i 7765/1343 - 561/1490ire = 5/4 305/329 7765/1343im = 7 2565/899 -561/1490ab = 2055/289 3 4305/743an = 283/203 142/113 -82/1261co = 5/4-7i 305/329- 2565/899i 7765/1343+561/1490i2.2 对于两个复常数之间进行乘法、除法运算及复方程求根在MATLAB中，两个

6、复数之间的乘法、除法可以使用“*”、“/”来实现，求复方程的解使用solve(f(x=0来实现.例2 (1 a= b=+ 计算a*b.(2 +5=0求所有根.解:在命令窗口中输入如下：>> a=2/(1+5i;>> b=3/5i+3i/(2+4i;>> c=a*bc=-0.0692 - 0.2538i>> solve('x3+5=0'ans=-5(1/35(1/3*(3(1/2*i/2 + 1/2-5(1/3*(3(1/2*i/2 - 1/23 泰勒级数的展开定理1 （泰勒展开定理） 设在区域内解析，,为到D的边界上各点的最短距离

7、当时，为在处的泰勒级数.其中：= =0，1，2，用函数taylor来实现泰勒级数的展开，taylortool可以进行泰勒级数逼近分析.例3 求函数在x=0的泰勒展开式的6次幂多项式和16次幂多项式，并分别进行泰勒级数逼近分析.解：在命令窗口中输入：>> clear>> syms x>> f=exp(-x;>> T1=taylor(f,7T1 =x6/720 - x5/120 + x4/24 - x3/6 + x2/2 - x + 1>> T2=taylor(f,17T2 =然后运用taylortool命令进行泰勒级数逼近分析，图（1）

8、为6次幂多项式泰勒级数逼近分析，图（2）为16次幂多项式泰勒级数逼近分析.图（1）图（2）由上图可知，泰勒级数展开的项数越多，函数值越接近原函数，4 留数计算及积分计算和有理函数的部分分式展开4.1 留数计算及积分计算定义1 设f(x在0<- 解析，即 a 是 f(x 的孤立奇点，则称积分值 为 f(x 关于点 a 的留数，记作 Resf(x,a. 其中 c 为在 0< - 内包含在点 a 的任何一条正向简单闭曲线 . 定理2 设函数f(x在区域D内除有限个孤立奇点z,z,z.,z处处解析，c是D内包围所有奇点的一条正向简单闭曲线，那么=f(x,z在MATLAB中可用如下方法：假设

9、以知奇点a和m重数，则用下面的MATLAB语句可求出相应的留数Limit(f*(x-a,x,a %返回x=a的一级极点的留数Limit(diff(f*(x-am,x,m-1/prod(1:m-1,z,a %返回x=a的m级极点的留数例4 求积分dx,C为正向圆周：=4.解：可知x=0为的一级极点,x=1为的二级极点在编辑器中建立M文件002.m如下：syms xf=exp(x/x*(x-12a1=limit(f*x,x,0a2=limit(diff(f*(x-12,x,1/prod(1:1,x,1结果如下：a1 =1 a2 =0所以dx=2iResf(x,0+Resf(x,1=4.2 有理函数

10、的部分分式展开在MATLAB的函数语言中有数值函数residue(来求取有理函数的部分分式展开表示，该函数的调用格式为：R,P,K=residue(B,A: %回车R为留数%回车P为极点 %回车K为两个多项式的比值B(s/A(s的部分分式展开的直接项 向量B是分子的系数、向量A是分母系数，向量R是返回的留数，向量P是返回的极点，向量K由的商的多项式系数组成，若，则为空向量，否则，.如果存在级极点则有，展开式包括以下形式： 例5 试求函数的部分分式展开解：在命令窗口中输入如下:>> B=1,0;>> A=1,0,-3,-2;>> R,P,K=residue(B

11、,A;>> B,A1=rat(R;>> B,A1,Pans =2.0000 9.0000 2.0000-2.0000 9.0000 -1.00001.0000 3.0000 -1.0000结果为: 5 Fourier变换及其逆变换我们平常所用到的积分变换,就是把函数乘上一个确定的二元函数,然后计算积分,即这样变成另一个函数.定理3 若f(t在上满足：（1）在任何的有限区间上满足Dirichlet条件；（2）在无限区间上绝对可积（即收敛）：则有=定义2 如果函数f(t满足定理3，由 (1设 (2则 (3(2式称为的傅里叶变换，记为称为的象函数，并且这样的积分运算称为取的F

12、ourier变换，式(3称作的傅里叶逆变换式，记为使用fourier函数来实现Fourier变换,格式为fourier(f,逆变换可用ifourier(F来实现.例6 求钟形脉冲函数的频谱函数,然后绘制频谱图.解：在命令窗口中输入如下：>> syms t w>> f=4*exp(-2*t2;>> F=fourier(fF =(2*2(1/2*pi(1/2/exp(w2/8>> ezplot(F,-6,6所得频谱图（2）所示图（2）6 Laplace变换及其逆变换定义3 如果函数当时有定义，并且广义积分 (4在s的某一区域内收敛，则由（4）式所确定

13、的参数为s的函数叫做函数的Laplace变换.使用laplace函数来实现Laplace变换，使用ilaplace函数来实现拉普拉斯逆变换.例7 证明Laplace变换的时移性质.其中，f为任意的一个函数，u是阶跃函数，L表示的是Laplace变换.解：在命令窗口中输入:>> syms t s>> syms t0 positive>> ft=heaviside(t-t0*sym('f(t-t0'>> FS=laplace(ft,t,s;>> FS_t=ilaplace(FS,s,t;>> ftft =f(t

14、 - t0*heaviside(t - t0>> FSFS =laplace(f(t, t, s/exp(s*t0>> FS_tFS_t =f(t - t0*heaviside(t - t0从运算过程中可以看出，ft=,FS为ft函数对应的Laplace变换的结果.最后,FS_t的结果为FS函数的Laplace逆变换,结果为ft=.7 复变函数图形绘制设有所学知识可知绘制复变函数的图形，需要四维空间才能满足.为了避免这一困难，借用两张复平面：z平面与w平面点集间的对应关系来来描述复变函数.例8 做圆周=5在映射下的象.解：在编辑器中建立M文件003.m如下：syms x

15、 y z tt=-pi:0.001:pix=5*cos(ty=5*sin(tz=x+i*yw=z+5./zsubplot(2,1,1plot(ztitle(z=5*cos(t+i*5*sin(taxis equalsubplot(2,1,2plot(wtitle(w=3*z+5./zaxis equal运行结果如图(3所示：图（3）参考文献：1葛美宝.利用MATLAB促进复变函数与积分变换的教学改革J.科技信息，2009，314（30）：36-40.2霍新霞,张世唯.复变函数与积分变换课程教学探究J.科教文汇(上旬刊），2012，223（11）：86-87.3曹海涛,张伟杰.工科“复变函数与积

16、分变换”教学改革J.中国电力教育，2013，260（01）：80-89.4温录亮.工科复变函数与积分变换课程的教学改革探究J.济南职业学院学报，2011，87（04）：65-67.5陈静,贠书杰.MATLAB软件在高等数学教学中的应用J.河南机电高等专科学校学报，2008,64（5）：64-666周德强.MATLAB在工程数学教学中的应用J.现代计算机（专业版），2007，263（07）：37-39美国科研出版社,2011:3.8孟品超,汤新昌,姜志侠.浅析复变函数与积分变换课程的建设J.教育教学论坛，2012，72（31）：281-282.9孙妍,刘向丽.复变函数与积分变换课程一体化改革之浅

17、见J.中国科教创新导刊，2011，613（29）：95.1012Application of MATLAB in Complex Variable Function and Integral TransformAbstract:The theory of complex variable function and integral transform is strong,and it is a powerful tool in solving practical problems.The course has been deep into all branches of mathematics

18、, such as differential equation, integral equation, probability theory and number theory.However,the course content is abstract,boring and difficult to learn.Here MATLAB is utilized to solve several problems in complex variable function and and integral transform,such as:the plural computing,the Taylor series expansion,residue,the expansion of rational functions,Fourier transform,Laplace transform and complex function graphing.This will not only improve and perfect the practical function o