数学物理方法仿真.docVIP

仿真实验 第一部分：计算机仿真在复变函数中的应用 仿真实验1： 绘出对数函数的图形. 【Matlab源程序】 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 view -75,30 仿真实验2：若对应为的根，其中且取整数.试用计算机仿真编程验证下列数学恒等式 成立. 【Matlab源程序】 n round 1000*random beta,1,1 +1 % n input please enter n su 1; sum 0; for s 1:n N s exp i*2*s*pi/n ; end for k 1:n for s 1:n if s~ k su 1/ N k -N s *su; end end sum sum+su; su 1; end sum %仿真验证结果为：n 735 sum 2.2335e-016 -5.1707e-016i 仿真实验3： 对于复变函数，用Mtlab仿真计算极限，. 【Matlab源程序】 clear syms z f sin z /z; limit f,z,0 %仿真结果为ans 1 limit f,z,1+i %仿真结果为 ans 1/2*sin 1 *cosh 1 -1/2*i*sin 1 *cosh 1 +1/2*i*cos 1 *sinh 1 +1/2*cos 1 *sinh 1 仿真实验5： 研究电偶极子 Diploe 所产生的电势和电场强度【2】.设在处有电荷,在处有电荷.则在电荷所在平面上任何一点的电势为,其中 【Matlab源程序】 根据解析函数理论中求复势的方法,可由等势线求出电力线方程. 下面给出计算机仿真方法求解：仿真 MATLAB 程序和仿真结果 clear ; clf;q 2e-6;k 9e9;a 1.5;b -1.5;x -6:0.6:6;y x; [X,Y] meshgrid x,y ; rp sqrt X-a .^2+ Y-b .^2 ; rm sqrt X+a .^2+ Y+b .^2 ; V q*k* 1./rp-1./rm ; % 计算电势 [Ex,Ey] gradient -V ; %计算电场强度 AE sqrt Ex.^2+Ey.^2 ;Ex Ex./AE;Ey Ey./AE; %场强归一化 cv linspace min min V ,max max V ,49 ; %用黑实线绘等势线 contour X,Y,V,cv,k- %axis square title \fontname 宋体 \fontsize 22 电偶极子的场和等势线 ,hold on quiver X,Y,Ex,Ey,0.7 plot a,b,wo,a,b,w+ plot -a,-b,wo,-a,-b,w- xlabel x ; ylabel y ,hold off 说明：图21.5中黑实线代表等势线，箭头构成电力线．根据题中电荷的位置，不难看出图中右下方为正电荷，左上方为负电荷． 仿真实验6： 求积分 . 【Matlab源程序】 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 仿真实验7：求函数在奇点处的留数． 【Matlab源程序】 [R,P,K] residue [1,0,1],[1,1] %仿真结果为： R 2 P -1 K 1 -1 仿真实验8：求函数在原点 0，0 ，以及（1，a）点处的Taylor展式． 【Matlab源程序】 syms x y; f x^2-2*x *exp -x^2-y^2-x*y ; F maple mtaylor,f,[x,y],8 %在 0,0 点处的泰勒级数展开 latex collect F,x syms a; F maple mtaylor,f,[x 1,y a],3 ; %在 1,a 点处的泰勒级数展开 F maple mtaylor,f,[x a],3 ; %对单变量展开, 且x a 仿真实验8：求的Fourier变换． 【Matlab源程序】 syms t v fourier 1/t %仿真结果为: ans i*pi* 1-2*heaviside w 仿真实验9：求的Fourier逆变换． 【Matlab源程序】 syms t u w ifourier w*exp -3