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

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

第一节：计算机仿真在复变函数中的应用

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