课堂授课专题3：复变函数的计算机仿真.pptVIP

数学物理建模与计算机辅助设计 专题3：复变函数的计算机仿真 本专题主要内容与参考资料 主要内容 复数和复数运算 复变函数的图形 复变函数的极限与导数 解析函数的图形 复变函数的积分 留数和级数的计算 复变函数的变换和逆变换 参考资料 杨华军, 数学物理方法, 电子工业出版社 薛定宇, 高等应用数学问题的Matlab求解, 清华大学出版社 复数和复数运算 数学上复数的表示方法 代数表示法 z = x + iy 几何表示法 z = (x, y) 三角表示法 z = r (cosθ+ sinθ) 指数表示法 z = r eiθ MATLAB中复数的表示方法 代数表示法 z=a+b*i 指数表示法 z= r*exp(i*theta) 复数和复数运算 创建复数矩阵 直接输入小矩阵，M文件输入大矩阵 用复数的实部矩阵和虚部矩阵构造 用复数的模矩阵和辐角矩阵构造 复数和复数运算 复数的实部和虚部的计算 real(x) 返回复数的实部 imag(x) 返回复数的虚部 复数和复数运算 共轭复数的计算 conj(x) 返回复数(矩阵)的共轭 复数的除法 * 和 / 注意：(…)/5i= (…)/(5*i) 复数的幂 ^ z = sym(a+i*b); collect(z^2) ans = a^2+2*i*a*b-b^2 复数的平方根 sqart(x) 返回x的平方根 z=1+2i; rst = sqrt(z) rst = 1.2720 + 0.7862i 复数和复数运算 复数的指数和对数运算 exp(x) 返回复数的以e为底的指数的值 log(x) 返回复数的以e为底的对数的值 复数的三角运算 计算方法为利用欧拉公式和复数的指数求法 形式上和实数的三角运算类似 复变函数图形的绘制步骤 第一步： 用CPLXGRID构建极坐标的复数数据网络 Z=CPLXGRID(m) 创建(m+1)×(2m+1)个网格点的复数极坐标系下的网格 复变函数图形的绘制步骤 第二步： 用CPLXMAP绘制对复变函数图形 CPLXMAP(z, f(z), (optional bound)) optional bound为选择变量的范围 复变函数图形的绘制步骤 CPLXROOT 画复数n次根的函数曲面 CPLXROOT(n, m) 以m×m网格画复数的n次方根 不指定m值则缺省为m=20 既不指定m值、也不指定n值，则n值缺省为3 常见复变函数的图形 绘制幂函数z, z3, z1/2 , z1/3的图形 常见复变函数的图形 绘制复变函数1/z的图像 常见复变函数的图形 绘制1/(1-z)的图像，并验证泰勒展开和罗朗展开函数 常见复变函数的图形 绘制指数函数ez的图形和其泰勒展开的图形 常见复变函数的图形 绘制对数函数lnz的图形 常见复变函数的图形 三角函数图形 sin cos 复变函数的极限与导数 复变函数的极限 极限存在条件：实部和虚部同时存在极限 求法： limit(F, x, a) 复变函数的极限与导数 复变函数的导数 diff(f, z, N) 解析函数的仿真 物理意义：静电场的复势 解析函数的仿真 复变函数的积分与留数定理 非闭合路径的积分计算 利用int积分函数求解 复变函数的积分与留数定理 闭合路径的积分计算 求孤立奇点的留数 R=limit(F*(z-z0),z,z0) %单奇点 R=limit(diff(F*(z-z0)^m,z,m-1)/prod(1:m-1),z,z0) %m重奇点 复变函数的积分与留数定理 闭合路径的积分计算 求有理函数的留数 [B, R, K]=residue(B, A) 复变函数的积分与留数定理 利用留数定理计算闭合路径积分 留数定理：设函数f(z)在区域D内除有限个孤立奇点z1, z2, …, zn外处处解析，C为D内包围所有奇点的一条正向简单闭曲线，则 复变函数的级数展开 复数列的收敛性： 用求极限的方法判断 幂级数的收敛半径：用比值法或根值法 复变函数的级数展开 单变量函数的泰勒级数展开 taylor(f,a,n) 傅里叶变换及其逆变换 常用函数 F=fourier(Fun,v,u); 将函数Fun做以v和u为共轭变量的变换 f=ifourier(Fun,u,v);将函数Fun做以v和u为共轭变量的反变换 本专题小结 复数和复数运算 复变函数的图形 复变函数的极限与导数 解析函数的图形 复变函数的积分 留数和级数的计算 复变函数的变换和逆变换 数学物理建模与计算机辅助设计 * z=1+2*i z = 1.0000 + 2.0000i z=1+2i z = 1.0000 + 2.0000i z=2*exp(i*pi/4) z = 1.4142