%16020130038?李璞初?使用MATLAB畫電力線與等勢面大作業?2018.4.7
clearnumber=input（‘請輸點電荷數:‘）;
charge=zeros（number3）;%創建元數組以存放電荷信息
for?i=1:number
????s1=[‘請輸入第‘num2str（i）‘個電荷的電荷量（μC）:‘];
????s2=[‘請輸入第‘num2str（i）‘個電荷的x坐標:‘];
????s3=[‘請輸入第‘num2str（i）‘個電荷的y坐標:‘];
????charge（i1）=input（s1）*1e-6;
????charge（i2）=input（s2）;
????charge（i3）=input（s3）;
end?????%得到所需的所有電荷的信息
k=9*10^9;%這里取介質為無限大真空，化簡得到
x=-2.2:0.01:2.2;
y=-2.2:0.01:2.2;
[XY]=meshgrid（xy）;?%建立X_Y坐標系
U=zeros（441441）;
clf
for?j=1:number
????U=U+k*charge（j1）./sqrt（（X-charge（j2））.^2+（Y-charge（j3））.^2）;
end%得到電勢的值?注意U與XY的對應關系是U[YX]
for?a=1:number
for?b=-10:10?%橫坐標變化
for?c=-10:10?%縱坐標變化
if?abs（U（charge（a3）*100+221+bcharge（a2）*100+221+c））>abs（U（charge（a3）*100+221+5charge（a2）*100+221+5））
U（charge（a3）*100+221+bcharge（a2）*100+221+c）=U（charge（a3）*100+221+5charge（a2）*100+221+5）;
end
end
end
end%由于在電荷處，電勢理論上無窮大，所以通過設置閾值對其進行“削頂”
contour（XYU100）%畫等勢面
hold?on
%畫電場線與畫等勢線完全沒有關系╮（╯_╰）╭
total=sum（abs（charge（:1）））;%電荷模值求和
px=[];py=[];%存放電場線數據
for?a=1:number%畫由第a個電荷發出的電場線
????for?b=1:round（abs（charge（a1））/total*50*number）%電場線數目由歸一化電荷決定
????????j=round（abs（charge（a1））/total*50*number）;i=1;
????????px（abi）=charge（a2）+0.05*2^（1/2）*cos（b*2*pi/j）;
????????py（abi）=charge（a3）+0.05*2^（1/2）*sin（b*2*pi/j）;%分配電場線起點坐標
????????t=1;
????????while?t==1
????????????Ex=0;Ey=0;
????????????for?c=1:number
????????????????Rx=px（abi）-charge（c2）;
????????????????Ry=py（abi）-charge（c3）;
????????????????E=k*charge（c1）./（Rx^2+Ry^2）;
??????