FACULTAD DE INGENIERA

ESCUELA ACADMICO PROFESIONAL

DE INGENIERA MECNICA

TEMA:

Introduccin a Matlab

DOCENTE:
Ing. Guayan Huaccha, Eli

CURSO:

Computacin Bsica

ALUMNO:

Asencio Palma, Sthefano Deckert

TRUJILLO - PER

2017
 NDICE
CARATULA .................................................................................................................... 1
INDICE ........................................................................................................................... 2
DESARROLLO DE PREGUNTA 1 ................................................................................. 3
DESARROLLO DE PREGUNTA 2 ................................................................................. 8
DESARROLLO DE PREGUNTA 3 ............................................................................... 10
DESARROLLO DE PREGUNTA 4 ............................................................................... 15
DESARROLLO DE PREGUNTA 5 ............................................................................... 17
COMENTARIOS .......................................................................................................... 21
CONCLUSIONES DEL TEMA Y APRECIACION CRITICA.22

BIBLIOGRAFIA...22

2
1. PREGUNTA 1

La densidad del aire, D (resultado de calcular la media de distintas medidas), a diferentes

alturas, h desde el nivel del mar hasta 33 km, viene dada por la tabla que se muestra a
continuacin:

H(Km) 0 3 6 9 12 15
D(Kg/m3) 1.2 0.91 0.66 0.47 0.31 0.19
H(Km) 18 21 24 27 30 33
D(Kg/m3) 0.12 0.075 0.046 0.029 0.018 0.011

1.1 DESARROLLO

Tabulamos los datos segn las escalas que nos piden:

A B C D
Lineal Logartmica Lineal Logartmica
Lineal Lineal Logartmica Logartmica

Luego creamos polinomios que se ajusten a los puntos para cada caso, con el comando
polyfit.
Tomamos el grafico que ms se aproxime a los datos y vemos que tipo de funcin es,
luego calculamos sus coeficientes de acuerdo a la grfica que escojamos, en este si es
lineal:
= (1) = (2)

Y graficamos.

1.2 PROGRAMA EN MATLAB (Archivo M-file)

%Pregunta 1
%Inciso 1
h=[0 3 6 9 12 15 18 21 24 27 30 33];
D=[1.2 0.91 0.66 0.47 0.31 0.19 0.12 0.075 0.046 0.029 0.018 0.011];
f1=polyfit(h,D,4);
%Grafica 1
figure
hf=0:0.01:33;
Df=polyval(f1,hf);
plot(h,D,'o');
title('Polinomio grado 4, Escala lineal');
xlabel('altura(Km)');
ylabel('Densidad(Kg/m^3)');

%Iniciso 2
h=[ 3 6 9 12 15 18 21 24 27 30 33];
D=[ 0.91 0.66 0.47 0.31 0.19 0.12 0.075 0.046 0.029 0.018 0.011];
r=log(h);

3
f2=polyfit(r,D,4);
%Grafica 2
figure
rf=1:0.01:log(33);
Df=polyval(f2,rf);
plot(r,D,'o');
title('Polinomio grado 4, Escala logaritmica vs lineal');
xlabel('logaritmo altura(u)');
ylabel('Densidad(Kg/m^3)');

%Inciso 3
h=[0 3 6 9 12 15 18 21 24 27 30 33];
D=[1.2 0.91 0.66 0.47 0.31 0.19 0.12 0.075 0.046 0.029 0.018 0.011];
q=log(D);
f3=polyfit(h,q,1);
%Grafica 3
figure
hf=0:0.01:33;
qf=polyval(f3,hf);
plot(h,q,'o');
title('Polinomio grado 1, Escala logaritmica vs lineal');
xlabel('altura(Km)');
ylabel('log(Densidad)(u)');

%Inciso 4
h=[3 6 9 12 15 18 21 24 27 30 33];
D=[0.91 0.66 0.47 0.31 0.19 0.12 0.075 0.046 0.029 0.018 0.011];
r=log(h);
q=log(D);
f4=polyfit(r,q,4);
%Grafica 4
figure
rf=1:0.01:log(33);
qf=polyval(f4,rf);
plot(r,q,'o');
title('Polinomio grado 4, escala logaritmica')
xlabel('log(Altura) (u)');
ylabel('log(Densidad) (u)');

%Inciso 5

%La grafica 4 tine tendencia lineal y por lo tanto la ecuacion potencial

%se ajusta mejor a los datos
%Hallamos las costantes b y m
h=[0 3 6 9 12 15 18 21 24 27 30 33];
D=[1.2 0.91 0.66 0.47 0.31 0.19 0.12 0.075 0.046 0.029 0.018 0.011];
s=polyfit(h,log(D),1);
%m
m=s(1)
%b
b=exp(s(2))
%Grafica
figure
hm=0:0.1:33;

4
Dm=b*exp(m*hm);
plot(h,D,'o',hm,Dm)

Graficas:

5
1.3 EJECUCIN:

1.5 DIAGRAMA DE FLUJO

INICIO

h=[0 3 6 9 12 15 18 21 24 27 30 33]
D=[1.2 0.91 0.66 0.47 0.31 0.19 0.12 0.075 0.046 0.029 0.018 0.011]

plot(h,D,'o')

6
 plot(log(h),D,'o')

plot(h,log(D),'o')

plot(log(h),log(D),'o')

s=polyfit(h,log(D),1)

m=s(1) y b=exp(s(2))

hm=0:0.1:33
Dm=b*exp(m*hm)
plot(h,D,'o',hm,Dm
)

Imprimir: m b FIN

7
2. PREGUNTA 2

Un peso W cuelga de un anillo que a su ves est sujeto por dos cables unidos a dos
bisagras, tal y como se muestra en la figura. La bisagra del punto A se encuentra fija,
mientras que la bisagra del punto B se puede desplazar (sin friccin) en la direccin
horizontal. La fuerza F y F depende de la posicin de la bisagra B (distancia X); y se
puede calcular mediante las ecuaciones:

FacCOS (alpha) = FacCOS(beta)

FacSEN(alpha) + FacSEN(beta) = W

a) Utilice Matlab para obtener las expresiones de las fuerzas Fac y Fbc en funcin de
x, W y la longitud de los cables, Lac y Lbc.
b) Utilice el comando subs para sustitur W= 2000 N, Lac =0,3 m y Lbc =0,5 m en la
expresin obtenida en el apartado a. Esto proporcionar la fuerza en los cables en
funcin de la distancia x.
c) Utilice el comando ezplot para representar las fuerzas Fac y Fbc (ambas en el
mismo grfico) en funcin de x, empezando en 0,4 m y acabando cerca de los 0,8
m. Qu sucede cuando x se aproxima a 0,8 m?

2.2 PROGRAMA EN MATLAB (Archivo M-file)

%Pregunta 2
syms x W la lb
a=(la^2-lb^2+x^2)/(2*x) ;
b=(lb^2-la^2+x^2)/(2*x);
h =(la^2 - (la^2 - lb^2 + x^2)^2/(4*x^2))^(1/2);
Fa =(W*la)/((la^2 - (la^2 - lb^2 + x^2)^2/(4*x^2))^(1/2)*((la^2 - lb^2 + x^2)/(- la^2 + lb^2 + x^2) + 1));
Fb =(W*lb)/((la^2 - (la^2 - lb^2 + x^2)^2/(4*x^2))^(1/2)*((- la^2 + lb^2 + x^2)/(la^2 - lb^2 + x^2) + 1));
FA=subs(Fa,{W,la,lb},{200,0.3,0.5});
FB=subs(Fb,{W,la,lb},{200,0.3,0.5});
ezplot(FA,[0.4,0.8])
grid
hold on
ezplot(FB,[0.4,0.8])
xlabel('longitud de la visagra "x"')
ylabel('Fuerza de la cuerda A y B')

8
GRAFICA:

2.3 EJECUCION

2.4 DIAGRAMA DE FUJO

INICIO

syms x W la lb
a=(la^2-lb^2+x^2)/(2*x) ;
b=(lb^2-la^2+x^2)/(2*x);
h =(la^2 - (la^2 - lb^2 + x^2)^2/(4*x^2))^(1/2);

Fa =(W*la)/((la^2 - (a)^2/(4*x^2))^(1/2)*((b)/(b) + 1));

Fb =(W*lb)/((la^2 - (a)^2/(4*x^2))^(1/2)*((b)/(a) + 1));

FA=subs(Fa,{W,la,lb},{200,0.3,0.5})
FB=subs(Fb,{W,la,lb},{200,0.3,0.5})

ezplot(FA,[0.4,0.8])
ezplot(FB,[0.4,0.8])

FIN

3.1PREGUNTA 3

Considere el siguiente conjunto de siete ecuaciones:

3x1 + 4x2 + 2x3 x4 + x5 + 7x6 + x7 =42

2x1 2x2 + 3x3 x4 + 5x5 + 2x6 + 8x7 = 32

x1 + 2x2 + 3x3 + x4 + 2x5 + 4x6 + 6x7 = 12

5x1 + 10x2 + 4x3 + 3x4 + 9x5 2x6 + x7 = -5

3x1 + 2x2 2x3 4x4 5x5 6x6 + x7 = 10

-2x1 + 9x2 + x3 + 3x4 3x5 + 5x6 + x7 = 18

x1 2x2 8x3 + 4x4 + 2x5 + 4x6 + 5x7 = 17

Defina una variable simblica para cada una de las ecuaciones y use la capacidad
simblica de Matlab para resolver cada una de las incgnitas.

10
3.2 DESARROLLO

Hacemos simblicas todas las variables para poder trabajar

Luego denotamos cada ecuacin con un valor diferente a, b, c, etc.

A=3*x1+4*x2+2*x3-1*x4+1*x5+7*x6+1*x7-42;
B=2*x1-2*x2+3*x3-1*x4+5*x5+2*x6+8*x7-32;
C=1*x1+2*x2+3*x3+1*x4+2*x5+4*x6+6*x7-12;
D=5*x1+10*x2+4*x3+3*x4+9*x5-2*x6+1*x7+5;
E=3*x1+2*x2-2*x3-4*x4-5*x5-6*x6+7*x7-10;
F=-2*x1+9*x2+1*x3+3*x4-3*x5+5*x6+1*x7-18;
G=1*x1-2*x2-8*x3+4*x4+2*x5+4*x6+5*x7-17;

Con la ayuda del comando solve creamos un vector con las soluciones del sistema de
ecuaciones.

= (, , , , , , )

3.2 PROGRAMA EN MATLAB (Archivo function)

%Pregunta 3
syms x1 x2 x3 x4 x5 x6 x7
%Definimos las ecuaciones
A=3*x1+4*x2+2*x3-1*x4+1*x5+7*x6+1*x7-42;
B=2*x1-2*x2+3*x3-1*x4+5*x5+2*x6+8*x7-32;
C=1*x1+2*x2+3*x3+1*x4+2*x5+4*x6+6*x7-12;
D=5*x1+10*x2+4*x3+3*x4+9*x5-2*x6+1*x7+5;
E=3*x1+2*x2-2*x3-4*x4-5*x5-6*x6+7*x7-10;
F=-2*x1+9*x2+1*x3+3*x4-3*x5+5*x6+1*x7-18;
G=1*x1-2*x2-8*x3+4*x4+2*x5+4*x6+5*x7-17;
%Vector solucion
[x1 x2 x3 x4 x5 x6 x7]=solve(A,B,C,D,E,F,G);
%Resulatos con 5 cifras significativas
X1=vpa(x1,5)
X2=vpa(x2,5)
X3=vpa(x3,5)
X4=vpa(x4,5)
X5=vpa(x5,5)
X6=vpa(x6,5)
X7=vpa(x7,5)

11
3.3 EJECUCION

3.4 DIAGRAMA DE FLUJO

INICIO

syms x1 x2 x3 x4 x5 x6 x7

A=3*x1+4*x2+2*x3-1*x4+1*x5+7*x6+1*x7-42;
B=2*x1-2*x2+3*x3-1*x4+5*x5+2*x6+8*x7-32;
C=1*x1+2*x2+3*x3+1*x4+2*x5+4*x6+6*x7-12;
D=5*x1+10*x2+4*x3+3*x4+9*x5-2*x6+1*x7+5;
E=3*x1+2*x2-2*x3-4*x4-5*x5-6*x6+7*x7-10;
F=-2*x1+9*x2+1*x3+3*x4-3*x5+5*x6+1*x7-18;
G=1*x1-2*x2-8*x3+4*x4+2*x5+4*x6+5*x7-17;

[x1 x2 x3 x4 x5 x6 x7]=solve(A,B,C,D,E,F,G)

IMPRIMIR: X1 X2 X3 X4 X5 X6 X7

12
 FIN

4 .PREGUNTA 4

Muchos procesos fsicos se pueden modelar mediante una ecuacin exponencial. Por
ejemplo, las tasas de reaccin qumica dependen de una constante de: tasa de reaccin
que es funcin de la temperatura y la energa de activacin.

= 0

Donde:

R = constante universal de los gases: 8,314 KJ/Kmol.K

Q = energa de activacin, en KJ/Kmol

T = Temperatura, en K

K0 = constante cuyas unidades dependen de caractersticas de la reaccin. Una

posibilidad es s-1.

Un enfoque para encontrar los valores de K0 y Q a partir de datos experimentales es

graficar el logaritmo natural de k en el eje y y 1/T en el eje x. Esto debe resultar en una
lnea recta con pendiente Q/R e intercepta ln(Ke); esto es,

1
ln() = ln(0 ) ( )( )

Pues la ecuacin de la recta tiene la forma:

Y = ax + b

Con y = ln (K), x = 1/T, a = -Q/R y B = ln (k)

Ahora considere los datos siguientes:

T, K K, S 1
200 1.46 x 107
400 0.0012
600 0.0244
800 0.1099
1000 0.2710

13
4.1 DESARROLLO

Sabemos la constante R =8,314 KJ/Kmol.K

Creamos un vector con la primera columna de la tabla:

T=[200 400 600 800 1000 ]

Y otro de la segunda:
K=[1.46*10^(-7) 0.0012 0.0244 0.1099 0.2710]

Y como nos dice que x = 1/T y Y=ln(K)

Creamos una funcin polinomica que es una recta asi que de grado 1.
Donde m=s(1) y b= es(2)

Adems como:
= (0 ) = /

Tenemos:
0 = (2) =

4.2 PROGRAMA EN MATLAB (Archivo function)

%Pregunta 4
R=8.314;
T=[200 400 600 800 1000 ];
K=[1.46*10^(-7) 0.0012 0.0244 0.1099 0.2710];
%Aplicamos logaritmo natural
t=1./T;
k=log(K);
s=polyfit(t,k,1);
m=s(1);
K0=exp(s(2));
Q=-m*R;
fprintf('El valor de K0 es: %4.5f\n y el valor de Q es: %4.5f\n',K0,Q);
%GRAFICA
x=1/1000:0.00001:1/200;
y=polyval(s,x);
plot(x,y,'r');
xlabel('1/T (1/K)');
ylabel('log(K)');
title('Constante de tasa de reaccin')
text(1.1*10^(-3),2.3010,'Intercepto')

14
GRAFICA :

4.3 EJECUCION

4.4 DIAGRAMA DE FLUJO

INICIO

R=8.314;
T=[200 400 600 800 1000 ];
K=[1.46*10^(-7) 0.0012 0.0244 0.1099 0.2710];

t=1./T;
k=log(K
)

15


s=polyfit(t,k,1);

m=s(1);
K0=exp(s(2));
Q=-m*R;

x=1/1000:0.00001:1/200;
y=polyval(s,x);
plot(x,y,'r');

IMPRIMIR: Q K0

FIN

16
5. PREGUNTA 5

Hacer un programa en Matlab, para determinar el rea de la figura comprendida entre la

parbola y= - x2 +4x -3 y las tangentes a sta en los puntos (0;-3) y (3;0)

Por el mtodo de trapecios de Simpson, en 10 tramos. Tambin desde el programa

deber graficar las funciones dadas.

a) Escribir la codificacin del programa en un (M-File) en forma completa, para la

grfica dar un mnimo y mximo de 0,8 sobre el dominio de x. La curva de color
magenta y azul y con cuadrcula. La respuesta debe salir en m2.

5.1 DESARROLLO

Hallar las ecuaciones de la recta

Ecuacin de la parbola :
2 + 4 3 =
2 4 + 3 =

1
( 2)2 = 4( )( 1)
4

Entonces:
p=-1/4, h=2, k=1

Ecuacion de la recta tangente 1

4(( + 0)
(0 )( ) = )
2
1(3 + )
2( 2) = 1
2
4 3 =

Ecuacion de la recta tangente 2

4(( + 0)
(0 )( ) = )
2

1()
1( 2) = 1
2

2 + 6 =

17
5.2 PROGRAMA EN MATLAB (Fichero M-file)

%Pregunta 2
syms x y
%Hallar las ecuaciones de la recta
%Ecuacion de la parabola
% -x^2+4*x-3 = y
% x^2-4*x+3=-y
% (x^2-4*x+4)-1=-y
% (x-2)^2=-1*(y-1)
% (x-2)^2=4*(-1/4)*(y-1)
%Entonces :
%p=-1/4;
%h=2;
%k=1;
%Primer punto (0;-3)
%x0=0;y0=-3;
%Ecuacion de la recta tangente
%(x0-h)*(x-h)=(4*p*(((y+y0)/2)-k))
%-2*(x-2)=-1*((-3+y)/2)-1
%-2*x+4=((3-y)/2)+1
%-2*x+3=(3-y)/2
%-4*x+6=3-y
%-4*x+3=-y
%4*x-3=y
%Segundo punto (3;0)
%Ecuacion de la recta tangente
%(x0-h)*(x-h)=(4*p*(((y+y0)/2)-k))
% 1*(x-2)=-1*(y/2)-1
% x-2=(-y/2)+1
% x-2-1=-y/2
% 2*(x-3)=-y
% 2*x-6=-y
% y=-2*x+6
%Hallamos el limite superior
%4*x-3=-2*x+6
%6*x=9
%x=1.5
%y=3
%AREA DEBAJO DE LAS RECTAS
%Recta 1 y=4*x-3
%A1=b*h/2
%tomando nivel de referencia horizontal y=-3
%A1=6*1.5/2=4.5 m^2
%Recta 2 y=-2*x+6
%A2=(B+b/2)*h
%A2=(6+3/2)*1.5
%A2=(9/2)*1.5=6.75 m^2
f=input('Ingrese la funcion f a integrar: ','s');
g=input('Ingrese la funcion g a integrar: ','s');
e=input('Ingrese la funcion h a integrar: ','s');
a=input('Ingrese limite inferior: ');
b=input('Ingrese limite superior: ');
disp('Ingrese el numero de trapecios a')
n=input('considerar en la integracion: ');

18
xmin=a-0.6;xmax=b+0.6;
h=(b-a)/n;
%Dos rectas
%area de la recta1
base1=6;
altura1=1.5;
A1=base1*altura1/2;
%Area de la recta 2
Base2=6;
base2=3;
altura2=1.5;
A2=((base2+Base2)*altura2)/2;
%Area debajo de la curva
%Acomodando la ecuacion de la parbola como nivel de referencia y=-3
fr='-x.^2+4*x';
x=a:h:b;
fx=eval(fr);y=abs(fx);
A=y(1)+y(n+1);
B=2*sum(y(2:n));
integral=(h/2)*(A+B);
INTEGRAL=(A1+A2)-integral;
fprintf('El area es: %3.4f m^2\n',INTEGRAL);
%Grafica
xd=xmin:0.1:xmax;
x=xd;
y1=eval(f);
y2=eval(g);
y3=eval(e);
hold on
plot(xd,y1,'r',xd,y2,'g',xd,y3,'b');
grid on
xlabel('EJE X');
ylabel('EJE Y');
title('Grafica del area');

5.3 EJECUCION

19
5.4 DIAGRAMA DE FLUJO

INICIO

Leer: f, g, e, a, b, n

xmin=a-0.8;xmax=b+0.8
h=(b-a)/n

A1=base1*altura1/2
A2=((base2+Base2)*altura2
)/2

fr=-x.^2+4*x

x=a:h:b;
fx=eval(fr);y=abs(fx);
A=y(1)+y(n+1);
B=2*sum(y(2:n));
integral=(h/2)*(A+B);

INTEGRAL=(A1+A2)-integral;

Imprimir: INTEGRAL FIN

20
COMENTARIOS

PREGUNTA 1:
En esta situacin, utilizamos conocimientos sobre funciones y sus grficas.
Al saber esto pudimos aplicar lo que es el ajuste de curvas.

PREGUNTA 2:
En este ejercicio utilizamos conceptos tanto de fsica como de trigonometra.
Al llegar a las frmulas que obtenemos podemos ingresarlas en matlab y al
reemplazarla por cualquier valor obtener el resultado deseado.
PREGUNTA 3:
Por algebra lineal sabemos que si tenemos un sistema de ecuaciones con siete
incgnitas y siete variables siempre es posible encontrar una solucin para cada
variable por eso usando la capacidad de clculo hallamos fcilmente esta
soluciones, usando el comando solve.

PREGUNTA 4:
Como tenemos una funcin exponencial y queremos hallar unas constantes que
estn el exponente procedemos a linealizarla mediante el logaritmo natural.
Luego con el tema de ajuste de curvas hallamos los coeficientes.
PREGUNTA 5:
Tomamos como referencia del eje Y la recta Y=-3 para facilitar el clculo.
Hallamos las reas debajo de las rectas y luego el rea debajo de la curva
mediante el mtodo de Simpson luego restamos y hallamos el rea comprendida
entre ambas

21
 CONCLUSIONES Y APRECIACION CRTICA

Conforme pasa el tiempo hemos visto los cambios innovadores dentro de la Tecnologa y
la ciencia. En este segundo ciclo me toco aprender a manejar MATLAB, la cual considero
una potente herramienta que permite efectuar clculos y grficos, aunque suene
fastidioso los comando la practica garantiza una muy buena manera de aprender a la
perfeccin esta herramienta que sin duda me servir para mi futura carrera como
Ingeniero.

BIBLIOGRAFIA

Matlab y sus aplicaciones en las Ciencias y la Ingeniera. Cesar Prez.

Prentice Hall, Madrid, 2002.
Mtodos numricos Teora, probelmas y prcticas con MATLAB. Infante
del Ro J-A. & Rey Cabezas J. M. 2da Edicin - Pirmide. 2002.
Mtodos Numricos con Matlab. Mathews J.H., & Fink K.D. 3ra Edicin -
Prentice Hall 2000.

Intriduccin a Matlab. Sigmon, K. Department of Mathematics-University of

Florida

22
23