二、实验步骤
1、打开 MATLAB 软件，新建命令文档并保存，并在文档中输入程序。
2、输入点电荷 q1 的坐标（q1x，q1y）, 以及 q1 所带的电量。调用 input 函数。
如果不知道该函数的使用方法可在 MATLAB 命令行处键入 doc input。
3、输入点电荷 q1 的坐标（q1x，q1y）, 以及 q1 所带的电量。
4、定义比例常系数
， 命令为 k=9e9。
5、定义研究的坐标系范围为
[ ] [ ]
5,5 , 5,5x yÎ - Î -
,步长值为 0.1。
6、将x,y两组向量转化为二维坐标的网点结构，函数为meshgrid。命令为
[X,Y]=meshgrid(x,y)，如果不知道该函数的使用方法可在MATLAB命令行处键入
doc meshgrid。
7、计算任意一点与点电荷之间的距离 r，公式为
2 2
1 1 1
( ) ( )r x q x y q y= - + -
，
2 2
2 2 2
( ) ( )r x q x y q y= - + -
8、计算由 q1,q2 两个点电荷共同产生的电势
9、注意，由于在 q1 和 q2 位置处计算电势函数为无穷大或者无穷小，因此要把
这两点去掉掉，以方便下面绘制等势线。具体命令可参考
Vinf1=find(V==inf);
V(Vinf1)=NaN;
Vinf2=find(V==-inf);
V(Vinf2)=NaN;
如果是可以解释这四句话的原理，可以有加分！
10、根据天长强度与电位函数的关系
，可直接计算 E，调用 gradient 函
数。如果不知道该函数的使用方法可在 MATLAB 命令行处键入 doc gradient。
参考命令为
[Ex,Ey]=gradient(-V)
11、计算 E 的模值
，注意在计算时运算要加点，Ex.^2
12、计算电场强度的单位矢量，
，
，注意在计算时运算要
加点，Ey=Ey./ Eq
13、生成你要绘制的等位线的数量与每条等位线上的电位值
cv=linspace(min(min(V)),max(max(V)),49)
该命令表示在最大电位与最小电位之间插入 49 个点，形成一个向量 cv
