有關matlab畫電力線

b10007133

小的最近在用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

ssamxd

学习一下，谢谢了！