Измерение средней частоты, степени, полосы пропускания

Сгенерировать 1024 выборки щебета, отобранного с частотой дискретизации 1024 кГц. Щебет имеет начальную частоту 50 кГц и достигает 100 кГц в конце дискретизации. Добавьте белый Гауссов шум таким образом, чтобы отношение сигнал/шум составляло 40 дБ.

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);

Оцените 99% занимаемой полосы сигнала и аннотируйте его на графике спектральной плотности степени (PSD).

obw(x,Fs);

Вычислите степень в полосе и проверьте, что она составляет 99% от общего количества.

[bw,flo,fhi,powr] = obw(x,Fs);

pcent = powr/bandpower(x)*100

pcent = 99.0000

Сгенерируйте еще один щебет. Задайте начальную частоту 200 кГц, конечную частоту 300 кГц и амплитуду, которая в два раза больше, чем у первого сигнала. Добавьте белый Гауссов шум.

x2 = 2*chirp(t,200e3,nSamp/Fs,300e3);
x2 = x2+randn(size(x2))*std(x2)/db2mag(SNR);

Добавьте два щебета, чтобы сформировать новый сигнал. Постройте график PSD сигнала и аннотируйте его медианную частоту.

medfreq([x+x2],Fs);

Постройте график PSD и аннотируйте среднюю частоту.

meanfreq([x+x2],Fs);

Теперь рассмотрим каждый щебет, чтобы представлять отдельный канал. Оцените среднюю частоту каждого канала. Аннотировать средние частоты на графике PSD.

meanfreq([x x2],Fs)

ans = 1×2
105 ×

    0.7503    2.4999

Оцените половину полуполосы пропускания мощности каждого канала. Аннотируйте 3-dB полосы пропускания на графике PSD.

powerbw([x x2],Fs)

ans = 1×2
104 ×

    4.4386    9.2208

См. также

bandpower | meanfreq | medfreq | obw | powerbw