1

Universidad Politécnica Salesiana. Teoría Electromagnética

SOLUCIÓN NUMÉRICA DE LA ECUACIÓN DE

LAPLACE, POTENCIAL ELECTROSTÁTICO
Guamán E. Telmo/ Quizhpi C. Mateo/ Velecela C. Juan

Abstract — this paper presents the solution the Laplace

equation using a numerical method for electric potential in a
certain region of space, knowing their behavior on the border En el caso bidimensional tratamos la placa como una
of the region. malla de puntos discretos.

1. Introducción

En este documento de tipo educativo se presenta la

solución de la ecuación de Laplace mediante un método
numérico conocido como diferencias finitas, para
potenciales eléctricos en cierta región del espacio,
conociendo su comportamiento o valor en la frontera de
dicha región.
De igual manera se procederá a graficar de las líneas de
flujo magnético que se producen en dicha región.

2. Objetivos

o Determinar el número de ecuaciones adecuadas

de potencial eléctrico. Luego aproximamos las derivadas parciales en cada
o Definir o plantear las ecuaciones punto de la malla transformando la ecuación diferencial
correspondientes a cada punto de la región. en una ecuación algebraica.
o Obtener las soluciones de los potenciales en Las diferencias centrales basadas en la malla de la
figura anterior son:
cada punto.
o Graficar las líneas de flujo. V(i 2V ( i , j ) V(i
V 1 j) 1, j )
2 2
3. Marco Teórico x x

V V(i, j 1)
2V ( i , j ) V(i, j 1)
3.1 Ecuación de Laplace 2 2
y y
La ecuación de Laplace, se utiliza para modelar diversos
Las cuales tienen errores de O x
2
y O y
2
:
problemas que tienen que ver con el potencial de una
variable desconocida. Sustituyendo en la ecuación de Laplace para el
Expresión de la ecuación de Laplace para potencial potencial
eléctrico.
2 2 V(i 2V ( i , j ) V( i V( i , j 2V ( i , j ) V( i , j
V V 1, j ) 1, j ) 1) 1)
0 2 2
0
2 2
x y x y
3.2 Solución numérica En la malla cuadrada de la figura anterior, x y
La solución numérica, se basa en el método de Reagrupando términos la ecuación nos queda:
diferencias finitas. V( i V( i V( i , j V( i , j 4V ( i , j ) 0
1, j ) 1, j ) 1) 1)
 2
Universidad Politécnica Salesiana. Teoría Electromagnética

Esta relación, que satisface para todos los puntos → Ahora procedemos a obtener las ecuaciones
interiores de la placa, se conoce como ecuación correspondientes a cada nodo, mediante la ecuación
Laplaciana en diferencias. definida con anterioridad:

3.2 Condiciones frontera V( i 1, j )

V( i 1, j )
V( i , j 1)
V( i , j 1)
4V ( i , j ) 0
Debemos además especificar las condiciones de frontera Nodo 1:
en los extremos de la placa para obtener una solución
única. V 2 V13 4V1 220
El caso más simple es aquel donde el potencial eléctrico Nodo 2:
en la frontera es un valor fijo, a este tipo de condición se V1 V 3 V14 4V 2 100
la conoce como “condición de frontera de Dirichlet”.
Otro tipo de condición es la “condición de frontera de Nodo 3:
Neumann” la cual tiene como dato la derivada en la V2 V4 V15 4V 3 100
frontera. [1]
Nodo 4:
4. Metodología V3 V5 V16 4V 4 100
4.1 Identificación Nodo 5:
→ Primero debemos definir el punto inicial y tratar a la V4 V6 V17 4V 5 100
figura como una malla de puntos discretos; para nuestra
Nodo 6:
figura hemos considerado adecuado tomar al punto
V5 V7 V18 4V 6 100
ubicado en el extremo superior izquierdo como el punto
de coordenadas (1,1) al cual denominaremos como V1, Nodo 7:
de igual manera se va dando denominaciones a cada V 6 V8 V19 4V 7 100
punto coordenado de forma esquemática o secuencial;
como se observa en la siguiente figura: Nodo 8:
V7 V9 V 20 4V 8 100

Nodo 9:
V 8 V10 V 21 4V 9 100
 3
Universidad Politécnica Salesiana. Teoría Electromagnética
Nodo 10: Nodo 28:
V 9 V11 V 22 4V10 100 V18 V 27 V 29 V 36 4V 28 0
Nodo 11: Nodo 29:
V10 V12 V 23 4V11 100 V19 V 28 V 30 V 37 4V 29 0
Nodo 12: Nodo 30:
V11 V 24 4V12 200 V 20 V 29 V 38 4V 30 20
Nodo 13: Nodo 31:
V1 V14 V 25 4V13 120 V 23 V 32 V 39 4V 31 20
Nodo 14: Nodo 32:
V 2 V13 V15 V 26 4V14 0 V 24 V 31 V 40 4V 32 100
Nodo 15: Nodo 33:
V 3 V14 V16 4V15 60 V 25 V 34 V 41 4V 33 120
Nodo 16: Nodo 34:
V 4 V15 V17 4V16 60 V 26 V 33 V 42 4V 34 60
Nodo 17: Nodo 35:
V 5 V16 V18 V 27 4V17 0 V 27 V 36 V 45 4V 35 60
Nodo 18: Nodo 36:
V 6 V17 V19 V 28 4V18 0 V 28 V 35 V 37 V 46 4V 36 0
Nodo 19: Nodo 37:
V 7 V18 V 20 V 29 4V19 0 V 29 V 36 V 38 V 47 4V 37 0
Nodo 20: Nodo 38:
V 8 V19 V 21 V 30 4V 20 0 V 30 V 37 V 48 4V 38 20
Nodo 21: Nodo 39:
V 9 V 20 V 22 4V 21 20 V 31 V 40 V 51 4V 39 20
Nodo 22: Nodo 40:
V10 V 21 V 23 4V 22 20 V 32 V 39 V 52 4V 40 100
Nodo 23: Nodo 41:
V11 V 22 V 24 V 31 4V 23 0 V 33 V 42 V 53 4V 41 120
Nodo 24: Nodo 42:
V12 V 23 V 32 4V 24 100 V 34 V 41 V 43 V 54 4V 42 0
Nodo 25: Nodo 43:
V13 V 26 V 33 4V 25 120 V 42 V 44 V 55 4V 43 60
Nodo 26: Nodo 44:
V14 V 25 V 34 4V 26 60 V 43 V 45 V 56 4V 44 60
Nodo 27: Nodo 45:
V17 V 28 V 35 4V 27 60 V 35 V 44 V 46 V 57 4V 45 0
 4
Universidad Politécnica Salesiana. Teoría Electromagnética

Nodo 46: Nodo 64:

V 36 V 45 V 47 V 58 4V 46 0 V 52 V 63 V 68 4V 64 100

Nodo 47: Nodo 65:

V 37 V 46 V 48 V 59 4V 47 0 V 59 V 66 V 69 4V 65 140

Nodo 48: Nodo 66:

V 38 V 47 V 49 V 60 4V 48 0 V 60 V 65 V 70 4V 66 60

Nodo 49: Nodo 67:

V 48 V 50 V 61 4V 49 20 V 63 V 68 V 71 4V 67 60

Nodo 50: Nodo 68:

V 49 V 51 V 62 4V 50 20 V 64 V 67 V 72 4V 68 100

Nodo 51: Nodo 69:

V 39 V 50 V 52 V 63 4V 51 0 V 65 V 70 V 73 4V 69 140

Nodo 52: Nodo 70:

V 40 V 51 V 64 4V 52 100 V 66 V 69 V 74 4V 70 60

Nodo 53: Nodo 71:

V 41 V 54 4V 53 260 V 67 V 72 V 77 4V 71 60

Nodo 54: Nodo 72:

V 42 V 53 V 55 4V 54 140 V 68 V 71 V 78 4V 72 100

Nodo 55: Nodo 73:

V 43 V 54 V 56 4V 55 140 V 69 V 74 V 79 4V 73 140

Nodo 56: Nodo 74:

V 44 V 55 V 57 4V 56 140 V 70 V 73 V 75 V 80 4V 74 0

Nodo 57: Nodo 75:

V 45 V 56 V 58 4V 57 140 V 74 V 76 V 81 4V 75 60

Nodo 58: Nodo 76:

V 46 V 57 V 59 4V 58 140 V 75 V 77 V 82 4V 76 60

Nodo 59: Nodo 77:

V 47 V 58 V 60 V 65 4V 59 0 V 71 V 76 V 78 V 83 4V 77 0

Nodo 60: Nodo 78:

V 48 V 59 V 61 V 66 4V 60 0 V 72 V 77 V 84 4V 78 100
Nodo 61: Nodo 79:
V 49 V 60 V 62 4V 61 60 V 73 V 80 4V 79 260
Nodo 62: Nodo 80:
V 50 V 61 V 63 4V 62 60 V 74 V 79 V 81 4V 80 120
Nodo 63: Nodo 81:
V 51 V 62 V 64 V 67 4V 63 0 V 75 V 80 V 82 4V 81 120
 5
Universidad Politécnica Salesiana. Teoría Electromagnética

Nodo 82: → Con lo cual obtuvimos los siguientes resultados:

V 76 V 81 V 83 4V 82 120
V1 = 103.1692 V29 = 66.5954 V57 = 113.9050
Nodo 83: V2 = 95.3165 V30 = 50.7229 V58 = 108.6436
V 77 V 82 V 84 4V 83 120 V3 = 90.8177 V31 = 56.8711 V59 = 92.0812
Nodo 84: V4 = 89.2299 V32 = 80.6902 V60 = 74.5548
V 78 V 83 4V 84 220 V5 = 89.3010 V33 = 98.6118 V61 = 60.5338
V6 = 88.7244 V34 = 84.5878 V62 = 59.8298
→ Entonces armamos el sistema de matrices Ax=B V7 = 86.6095 V35 = 74.9719 V63 = 70.8871
Como el sistema de matriz es de 84x84; realizamos un V8 = 83.1525 V36 = 76.1645 V64 = 85.1911
programa en matlab que nos permita obtener las V9 = 79.9485 V37 = 67.6837 V65 = 107.3954
soluciones de forma inmediata empleado el método de V10 = 80.8699 V38 = 50.2442 V66 = 83.0343
Eliminación de Gauss. V11 = 86.4455 V39 = 55.0933 V67 = 74.5144
Seudocódigo del programa en matlab.
V12 = 93.2115 V40 = 79.4892 V68 = 87.7041
disp('Solucion Numerica de la Ecuacion de V13 = 97.3603 V41 = 110.8664 V69 = 114.4660
Laplace');
ls; V14 = 87.2791 V42 = 102.0243 V70 = 90.1870
A = xlsread('Laplace.xlsx'); %Importamos V15 = 78.7244 V43 = 92.1919 V71 = 79.4665
la matriz de Coeficientes desde Excel V16 = 76.8010 V44 = 89.7921 V72 = 91.1109
B = xlsread('Laplace.xlsx',2);
%Importamos la matriz de Terminos V17 = 79.2495 V45 = 91.8143 V73 = 120.2817
Independientes desde Excel, Hoja 2 V18 = 78.9872 V46 = 88.5883 V74 = 103.2477
[n,m] = size(A); V19 = 74.5610 V47 = 77.7308 V75 = 89.1517
C = [A,B];
V20 = 66.0520 V48 = 62.5701 V76 = 86.5384
disp('La Matriz resultante es: '); V21 = 55.7715 V49 = 47.7506 V77 = 92.2408
%Matriz Aumentada V22 = 57.0855 V50 = 47.8984 V78 = 97.2728
disp(C); V23 = 71.7007 V51 = 64.0130 V79 = 123.4130
V24 = 86.4006 V52 = 82.1734 V80 = 113.3704
for k=1:(n-1) %Eliminacion hacia Abajo
for i=(k+1):n V25 = 78.9930 V53 = 122.8294 V81 = 106.8208
m(i,k)= C(i,k)/C(k,k); V26 = 77.7150 V54 = 120.4513 V82 = 104.7611
for j=k:(n+1) V27 = 71.9089 V55 = 116.9513 V83 = 105.6854
C(i,j)=(C(i,j)-(m(i,k)*C(k,j)));
end V28 = 73.4140 V56 = 115.1621 V84 = 105.7395
end
end → Ahora procedemos a obtener las gráficas de las líneas
de flujo y de las superficies equipotenciales, de igual
for i=n:-1:1 %Sustitucion hacia Arriba manera desarrollaremos un programa en matlab, el cual
s = 0; nos permita obtener dichas gráficas.
for b = (i+1):n
s = s + (C(i,b)*X(b));
end
X(i) = (C(i,n+1)-s)/(C(i,i));
end
disp('Matriz de Soluciones');
disp(X); %Soluciones

disp('Las soluciones son: ');

%Imprimimos las Soluciones V
for i=1:n
Xi = X(1,n);
fprintf('\nV%g', i);
disp(X(i));
end
 6
Universidad Politécnica Salesiana. Teoría Electromagnética

 Curvas de nivel

 Líneas de Flujo
 7
Universidad Politécnica Salesiana. Teoría Electromagnética

 Superficies Equipotenciales
 8
Universidad Politécnica Salesiana. Teoría Electromagnética

Tabla de potenciales en porcentaje según la posición cartesiana.

y/x 0 1 2 3 4 5 6 7 8 9 10 11 12 13
0 0 0 0 0 0 0 120 120 120 120 120 120 120 120
1 0 0 0 0 0 0 140 123,41 113,37 106,82 104,76 105,69 105,74 100
2 0 0 0 0 0 0 140 120,28 103,25 89,15 86,54 92,24 97,27 100
3 0 0 0 0 0 0 140 114,47 90,19 60,00 60,00 79,47 91,11 100
4 160 140 140 140 140 140 140 107,40 83,03 60,00 60,00 74,51 87,70 100
5 120 122,8 120,5 117,0 115,2 113,9 108,6 92,1 74,6 60,5 59,8 70,9 85,2 100
6 120 110,9 102,0 92,2 89,8 91,8 88,6 77,7 62,6 47,8 47,9 64,0 82,2 100
7 120 98,6 84,6 60,0 60,0 75,0 76,2 67,7 50,2 20,0 20,0 55,1 79,5 100
8 120 79,0 77,7 60,0 60,0 71,9 73,4 66,6 50,7 20,0 20,0 56,9 80,7 100
9 120 97,4 87,3 78,7 76,8 79,2 79,0 74,6 66,1 55,8 57,1 71,7 86,4 100
10 120 103,2 95,3 90,8 89,2 89,3 88,7 86,6 83,2 79,9 80,9 86,4 93,2 100
11 220 100 100 100 100 100 100 100 100 100 100 100 100 100

Potenciales representación porcentual

250

200

200-250
150
150-200
100-150
100
50-100
0-50
50
9
6
0
0 3
1 2 3 4 5 6 7 8 0
9 10 11 12 13
 9
Universidad Politécnica Salesiana. Teoría Electromagnética

5. Conclusiones

Luego de la realización de este proyecto podemos

concluir que:
El proyecto en sí fue algo tedioso,
específicamente al momento de plantear las
ecuaciones para los puntos, lo cual fue un
proceso sencillo, pro que tomó tiempo, también
al momento de armar la matriz en Excel.
Es bastante conveniente pasar las ecuaciones a
una hoja de cálculo en Excel para de allí poder
importar al programa en matlab directamente y
mandar a ejecutar para una rápida obtención de
las soluciones.
Podemos constatar que los potenciales
obtenidos nos ayudan a identificar la manera en
que está distribuido el mismo sobre la región.
Si se hubiese tomado mayor número de
divisiones, el cálculo hubiera resultado mucho
mejor, pero también más tedioso, sin embargo
los resultados obtenidos fueron suficientes para
establecer las líneas de flujo y las superficies
equipotenciales, lo cual era el objetivo d este
proyecto.

6. Referencias

[1] APPLIED NUMERICAL METHODS USING

MATLAB, Won Young Yang.

[2] CLASSICAL ELECTRODYNAMICS, John

Davis Jackson, John Wiley & Sons, Inc., Publication