复变函数在计算机领域的应用,第一节计算机仿真在复变函数中的应用.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(-1)

[解]：对的单极点1求留数

>> clear

syms z

f=z/((z-1)*(z-2)^2);

>> R=limit(f*(z-1),z,1)

(仿真结果)R =1

[解]：对的二阶极点2求留数

>>clear

syms z

f=z/((z-1)*(z-2)^2);

>> R=limit(diff(f*(z-2)^2,z,1)/prod(1:1),z,2)

(仿真结果)R =-1

[解]：是有理分式，可用专门的函数residue()求各极点的留数

>>clear

[R,P,K]= residue ([1,0],[1,0,0,0,-1])

(仿真结果)(各点留数)R =

0.2500

0.2500

-0.2500 + 0.0000i

-0.2500 - 0.0000i

(奇点取值)P =

-1.0000

1.0000

0.0000 + 1.0000i

0.0000 - 1.0000i

(无奇点项)K =[]

9用留数定理求闭和环路积分

[解]：积分一

>> clear

syms t z

z=2*cos(t)+i*2*sin(t);

f=1/(z+i)^10/(z-1)/(z-3);

>> inc=int(f*diff(z),t,0,2*pi)

(仿真结果)inc =779/78125000*i*pi+237/312500000*pi

>> vpa(inc,3)

(输出三位有效数字)ans =.238e-5+.313e-4*i

[解]：积分二

>>clear

[R,P,K]= residue ([1,0],[1,0,0,0,-1]);

>> I=2*pi*i*sum(R)

(仿真结果)I =0

10 求下列函数的傅立叶变换和拉普拉斯变换

[解]：对

>>syms t v

>> fourier(1/t)

(仿真结果)ans =i*pi*(Heaviside(-w)-Heaviside(w))

[解]：对 (a>0)

>>syms t w

syms a positive

f=1/(t^2+a^2);

>> F=fourier(f,t,w)

(仿真结果)

F = 1/a*exp(a*w)*pi*Heaviside(-w)+1/a*exp(-a*w)*pi*Heaviside(w)

11 求下列函数的拉普拉斯变换

[解]：对 [解]：对

>>syms s t w >>syms s x;

>> F=ilaplace(1/(s-1)) >> F=ilaplace(1/(s^2+1))

(仿真结果) (仿真结果)

F =exp(t) F =sin(t)

展开阅读全文