Universidad

Cálculos numéricos Nacional

de Trujillo

UNIVERSIDAD NACIONAL DE TRUJILLO

ESCUELA DE FÍSICA

Universidad
Nacional de Trujillo
La única

“MÉTODO DE IMÁGENES CON MATLAB”

Curso
ELECTROMAGNETISMO

Autor
Prof. Kelman W. Marín Rengifo

07 de octubre, 2019
 Universidad
Cálculos numéricos Nacional
de Trujillo

Contenidos

I. Carga inducida en un conductor plano a Tierra.

II. Carga inducida en un conductor esférico a Tierra.
III. Esfera cargada cercana a un plano conductor a Tierra.
IV. Dos esferas conductoras con potenciales 0 y V.
V. Fuerza entre dos esferas conductoras.
 Universidad
Carga inducida en un plano conductor a Tierra Nacional
de Trujillo

El caso más sencillo es el de una carga q situada a una distancia d de

una placa conductora conectada a Tierra. La placa puede
reemplazarse por una carga imagen -q, tal como se muestra en la
Figura 1.
Consideremos el sistema de cargas de la figura,
formada por dos cargas +q y -q separadas una
distancia 2d. El potencial en el punto P de
coordenadas (x,y) producido por las dos cargas es

Por simetría, el potencial en el eje Y y en cualquier

punto del plano perpendicular a la línea que une
ambas cargas que contiene a dicho eje, es cero.
El campo eléctrico E en las proximidades de dicho
plano conductor es perpendicular al plano, (tiene
componente x: Ex) tal como se ve en las líneas de
fuerza de color blanco de la Figura 1.
 Graficar en matlab 𝑬𝒙 .
Carga inducida en el plano conductor
En los puntos del plano conductor, x=0.

Figura 02:
Su respectivo código en Matlab, es:

Al compilarlo, obtenemos:
El campo eléctrico en las proximidades de un
conductor, tiene por módulo: Ex=σ/ε0. Donde σ es
la densidad de carga. Como el sentido del campo
eléctrico es hacia el conductor, la densidad de
carga σ es negativa σ=-Ex·ε0

Su código en Matlab, será:

 -9 Densidad de carga en plano conductor a Tierra
x 10
4

3.5

2.5
s (C/m )
2

1.5

0.5

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
y (m)
Obs: La distribución de carga inducida en el plano
conductor tiene simetría cilíndrica alrededor del eje
perpendicular al plano y que pasa por la carga. La
carga dq inducida en al anillo de radio y, de
anchura dy es dq=σ·2πy·dy (ver Figura 02).

Carga total inducida en el plano conductor:

Se calcula integrando respecto de y entre 0 e ∞:

Obs: La carga inducida en el plano conductor es

tanto mayor cuanto mayor sea la intensidad del
campo eléctrico Ex, en el origen y va disminuyendo a
medida que nos alejamos radialmente a lo largo del
plano.
 Universidad
Carga inducida en un conductor esférico a Tierra Nacional
de Trujillo

Supongamos un sistema formado por una carga puntual Q en las

proximidades de una esfera conductora a potencial cero a una
distancia d de su centro. El método de las imágenes nos permite
sustituir el conductor por una carga "imagen" q que anulará el
potencial sobre la superficie esférica de radio R.

Figura 03.
El potencial en el punto P1 de la superficie esférica deberá ser
cero:
El potencial en el punto P3 diametralmente opuesto
deberá ser cero.

Tenemos un sistema de dos ecuaciones con dos

incógnitas de las que despejamos q y b:

Potencial electrostático:
Vamos a calcular el campo en cualquier punto P
exterior a la esfera conductora. Primero calculamos
el potencial
donde:

Carga inducida en el conductor esférico

La componente radial del campo es:

El campo en las proximidades de una

superficie conductora es
perpendicular a dicha superficie: Er.
Evaluamos Er para r=R y luego la 𝜎
en la esfera, considerando que el
campo en las proximidades de la
superficie de un conductor
es Er=σ /ε0
Obs: Si integramos esta densidad de carga sobre la
superficie de la esfera, obtenemos la carga total
inducida q, es decir:

La carga total inducida sobre la

esfera conductora es igual a la
carga imagen que sustituye a la
esfera.
Obs: Este resultado es consecuencia de la ley de
Gauss: si dibujamos una superficie cerrada que abarque a la
esfera pero que esté muy próxima a ésta, el flujo del campo
eléctrico a través de esta superficie debe ser la misma,
independientemente, de que la carga esté distribuida sobre la
esfera conductora o haya sido sustituida por la carga puntual
imagen situada en su interior.
Energía electrostática
Cuando la carga Q está a una distancia x del centro
de la esfera conductora conectada a tierra, la fuerza
de atracción F (flecha de color rojo) entre ambas es
igual a la fuerza de atracción entre la carga Q y su
imagen q situada a una distancia x-b

Calculamos el trabajo de la fuerza Fm=F que hay

que hacer para traer la carga Q desde el infinito
hasta una distancia d del centro de la esfera
conductora.
El código, es:
Al compilarlo, obtenemos:
 Universidad
Esfera cargada cercana a un plano conductor a Tierra Nacional
de Trujillo

Sustituiremos la esfera y el plano por una sucesión de cargas

puntuales de signos contrarios que converge a cero rápidamente
y que hacen que las dos superficies (esfera y plano) sean
equipotenciales.
Supongamos que la esfera de radio r, está a un potencial V y su
centro dista d>r del plano a potencial cero.
1. Colocamos una
carga q0 en el centro de
la esfera. Esto hace que
la superficie esférica de
radio r esté a un
potencial V.

2. Colocamos una carga –q0 a

una distancia 2d del centro
de la esfera. Esto hace que
el plano sea una superficie
equipotencial, pero ya no lo
es la esfera.
3. Colocamos una carga q1 en
el interior de la esfera.
Calculamos el valor de q1 y
su posición x1 para que la
esfera sea una superficie
equipotencial, aunque deje
de serlo el plano.
• ΦA (-r, 0) = 0: debido a las cargas –q0 y q1:

• ΦB (r, 0) = 0: debido a las cargas –q0 y q1:

Luego, despejamos q1 y x1 del sistema de ecuaciones:

4. Colocamos una carga –

q1 simétrica a q1 para el
plano sea equipotencial,
pero deja de serlo la
superficie esférica
5. Colocamos una carga q2 en el interior de la esfera, para que
esta sea equipotencial, aunque el plano deje de serlo

• ΦA (-r, 0) = 0: debido a las cargas –q1 y q2:

• ΦB (r, 0) = 0: debido a las cargas –q1 y q2:

Luego, despejamos q2 y x2 del sistema de ecuaciones:

Relaciones recursivas
Calculamos la sucesión de cargas qi y sus posiciones xi mediante
las relaciones recursivas:

Campo y potencial eléctrico:

El campo y el potencial
producido en el punto
P (x, y) por el par de
cargas qi situada en el
punto xi y su simétrica
–qi en la posición 2d-
xi:
El módulo del campo E1 producido por la carga qi es:

El módulo del campo E2 producido por la carga simétrica -

qi es:

Como las componentes del campo eléctrico total Ei son:

Entonces:
• El potencial Vi en el punto P debido a las dos cargas es:

Sabemos que el campo y el potencial total es la suma de

todos los campos y potenciales producidos por los pares de
cargas dispuestas simétricamente al plano, es decir:
Al compilarlo, obtenemos:
 Universidad
Dos esferas conductoras de potenciales 0 y V Nacional
de Trujillo

Aplicaremos el método de las imágenes al cálculo del campo

eléctrico producido por dos esferas conductoras, con una de las
cuales está a potencial cero.

Consideremos el caso de dos esferas de radios r y R cuyos

centros están separados una distancia d>r+R. La primera
esfera está a una potencial V y la segunda esfera está
conectada a Tierra, V=0.

Sustituiremos las dos esferas por dos sucesiones de cargas

puntuales que convergen rápidamente a cero y que hacen que
las dos superficies esféricas sean equipotenciales.
1. Colocamos en el centro de la
primera esfera de radio r una
carga q0, de modo que, el
potencial de la esfera sea V.

2. La superficie esférica de
radio R deja de ser
equipotencial, por lo que
ponemos una carga Q1 en el
interior de la segunda esfera a
una distancia X1 de su centro
Calculamos el valor de Q1 y su posición X1 para que la segunda
esfera de radio R sea una superficie equipotencial, aunque
deje de serlo la primera esfera de radio r.
• ΦC (d-R, 0) = 0 debido a las cargas q0 y Q1:
 • ΦD (d+R, 0) = 0, debido a las cargas q0 y Q1:

Despejamos Q1 y X1 de este sistema de dos ecuaciones:

3. La superficie esférica de radio r deja de ser equipotencial, por

lo que ponemos una carga q1 en el interior de primera esfera
a una distancia x1 de su centro

Calculamos el valor de q1 y su posición x1 para que la primera

esfera de radio r sea una superficie equipotencial, aunque
deje de serlo la segunda esfera de radio R.
• ΦA (-r, 0) = 0, debido a las cargas q1 y Q1:

• ΦB (r, 0) = 0, debido a las cargas q1 y Q1:

Despejamos q1 y x1 en este sistema de dos ecuaciones:

4. La superficie esférica de radio R deja de ser equipotencial,

por lo que ponemos una carga Q2 en el interior de la
segunda esfera a una distancia X2 de su centro
Calculamos el valor de Q2 y su posición X2 para que la segunda
esfera de radio R sea una superficie equipotencial, aunque
deje de serlo la primera esfera de radio r.

• ΦC (d-r, 0) = 0, debido a las cargas Q2 y q1:

• ΦD (d+r, 0) = 0, debido a las cargas Q2 y q1:

Despejamos Q2 y X2 en este sistema de dos ecuaciones:

5. La superficie esférica de radio r deja de ser equipotencial,
por lo que ponemos una carga q2 en el interior de primera
esfera a una distancia x2 de su centro

Calculamos el valor de q2 y su posición x2 para que la primera

esfera de radio r sea una superficie equipotencial, aunque deje
de serlo la segunda esfera de radio R.

• ΦA (-r, 0) = 0, debido a las cargas q2 y Q2:

• ΦA (-r, 0) = 0, debido a las cargas q2 y Q2:

Despejamos q2 y x2 en este sistema de dos ecuaciones:

y así, sucesivamente.

Relaciones recursivas
Calculamos la sucesión de cargas qi, Qi y sus posiciones xi,
Xi mediante las relaciones recursivas:
Campo y potencial producido por la sucesión de cargas
Calculamos el campo y el potencial producido en el punto P
(x, y) por el par de cargas qi situada en el punto xi y Qi en la
posición d-Xi

El campo E1 producido
por la carga qi es:

El campo E2 producido
por la carga Qi es

Las componentes del campo total Ei son:

Reemplazando:
• El potencial Vi en el punto P debido a las dos cargas es:

• El campo y el potencial total es la suma de todos los

campos y potenciales producidos por los pares de cargas
dispuestas simétricamente al plano:
• El correspondiente código en Matlab, es:

Continúa:
Al compilarlo, obtenemos:
• También, obtenemos:
 Universidad
Fuerza entre dos esferas conductoras Nacional
de Trujillo

Aplicaremos el método de las imágenes a dos esferas igual

radio R e igual carga Q, separadas entre sus centros una
distancia d>2R. Situaremos el origen en el centro de la primera
esfera.
La simetría del problema nos permite reemplazar las esferas
cargadas por dos sucesiones idénticas de cargas puntuales q0,
q1…qn… situadas en la línea que une los centros en las
posiciones x0, x1…xn … y d-x0, d-x1…d-xn …, respectivamente.
Nos fijaremos en la primera esfera cuyo centro está en el
origen.
1. Colocamos una carga q0 en el centro de la primera esfera,
en el origen x0=0. La superficie de la primera esfera de
radio R es equipotencial. Colocamos una carga
idéntica q0 en el centro de la segunda esfera. La superficie
de la primera esfera deja de ser equipotencial.

2. Colocamos una carga q1 a una distancia x1 del centro de la

primera esfera
El potencial producido por la carga q0 situada en el centro de la
segunda esfera y por la carga adicional q1 debe cancelarse en
la superficie de la primera esfera, para que siga siendo
equipotencial. El potencial en el punto A (R, 0) y B(-R, 0)
producido por ambas cargas debe ser cero, es decir:

Solucionando el sistema de dos ecuaciones con dos incógnitas:

3. Si colocamos una carga q1 en la posición d-x1, la superficie

de la primera esfera deja de ser equipotencial.
4. Colocamos una carga q2 en la posición x2, tal como se
muestra en la figura.

El potencial producido por la carga q1 situada en la segunda

esfera y por la carga adicional q2 debe cancelarse en la
superficie de la primera esfera, para que siga siendo
equipotencial. El potencial en el punto A (R, 0) y B(-R, 0)
producido por ambas cargas debe ser cero, es decir:
Para solucionar el sistema de ecs, eliminamos q1 y q2 y
despejamos x2. Luego, despejamos q2 de la primera o de la
segunda ecuación, obteniendo:

Ó también:

5. Si colocamos una carga q2 en la posición d-x2 la superficie

de la primera esfera deja de ser equipotencial.

6. Colocamos una carga q3 en la posición x3, tal como se

muestra:
El potencial producido por la carga q2 situada en la segunda
esfera y por la carga adicional q3 debe cancelarse en la
superficie de la primera esfera, para que siga siendo
equipotencial. El potencial en el punto A (R, 0) y B(-R, 0)
producido por ambas cargas debe ser cero, es decir:

Para solucionar el sistemas de ecs, eliminamos q2 y q3 y

despejamos x3, luego, despejamos q3 de la primera o de la
segunda ecuación, obteniendo:

Ó también:
7. Continuamos es proceso iterativo, obteniendo el valor de la
carga imagen q4 en la posición x4.

• Obs: La carga total de cada esfera es:

que es proporcional a q0.

Cálculo aproximado de la fuerza de repulsión entre dos

esferas
Observar que las cargas imágenes q0…qn y sus
posiciones x0…..xncorrespondientes a la primera esfera, que
son las mismas que para la segunda esfera pero en las
posiciones d-x0…..d-xn, mediante las relaciones de recurrencia:
Se calcula la carga inicial q0 teniendo en cuenta que la carga de
cada esfera es Q.

Se calcula la fuerza entre los dos conjuntos de cargas imágenes

Luego se compara con la fuerza repulsiva de Coulomb para 2

cargas puntuales.

Continúa….
Al compilarlo, obtenemos: