Introducción - procesado digital 1

Para la primera entrada de este nuevo tutorial solo pondré como mostrar una señal sinusoidal muestreada cada intervalo Tm usando matlab u octave:

fsin=500;               % frecuencia de la señal en Hz

tiempo=0.05;         % tiempo total

Tm=1/20000;        % tiempo cada cuanto se toma un valor(muestreo)*

%vector de tiempos, desde la primera muestra : sumando muestras : hasta tiempo total:

vt=(Tm:Tm:tiempo);

fase=pi/8;                %añadimos una fase inicial

x=sin(2*pi*fsin*vt+fase);   %definimos una señal sinusoidal con los parámetros anteriores

plot(vt,x,'.');             %mostramos por pantalla la señal con respecto a los valores del vector tiempo

xlabel('tiempo');       % Nombramos al eje x como eje de "tiempo"

ylabel('amplitud');    % Nombramos al eje y como eje de "Amplitud"

*La frecuencia de muestreo(es decir 1/Tm) tiene que cumplir el teorema del muestro de Nyquist para que la señal original pueda ser recuperada de la señal muestreada, para eso tiene que cumplir la siguiente ecuación:

frecuencia de muestreo > 2*frecuencia máxima de la señal original

Nota: Los '%' indican que empieza un comentario.