1.5.
Considere gráficamente la señal X a (t)=3 sen(100 πt)
a) Represente gráficamente la señal X a (t) para 0≤ t ≤ 30 ms
Código:
clc
close all
clear all
ts=0:0.001:0.03;
xt=3*sin(100*pi*ts);
plot(ts,xt,'b');
title('Ejercicio 5 Literal a)');
grid on
axis([-0.001 0.031 -3 3])
Gráfica
b) La señal X a (t) se muestra con una tasa de muestreo de F s= 300 muestras/s. Determine la
frecuencia de la señal discreta en el tiempo x ( n )=x a ( nT ) , T =1/F s , y demuestre que es
periódica.
x ( n )=X a ( nT )
( F1 )
x ( n )=X a n
s
x ( n )=3 sin n( 100300π )
πn
x ( n )=3 sin ( )
3
πn
f =(
o
6 πn )
f o= ( 16 )
T =6
(Es una señal periódica no hay números racionales)
c) Calcule los valores de las muestras en un período de x(n). Dibuje en la misma gráfica que
X a (t) .
¿Cuál es el período de la señal discreta en el tiempo en milisegundos?
N x(n) πn
0 0
x ( n )=3 sin ( 100 )
1 2.598
2 2.598 100 π =2 πf
3 0
4 -2.598 100 π
f=
5 -2.598 2π
6 0
f =50
1
T= =20 ms
50
Fig. 1.1. Señal tiempo discreto
1.15 Muestreo de señales sinusoidales: aliasing. Considere la siguiente señal sinusoidal
continua en el tiempo:
x a ( t )=sen ( 2 πt F0 ) −∞ <t ←∞
Puesto que X a ( t ) está descrita en forma matemática, su versión muestreada puede describirse
mediante sus valores cada T segundos. La señal muestreada se describe mediante la fórmula
siguiente
F0
(
x a ( n )=x a ( nT )=sen 2 πt
Fs)n , −∞<t ←∞
donde Fs = 1/T es la frecuencia de muestreo.
(a) Represente gráficamente la señal x(n), 0 ≤ n ≤ 99 para Fs = 5 kHz y F0 = 0.5, 2, 3 y 4.5
kHz.
F S=5 kHz y F o=500 Hz
Fig. 1.2. Señal tiempo discreto a Fs=5kHz y Fo=500Hz
F S=5 kHz y F o=2000 Hz
Fig. 1.3. Señal tiempo discreto a Fs=5kHz y Fo=2000Hz
F S=5 kHz y F o=3000 Hz
Fig. 1.4. Señal tiempo discreto a
Fs=5kHz y Fo=3000Hz
F S=5 kHz y F o=4500 Hz
Fig. 1.5. Señal tiempo discreto a Fs=5kHz
y Fo=4500Hz
Explique las similitudes y diferencias
entre las distintas representaciones.
- Como podemos observar se muestreo la señal x a ( t ) en el dominio de tiempo
discreto en la Fig. 1.2 discretizamos la señal con una frecuencia de 5000Hz a
una frecuencia de 500Hz, se observa que se cumple el teorema Nysquit
F s ≥ 2 Fo para evitar el aliasing.
- En la Fig. 3. Aumentamos F oa 2000Hz, con esta frecuencia todavía se cumple
el teorema de Nysquit, por lo tanto aun no se produce el aliasing, además
que podemos interpolar la señal para obtener toda la información de la
misma.
- En la Fig. 4. Trabajamos con F o=3000Hz, aquí observamos que ya se esta
produciendo un aliasing debido que no estamos cumpliendo con el teorema
correspondiente de la señal.
- En la Fig. 5 Trabajamos con F o=4500Hz, podemos determinar que no se
cumple para nada el teorema de Nysquit, por lo tanto se puede observar el
aliasing en la señal producida.
(b) Suponga que F0 = 2 kHz y Fs = 50 kHz.
1. Represente gráficamente la señal x(n). ¿Cuál es la frecuencia f0 de la señal x(n)?
πn
x ( n )=3 sin ( 100 )
100 π =2 πf
100 π
f= kHz
2π
f =25 kHz
2. Represente gráficamente la señal y(n) obtenida tomando las muestras pares de x(n). ¿Es una
señal sinusoidal? ¿Por qué? En caso afirmativo, ¿cuál es su frecuencia?
1.16 Error de cuantificación en la conversión A/D de una señal sinusoidal Sea x q(n) la señal
obtenida al
cuantificar la señal x(n)=sen2π f0n. La potencia del error de cuantificación Pq se define como
sigue:
La “calidad” de la señal cuantificada puede medirse mediante la relación señal-ruido de
cuantificación (SQNR) definida como sigue:
donde Px es la potencia de la señal sin cuantificar x(n).
(a) Para f0 = 1/50 y N = 200, escriba un programa para cuantificar la señal x(n), utilizando
truncamiento
para 64, 128 y 256 niveles de cuantificación. En cada caso, represente gráficamente las
señales x(n), xq(n) y e(n), y calcule la relación SQNR correspondiente.
X(n) señal de 64 niveles
Codigo:
A=input('Ingrese la amplitud: ');
fo=input('Ingrese la frecuencia: ');
p=input('Ingrese la fase: ');
n=input('Ingrese las n muestras ');
b=input('Ingrese los n Bits: ');
disp('ARCHIVO EXCEL');
file=input('\n\nNombre del archivo: ','s');
hoja=input('Hoja en la que deseas trabajar:','s');
desc=[{'N de muestra'},{'Tiempo(seg)'},{'Voltaje'},{'Valor binario'}];
T=1/fo; %Periodo
Tm=T/n;%Tiempo de muestreo
nb=b;
switch b
case b
for x=1:1:n
t=Tm*x;
y=sin(2*pi*fo*t+p);
if y>(0) && y<=(A/nb);
f=dec2bin(x,nb);
fprintf('%g %f %f %s \n\n',x,t,y,f)
elseif y>(x*A/64) && y<=((x+1)*A/64)
f=dec2bin(1,64);
fprintf('%g %f %f %s \n\n',x,t,y,f)
end
v(x)={f};
end
disp(v);
end
t=0:T/n:2*T;
y=(A/2)*sin(2*pi*fo*t+p)+A/2;
z=A/2*sin(2*pi*fo*t+p)+A/2;
stem(t,z);
hold on;
plot(t,y,'r');
grid on;
title('Muestreo (Señal Seno)')%titulo
xlabel('Tiempo(seg)'); %eje horizontal
ylabel('Amplitud(Volt.)')%eje vertical
%axis([-0.01 0.02 0 2])
X(n) señal de 128 niveles
X(n) señal de 256 niveles
