上机练习专题数学物理方程的计算机求解和可视化.pptVIP

数学物理建模与计算机辅助设计 数学物理建模与计算机辅助设计 上机练习专题4：数学物理方程的计算机求解和可视化 Page * 本专题主要内容与参考资料 主要内容 直接对代数表达式求解并进行可视化 利用差分法求解并进行可视化 利用Matlab的PDE工具箱求解并进行可视化 参考资料 杨华军, 数学物理方法, 电子工业出版社 彭芳麟, 数学物理方程的Matlab解法与可视化, 清华大学出版社 Page * 对代数表达式进行可视化 三维拉普拉斯方程 环形电流磁感应强度 圆线圈的半径为a，通有电流I，求空间任意一点的磁感应强度B Page * 对代数表达式进行可视化 三维拉普拉斯方程程序和可视化%P104_1.m a=0.35;[X,Y]=meshgrid(0:0.1:1);[theta,r]=cart2pol(Y,X); rout=r; rin=r;rin(find(rina))=NaN;rout(find(routa))=NaN; Boutr=0; Binr=0; BoutQ=0; BinQ=0; for k=1:2:43 fun10=legendre(k,0);fun0=fun10(2);fun=legendre(k,cos(theta)); rfun=squeeze(fun(1,:,:));Qfun=squeeze(fun(2,:,:)); % Cr=1/a^k*rin.^(k-1);Bkinr=-1/2*fun0*Cr.*rfun;Binr=Binr+Bkinr; CQ=1/k/a^k*rin.^(k-1);BkinQ=-1/2*fun0*CQ.*Qfun;BinQ=BinQ+BkinQ; % Dr=a.^(k+1)./rout.^(k+2);Bkoutr=-1/2*fun0*Dr.*rfun;Boutr=Boutr+Bkoutr; DQ=1/(k+1)*a^(k+1)./rout.^(k+2);BkoutQ=1/2*fun0*DQ.*Qfun;BoutQ=BoutQ+BkoutQ; end Page * 对代数表达式进行可视化 L1=isnan(Binr); Binr(find(L1==1))=0;BinQ(find(L1==1))=0; L2=isnan(Boutr); Boutr(find(L2==1))=0;BoutQ(find(L2==1))=0; Br=Binr+Boutr;BQ=BinQ+BoutQ; By=Br.*cos(theta)-BQ.*sin(theta); Bx=Br.*sin(theta)+BQ.*cos(theta); vy=0; vx=[0:0.01:0.35]; [Vx,Vy]=meshgrid(vx,vy); cla streamline(X,Y,Bx,By,Vx,Vy); %streamline绘制二维、三维流线图 hold on streamline(-X,Y,-Bx,By,-Vx,Vy); streamline(-X,-Y,-Bx,-By,-Vx,-Vy); streamline(X,-Y,Bx,-By,Vx,-Vy); title(B 的场线) 对代数表达式进行可视化 Page * Page * 对代数表达式进行可视化 泊松方程和格林函数 球域的格林函数 在半径为a的导体球内(或导体球外)，距球心r0处放置电量为 的点电荷，求它形成的静电场 定解问题是 解的代数表达式是 Page * 对代数表达式进行可视化 泊松方程和格林函数 球域的格林函数的程序及可视化：(1)球外有点电荷 %P133.m a=1; r1=2; q=1; [X,Y]=meshgrid(-3:0.01:3,0:0.01:4); [Q,R]=cart2pol(Y,X); R(R1)=NaN; u1=1./sqrt(r1^2-2*r1*R.*cos(Q)+R.^2); ar=a/r1; u2=-ar./sqrt((a*ar)^2-2*a*ar.*R.*cos(Q)+R.^2); u=u1+u2; contour(X,Y,u,[-1:20,3000]); hold on [ex,ey]=gradient(-u); t=0:pi/10:2*pi; sx=0.1*cos(t); sy=2+0.1*sin(t); streamline(X,Y,ex,ey,sx,sy); axis equal tt=0:pi/30:2*pi; plot(exp(i*tt),r) 对代数表达式进行可视化 Page * Page * 泊松方程和格林函数 球域的格林函数的程序及可视化：(2)球内有点电荷 %P135_1.m a=2;