Laboratorio 09

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 10

Matemáticas Aplicadas a la

Electrónica

Laboratorio N°9
Ecuación Diferencial de Bernoulli
Alumno: CABALLERO CARHUARICRA
JOSUE
Ciclo: III Sección: A

Profesor: Godinez De La Cruz Ernesto


Fecha de realización: 30/04/19
Fecha de presentación: 06/05/19
2019-1
1
Objetivo

El presente laboratorio tiene como objetivo utilizar los comandos del MatLab para resolver una
ecuación diferencial ordinaria de Bernoulli.

Fundamento Teórico

La ecuación diferencial ordinaria de Bernoulli tiene la siguiente forma

dy
+ f ( t ) y = g (t ) y n
dt

Donde f ( t ) y g ( t ) son funciones de t o constantes y n es un numero real

Procedimiento

dy
1. Resolver la siguiente ecuación diferencial + 5 y = 4 y1/2
dt

a) Obtener la solución general

b) Obtener la solución particular considerando la condición inicial y(0) = 0

2
c) Graficar la solución particular de t=0 a t=5s

d) Calcular y(1,5)

3
e) Calcular aproximadamente el valor de t para y=0.4

2. Se tiene un circuito serie RC en el cual R=10Ω y C=0.02F. La tensión de la


fuente es una onda cuadrada de amplitud 12V y un periodo de 5s. Realizar
el código MatLab para graficar de 0 a 20s la tensión de la fuente, la tensión
y la corriente en el condensador.

CODIGO MATLAB:

clc
clear all
close all

%DATOS
R=10;
C=0.02;
A=12;
T=5;
f=1/T;
h=0.001;
t=0:h:4*T;

%SEÑAL CUADRADA
u= A*square(2*pi*f*t);

%NUMERO DE MUESTRA
n=length (u);
vc=zeros(1,n); %VOLTAJE CONDENSADOR
ic=zeros(1,n); %CORRIENTE CONDENSADOR
dvc=zeros(1,n); % DERIVADA DEL VOLTAJE

4
for k=1:n-1;
dvc(k)=(u(k)-vc(k))/(R*C);
vc(k+1)= vc(k)+h*dvc(k);
ic(k)= C*dvc(k);
end

subplot(311)
plot(t,u,'b','linewidth',3)
xlabel('t')
ylabel('u')
title('GRAFICA DE SOLUCION PARTICULAR')
grid on

subplot(312)
plot(t,vc,'r','linewidth',3)
xlabel('t')
ylabel('vc')
title('GRAFICA DE VOLTAJE')
grid on

subplot(313)
plot(t(1:n-1),ic(1:n-1),'g','linewidth',3)
xlabel('t')
ylabel('ic')
title('GRAFICA DE CORRIENTE')
grid on

5
3. Se tiene un circuito serie RL en el cual R=5Ω y L=0.05H. La tensión de la
fuente es una onda senoidal de amplitud 12V y un periodo de 5s. Realizar el
código MatLab para graficar de 0 a 20s la tensión de la fuente, la tensión y la
corriente de la bobina.

CODIGO MATLAB:
clc
clear all
close all

%DATOS
R=5;
L=0.05;
A=12;
T=5;

6
f=1/T;
h=0.001;
t=0:h:4*T;
%SEÑAL SENOIDAL
u= A*sin(2*pi*f*t);

%NUMERO DE MUESTRA
n=length (u);
vl=zeros(1,n); %VOLTAJE BOBINA
il=zeros(1,n); %CORRIENTE BOBINA
dil=zeros(1,n); %DERIVADA DEL VOLTAJE

for k=1:n-1;
dil(k)=(u(k)-R*il(k))/(L);
il(k+1)= il(k)+h*dil(k);
vl(k)=L*dil(k);
end

subplot(311)
plot(t,u,'b','linewidth',3)
xlabel('t')
ylabel('u')
title('GRAFICA DE SOLUCION PARTICULAR')
grid on

subplot(312)
plot(t,vl,'r','linewidth',3)
xlabel('t')
ylabel('vl')
title('GRAFICA DE VOLTAJE')
grid on

subplot(313)
plot(t(1:n-1),il(1:n-1),'g','linewidth',3)
xlabel('t')
ylabel('il')
title('GRAFICA DE CORRIENTE')
grid on

7
Aplicación

8
CODIGO MATLAB:
clear all
close all
syms t
%SOLUCION PARTICULAR
vp=dsolve('Dv+(30*v^2*(t+1)^3)/1500=0,v(0)=90');
disp('vp=')
disp(vp)
%GRAFICA VP
t=0:0.001:3;
v=1./((t + 1).^4/200 + 11/1800);
plot(t,v,'b','linewidth',3)
xlabel('x')
ylabel('v')
title('Gráfica de V')
grid on

9
10

También podría gustarte