清空筛选条件
清空筛选条件

How to apply sum function for simulating charges around a circle?

1 次查看（过去 30 天）
Rodrigo Ayala
Rodrigo Ayala 2021-3-9
评论： Rodrigo Ayala 2021-3-9
I'm trying to model electric potential energy of protons arranged in a circle, graphically is(the next plot was made manually in geogebra 3d),
For :
Asumming k=1, q1=1, q2=2 , mathematically can been described as:
The thing is that i want to generalize this into a summation function in matlab for n charges separated in the circle, mathematically it will be:
I try the next code, but it doesn't work. Any suggestions? I can't understand how the sum function works exactly:
[x,y,t]=meshgrid(-15:0.5:15,-15:0.5:15, 0:pi/8:pi/2);
r=5;
z=sum(1./sqrt((x-r*cos(t)).^2+(y-r*sin(t).^2)), 3);
surf(x(:,:,1),y(:,:,1),z);

请先登录，再回答此问题。

回答（2 个）

Matt J
Matt J 2021-3-9
编辑：Matt J 2021-3-9
Ran in:
xs=linspace(0,10,500);
[x,y,t]=meshgrid(xs,xs, 0:pi/8:pi/2);
r=5;
z=sum(1./ ( (x-r*cos(t)).^2 + (y-r*sin(t)).^2 ), 3);
surf(x(:,:,1),y(:,:,1),z,'EdgeColor','none','FaceAlpha',0.3,'FaceColor','g');
  6 个评论
显示 4更早的评论
Matt J
Matt J 2021-3-9
编辑：Matt J 2021-3-9
Ran in:
Is the same r being used in the geogebra drawing? Also, you've presented different versions of your equations in different places, some with sqrt()'s and some without. Including or removing sqrt would affect the sharpness of the peaks.
Here's what we get with r=1 and including the sqrt():
r=1;
xs=linspace(0,1.5*r,500);
[y,x,t]=ndgrid(xs,xs, 0:pi/8:pi/2);
z=sum(1./ sqrt( (x-r*cos(t)).^2 + (y-r*sin(t)).^2 ), 3);
surf(x(:,:,1),y(:,:,1),z,'EdgeColor','none','FaceAlpha',0.3);
zlim([0,30]); view([-45,30])
axis vis3d
grid off

请先登录，再进行评论。

Walter Roberson
Walter Roberson 2021-3-9
I showed the strategy in https://www.mathworks.com/matlabcentral/answers/766426-how-can-i-handle-division-by-zero-inside-the-symsum-function#answer_642441
It would only take a small modification to the z = formula I showed there.

类别

在 Help Center 和 File Exchange 中查找有关 Axis Labels 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by