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

本科毕业论文 题目： MATLAB在复变函数与积分变换的应用 学院： 数学与计算机科学学院 班级： 数学与应用数学2009级班 姓名： 指导教师： 职称： 副教授 完成日期： 2013 年 05 月 10 日MATLAB在复变函数与积分变换的应用摘要：复变函数与积分变换理论性较强,又是解决实际问题的强有力的工具.该课程已深入到数学的各个分支,如微分方程、积分方程、概率论和数论等多个学科.然而该课程的很多内容比较抽象,学起来比较枯燥且难学.本文利用MATLAB讨论了复变函数与积分变换中的复数运算、泰勒级数的展开、留数、有理函数展开、Fourier变换、Laplace变换和复变函数图形绘制等几个问题.这样不仅提高和完善复变函数与积分变换方法的实用性,同时可以培养学习者运用MATLAB语言编程的能力,对学习者以后的专业课及工作中使用数学软件进行数据处理有很大帮助.关键词：MATLAB; 复变函数; 积分变换目 录1 引言(1)2 复常数的运算(1) 2.1 求复数的实部、虚部、模、幅角、共轭复数(1) 2.2 对于两个复常数之间进行乘法、除法运算及复方程求根(2)3 泰勒级数的展开(3)4 留数计算及积分计算和有理函数的部分分式展开(4) 4.1 留数计算及积分计算(4) 4.2 有理函数的部分分式展开(5)5 Fourier变换及其逆变换(6)6 Laplace换变换及其逆变换(8)7 复变函数图形绘制(9)参考文献（10）1 引言复变函数与积分变换是电力工程、控制领域和通讯等理工科必备的重要课程,同时在解决实际问题中也有十分重要的作用.但是大多数人在学习这门课程时都会感觉内容抽象,学起来感觉枯燥且难学.如何应用现代高科技信息技术,让比较难理解的理论与繁杂枯燥的内容变得生动有趣,激发学习的兴趣,以及可以提高计算能力、实践能力就相当重要.在国际学术界,MATLAB已经被接受为一种准确、可靠的标准计算软件.用户可以直接在Command Window内输入执行命令，或者可以建立一个M文件，输入较大应用程序，编译完成后一起运行.现在常用的MATLAB语言是基于最为流行的C+语言基础之上的，因此语法与C+语言有很大的相识，而且较C+语言更加简单，更符合研究人员对数学表达式的书写格式.使之更便利与非专业人员的使用.并且这种语言可拓展性极强，具有良好的可移植性，这也是在各个领域流行MATLAB的重要原因.本文把复变函数与积分变换的学习过程和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)+5re=real(X)im=imag(X)ab=abs(X)an=angle(X)co=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中，两个复数之间的乘法、除法可以使用“*”、“/”来实现，求复方程的解使用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/3)5(1/3)*(3(1/2)*i)/2 + 1/2)-5(1/3)*(3(1/2)*i)/2 - 1/2)3 泰勒级数的展开 定理1 （泰勒展开定理） 设在区域内解析，,为到D的边界上各点的最短距离当时，为在处的泰勒级数.其中：= =0，1，2， 用函数taylor来实现泰勒级数的展开，taylortool可以进行泰勒级数逼近分析. 例3 求函数在x=0的泰勒展开式的6次幂多项式和16次幂多项式，并分别进行泰勒级数逼近分析. 解：在命令窗口中输入： clear syms x f=exp(-x); T1=taylor(f,7) T1 =x6/720 - x5/120 + x4/24 - x3/6 + x2/2 - x + 1 T2=taylor(f,17) T2 =x16/20922789888000 - x15/1307674368000 + x14/87178291200 - x13/6227020800 + x12/479001600 - x11/39916800 + x10/3628800 - x9/362880 + x8/40320 - x7/5040 + x6/720 - x5/120 + x4/24 - x3/6 + x2/2 - x + 1 然后运用taylortool命令进行泰勒级数逼近分析，图（1）为6次幂多项式泰勒级数逼近分析，图（2）为16次幂多项式泰勒级数逼近分析.图（1）图（2） 由上图可知，泰勒级数展开的项数越多，函数值越接近原函数，4 留数计算及积分计算和有理函数的部分分式展开4.1 留数计算及积分计算 定义1 设f(x)在0-R解析，即a是f(x)的孤立奇点，则称积分值为f(x)关于点a的留数，记作Resf(x),a.其中c为在0- B=1,0; A=1,0,-3,-2; R,P,K=residue(B,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) 式称为的傅里叶变换，记为 称为的象函数，并且这样的积分运算称为取的Fourier变换，式(3)称作的傅里叶逆变换式，记为 使用fourier函数来实现Fourier变换,格式为fourier(f),逆变换可用ifourier(F)来实现. 例6 求钟形脉冲函数的频谱函数,然后绘制频谱图. 解：在命令窗口中输入如下： syms t w f=4*exp(-2*t2); F=fourier(f) F =(2*2(1/2)*pi(1/2)/exp(w2/8) ezplot(F,-6,6)所得频谱图（2）所示图（2）6 Laplace变换及其逆变换 定义3 如果函数当时有定义，并且广义积分 (4)在s的某一区域内收敛，则由（4）式所确定的参数为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); ft ft =f(t - t0)*heaviside(t - t0) FS FS =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 y z tt=-pi:0.001:pix=5*cos(t)y=5*sin(t)z=x+i*yw=z+5./zsubplot(2,1,1)plot(z)title(z=5*cos(t)+i*5*sin(t)axis equalsubplot(2,1,2)plot(w)title(w=3*z+5./z)axis equal运行结果如图(3)所示：图（3）参考文献：1葛美宝.利用MATLAB促进复变函数与积分变换的教学改革J.科技信息，2009，314（30）：36-40.2霍新霞,张世唯.复变函数与积分变换课程教学探究J.科教文汇(上旬刊），2012，223（11）：86-87.3曹海涛,张伟杰.工科“复变函数与积分变换”教学改革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-397金琦,宋明福,董玉才.军队学历教育院校复变函数与积分变换教学方法与手段改革的探索C./中国教育技术协会实践教学委员会、上海高职电子信息类职业教学指导委员会.2011年全国电子信息技术与应用学术会议论文集.上海：美国科研出版社,2011:3.8孟品超,汤新昌,姜志侠.浅析复变函数与积分变换课程的建设J.教育教学论坛，2012，72（31）：281-282.9孙妍,刘向丽.复变函数与积分变换课程一体化改革之浅见J.中国科教创新导刊，2011，613（29）：95.10S.K.Sen,Gholam Ali Shaykhian.MatLab tutorial for scientific and engineering computationsJ. Nonlinear Analysis,2009,71(12)：46-4711V.A.Abilov,F.V.Abilova,M.K.Kerimov.Sharp estimates for the convergence rate of Fourier series of complex variable functions in L2(D, p(z)J. Computational Mathematics and Mathematical Physics,2010,50(6)：51-5212V.A.Abilov,F.V.Abilova,M.K.Kerimov.Sharp estimates for the convergence rate of Fourier series in terms of orthogonal polynomials in L2(a, b), p(x)J. Computational Mathematics and Mathematical Physics,2009,49(6):31-33.Application 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, 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