Intemerdio > Matematica aplicada > Ecuaciones diferenciales

EDO

Objetivo

Utilizar las funciones de MatLab para obtener el valor eficaz y medio de una señal periodica

Fundamento teórico

Software MatLab

Es una herramienta de software matemático que ofrece un entorno de desarrollo integrado (IDE) con un lenguaje de programación propio (lenguaje M).

Entre sus prestaciones básicas se hallan: la manipulación de matrices, la representación de datos y funciones, la implementación de algoritmos, la creación de interfaces de usuario y la comunicación con programas en otros lenguajes y con otros dispositivos hardware. El paquete MATLAB dispone de dos herramientas adicionales que expanden sus prestaciones, a saber, Simulink (plataforma de simulación multidominio) y GUIDE (editor de interfaces de usuario - GUI). Además, se pueden ampliar las capacidades de MATLAB con las cajas de herramientas (toolboxes); y las de Simulink con los paquetes de bloques.

Es un software muy usado en universidades y centros de investigación y desarrollo. En los últimos años ha aumentado el número de prestaciones, como la de programar directamente procesadores digitales de señal o crear código VHDL.

Valor medio

Se llama valor medio de una tensión (o corriente) alterna a la media aritmética de todos los valores instantáneos de tensión (o corriente), medidos en un cierto intervalo de tiempo.

En una corriente alterna sinusoidal, el valor medio durante un período es nulo: en efecto, los valores positivos se compensan con los negativos. Vm = 0

Periodo de la señal

Valor eficaz

Se define el valor eficaz de una corriente alterna, como aquel valor que llevado a corriente continua nos produce los mismos efectos caloríficos. Es un valor característico, que por otra parte es el que proporcionan los instrumentos de medida, ya sean analógicos o digitales. Aunque en la actualidad ya existen instrumentos digitales que proporcionan otros parámetros de la señal alterna.

Periodo de la señal:

Desarrollo

Ecuaciones diferenciales

Grafica, valor medio y valor eficaz
Código de MatLab:

Solución

Primero veamos la función
- %valor medio y eficaz de una señal trapecial-triangular
- clc
- close all
- clear all
- %señal trapecial-triangular
- T=14;
- fr=1/T;
- t=0:0.01:T;
- f1=2*t/3;
- f2=2;
- f3=t-3;
- f4=-t+11;
- f5=2;
- f6=(-2*t/3)+9.3;
- f7=0;
- ft=f1+(f2-f1).*heaviside(t-(3*T/14))+(f3-f2).*heaviside(t-(5*T/14))+(f4-f3).*heaviside(t-T/2)+(f5-f4).*heaviside(t-9*T/14)+(f6-f5).*heaviside(t-11*T/14)+(f7-f6).*heaviside(t-T);
- plot(t,ft,'LineWidth',4)
- axis([0 16 -1 7])
- grid on
- syms t
- %intervalo de integración
- d= [0 3*T/14 5*T/14 T/2 9*T/14 11*T/14 T];
- %función trapecial-triangular
- f1=2*t/3;
- f2=2;
- f3=t-3;
- f4=-t+11;
- f5=2;
- f6=(-2*t/3)+9.3;
- f=[f1 f2 f3 f4 f5 f6];
- f=syms(f);
- %valor medio inicial
- vm0=0;
- for i=1:length(f)
- vm0=vm0+int(f(i),'t',d(i),d(i+1));
- end
- %valor medio final
- vm=(1/T)*vm0;
- %valor eficaz inicial
- vef0=0;
- for i=1:length(f)
- vef0=vef0+int(f(i)^2,'t',d(i),d(i+1));
- end
- %valor eficaz final/
- vef=((1/T)*vef0)^0.5;
Gráfica de la función:

Resultados
- vm = 37/20 = 1.85
- vef= (4200^(1/2)*18341^(1/2))/4200 = 2.0897
Grafica, valor medio y valor eficaz:

función

Código de MatLab:
- %valor medio y eficaz de una señal escalonada
- clc
- close all
- clear all
- %señal escalonada
- T=10;
- fr=1/T;
- t=0:0.01:T;
- f1=3;
- f2=5;
- f3=7;
- f4=5;
- f5=3;
- f6=0;
- ft=f1+(f2-f1).*heaviside(t-T/5)+(f3-f2).*heaviside(t-2*T/5)+(f4-f3).*heaviside(t-3*T/5)+(f5-f4).*heaviside(t-4*T/5)+(f6-f5).*heaviside(t-T);
- plot(t,ft,'LineWidth',4)
- axis([0 11 0 8])
- grid on
- syms t
- %intervalo de integración;
- d=[0 T/5 2*T/5 3*T/5 4*T/5 T];
- %función escalonada
- f1=3;
- f2=5;
- f3=7;
- f4=5;
- f5=3;
- f6=0;
- f=[f1 f2 f2 f3 f4 f5];
- f=sym(f);
- %valor medio inicial
- vm0=0;
- for i=1:length(F)
- vm0=vm0+int(f(i),'t',d(i),d(i+1));
- end
- %valor medio final
- vm=(1/T)*vm0;
- %valor eficaz inicial
- vef0=0;
- for i=1:length(f)
- vef0=vef0+int(f(i)^2,'t',d(i),d(i+1));
- end
- %valor eficaz final
- vef=((1/T)*vef0)^0.5;
Gráfica de la función:

Resultados:
- vm = 23/5 = 4.6
- vef = (5^(1/2)*117^(1/2))/5 = 4.8374

Circuito RC

Función:

Código para MatLab
- %Analisis de circuito rc usando ecuaciones diferenciales ordinarias (EDO)
- clc
- clear all
- close all
- syms t v R C
- vc1= dsolve ('R*C*Dvc+vc=V','vc(0)=0','t');
- vc=v - v*exp(-t/(C*R));
- i1=C*diff(vc,'t');
- i=(v/R)*exp(-t/(R*C));
- vr=R*i;
- R=100;C=0.001;v=15;
- t=0:0.01:100;
- vc2=v - v.*exp(-t/(C*R));
- i2=(v/R).*exp(-t/(R*C));
- vr2=R*i2;
- plot(t,vc2,'LineWidth',4)
- ylabel('vc(V)')
- title('voltaje en el condensador')
- axis([0 16 -1 20])
- grid on
Gráfica de la función:

Resultado:
- vc= 15 - 15*exp(-t/(0.001*100))

Observaciones:

Para determinar el valor medio y eficaz se tiene que declarar los parámetros de integración y variables de integración después de el comando “plot” que sirve para graficar la señal.
Se debe tener en cuenta que en ciertas señales no son periódicas puesto se utiliza solo un periodo para determinar el valor medio y eficaz.
Para la gráfica de señal senoidal de onda completa, se integró la primera parte de la onda y se multiplico por 2 ya que la onda es igual.
Para poder graficar las señales se debe remplazar la frecuencia (f) por la frecuencia (t), además de comentar la línea “syms” y darles valores a las variables. Como recomendación añadir el comando axis para centrar la señal en mitad de la gráfica.

Conclusiones:

Al conocer el funcionamiento de los comandos para generar funciones, se pueden modificar los parámetros para dar formas a las señales y de esta forma convertirlas en la señal deseada.
Poder generar estas funciones digitalmente es útil para simular algunos procesos y así estudiar los posibles fenómenos a estudiar.
Matlab es un entorno amigable para entender y aprender código de una manera fácil y rápida.
Es útil poder hallar el valor medio y eficaz de señales comunes en electrónica de una manera tan fácil y visual.

EDO