Guia3COE-Ley de Gauss

Solucion de leyes Maxwell con Matlab
Jorge Benjamín Torres López. Wilber Alejandro Reyes Álvarez. Luis Alexander Canjura Sánchez.
Estudiante de Ingeniería en Telecomunicaciones. Estudiante de Ingeniería en Automatización. Estudiante de Ingeniería Eléctrica.
Universidad Don Bosco. Universidad Don Bosco. Universidad Don Bosco.
San Salvador, El Salvador. San Salvador, El Salvador. San Salvador, El Salvador.
benjaykris06@gmail.com reyesalvarezwilberalejandro11@gmail.com luis_10_sanchez@hotmail.com
Summary— Este documento demostramos la utilización del un campo vectorial se alejan o se aproximan a un vector normal
software Matlab para graficar superficies Gaussianas e integrar a esta superficie. El rotacional de un campo vectorial es,
funciones para aplicarlo a la resolución de problemas de flujo también, un campo vectorial, o sea a cada punto del espacio
eléctrico mediante la Ley de Gauss. donde definimos el rotacional será dado por un vector.
Keywords— divergencia, Rotacional, vector, integrales, producto
punto, producto cruz, ángulo. La Ley de Gauss establece que el flujo de campo eléctrico
que atraviesa una superficie cerrada es igual a la carga neta
situada en su interior dividida por la constante dieléctrica del
medio.
I. NOMENCLATURA
ΦE=∮SE→⋅dS→ = Qε
Símbolo Nombre donde:
𝛁. F Divergencia ΦE es el flujo neto de carga
𝛁xF Rotacional
D Vector de densidad de flujo eléctrico E→ es la intensidad de campo eléctrico
𝜌𝑣 Densidad de carga volumétrica dS → es un diferencial del vector de superficie (trozo
ʃ Integral lineal elemental de superficie)
Q es la carga contenida en la superficie
II. INTRODUCCIÓN ε es la constante dieléctrica del medio.
Las integrales dobles son una manera de integrar sobre una

En esta práctica demostraremos la aplicación correcta de los region bidimensional. Entre otras cosas, nos permiten calcular el
comandos para integrales en Matlab para el desarrollo de volumen bajo una superficie.
integrales dobles y triples.
También demostraremos los ejercicios de aplicación de Gauss.
Las integrales triples son el análogo de las integrales dobles
Estos conocimientos y las competencias se lograron conociendo para tres dimensiones. Son una herramienta para sumar infinitas
los conceptos esenciales tales como magnitudes físicas y sus cantidades infinitesimales asociadas con puntos de una región
respectivas clasificaciones, Producto Cruz, Producto punto, tridimensional
Divergencia, Rotacional, Ley de Gauss y también el En Las leyes de Maxwell se muestra como los campos
conocimiento sobre integrales: DOBLES Y TRIPLES. magnéticos y los campos eléctricos están relacionados entre sí,
tal que cualquier variación del campo eléctrico produce un
campo magnético y cualquier variación del campo magnético
III. SOPORTE TEÓRICO DE DIVERGENCIA, ROTACIONAL, produce un campo eléctrico.
LEY DE GAUSS, INTEGRALES DOBLES Y TRIPLES Y
LEYES DE MAXWELL IV. PASOS PARA GRAFIGAR LAS FIGURAS EN EL
ESPACIO 3D
La Divergencia es un operador que toma una función
vectorial que define a este campo vectorial y arroja como valor Procedimos a inicializar el script para la primera grafica de
de salida una función escalar que mide el cambio de la densidad superficie que en este caso es una esfera.
del fluido en cada punto.
.
El Rotacional En el cálculo vectorial es un operador que
calcula por una superficie infinitesimal cuanto los vectores de
Fig. 1. Comandos para grafica de esfera
Fig. 4. Grafica de cilindro 3D.
V. REALIZACION DE INTEGRALES DEFINIDAD DOBLES

Y TRIPLES .
Creamos un nuevo script para resolver integrales definidas. La

primera a realizar es la función: Determine la integral de x 2 con
límites de 0 a 5
Fig. 5. Comando de la integral definida.

Fig. 2. Grafica de esfera 3D.
Fig. 6. Respuesta de la integral definida
Determine la integral de x 2 con límites de 0 a 5
Fig. 3. Comandos para grafica de cilindro
Fig. 7. Comando de integral definida-

Encuentre el volumen de un cilindro cerrado de altura z=4m
Fig. 8. Respuesta integral al infinito.
Encuentre la Integral doble de x2 y con límites en “x” de 0 a 5 y

en “y” de 0 a 3
Fig. 9. Integral doble,

Fig. 12. Volumen de un cilindro con integrales triple.
Encuentre la Integral doble de x2 y con límites en “x” de 0 a 5
y en “y” de 0 a 3 y de “z” de -1 a 5
VII. DIVERGENCIA Y ROTACIONAL DE UN

VECTOR
Creamos un nuevo archivo .m Calcule el rotacional y

divergencia de los siguientes campos vectoriales: 𝑥2 𝑦𝑧 𝑎𝑥 +
𝑥𝑦2 𝑧 𝑎𝑦 + 𝑥𝑦𝑧2 𝑎𝑧
Fig. 10. Integral triple.
PARTE A
VI. RESOLUCION DE AREAS Y VOLUMENES A TRAVES

DE INTEGRALES
Encontramos el área de la superficie de una esfera con radio r

= 3m aplicando integrales.
Fig. 13. Resolución de Rotacional.

Fig. 11. Área de una esfera, resuelta por integrales dobles.
VIII. LEY DE GAUSS
a) Determinamos el flujo eléctrico total que pasa por un

cubo cuyas aristas mide d=0.20 m, en donde pasa un
campo eléctrico de 𝑬 = 1500𝑥 𝑎𝑥 [𝑉/𝑚]
(n=por cada cara del cubo)
𝐹𝑙𝑢𝑗𝑜 𝑒𝑙𝑒𝑐𝑡𝑟𝑖𝑐𝑜 𝑡𝑜𝑡𝑎𝑙 = ∮ 𝐸⃗ 𝑑𝐴
𝐹𝑙𝑢𝑗𝑜 = 𝐸⃗ ∮ 𝑑𝐴
Fig. 14. Resolución de Divergencia. ∮ 𝑑𝐴 = 𝐿2 𝐸⃗ = 1500𝑥ar
2. 𝑥 cos 𝑦 𝒂𝑥 + 𝑥𝑦 2 𝒂y 𝐹𝑙𝑢𝑗𝑜 = ∮ 𝐸⃗ 𝑑𝑛1 + ∮ 𝐸⃗ 𝑑𝑛2 + ∮ 𝐸⃗ 𝑑𝑛3 + ∮ 𝐸⃗ 𝑑𝑛4

+ ∮ 𝐸⃗ 𝑑𝑛5 + ∮ 𝐸⃗ 𝑑𝑛6
𝐹𝑙𝑢𝑗𝑜 = ∮ 𝐸⃗ 𝑑𝑛1 + ∮ 𝐸⃗ 𝑑𝑛2
𝐹𝑙𝑢𝑗𝑜 = −𝐸⃗ 𝐿2 + 𝐸⃗ 𝐿2
𝐹𝑙𝑢𝑗𝑜 𝑒𝑙𝑒𝑐𝑡𝑟𝑖𝑐𝑜 𝑡𝑜𝑡𝑎𝑙 = 0
Fig. 15. Resolución de Rotacional
Fig. 16. Resolución de Divergencia.
b) Determinamos la carga total En una esfera de radio

de 5cm con un campo eléctrico de 𝑬 = 500𝑟 𝑎𝑟 [𝑉/𝑚]
IX. TRABAJO AUTONOMO
Realice la parte 4 en matlab utilizando ∮ 𝜺𝟎𝑬. 𝒅𝒔 ∫ 𝛁. 𝑫. 𝒅𝑽

para los literales a) y b).
1) Análisis y discusión: ¿Los resultados obtenidos en la
parte 4 para los literales a) y b) son los mismos?
Si, los resultados son los mismos para el literal a y para el

b. Para el literal a El flujo de campo eléctrico estará entrando
por una cara del cubo y saliendo por otra cara por ellos el flujo
eléctrico total tendrá un valor de cero, mientras las otras 4
caras no tendrán un flujo de campo debido a la dirección que
tiene el campo eléctrico y comprobamos los resultados en
Matlab. Así también para el literal b los resultados son los
mismos ya que la figura contiene una carga dentro de ella.
Evaluación: Investigue sobre las funciones curl y
divergence de Matlab, además del comando Jacobian para
determinar el rotacional.
a. Curl: La función curl (V,X) devuelve el rizo del

campo vectorial V con respecto al vector X. El campo
vectorial V y el vector X son tridimensionales.
curl(V)devuelve el rizo del campo vectorial Vcon
respecto al vector de variables devuelto
por symvar(V,3).
c) Sea 𝑬 = 𝑥𝑦𝒂𝑥 + 𝑥 2𝒂𝑦, halle la densidad de carga b. Divergence: Está función calcula la divergencia del
volumétrica ρv. campo vectorial; divergence(X,Y,Z,U,V,W)calcula la
divergencia de un componentes del vector campo
vectorial 3-D que tienen U, V, W. Las matrices X, Y
y Z, que definen las coordenadas de los componentes
del vector U, V y W, debe ser monótona, pero no
necesitan estar espaciados de manera uniforme. X, Y
y Z debe tener el mismo número de elementos.
c. Jacobian: El comando jacobian(f,v) calcula la matriz

jacobiana de fcon respecto a v. El elemento ( i , j ) del
resultado es ∂ f(yo)∂ v(j).
CONCLUSIÓN
Matlab es una herramienta que nos facilita muchas cosas

en esta guía práctica se demostró que es muy útil para
graficar superficies, resolver integrales, gradientes y
rotacional de una manera relativamente sencilla. Y todo esto
nos puede ayudar a resolver problemas más complejos de
una manera más rápida y segura como es el caso de los
ejercicios de la Ley de Gauss.
RECONOMICIENTOS
Luis Canjura, nació en Santa Tecla, El
Al Ing. José F. Sarmiento, por la inducción brindada en la Salvador, el 02 de Julio de 1997.
utilización del software Matlab, el cual nos permitió realizar la Obtuvo su título de bachiller en la
especialidad de Técnico Vocacional
guía práctica de manera exitosa.
Electrotecnia en el Instituto Nacional
Técnico Industrial (2015).
Actualmente continua sus estudios
REFERENCIAS universitarios en la Universidad Don
Bosco, en la carrera de ingeniería
[1] Matthew N. O. Sadiku, Elementos de electromagnetismo, eléctrica.
3ra. ed. Oxford University Press México, 2002, ISBN 970-613-
672-X.
[2] William H. Hayt Jr., John A. Buck, Teoría Benjamín Torres, nació en
electromágnetica, 7ma .ed., McGraw-Hill, 2006, México D.F, Cuscatancingo, El Salvador, el 18 de
C.P. 01376. enero 1997. Obtuvo su titulo de
bachiller general en Colegio Hebrón
de El Salvador (2014). Actualmente
continua sus estudios universitarios
BIOGRAFÍAS en la Universidad Don Bosco, en la
carrera de ingeniería en
Telecomunicaciones.
Wilber Reyes, nació en San Salvador, El
Salvador, el 11 de Diciembre de 1997.
Obtuvo su título de bachiller en la
especialidad de Técnico en
Mantenimiento Automotriz en el
Colegio Don Bosco (2016).
Actualmente continúa sus estudios
universitarios en la Universidad Don
Bosco, en la carrera de Ing. en
Automatización.
Actualmente es voluntario en el Programa de Orientación de

Carrera (POC) de la Universidad Don Bosco.

Guia3COE-Ley de Gauss

Cargado por

Copyright:

Formatos disponibles

Guia3COE-Ley de Gauss

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Guia3COE-Ley de Gauss

Cargado por

Copyright:

Formatos disponibles

Solucion de leyes Maxwell con Matlab

II. INTRODUCCIÓN ε es la constante dieléctrica del medio.

Las integrales dobles son una manera de integrar sobre una

Fig. 4. Grafica de cilindro 3D.

V. REALIZACION DE INTEGRALES DEFINIDAD DOBLES

Creamos un nuevo script para resolver integrales definidas. La

Fig. 5. Comando de la integral definida.

Fig. 6. Respuesta de la integral definida

Determine la integral de x 2 con límites de 0 a 5

Fig. 3. Comandos para grafica de cilindro

Fig. 7. Comando de integral definida-

Fig. 8. Respuesta integral al infinito.

Encuentre la Integral doble de x2 y con límites en “x” de 0 a 5 y

Fig. 9. Integral doble,

VII. DIVERGENCIA Y ROTACIONAL DE UN

Creamos un nuevo archivo .m Calcule el rotacional y

VI. RESOLUCION DE AREAS Y VOLUMENES A TRAVES

Encontramos el área de la superficie de una esfera con radio r

Fig. 13. Resolución de Rotacional.

a) Determinamos el flujo eléctrico total que pasa por un

Fig. 14. Resolución de Divergencia. ∮ 𝑑𝐴 = 𝐿2 𝐸⃗ = 1500𝑥ar

2. 𝑥 cos 𝑦 𝒂𝑥 + 𝑥𝑦 2 𝒂y 𝐹𝑙𝑢𝑗𝑜 = ∮ 𝐸⃗ 𝑑𝑛1 + ∮ 𝐸⃗ 𝑑𝑛2 + ∮ 𝐸⃗ 𝑑𝑛3 + ∮ 𝐸⃗ 𝑑𝑛4

𝐹𝑙𝑢𝑗𝑜 = ∮ 𝐸⃗ 𝑑𝑛1 + ∮ 𝐸⃗ 𝑑𝑛2

𝐹𝑙𝑢𝑗𝑜 𝑒𝑙𝑒𝑐𝑡𝑟𝑖𝑐𝑜 𝑡𝑜𝑡𝑎𝑙 = 0

Fig. 15. Resolución de Rotacional

Fig. 16. Resolución de Divergencia.

b) Determinamos la carga total En una esfera de radio

Realice la parte 4 en matlab utilizando ∮ 𝜺𝟎𝑬. 𝒅𝒔 ∫ 𝛁. 𝑫. 𝒅𝑽

Si, los resultados son los mismos para el literal a y para el

a. Curl: La función curl (V,X) devuelve el rizo del

c. Jacobian: El comando jacobian(f,v) calcula la matriz

Matlab es una herramienta que nos facilita muchas cosas

Actualmente es voluntario en el Programa de Orientación de

También podría gustarte