Summing over indices using symsum

3 次查看（过去 30 天）
RH
RH 2023-3-17
编辑： Torsten 2023-3-17
Ran in:
I am trying to calculate electric field at position due to charges
, at positions .
The equation is .
If I can figure out how to for example just get this sum , I would be able to figure out the equation for the electric field. I tried
u=symunit;
syms k
e0=8.854E-12*(u.s)^4*(u.A)^2/((u.kg)*(u.m)^3);
Q1=10E-9*(u.C);
Q2=15-9*(u.C);
Q3=7E-9*(u.C);
r1=[1 4 2]*1E-3*(u.m);
r2=[-2 0 5]*1E-3*(u.m);
r3=[-1 -2 3]*1E-3*(u.m);
r=[5 5 5]*1E-3*(u.m);
Q=symsum(Qk,k,1,3)
Unrecognized function or variable 'Qk'.
That obviously didn't work. Is there a way to make it work with 'symsum'?
  5 个评论
显示 3更早的评论
Walter Roberson
Walter Roberson 2023-3-17
If you have somehow created a list of numbered variables that is too long to make it practical to put them into a list by name, then what you should do is...
Not create a list of numbered variables that is too long to make it practical to put them into a list by name.
Dyuman Joshi
Dyuman Joshi 2023-3-17
编辑：Dyuman Joshi 2023-3-17
Ran in:
Note that the sum isn't just over Q, but rk as well
(Without syms and symunit but values according SI units)
e0=8.854E-12;
Q=[10E-9 15E-9 7E-9];
R=[1 4 2;-2 0 5;-1 -2 3]*1E-3;
r=[5 5 5]*1E-3;
r0=vecnorm(r-R,2,2).^3;
out=sum(Q.*(r-R)./(4*pi*e0*r0),1)
out = 1×3
1.0e+06 * 4.3424 3.1998 1.5735

请先登录，再进行评论。

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

采纳的回答

Torsten
Torsten 2023-3-17
Ran in:
u=symunit;
Q(1)=10E-9*(u.C);
Q(2)=15E-9*(u.C);
Q(3)=7E-9*(u.C);
Qsum = sum(Q)
Qsum = 
  1 个评论
Torsten
Torsten 2023-3-17
编辑：Torsten 2023-3-17
Ran in:
u=symunit;
e0=8.854E-12*(u.s)^4*(u.A)^2/((u.kg)*(u.m)^3);
Q(1)=10E-9*(u.C);
Q(2)=15E-9*(u.C);
Q(3)=7E-9*(u.C);
rk(1,:)=[1 4 2]*1E-3*(u.m);
rk(2,:)=[-2 0 5]*1E-3*(u.m);
rk(3,:)=[-1 -2 3]*1E-3*(u.m);
r=[5 5 5]*1E-3*(u.m);
Qsum = zeros(1,3);
for k = 1:3
Qsum = Qsum + Q(k)*(r-rk(k,:))/(norm(r-rk(k,:)))^3;
end
Qsum = Qsum/(4*pi*e0);
simplify(Qsum.')
ans = 
Data = separateUnits(Qsum.');
double(Data)
ans = 3×1
1.0e+06 * 4.6438 2.2614 2.1837

请先登录，再进行评论。

更多回答（1 个）

Walter Roberson
Walter Roberson 2023-3-17
You cannot use a symbolic variable as an index.
I think it is unlikely that it will ever be supported to use symbolic variables as indices.
Create the definite values and sum() the definite values.

类别

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

标签

产品


版本

R2022a

Community Treasure Hunt

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

Start Hunting!

Translated by