计算机仿真在复变函数中的应用.docVIP

第一节：计算机仿真在复变函数中的应用 1 求下列复数的实部，虚部、模，辐角。 [解]将要求的复数建立复数矩阵 a=[i 1+i*sqrt(3) (3+4i)^3 exp(1+i) (1-i)/(1+i)]； real(a) （仿真结果）ans = 0 1.0000 -117.0000 1.4687 0 imag(a) （仿真结果）ans =1.0000 1.7321 44.0000 2.2874 -1.0000 abs(a) （仿真结果）ans =1.0000 2.0000 125.0000 2.7183 1.0000 angle(a) （仿真结果）ans =1.5708 1.0472 2.7819 1.0000 -1.5708 2 求下列各式的值。 [解] sqrt(3+4i) （仿真结果）ans = 2.0000 + 1.0000i [解]： i^1/3 （仿真结果）ans = 0 + 0.3333i [解]： i^1/i （仿真结果）ans =1 3 求解下列方程的根 [解] solve(z^3=8) （仿真结果）ans =[ 2]，[ -1+i*3^(1/2)]，[ -1-i*3^(1/2)] solve(z^8=1) （仿真结果）ans =[ 1]，[-1]，[i]，[-i]， [ 1/2*2^(1/2)+1/2*i*2^(1/2)]，[ -1/2*2^(1/2)-1/2*i*2^(1/2)]，[ 1/2*2^(1/2)-1/2*i*2^(1/2)]，[ -1/2*2^(1/2)+1/2*i*2^(1/2)] 4 作出以下复变函数的图形 [解]： z=cplxgrid(30);cplxroot(3);colorbar(vert);title(z^[1/3]) （仿真结果） [解]： z=5*cplxgrid(30); cplxmap(z,sin(z));colorbar(vert);title(sin(z)) （仿真结果） [解]： z=cplxgrid(20);w=log(z); for k=0:3 w=w+i*2*pi; surf(real(z),imag(z),imag(w),real(w)); hold on title(Lnz) end （仿真结果） [解]： z=cplxgrid(20);w=exp(z);surf(real(z),imag(z),real(w),imag(w)); （仿真结果） 5 求下列幂级数的收敛半径 [解]：级数一 clear syms n C1=(-1)^n C1 =(-1)^n R1=abs(limit(C1^(-1/n),n,inf)) （仿真结果）R1 =1 [解]：级数二 clear syms k C1=(1/k)^k R1=abs(limit(C1^(-1/k),k,inf)) （仿真结果）R1 =∞ 6 求下列函数的泰勒级数（8阶）展开式及图形 [解]： clear syms z z0=-20:0.05:20;y=exp(z); p1=taylor(y,z,8),y1=subs(p1,z,z0); plot(z0,y1);title(exp(z)的8阶图形) （仿真结果）p1 =1+z+1/2*z^2+1/6*z^3+1/24*z^4+1/120*z^5+1/720*z^6+1/5040*z^7 （仿真图形） [解]： clear syms z z0=-20:0.05:20;y=sin(z); p1=taylor(y,z,8),y1=subs(p1,z,z0); plot(z0,y1);title(sin(z)的8阶图形) （仿真结果p1 =z-1/6*z^3+1/120*z^5） （仿真图形） 7 求解析函数的非闭和回路的积分 [解]：积分一 clear syms z x1=int(cosh(3*z),z,pi/6*i,0) （仿真结果）x1 =-1/3*i [解]：积分二 clear syms z x1=int(exp(z),z,pi/2*i,0) （仿真结果） x1 =1-i 8 求以下函数的孤立奇点处的留数 [解]：对的单极点-1求留数 clear syms z f=exp(z)/(1+z)； R=limit(f*(z+1),z,-1) （仿真结果）R =exp(-1) [解]：对的单极点1求留数 clear syms z f=z/((z-1)*(z-2)^2); R=limit(f*(z-1