Scribd
Cargar
173 vistas5 páginas

Métodos Numéricos en Ingeniería Electrónica

Este documento contiene 5 problemas de métodos numéricos resueltos con Matlab. El primer problema integra numéricamente las ecuaciones de una reacción química en un reactor batch no isotérmico usando el método de Euler. El segundo problema resuelve la ecuación de difusión unidimensional usando diferencias finitas. El tercer problema determina los coeficientes de un spline cúbico que interpola 3 puntos dados. El cuarto problema calcula el volumen requerido para rellenar un pantano después de drenarlo usando integrales numéricas

Cargado por

Quino B. Jeffry
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
173 vistas5 páginas

Métodos Numéricos en Ingeniería Electrónica

Este documento contiene 5 problemas de métodos numéricos resueltos con Matlab. El primer problema integra numéricamente las ecuaciones de una reacción química en un reactor batch no isotérmico usando el método de Euler. El segundo problema resuelve la ecuación de difusión unidimensional usando diferencias finitas. El tercer problema determina los coeficientes de un spline cúbico que interpola 3 puntos dados. El cuarto problema calcula el volumen requerido para rellenar un pantano después de drenarlo usando integrales numéricas

Cargado por

Quino B. Jeffry
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

METODOS NUMERICOS

CICLO 2019-A

TEMA: EXAMEN FINAL


ESTUDIANTES: QUINO BRICEÑO JEFFRY.

CODIGO: 1713210041

ESCUELA PROFESIONAL: INGENIERIA ELECTRONICA.


GRUPO HORARIO: 01L

LIMA-CALLAO – 2019
1. Un reactor Batch no isotérmico se puede describir mediante las siguientes
ecuaciones.
dC
 e( 10/(T  273))C
dt
dT
 1000e( 10/(T  273))C  10(T  20)
dt
Donde C es la concentración del reactante y T es la temperatura del reactor. Inicialmente,
el reactor se encuentra a 25 °C y tiene una concentración de reactante de C de 1.0 gmol/L.
Encuentre la concentración y temperatura del reactor como función del tiempo.
clc
clear
disp('INTEGRACION POR METODO DE EULER')
%Concentracion del reactante
f = @(t,C)(-C*exp(-10*t*(T+273)));
%Temperaura del reactor
g = @(t,T)(1000*C*exp(-10*t*(T+273))- 10*(T-20));
%se inicializa la suma----------------------------
h = 0.02;
N = t/h;
C(1)=1;
T(1)=25;
t=0:h:2;
%para la funcion f__________
for i=1:N
dC=-C(i)*exp(-10/(T(i)+273));
C(i+1)= C(i)+h*dC;
%para la funcion g__________
dT=1000*C(i)*exp(-10/(T(i)+273))- 10*(T(i)-20);
T(i+1)= T(i)+h*dT;
end
plot (t,C,'r-')
hold on
plot(t,T,'b.')
grid on
legend('Concentracion','Temperatura')

2. Las ecuaciones parabólicas se conocen como la ecuación de difusión y expresa:


𝝏𝟐 𝑻 𝝏𝑻
𝒌 𝝏𝒙𝟐 = 𝝏𝒕
Utilizando diferencias hacia adelante para la parte temporal y diferencias
centradas para ecuación que permite encontrar todos los nodos inferiores de una
barra unidimensional.
K= cte que depende del fenómeno

∂2 T ∂T
k 2= … … … … … … . . (1)
∂x ∂t
𝑗+1 𝑗
𝜕𝑇 𝑇𝑖 − 𝑇𝑖
≅ … … … … … (2)
𝜕𝑡 𝑚
𝑗 𝑗 𝑗
𝜕𝑇 2 𝑇𝑖+1 − 2 ∗ 𝑇𝑖 + 𝑇𝑖−1
≅ … … … … … (3)
𝜕𝑥 2 ℎ2
Reemplazando (2) y (3) en (1)
𝑗+1 𝑗 𝑗 𝑗 𝑗
𝑇𝑖 − 𝑇𝑖 𝑇 − 2 ∗ 𝑇𝑖 + 𝑇𝑖−1
= 𝑘 ∗ ( 𝑖+1 )
𝑚 ℎ2
𝑗 𝑗 𝑗
𝑗+1 𝑗 𝑚𝑘 𝑇 − 2 ∗ 𝑇𝑖 + 𝑇𝑖−1
𝑇𝑖 − 𝑇𝑖 = ( 2 ) ∗ ( 𝑖+1 )
ℎ ℎ2
𝑚
Donde 𝜆 = ( 2)𝑘
ℎ
𝑗+1 𝑗 𝑗 𝑗
𝑇𝑖 = (1 − 2𝜆) ∗ 𝑇𝑖 + 𝜆 ∗ (𝑇𝑖+1 + 𝑇𝑖−1
3. Un Spline Cubico S es definido por:
𝑆0 (𝑥) = 1 + 𝐵(𝑥 − 1) − 𝐷(𝑥 − 1)3 𝑠𝑖 1 ≤ 𝑥 ≤ 2
𝑆(𝑥) = { 3
𝑆1 (𝑥) = 1 + 𝑏(𝑥 − 2) − (𝑥 − 2)2 + 𝑑(𝑥 − 2)3 𝑠𝑖 2 ≤ 𝑥 ≤ 3
4
Si S es la interpolación de datos (1,1), (2,1) y (3,0), encuentra B, D, b y d.
*propiedad para determinar los valores:

𝑆0 (𝑥1 ) = 𝑆0 (2) = 1 + 𝐵 − 𝐷
𝑆1 (𝑥1 ) = 𝑆1 (2) = 1
𝑆1 (2) = 𝑆0 (2) 𝐵=𝐷
* Propiedad para determinar los valores:

𝑆1 ′(𝑥1 ) = 𝑆1 ′(2) = 𝑏
𝑆0 ′(𝑥1 ) = 𝑆0 ′(2) = 𝐵 − 3𝐷
𝑆1 ′(2) = 𝑆0 ′(2) 𝑏 = −2𝐷
* Propiedad para determinar los valores:
6
𝑆1 ′′(𝑥1 ) = 𝑆1 ′′(2) = −
4
𝑆0 ′′(𝑥1 ) = 𝑆0 ′′(2) = −6𝐷
1
𝑆1 ′′(2) = 𝑆0 ′′(2) 𝐷=
4
* Propiedad para determinar los valores:

𝑆1 ′′(𝑥𝑁 ) = 𝑆1 ′′(3) = 0
6
− + 6𝑑 = 0
4
Hallamos los valores:
1
𝐷=
4
1
𝐵=
4
1
𝑏=−
2
1
d=
4
4. Una ciudad desea drenar y rellenar el pequeño pantano mostrado en la figura. El
pantano tiene una profundidad media de 5m ¿Cuántos metros cúbicos de basura
se requerirán para llenar esa ara después de drenar el pantano?

Integral del trapecio

h N 1

IT  
2
f ( a )  f (b )  2 
i 1
f ( x) 

20
IT  146  14  2(14  30  40  54  76  122)  8320
2
Integral de Simpson

 N N
1 
h
f (a)  f (b)  4 f ( x2i 1 )  2  f ( x2i ) 
2 2
IS 
3  i 1 i 1

 

20
IS  146  14  4(122  54  30)  2(76  40  14)  8293.3
3
5. Para realizar un ajuste lineal de los datos: Como realizaría un script en Matlab
11 12 15 17 19
8 7 10 12 12

clear all
x=[11,12,15,17,19];
y=[8,7,10,12,12];
c=polyfit(x,y,1)
plot(x,y,'ro')

También podría gustarte