电磁场 matlab 仿真实验一
实验一：[例 7-5]试分析一对等量异号的电荷周围空间上的电位和电场分布情况。
分析：将等量异号的电荷的几何中心放置于坐标原点位置，则它们在空间某点 p 处产
生的点位为：
其中 G 为格林函数
将 G 用片面积坐标表示为
在编程时，将 G 当作点位函数处理，并利用梯度求出唱腔 E=-▽φ。用 matlab 的 m 语言
编写的程序如下：
[x,y]=meshgrid(-10:0.1:10);
[Q,R]=cart2pol(x,y);
R(R<=1)=NaN;
q=input('请输入电偶极子的电量 q＝') %原程序有误，以此为准
d=input('请输入电偶极子的间距 d＝') %原程序有误，以此为准
E0=8.85*1e-12;
K0=q/4/pi/E0;
g1=sqrt((d./2).^2-d.*R.*cos(Q)+R.^2); %原程序有误，以此为准
g2=sqrt((d./2).^2+d.*R.*cos(Q)+R.^2); %原程序有误，以此为准
G=log(K0*g2./g1);
contour(x,y,G,17,'g');
hold on
[ex,ey]=gradient(-G);
tt=0:pi/10:2*pi; %原程序未定义 tt，以此为准
sx=5*sin(tt);sy=5*cos(tt);
streamline(x,y,ex,ey,sx,sy);
xlabel('x');ylabel('y');
hold off;
当运行此程序后，按提示输入电偶极子电量和嗲耨集子间距如下：
请输入电偶极子的电量 q＝0.5*1e-10
请输入电偶极子的间距 d＝0.01
即可汇出入图说使得嗲耨集资周围的长的分布图。
资源评论
