matlab处理数学物理方法,MATLAB在《数学物理方法》中的应用

《数学物理方法》是许多理工专业的必修课和重要基础课,也是一门公认的难度大的课程。因课程内容抽象,数学推导繁琐,学生学习起来感到非常枯燥。MATLAB是高性能的数值计算型数学类科技应用软件,具有优秀的数值计算功能和强大的数据可视化能力。在《数学物理方法》中应用MATLAB进行习题求解和计算机仿真,一方面可以提高解题的速度,另一方面可将抽象的解和一些特殊函数以图形形式显示出来,直观明了,物理意义明确。1.MATLAB在复变函数中的应用[1,2]1.1复数计算例1:求下列复数的实部、虚部、共轭复数、模与辐角。(1)(1+i)6(2)i1+i+1-ii定义复数后,利用简单的MATLAB语句:real)(、imag()、conj()、abs()、angle()可直接得出复数的的实部、虚部、共轭复数、模与辐角。解:>>a=[(1+i)^6i/(1+i)+(1-i)/i]a=0-8.0000i-0.5000-0.5000i>>real(a)ans=0-0.5000>>imag(a)ans=-8.0000-0.5000>>conj(a)ans=0+8.0000i-0.5000+0.5000i>>abs(a)ans=8.00000.7071>>angle(a)ans=-1.5708-2.3562例2:求方程z3+8=0的所有的根。调用内部函数Solve进行求解。解:>>Solve(‘z^3+8=0’)ans=-21+i*3^(1/2)1-i*3^(1/2)1.2Taylor级数展开Taylor级数的展开可以用符号工具箱中的Taylor()函数直接导出。例:将函数1(1-z)2在z=0处展成幂级数。解:>>symsz;>>Taylor(1/(1-z)^2,0)ans=1+2*z+3*z^2+4*z^3+5*z^4+6*z^51.3留数计算通过求极限的方法计算留数。例:确定函数f(z)=z+2iz5+4z3的极点,并求出函数在极点处的留数。解:通过对函数的分析,可知z=0是三阶极点,z=2i,-2i是单极点。>>symsz;>>f=(z+2*i)/(z^5+4*z^3);>>limit(diff(f*z^3,z,2)/prod(1:2),z,0)ans=-1/8*i>>limit(f*(z-2i),z,2i)ans=1/8*i>>limit(f*(z+2i),z,-2i)ans=0计算留数,还可直接调用MATLAB函数库中的Residue函数直接计算,见下题。1.4积分计算积分计算可以转化为留数计算,先求被积函数的留数,利用留数定理求复变函数的积分。例:求c矣(z2+1d)z(z2+4)(其中C是逆时针方向的圆周:z=23)解:>>[R,P,K]=residue([1],[1,0,5,0,4])R=%求被积函数的留数0+0.0833i0-0.0833i0-0.1667i0+0.1667iP=%四个极点0+2.0000i0-2.0000i0+1.0000i0-1.0000iK=%(展开式中无此解析项)][可见在圆周z=32内有四个极点,所以积分值S为:>>S=2*pi*i*sum(R)%计算出积分的值S=0利用MATLAB进行复数运算、级数展开、计算留数、计算积分,操作简单,运算速度快,比起手工计算省时、方便。2.MATLAB在积分变换中的应用[3,4]利用MATLAB语句:Fourier(f,x,w)、laplace(f,x,w),可将函数f(x)进行傅里叶变换和拉普拉斯变换。例1:写出函数f(x)=1x2+a2(a>0)的Fourier变换式。解:>>symsxw;