domingo, 12 de febrero de 2012

Operando una función infinita...

Respondiendo al ejercicio planteado de graficar tres diferentes situaciones de una función infinita, he elegido como mi función:     x(t) = 3cos(2t).
A través de la ejecución del siguiente script de Matlab se obtienen dos figuras cada una con dos gráficas para un total de 4 gráficas.

1 - Gráfica de x(t)
2 - Gráfica de y(t) = Una constante A + x(t)
3 - Gráfica de z(t) = x(t) + v(t) ; donde v(t) es una función infinita con un periodo igual a N veces enteras el periodo de x(t)
4 - Gráfica de q(t) = x(t) + u(t) ; donde u(t) es una función infinita con un periodo diferente a x(t) pero no necesariamente un múltiplo de este.

Se usan 3 funciones que deben ser previamente guardadas en el directorio de Matlab para la correcta ejecución del programa y son descritas al final del script.


%SCRIPT:

t=0:5:100;  %Vector tiempo

%Grafico funcion x(t)
x=f_cost(t);
figure
subplot(2,1,1)
plot(t,x)
grid on
title('Funcion x(t) = 3cos(2t)')
xlabel('t')
ylabel('x')


a=5; %Constante para sumar

%Grafico funcion y(t)
y=a+f_cost(t);
subplot(2,1,2)
plot(t,y)
grid on
xlabel('t')
ylabel('y')
title('Funcion y(t) = a + 3cos(2t) // a = 5 ')

w=2;
%Frecuencia angular w=2=2PI/T  -->  T=2PI/2=PI
%Para V(t) --> To=nT = 4PI ---> w=2PI/4PI = 0.5

%Grafico funcion z(t)
v=fv_cost(t);
z=x+v;
figure
subplot(2,1,1)
plot(t,z)
grid on
xlabel('t')
ylabel('z')
title('Funcion z(t) = x(t) + v(t) // v(t) = 3cos(0.5t)')

%Grafico funcion q(t)
u=fu_cost(t);
q=x+u;
subplot(2,1,2)
plot(t,q)
grid on
title('Funcion q(t) = x(t) + u(t) // u(t) = 3cos(3t)')
xlabel('t')
ylabel('q')

%Se obtiene como resultado dos ventanas (usualmente aparecen una sobre otra)
%cada una contiene dos graficas.

 ----------------------------------//

%FUNCIONES PREVIAS:

FUNCION f_cost.m

function [ x ] = f_cost( t )
x=3*cos(t.*2);
end


FUNCION fv_cost.m

function [ v ] = fv_cost( t )
v=3*cos(t.*0.5);
end


FUNCION fu_cost.m

function [ u ] = fu_cost( t )
u=3*cos(t.*3);
end

-----------------------------//

Se obtienen las siguientes imágenes:



 Eso es todo por ahora, habrá mas cada vez que aprenda algo nuevo....

1 comentario:

  1. Para que la función tenga una forma más senoidal y curva, se sugiere modificar los intervalos de tiempo en el eje x.

    ResponderEliminar