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

Генерируют 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