3、MATLAB求复变函数积分(1)用函数int求解非闭合路径的积分【Matlab源程序】syms z a bf=返回积分结果int(f,z,a,b)例 4 求积分 xl=[°ch3zdz; x2=["(z-1)e"dzA解 syms zx1=int(cosh(3*z),z,pi/6*i,0)x2=int(z-1)*exp(-z),z,0,i)结果为:x1 = -1/3*ix2 = -i/exp(i)
3、 MATLAB求复变函数积分 (1)用函数int求解非闭合路径的积分. 【Matlab源程序】 syms z a b f= int(f,z,a,b) 返回积分结果 例 4 求积分 π 6 0 i i 0 x1= ch3zdz; x2 ( 1) dz z e z − = − 解 syms z x1=int(cosh(3*z),z,pi/6*i,0) x2=int((z-1)*exp(-z),z,0,i) 结果为: x1 = -1/3*i x2 = -i/exp(i)
(2)用函数int求解闭合路径的积分117例5计算积分z的值.9 =-2 (2z +i)°(z - 1)(z - 3)解【Matlab源程序】symstzz=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若只输出6位有效数值,使用语句vpa(inc,6)结果为ans=.238258e-5+.313254e-4*i
(2) 用函数int 求解闭合路径的积分. 例5 计算积分 10 | | 2 1 d ( i) ( 1)( 3) z z z z z = + − − 的值. 解 【Matlab源程序】 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 若只输出6位有效数值,使用语句 vpa(inc,6) 结果为 ans =.238258e-5+.313254e-4*i