Generacion de una señal alias tiempo discreto

Generación de una señal alías
En esta práctica consistió en observar la generación de una señal alías cuando aumentamos la frecuancia de la señal a muestrear, pero la frecuencia de muestraes la misma.
El tiempo de muestra se base en el teorema de Nyquist, lo cual en pocas palabras dice que la frecuencia de muestra tiene que ser por mínimo dos veces la frecuencia de laseñal.
ωs≥2ω0
Lo que implica
ts=2πωs=tiempo de muestra
Otro punto de la practica era ver que cuando aumentabamos la frecuencia, había en una frecuencia en especial donde las señales eran lasmismas.
La frecuencia en la que las señales eran iguales es
ωs-ωaumentada=ω0
El código en MATLAB es el siguiente:

clc,clear all, close all;

for i=1:60

%% funcion continuaoriginal
subplot(2,2,1)
t=-10:0.01:10;
w0=pi;
ws=3*w0;
ts=(2*pi)/ws;
f=cos(w0*t);
tM=-10:ts:10;
fM=cos(w0*tM);
plot(t,f)
hold on
stem(tM,fM,'r')
hold off
%espectrooriginal
subplot(2,2,2)

frec_ang=[-(ws+w0),-(ws-w0),-w0,w0,(ws-w0),(ws+w0)];
F=ones(1,6);
stem(frec_ang,F)

%% Funcion continua animada
subplot(2,2,3)
w0v=i;
%wsv=10*w0v;%tsv=(2*pi)/wsv;
fv=cos(w0v*t)
%tMv=-10:tsv:10;
fMv=cos(w0v*tM);
plot(t,fv)
hold on
stem(tM,fM,'r')
hold off

%espectro animado
subplot(2,2,4)frec_angv=[-(3*ws+w0v),-(3*ws-w0v),-(2*ws+w0v),-(2*ws-w0v),-(ws+w0v),-(ws-w0v),-w0v,w0v,(ws-w0v),(ws+w0v),(2*ws-w0v),(2*ws+w0v),(3*ws-w0v),(3*ws+w0v)];
F=ones(1,14);
stem(frec_angv,F)
pause()
if (ws-w0v)==w0break
end
end

Conclusión: en mi caso no pude acabar bien la práctica, pero el sentido teorico si lo pude comprender. Y mi experiencia fue que si no cumple con el teorema de Nyquist,la señal no puede ser la misma, generandose una señal alías.
Por otra parte pude observar que hay una frecuencia en la cual las señales se repite tanto en el tiempo como en la frecuencia.