Medir frecuencia media, potencia, ancho de banda
Genere 1024 muestras de un chirp muestreado a 1024 kHz. El chirp presenta una frecuencia inicial de 50 kHz y alcanza los 100 kHz al finalizar el muestreo. Añada ruido blanco gaussiano hasta que la relación señal-ruido sea de 40 dB.
nSamp = 1024; Fs = 1024e3; SNR = 40; t = (0:nSamp-1)'/Fs; x = chirp(t,50e3,nSamp/Fs,100e3); x = x+randn(size(x))*std(x)/db2mag(SNR);
Calcule el 99% del ancho de banda ocupado de la señal y anótelo en una gráfica de densidad del espectro de potencia (PSD).
obw(x,Fs);
Calcule la potencia en la banda y verifique que es el 99% del total.
[bw,flo,fhi,powr] = obw(x,Fs); pcent = powr/bandpower(x)*100
pcent = 99.0000
Genere otro chirp. Especifique una frecuencia inicial de 200 kHz, una frecuencia final de 300 kHz y una amplitud que sea el doble que la de la primera señal. Añada ruido blanco gaussiano.
x2 = 2*chirp(t,200e3,nSamp/Fs,300e3); x2 = x2+randn(size(x2))*std(x2)/db2mag(SNR);
Añada los dos chirps para formar una nueva señal. Represente la PSD de la señal y anote su frecuencia mediana.
medfreq([x+x2],Fs);
Represente la PSD y anote la frecuencia media.
meanfreq([x+x2],Fs);
Ahora, piense que cada chirp representa un canal independiente. Calcule la frecuencia media de cada canal. Anote las frecuencias medias en una gráfica de las PSD.
meanfreq([x x2],Fs)
ans = 1×2
105 ×
0.7503 2.4999
Calcule el ancho de banda de media potencia de cada canal. Anote los anchos de banda de 3 dB en una gráfica de las PSD.
powerbw([x x2],Fs)
ans = 1×2
104 ×
4.4386 9.2208