作者siang820420 (大德)
看板comm_and_RF
標題[請益] 有關MATLAB畫電力線
時間Mon Mar 25 00:35:13 2013
小的最近在用matlab畫電荷的電力線但是做出來為什麼r>1的地方沒有電力線呢?
能否請高手教導我程式要怎麼改
謝謝
我的程式如下:
clear %清除變量
E0=8.85e-12; %真空介電係數
C0=1/4/pi/E0;
q=1.6*10^(-19); %電荷電量
q=-2; %電量比
xm=2.5; %橫坐標範圍
ym=2; %縱坐標範圍
x=linspace(-xm,xm); %橫坐標向量
y=linspace(-ym,ym); %縱坐標向量
[X,Y]=meshgrid(x,y);
R1=sqrt((X+1).^2+Y.^2);
R2=sqrt((X-1).^2+Y.^2);
U=1./R1+q./R2;
u=1:0.5:4;
figure %創圖形
grid on %加網格
legend(num2str(u')) %圖例
hold on %保持圖像
plot([-xm;xm],[0;0]) %水平線
plot([0;0],[-ym;ym]) %垂直線
plot(-1,0,'o','MarkerSize',12) %畫第一電荷
plot(1,0,'o','MarkerSize',12) %畫第二電荷
[Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%梯度
dth1=20;
th1=(dth1:dth1:180-dth1)*pi/180;
r0=0.1; %電力線起點半徑
x1=r0*cos(th1)-1; %電力線起點橫坐標
y1=r0*sin(th1); %電力線起點縱坐標
streamline(X,Y,Ex,Ey,x1,y1) %畫左上電力線
streamline(X,-Y,Ex,-Ey,x1,-y1) %畫左下電力線
dth2=dth1/q;
th2=(180-dth2:-dth2:dth2)*pi/180;
x2=r0*cos(th2)+1;
y2=r0*sin(th2);
streamline(X,Y,Ex,Ey,x2,y2)
streamline(X,-Y,Ex,-Ey,x2,-y2)
axis equal tight %刻度相等
title('電力線分佈','fontsize',20)
xlabel('r','fontsize',16) %顯示橫坐標
ylabel('E(U)','fontsize',16) %顯示橫坐標
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 1.160.39.177