señales
ESCUELA POLITÉCNICA
DEL EJÉRCITO
CONTROL DIGITAL
DEBER #1
REGLA TRAPEZOIDAL INTEGRACION NUMERICA
DIEGO M. FLORES F.
Gabriel Moncayo
17/03/2013TEMA: Regla Trapezoidal.
OBJETIVOS:
PRINCIPAL:
Encontrar tres soluciones computacionales para resolver la regla trapezoidal, la solución debe ser por programación directa, por diagrama de flujoy mediante la función FILTER.
SECUNDARIO:
Encontrar la ecuación a diferencias.
Encontrar la función de transferencia de esta regla de integración.
Calcule el valor aproximado.
REGLATRAPEZOIDAL
CALCULO DEL ÁREA
ÁREA TOTAL.
FUNCION DE TRANSFERENCIA
MATLAB
FORMA CANONICA.%FORMA POR DIAGRAMA DE FLUJOS O CANONICA
xk1=0; %LIMPIAR MEMORIA ENTRADA
yk1=0; %LIMPIAR MEMORIA SALIDA
T=0.01; %INTERVALO DE TIEMPO PARA LA MUESTRA
k=100;%VALOR DESPEJADO DE LOS DATOS DEL EJERCICIO
td=[0:T:k*T]; %VECTOR QUE CONTIENE LOS INTERVALOS DE TIEMPO
x= exp(-td); %ENTRADA DEL EJERCICIO
for l=1:length(td)
yk2=yk1+T/2*(xk1+x(l));%ECUACION A DIFERENCIAS
xk1=x(l); %ACTUALIZACION DE MEMORIA ENTRADA
yk1=yk2; %ACTUALIZACION DE MEMORIA SALIDA
end
yk2 %PRESENTACION DE RESULTADO
Ventana en commandWindows
T=0.01
T=0.1
FORMA DIRECTA.
% PROGRAMA FORMA DIRECTA
T=0.1; %PERIODO DE TIEMPO PARA LA MUESTRA
k=10; %VALOR OBTENIDO DE LOS DATOS DEL EJERCICIOtd=[0:T:k*T]; %VECTOR DE INTERVALOS DE TIEMPO
x=exp(-td); %ENTRADA
a=0;
for l=1:length (td)
z=x(l)+a; %CON VARIABLE AUXILIAR
w=T/2*(z+a); %ECUACION DE LA SALIDA
a=z; %ACTUALIZACIONDE LA MEMORIA
end
w %SALIDA
Ventana en command Windows
T=0.1
T=0.01
FUNCION FILTER.
%FUNCION FILTER
y=filter(T/2.*[1 1],[1 -1],x);
y
Ventana en command...
Regístrate para leer el documento completo.