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

Page * 解析函数的仿真 物理意义：静电场的复势 例1：已知平面电场的复电势是 ，并设电势对应复电势的实部，试作出它的电力线和等势线。 抛物线族 抛物线族 Page * 解析函数的仿真 解法2：计算机仿真 z=cplxgrid(20); f=i*sqrt(z); u=real(f);v=imag(f); x=real(z);y=imag(z); contour(x,y,u,‘-’); %画电势 hold on contour(x,y,v,:) hold off 绘制复电势的实部和虚部 Page * 解析函数的仿真 绘制电力线和等势线 [x,y]=meshgrid(-1:0.1:1); z=x+i*y; f=i*sqrt(z); u=real(f);[Ex,Ey]=gradient(-u); %求电势的梯度—电力线 x=real(z);y=imag(z); contour(x,y,u,-); hold on quiver(x,y,Ex,Ey,1.0) %画电力线的箭头 hold off Page * Page * 解析函数的仿真 例2：研究电偶极子所产生的电势和电场强度。设在(a,b)处有电荷+q，在(-a,-b)处有电荷-q，则在电荷所在平面上任何一点的电势为 其中： clear;clf; q=2e-6;k=9e9;a=1.5;b=-1.5;x=-6:0.6:6; [X,Y]=meshgrid(x); rp=sqrt((X-a).^2+(Y-b).^2); rm=sqrt((X+a).^2+(Y+b).^2); V=q*k*(1./rp-1./rm); %计算电势 [Ex,Ey]=gradient(-V); AE=sqrt(Ex.^2+Ey.^2); Ex=Ex./AE; Ey=Ey./AE; cv=linspace(min(V(:)),max(V(:)),49); contour(X,Y,V,cv,‘k-’) %绘制电势 axis square title(\fontname{宋体}\fontsize{18}电偶极子的场和等势线); hold on Page * 解析函数的仿真 quiver(X,Y,Ex,Ey,0.7) plot(a,b,wo,a,b,w+) plot(-a,-b,wo,-a,-b,w-) xlabel(x); ylabel(y); hold off Page * 复变函数的积分与留数定理 非闭合路径的积分计算 利用int积分函数求解 syms z; x1=int(cosh(3*z),z,pi/6*i,0) x1 = -1/3*i x2=int((z-1)*exp(-z),z,0,i) x2 = -i/exp(i) Page * 复变函数的积分与留数定理 闭合路径的积分计算 留数的定义：设z0为函数f(z)的孤立奇点， f(z)在z0处的留数为： 其中，C为去心领域0|z-z0|δ内任意一条绕z的正方向简单曲线。 （1）若z0为函数f(z)的一阶极点，则 （2）若z0为函数f(z)的m阶极点，则 Page * 复变函数的积分与留数定理 求孤立奇点的留数的MATLAB程序： R=limit(F*(z-z0),z,z0) %单奇点R=limit(diff(F*(z-z0)^m,z,m-1)/prod(1:m-1),z,z0) %m重奇点 syms z f=sin(z+pi/3)*exp(-2*z)/(z^3*(z-1)); R1=limit(diff(f*z^3,z,2)/prod(1:2),z,0) %三重奇点的仿真结果为：R1=1/2 - 3^(1/2)/4 R2=limit(f*(z-1),z,1) %单奇点的仿真结果为：R2= sin(pi/3+1)/exp(2) z=0是三重奇点 z=1是单奇点 例：求函数 的孤立奇点处的留数 Page * 复变函数的积分与留数定理 对于有理函数的留数计算，在MATLB中调用格式： [R, P, K]=residue(B, A) 向量B为f(z)的分子系数；向量A为f(z)的分母系数；向量R为留数；向量P为奇点位置；向量K为直接项，即为没有奇点的项，如： Page * 复变函数的积分与留数定理 [R, P, K]=residue([1 0 1], [1 1]) R = 2 P = -1 K = 1 -1 方法1：分子系数B=[1 0 1]，分母的系数A= [1 1] 故在奇点z=-1处的