Matlab在《复变函数与积分变换》中的应用.docVIP

Matlab在《复变函数与积分变换》中的应用 摘要：复变函数与积分变换工程是应用必备的基础课，在解决实际问题中也有十分重要的意义。MATLAB是集数值计算、图形处理、图像处理、符号计算文字处理、数字建模、实时控制、动态仿真、信号处理等功能为一体的数学应用软件【1】。将MATLAB应用到复变函数与积分变换的学习中，可以为复变函数与积分变换的计算和应用带来极大的方便。 关键字：复变函数 积分变换 MATLAB 应用 前言：把复变函数与积分变换的学习和MATLAB结合起来，就可以把复杂繁琐的计算交于计算机,而把主要精力集中在建立和优化数学模型上。可以利用MATLAB可以对一个复常数进行基本的求模，求幅角，求实部、虚部的运算。更进一步地，还可以求复数的指数、对数，对复数进行三角运算。在对于复变函数的研究中，可以求解复变函数的留数，并用来求复变函数的积分，对复变函数进行泰勒级数展开。在积分变换方面，可以对函数进行傅里叶变换、逆变换，进行拉普拉斯变换、逆变换。 复数的基本运算 1.1复数的实部和虚部 real (z) 返回复数z 的实部 imag (z) 返回复数z 的虚部 1.2共轭复数 conj (z) 返回复数z 的共轭复数 1.3复数的模和辐角 abs (z) 返回复数z 的模 angle (z) 返回复数z 的辐角 1.4 复数的乘除法 * 乘法：模相乘，辐角相加 / 除法：模相除，辐角相减 1.5 复数的平方根 sqrt (z) 返回复数z 的平方根值 1.7复数的幂运算 z^n 返回复数z 的n 次幂 1.8 复数的指数运算和对数运算 exp(z) 返回复数z的以e为底的指数值 log(z) 返回复数z的以e为底的对数值 1.9 复数的三角运算 sin (z)、cos (z)、tan (z)、cot (z)、sec (z)、asin (z)、…等函数，返回复数z 的函数值。 1.10复数方程求根 solve (‘f (x) = 0’) 求方程f (x) = 0 的根 例1 求下列复数的实部、虚部、共轭复数、模、辐角【2】。 (1) 1-i1+i7 (2) i1-i+1-i 解：编写程序 Z=[((1-i)/(1+i))^7,i/(1-i)+(1-i)/i ,i^5+3*i^18+i] re=real(Z) %求实部 im=imag(Z) %求虚部 Z1=conj(Z) %求共轭复数 r=abs(Z) %求模 theta=angle(Z) %求辐角 运行结果为： 例2 求z4 解：编写程序 solve (‘z^4+1 = 0’) 运行结果为： 2 计算复变函数的极限、微分与积分 解析函数是《复变函数与积分变换》中最主要的研究对象，它在理论和实际问题中有着广泛的作用。 2.1 MATLAB求复变函数极限 例3 解：编写程序 syms z; f=z/(sin(z)) limit(f,z,0) 即f(z)在z=0的极限是1 2.2 MATLAB求复变函数微分 例4 解：编写程序 syms z f=z/((1+z)*(sin(z))); diff(f) 2.3 MATLAB求复变函数积分 例5 解 syms t z z=2*cos(t)+i*2*sin(t); f=1/(z+i)^10/(z-2)/(z-5); inc=int(f*diff(z),t,0,2*pi) 3 用MATLAB计算留数 根据留数的定义，设Z0是解析函数f(z)的孤立奇点，我们把f(z)在Z0处的洛朗展开式中负一次幂项的系数C-1称为f(z)在Z0出的留数，记作Res[f(z),z]，即Res[f(z),z]=C-1。【2】 在matlab中利用residue（）函数可以求出任意复变函数的孤立奇点以及在该点处的留数值。格式如下：[r,p,k]=residue[b,a] 【3】。 例6求下列函数在奇点处的留数 f(z)=z+1 解：[r,p,k]=residue([1,1],[1,-2,0]) 所以 Res[f(z),2]= 1.5 Res[f(z),0]=-0.5 。 例7 计算积分 Czdz 解：先求被积函数