INGENIERIA MECANICA AUTOMOTRIZ

MATEMATICAS AVANZADAS

TEMA: CAPITULO 10 y 11

TRABAJO EXTRACLASE

Docente:
Ing. Juan Pablo Montero

Integrantes:
Juan Quizhpe
Pedro Troya

Ciclo:
5to - Grupo 4

Fecha:
27 de Enero del 2016
 EJERCICIOS DE LA SECCION 10.9
EJERCICIO 1
CODIFICACION
clc
clear all
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 10.9
%Integral de Fourier
%Ejercicio 1
format short
%variables
syms u w b x

display('Aplique la integral de fourier y demostrar que')

display('int((cosxw+wsenxw)/(1+w^2)),0,inf,w')
t=linspace(-2*pi,2*pi,1000);
f=(t<=0).*(0)+(t>0).*(pi*exp(-t));
plot(t,f)
grid on

display('Aw')
A1=int(pi*exp(-u)*cos(w*u),u,0,b);
A=limit(A1,b,inf);
pretty(A)
display('Bw')
B1=int(pi*exp(-u)*sin(w*u),u,0,b);
B=limit(B1,b,inf);
pretty(B)

display('se forma la representacion')

f=pi/(1+w^2)*cos(w*x)+(pi*w)/(1+w^2)*sin(w*x);
f1=1/(1+w^2)*cos(w*x)+(w)/(1+w^2)*sin(w*x);
simplify(f1);
pretty(f1)
RESULTADOS
EJERCICIO 3

CODIFICACION
clc
clear all
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 10.9
%Integral de Fourier
%Ejercicio 3
format short
%variables
syms u w b x

display('Aplique la integral de fourier y demuestre que')

display('int((sin(pi*w)+sen(x*w))/(1-w^2)),0,inf,w')
t=linspace(-2*pi,2*pi,1000);
f=((t<0).*(0)+((t>=0)&(t<=1)).*(1)+(t>=1).*(0));
plot(t,f)
grid on
display('Bw')
B1=2*int(pi*sin(x)*sin(w*x),x,0,pi);
pretty(B1)
f=(2*sin(pi*w))/(1-w^2);
pretty(f)

RESULTADOS
EJERCICIO 8

CODIFICACION
clc
clear all
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 10.9
%Integral de Fourier
%Ejercicio 8
format short
syms u w b x
display('1 cuando 0<x<1')
display('0 cuando x>1')
t=linspace(-2*pi,2*pi,1000);
f=((t<0).*(0)+((t>=0)&(t<=1)).*(1)+(t>=1).*(0));
plot(t,f)
grid on
display('Aw')
A1=2*int(cos(w*u),u,0,1);
pretty(A1)

display('representacion')
f=A1*cos(w*x)*1/pi;
pretty(f)

RESULTADOS
EJERCICIO 10

CODIFICACION
clc
clear all
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 10.9
%Integral de Fourier
%Ejercicio 10
format short

syms u w b x

display('x cuando 0<x<1')

display('2-x cuando 1<x<2')
display('0 cuando x>2')
t=linspace(-2*pi,2*pi,1000);
f=(((t>=0)&(t<=1)).*(t)+((t>=1)&(t<=2)).*(2-t)+(t>=2).*(0));
plot(t,f)
grid on
display('Aw')
A1=2*(int(x*cos(w*x),x,0,1)+int((2-x)*cos(w*x),x,1,2));
pretty(A1)
display('representacion')
f=A1*cos(w*x)*1/pi;
pretty(f)

RESULTADOS
EJERCICIO 11

CODIFICACION
clc
clear all
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 10.9
%Integral de Fourier
%Ejercicio 11
format short
syms u w b x a
display('x cuando 0<x<a')
display('0 cuando x>a')
t=linspace(-3,3,1000);
f=(((t>=0)&(t<=1)).*(t)+(t>=1).*(0));
plot(t,f)
grid on

display('Aw')
A1=2*(int(x*cos(w*x),x,0,a));
pretty(A1)

display('representacion')
f=A1*cos(w*x)*1/pi;
pretty(f)
RESULTADOS
EJERCICIO 13

CODIFICACION
clc
clear all
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 10.9
%Integral de Fourier
%Ejercicio 13
format short
syms u w b x a

display('x cuando 0<x<a')

display('0 cuando x>a')
t=linspace(-3,3,1000);
f=((t>=0).*(0)+(t>=0).*(exp(-t)+exp(-2*t)));
plot(t,f)
grid on
display('Aw')
A1=2*(int(((exp(-x)+exp(-2*x)*cos(w*x))),x,0,b));
A=limit(A1,b,inf);
pretty(A)
display('funcion')
f=(2/pi)*((3*(w^2+2))/((w^2+1)*(w^2+4)));
pretty(f)
RESULTADOS
 EJERCICIOS DE LA SECCION 11.1
EJERCICIO 6
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 6
syms x y
disp('Funcin')
u = x^3 - 3*x*y^2
disp('2derivada respecto a "x"')

dux2 = diff(u,'x',2)
disp('2da derivada respecto a "y"')
duy2 = diff(u,'y',2)
R = dux2 + duy2

RESULTADOS
EJERCICIO 7
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 7
syms x y

disp('Funcin')
u = sin(x) * sinh (y)
disp('2derivada respecto a "x"')
dux2 = diff(u,'x',2)
disp('2derivada respecto a "y"')
duy2 = diff(u,'y',2)
R = dux2 + duy2
RESULTADOS
EJERCICIO 8
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 8
syms x y
disp('Funcin')
u = x^3 - 3*x*y^2
disp('2derivada respecto a "x"')
dux2 = diff(u,'x',2)
disp('2derivada respecto a "y"')
duy2 = diff(u,'y',2)
R = dux2 + duy2

RESULTADOS
EJERCICIO 13
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 13
syms x t
disp('Funcin')
u = x^3 + 3*x*t^2
disp('2derivada respecto a "x"')
dux2 = diff(u,'x',2)
disp('2derivada respecto a "t"')
dut2 = diff(u,'t',2)
R = (dut2 / dux2)^(1/2)

RESULTADOS
EJERCICIO 14
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 14
syms x t w c
disp('Funcin')
u = sin(w*c*t) * sin(w*x)
disp('2derivada respecto a "x"')
dux2 = diff(u,'x',2)

disp('2derivada respecto a "t"')

dut2 = diff(u,'t',2)
r = (dut2 / dux2)

RESULTADOS
EJERCICIO 15
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 15
syms x y z
disp('Funcin')
u = 1 / (x^2 + y^2 + z^2)^(1/2)
disp('2derivada respecto a "x"')
dx2 = diff(u,'x',2)
disp('2derivada respecto a "y"')
dy2 = diff(u,'y',2)
disp('2derivada respecto a "z"')
dz2 = diff(u,'z',2)
disp('dx2 + dy2 + dz2 = 0')
r = dx2 + dy2 + dz2
r = factor(r)
RESULTADOS
EJERCICIO 20
CODIFICACION
clear all
clc
disp('MATEMATICA AVANZADA ')
disp('ALUMNOS: JUAN QUIZHPE - PEDRO TROYA')
disp('SECCION 11.1_EJERCICIO 20')

disp('du/dx = 2*x*y*u')
%Cambiando "x" por "t"
u=dsolve('Du=2*t*y*u'); %resolucin de la ec. diferencial
syms x t
u=subs(u,{t},{x})
EJERCICIO 24
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.1
%Ejercicio 24

disp('((d(du/dx))/dy) + du/dx + x + y + 1 = 0')

r=dsolve('Dr = - r - t - y - 1');
u=dsolve('Du = C3*exp(-t) - y - t')
syms x t
u=subs(u,{t},{x})
RESULTADOS
 EJERCICIOS DE LA SECCION 11.5
EJERCICIO 1
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.5
%Flujo Unidimensional de Calor
%Ejercicio 1
bn=1
c=1
l=pi
disp('Clculo de lambda')
syms n
lambda=(c*n*pi) / l
syms t x
disp('Funcin u(x,t)')
u = bn * sin( (n*pi*x) / l) * exp(-( (lambda)^2 )*t) %funcin u(x,t)
disp('Con n = 1')
u1 = subs(u,{n},{1})
disp('u1 con t = 0')
u10 = subs(u1,{t},{0}) %t = 0
subplot(3,1,1)
ezplot(u10,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on

disp('u1 con t = 1')

u11 = subs(u1,{t},{1}) %t = 1
ezplot(u11,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp(' ')
disp('u1 con t = 2')
u12 = subs(u1,{t},{2}) %t = 2
ezplot(u12,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u1 con t = 3')
u13 = subs(u1,{t},{3}) %t = 3
ezplot(u13,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
title('u1(x,t)')
u2 = subs(u,{n},{2}) %n = 2
disp('u2 con t = 0')
u20 = subs(u2,{t},{0}) %t = 0
subplot(3,1,2)
ezplot(u20,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u2 con t = 1')
u21 = subs(u2,{t},{1}) %t = 1
subplot(3,1,2)
ezplot(u21,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u2 con t = 2')
u22 = subs(u2,{t},{2}) %t = 2
subplot(3,1,2)
ezplot(u22,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u2 con t = 3')
u23 = subs(u2,{t},{3}) %t = 3
subplot(3,1,2)
ezplot(u23,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
title('u2(x,t)')
u3 = subs(u,{n},{3}) %n = 3
disp('u3 con t = 0')
u30 = subs(u3,{t},{0}) %t = 0
subplot(3,1,3)
ezplot(u30,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u3 con t = 1')
u31 = subs(u3,{t},{1}) %t = 1

subplot(3,1,3)
ezplot(u31,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u3 con t = 2')
u32 = subs(u3,{t},{2}) %t = 2
subplot(3,1,3)
ezplot(u32,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
disp('u3 con t = 3')
u33 = subs(u3,{t},{3}) %t = 3
subplot(3,1,3)
ezplot(u33,[0,pi])
grid on
axis([0 pi 0 1.1])
hold on
title('u3(x,t)')
RESULTADOS
GRAFICA
EJERCICIO 2
CODIFICACION
clear all
clc
%Trabajo realizado por Juan Quizhpe Pedro Troya
%Ejercicios de la seccion 11.5
%Flujo Unidimensional de Calor
%Ejercicio 2
x=0:.1:pi
t=0:.1:3
[X,T]=meshgrid(x,t)
u= sin(X).*exp(-T)
subplot(2,2,1)
mesh(u)
grid on
axis on

x=0:.1:pi;
t=0:.1:3;
[X,T]=meshgrid(x,t);
u= sin(2*X).*exp(-4*T);
subplot(2,2,2)
mesh(u)
grid on
axis on

x=0:.1:pi;
t=0:.1:3;
[X,T]=meshgrid(x,t);
u= sin(3*X).*exp(-9*T);
subplot(2,2,3)
mesh(u)
grid on
axis on
GRFICA
EJERCICIO 3

Comprobar las figuras y describir, el comportamiento diferente de las soluciones de las 2

ecuaciones

Podemos describir que las soluciones de la ecuacin de onda son peridicas en T

EJERCICIO 4

De qu manera depende la rapidez del decremento de la ecuacin para n fijo, del calor
especifico, la densidad y la conductividad trmica del material

CODIFICACION

% Matematica Avanzada
% Realizado por: Juan Quizhpe Pedro Troya
% Seccion 11.5;
% Ejercicio 4
syms Fn(x) Gn(t)
Un(x,t)=Fn(x)*Gn(t)

Un(x, t) = Gn(t)*Fn(x)

syms Bn n x L Yn
Gn(t)Fn(x)=Bn*sin(n*pi*x/L)*exp(t*(-Yn^2))

% K=Conductividad termica
% d= Calor especifico del cuerpo
% p= Densidad del material del cuerpo
% La rapidez del decremento depende del valor de c que es una
% constante que depende de cada uno de los valores de Yn y tambin de
la longitud L y considerando tambin si el espesor de la varilla es
gruesa o delgada, con seccin transversal constante, material homogneo
y perfectamente aislado
EJERCICIO 7
CODIFICACION
function graf
clear all
clc
x=0:0.005:10;
%definir la funcin
f=(((x>0)&(x<=5)).*(x) + ((x>5)&(x<10)).*(10-x));
%Grfica funcin inicial
subplot(2,1,1)
plot(x,f,'r')
grid on
title('Temperatura inicial')
axis([0 10 0 6])
clear all
clc
disp('INTEGRANTES')
disp('Juan Quizhpe Pedro Troya')
disp('Seccion 11.5')
disp('Ejercicio 7')

disp('Clculo de "c"')
disp('Conductividad trmica')
K=1.04
disp('Calor especfico del cuerpo')
s=0.056
disp('Densidad')
p=10.6
c=(K/(s*p))^(1/2)

disp('Clculo de "lambda"')
syms k;
l=10;
lambda=(c*k)/l
disp('Clculo de coeficiente')
disp(' ak')
syms x %declarar variable
f1=x*sin(k*pi*x/l);
i1=0; %lmites
s1=5;
f2=(10-x)*sin(k*pi*x/l);
i2=5;
s2=10;
ak=(2/l) * ((int(f1,'x',i1,s1)) + (int(f2,'x',i2,s2)) ) %clculo de ak

disp('Haciendo t = 5 para la grfica')

t=5
u=(ak * (sin( (k*pi*x)/l )) * exp(-( (lambda)^2 )*t))
u1=subs(u,{k},{1});
u3=subs(u,{k},{3});
u5=subs(u,{k},{5});
u7=subs(u,{k},{7});
s4= u1 + u3 + u5 + u7

subplot(2,1,2)
ezplot(s4,[0,10])
title('u(x,5)')
grid on
axis([0 10 0 5])
Clculo de "lambda"

lambda =(26^(1/2)*371^(1/2)*k)/742

Clculo de coeficiente ak
ak =(100*sin((pi*k)/2) - 200*cos((pi*k)/2)*sin((pi*k)/2) +
50*pi*k*cos((pi*k)/2))/(5*pi^2*k^2) + (2*(10*sin((pi*k)/2) -
5*pi*k*cos((pi*k)/2)))/(pi^2*k^2)
Haciendo t = 5 para la grfica
t =5
u =sin((pi*k*x)/10)*exp(-(65*k^2)/742)*((100*sin((pi*k)/2) -
200*cos((pi*k)/2)*sin((pi*k)/2) + 50*pi*k*cos((pi*k)/2))/(5*pi^2*k^2)
+ (2*(10*sin((pi*k)/2) - 5*pi*k*cos((pi*k)/2)))/(pi^2*k^2))
s4 =(40*exp(-65/742)*sin((pi*x)/10))/pi^2 - (40*exp(-
455/106)*sin((7*pi*x)/10))/(49*pi^2) - (40*exp(-
585/742)*sin((3*pi*x)/10))/(9*pi^2) + (8*exp(-
1625/742)*sin((pi*x)/2))/(5*pi^2)

GRAFICA
EJERCICIO 9
CODIFICACION
function graf
clear all
clc
disp('INTEGRANTES')
disp('Juan Quizhpe Pedro Troya)
disp('Seccion 11.5')
disp('Ejercicio 9')
disp('Clculo de "c"')
disp('Conductividad trmica')
K=1.04
disp('Calor especfico del cuerpo')
s=0.056
disp('Densidad')
p=10.6
c=(K/(s*p))^(1/2)

disp('Clculo de "lambda"')
syms k;
l=10;
lambda=(c*k)/l

disp('Clculo de coeficiente')
disp(' ak')
syms x %declarar variable
f1= (x *(10-x)) * sin(k*pi*x/l);
i1=0; %lmites
s1=10;
ak=(2/l) * (int(f1,'x',i1,s1)) %clculo de ak

disp('Haciendo t = 5 para la grfica')

t=5
u=(ak * (sin( (k*pi*x)/l )) * exp(-( (lambda)^2 )*t))
u1=subs(u,{k},{1});
u3=subs(u,{k},{3});
u5=subs(u,{k},{5});
u7=subs(u,{k},{7});
s4= u1 + u3 + u5 + u7

u=x*(10-x);
subplot(2,1,1)
ezplot(u,[0,10])
title('u(x,t)')
grid on
axis([0 10 0 28])

subplot(2,1,2)
ezplot(s4,[0,10])
title('u(x,5)')
grid on
axis([0 10 0 28])
Clculo de "lambda"
lambda =(26^(1/2)*371^(1/2)*k)/742
Clculo de coeficiente ak
ak =(800*sin((pi*k)/2)^2)/(pi^3*k^3) - (200*sin(pi*k))/(pi^2*k^2)
Haciendo t = 5 para la grfica
t =5
u =-sin((pi*k*x)/10)*exp(-(65*k^2)/742)*((200*sin(pi*k))/(pi^2*k^2) -
(800*sin((pi*k)/2)^2)/(pi^3*k^3))
s4 =(800*exp(-455/106)*sin((7*pi*x)/10))/(343*pi^3) + (800*exp(-
65/742)*sin((pi*x)/10))/pi^3 + (800*exp(-
585/742)*sin((3*pi*x)/10))/(27*pi^3) + (32*exp(-
1625/742)*sin((pi*x)/2))/(5*pi^3)

GRAFICA