UNIVERSIDAD DE LAS FUERZAS ARMADAS (E.S.P.

E)
EXTENCION LATACUNGA

DEPARTAMENTO DE CIENCIAS EXACTAS

CALCULO VECTORIAL
TEMA:

INTEGRANTES:

Richard Stalin Pila Caiza

Hugo Manuel Aucapiña

DOCENTE:

Ing. Jorge Sanchez

CURSO:

3° Electronica

Latacunga – Ecuador

2018
 INTRODUCCIÓN

La importancia del estudio del Cálculo Vectorial radica principalmente en que, en

diversas aplicaciones de la ingeniería, la concurrencia de variables espaciales y
temporales, hace necesario el análisis de fenómenos naturales cuyos modelos utilizan
funciones vectoriales o escalares de varias variables.

La asignatura está diseñada de manera que el estudiante pueda representar conceptos, que
aparecen en el campo de la ingeniería por medio de vectores; resolver problemas en los
que intervienen variaciones continuas; resolver problemas geométricos en forma
vectorial; graficar funciones de varias variables; calcular derivadas parciales; representar
campos vectoriales que provengan del gradiente de un campo escalar, así como su
divergencia y rotacional; resolver integrales dobles y triples; aplicar las integrales en el
cálculo de áreas y volúmenes. (Daniel, 2017)

Por esta razón el presente trabajo trata de implementar todo lo aprendido en clase de
Cálculo Vectorial para realizar una simulación o programación que ayude a encontrar la
longitud de arco, centro de masa entre otras temas tratados de cualquier función ingresada
por el teclado , y previamente poderlos visualizar respectivamente en una ventana gráfica,
para esto utilizaremos el software de Matlab que es uno de los programas mas factibles
para realizarlo ya que cuenta con un sistema de programación propio y representativo,
por su gran rendimiento al momento de realizar representación de datos y funciones

Esto nos ayudara a comprender de mejor forma la asignatura de calculo vectorial ya que
mediante este trabajo podremos comprobar los datos calculados en clase de una forma
más comprensible mediante graficas reales de las funciones que deseemos ingresar en la
plataforma de Matlab, también ayudara a fomentar el uso de plataformas o softwares que
ayudan a resolver dudas creadas al momento de realizar cálculos y no entender el
comportamiento de los resultados obtenidos
 DESARROLLO

1. Matlab

Es una herramienta interactiva basada en matrices para cálculos científicos y de ingeniería

(de hecho, el termino Matlab procede de matrix laboratory). Desde el punto de vista del
control, matlab se puede considerar un entorno matemático de simulación que puede
utilizarse para modelar y analizar sistemas. Permitir a el estudio de sistemas continuos,
discretos, lineales y no lineales, mediante descripción interna y externa, en el dominio
temporal y frecuencia. Matlab constituye un entorno abierto, para el cual numerosos
paquetes específicos adicionales (toolboxes) han sido desarrollados. En el caso que nos
ocupa se utilizara fundamentalmente el Control System Toolbox. Estos paquetes
específicos adicionales están constituidos por un conjunto de funciones que pueden ser
llamadas desde el programa y mediante las cuales se pueden realizar multitud de
operaciones. (Carlos, 2018)

2. Longitud de Arco

Longitud de arco es la medida de la distancia o camino recorrido a lo largo de una curva o

dimensión lineal. Históricamente, ha sido difícil determinar esta longitud
en segmentos irregulares; aunque fueron usados varios métodos para curvas específicas,
la llegada del cálculo trajo consigo la fórmula general para obtener soluciones cerradas
para algunos casos. (Perez, s.f.)

Figura 1. Longitud de arco

La forma esta dada por:

𝒃
𝑳 = ∫ √𝟏 + [𝒇′ (𝒙)]2 𝒅𝒙
𝒂
 2.1. Longitud de arco en coordenadas polares

Sea la curva en coordenadas polares definida por la función r = r(θ) donde θ∈[α, β] . La
región D encerrada por la curva cuando θ∈[α, β] es el área que queda encerrada entre las
semirrectas θ = α y θ = β, y la propia curva. Esto es, el conjunto de puntos en coordenadas
polares (Perez, s.f.)

Figura 2. Región encerrada en coordenadas polares

La fórmula está dada por:

𝒃
𝑳 = ∫ √[𝒓(𝜽)]𝟐 + [𝒓´(𝜽)]𝟐 𝒅𝜽
𝒂

2.2. Longitud de arco en coordenadas paramétricas

Si una curva suave S está dada por x=f(t) y y=g(t) y C no se corta a sí misma en el intervalo 𝑎 ≥
𝑡 ≥ 𝑏 (excepto quizá en los puntos terminales), entonces la longitud de arco de S en ese
intervalo está dada por:

𝒃
𝑳 = ∫ √[𝒇´(𝒕)]𝟐 + [𝒈´(𝜽)]𝟐 𝒅𝒕
𝒂

3. Centro de Masa
La región plana se va a tomar como una lámina bidimensional de densidad P(x) en
una región tiene un ejes de simetría, el centro de masa (si la densidad es uniforme )
estará sobre el o los ejes de simetría: Así un circulo tendrá su centro de masa en el
centro que es el punto de intersección de los diámetros, un rectángulo en el punto de
corte de sus diagonales, o en el punto de intersección de las rectas que bisecan sus
lados. (Scribd, 2011)
Por lo tanto, las fórmulas nos quedan:

Masa:

𝐛
𝐦 = ∫ 𝜹(𝐱)𝐟(𝐱)𝐝𝐱
𝐚

Momento en x, y

𝟏 𝒃
𝐌𝐱 = ∫ 𝜹(𝒙)𝒇(𝒙)𝟐 𝒅𝒙
𝟐 𝒂

𝒃
𝐌𝐲 ∫ 𝐱𝜹(𝒙)𝒇(𝒙)𝒅𝒙
𝒂

Por lo tanto centro de masa en x,y son

𝐌𝐲 𝑴𝐱
𝐱̅ = y 𝐲̅ =
𝐦 𝒎
 ANALISIS DE RESULTADOS

CODIGOS DE PROGRAMACIÓN

clc;
clear;
disp(' UNIVERSIDAD DE LAS FUERZAS ARMADAS ESPE' )
disp(' CALCULO VECTORIAL - UNIDAD I')
disp('INTEGRANTES:' )
disp('*RICHARD PILA' )
disp('*HUGO AUCAPIÑA' )
%disp('*____________' )
pause
clc;
disp(' -----SELECCIONE LA OPCION A
REALIZAR------')
disp('LONGITUD DE ARCO' )
disp('1) Cordenadas Rectangulares' )
disp('2) Cordenadas Polares' )
disp('3) Cordenadas Parametricas' )
disp(' ')
disp('SUPERFICIE DE REVOLUCION' )
disp('4) Respecto al Eje X' )
disp('5) Respecto al Eje Y' )
disp(' ')
disp('CENTROS DE MASA Y MOMENTOS' )
disp('6) En una Recta' )
disp('7) En una Superfice' )
disp('8) Entre dos funciones' )
disp(' ' )
opcion = input('Seleccione la opcion: ');
switch opcion
case 1
clc;
disp('COORDENADAS RECTANGULARES')
disp(' ')
syms x
y(x)=input('Ingrese la funcion y(x): ');
a=input('Ingrese el limite inferior a: ' );
b=input('Ingrese el limite inferior b: ' );
dy=diff(y);
l=int(sqrt(1+(dy^2)),a,b);
Longitud=vpa(l,5)
ezplot(y,[a,b])
title ('Grafica de la funcion F(x)');
ylabel('f(x)');
xlabel('x');
grid on
case 2
clc;
disp('COORDENADAS POLARES')
disp(' ')
syms theta
r(theta)=input('Ingrese la funcion r(theta): ' );
theta1=input('Ingrese el limite inferior theta 1: ' );
theta2=input('Ingrese el limite inferior theta 2: ' );
dr=diff(r);
l=int(sqrt(r^2+(dr^2)),theta1,theta2);
Longitud=vpa(l,5)
thetag=[theta1:0.1:theta2];
rr(theta)=r(thetag);
polarplot(thetag,rr,'r')
title ('Grafica de la funcion r(theta)');
grid on
 case 3
clc;
disp('COORDENADAS PARAMETRICAS')
disp(' ')
syms t
x=input('Ingrese la funcion x(t): ' );
y=input('Ingrese la funcion y(t): ' );
a=input('Ingrese el limite inferior a: ' );
b=input('Ingrese el limite inferior b: ' );
disp('x(t)')
pretty(x)
disp('y(t)')
pretty(y)
dx=diff(x);
dy=diff(y);
l=int(sqrt((dx^2)+(dy^2)),a,b);
Longitud=vpa(l,5)
ezplot(x,[a,b])
hold on;
ezplot(y,[a,b])
grid on
title ('Grafica de la funcion x(t) - y(t)');
case 4
clc;
disp(' EJE X')
disp(' ')
syms x
F(x)=input('Ingrese la funcion F(x): ' );
a=input('Ingrese el limite inferior a: ' );
b=input('Ingrese el limite superior b: ' );
dFx=diff(F(x));
funcion=F(x)*sqrt(1+(dFx^2));
A=2*pi*int(funcion,x,a,b);
AREA=vpa(A,5)
ezplot(F(x),[a,b])
figure(1)
title ('Grafica de la funcion F(x)');
ylabel('f(x)');
xlabel('x');
grid on
x=a:.1:b;
Fx=input('Ingrese la funcion F(x): ' );
figure(1)
[x,y,z]=cylinder(Fx);
surf(z,x,y),axis vis3d
title('Superficie')
case 5
clc;
disp(' EJE Y ')
disp(' ')
syms y
G(y)=input('Ingrese la funcion G(y): ' );
a=input('Ingrese el limite inferior a: ' );
b=input('Ingrese el limite superior b: ' );
dGy=diff(G(y));
funcion=G(y)*sqrt(1+(dGy^2));
A=2*pi*int(funcion,y,a,b);
AREA=vpa(A,5)
y = linspace(a,b);
x = G(y);
plot(x,y)
figure(1)
title ('Grafica de la funcion G(y)');
ylabel('y');
xlabel('G(y)');
grid on
y=a:.1:b;
Gy=input('Ingrese la funcion G(y): ' );
figure(1)
[x,y,z]=cylinder(Gy);
surf(x,y,z),axis vis3d
title('Superficie')
case 6
clc;
disp(' EN UNA RECTA')
disp(' ')
syms x d
dx=input('Ingrese la funcion densidad d(x): ' );
a=input('Ingrese el limite inicial: ' );
L=input('Ingrese la longitud L: ' );

mt=int(dx,x,a,L);
MasaTotal=vpa(mt,5)

Mt=int(dx*x,x,a,L);
Momento=vpa(Mt,5)

xt=Mt/mt;
Centro=vpa(xt,5)

X=[a L];
Y=[a a];

line(X,Y,'LineWidth',5)
hold on
plot(X,Y,'ro','LineWidth',5)
grid on
title ('Grafica de la funcion F(x)');
ylabel('f(x)');
xlabel('x');
case 7
clc;
disp(' EN UNA SUPERFICIE')
disp(' ')
%dx=1; Fx=sqrt(4*x); a=1; b=4;

syms x d
dx=input('Ingrese la funcion densidad d(x): ' );
Fx=input('Ingrese la funcion F(x): ' );
a=input('Ingrese el limite inferior a: ' );
b=input('Ingrese el limite superior b: ' );

m=int(dx*Fx,x,a,b);
Masa=vpa(m,5)

mx=(1/2)*int((dx*(Fx^2)),x,a,b);
MomentoX=vpa(mx,5)

my=int(x*dx*Fx,x,a,b);
MomentoY=vpa(my,5)
xc=my/m;
yc=mx/m;
Centro=[xc,yc]

ezplot(Fx,[a,b])
hold on
grid on
title ('Grafica de la funcion F(x)');
ylabel('f(x)');
xlabel('x');
case 8
clc;
disp(' ENTRE DOS FUNCIONES')
disp(' ')
syms x d
%dx=1;
%Fx=sqrt(3*x)+1;
%Gx=x+1;
%a=0;
%b=3;
dx=input('Ingrese la funcion densidad d(x): ');
Fx=input('Ingrese la funcion F(x): ');
Gx=input('Ingrese la funcion G(x): ');
a=input('Ingrese el limite inferior a: ' );
b=input('Ingrese el limite superior b: ' );
Ejemplo

Determine el centroide de la region del primer cuadrante limitada por la curva y 2 =

4x el eje de las x y las tectas x=1 , x=4

Donde

y = √4x

Primero escogemos en el menú la opción 7

Ponemos la funcion que desemos

Obtenemos el resultado y la grafica correspondiente
CONCLUSIONES

 Los software o plataformas digitales son muy importantes para la

comprobación de trabajos y estudios realizados en clase
 Ayuda a comprender el comportamiento de las funciones ingresadas ya sean
en coordenadas rectangulares, polares o paramétricas
 Podremos comprobar los cálculos matemáticos con el software para un óptimo
desempeño académico

RECOMENDACIONES

 Al momento de realizar la programación respectiva tener en cuenta todos los

comandos para no tener errores al momento de compilar el programa
 Verificar todos los posibles errores al momento de ingresar una función para
no tener problemas al momento de presentar el trabajo

BIBLIOGRAFÍA

Carlos, J. (20 de Febrero de 2018). Matlab. Obtenido de

https://juancarlosusomatlab2015.weebly.com/definicion-matlab.html
Daniel, V. (18 de Septiembre de 2017). Moodle_ITTUX. Obtenido de
http://moodle.ittux.edu.mx/course/info.php?id=5
Perez, R. (s.f.). Calculo II. Obtenido de
https://sites.google.com/site/calculovrubiop/temas/solidos-de-
revolucion/longitud-de-arco
Scribd. (03 de Julio de 2011). Obtenido de Sebastian Murillo:
https://es.scribd.com/doc/61658738/CENTROS-de-MASA-Ejercicios-Resueltos