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....
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