Сгенерируйте 1 024 выборки щебета, произведенного на уровне 1 024 кГц. Задайте щебет так, чтобы он имел начальную частоту 50 кГц и достиг 100 кГц в конце выборки. Добавьте белый Гауссов шум, таким образом, что отношение сигнал-шум составляет 40 дБ. Сбросьте генератор случайных чисел для восстанавливаемых результатов.

nSamp = 1024;
Fs = 1024e3;
SNR = 40;
rng default

t = (0:nSamp-1)'/Fs;

x = chirp(t,50e3,nSamp/Fs,100e3);
x = x+randn(size(x))*std(x)/db2mag(SNR);

Оцените медианную частоту щебета. Постройте спектральную плотность мощности (PSD) и аннотируйте медианную частоту.

medfreq(x,Fs)

ans = 7.4998e+04

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

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

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

y = medfreq([x x2],Fs)

y = 1×2
105 ×

    0.7500    2.4999

Постройте PSDs двух каналов и аннотируйте их медианные частоты.

medfreq([x x2],Fs);

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

medfreq(x+x2,Fs)

ans = 2.3756e+05

Сгенерируйте 1 024 выборки синусоиды на 100,123 кГц, произведенной на уровне 1 024 кГц. Добавьте белый Гауссов шум, таким образом, что отношение сигнал-шум составляет 40 дБ. Сбросьте генератор случайных чисел для восстанавливаемых результатов.

nSamp = 1024;
Fs = 1024e3;
SNR = 40;
rng default

t = (0:nSamp-1)'/Fs;

x = sin(2*pi*t*100.123e3);
x = x + randn(size(x))*std(x)/db2mag(SNR);

Используйте periodogram функция, чтобы вычислить спектральную плотность мощности (PSD) сигнала. Задайте окно Кайзера с той же длиной как сигнал и масштабный фактор 38. Оцените медианную частоту сигнала и аннотируйте его на графике PSD.

[Pxx,f] = periodogram(x,kaiser(nSamp,38),[],Fs);

medfreq(Pxx,f);

Сгенерируйте другую синусоиду, этого с частотой 257,321 кГц и амплитудой, которая является дважды больше чем это первой синусоиды. Добавьте белый шум.

x2 = 2*sin(2*pi*t*257.321e3);
x2 = x2 + randn(size(x2))*std(x2)/db2mag(SNR);

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

[Pyy,f] = periodogram([x x2],kaiser(nSamp,38),[],Fs);

y = medfreq(Pyy,f)

y = 1×2
105 ×

    1.0012    2.5731

Аннотируйте медианные частоты двух каналов на графике PSDs.

medfreq(Pyy,f);

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

[Pzz,f] = periodogram(x+x2,kaiser(nSamp,38),[],Fs);

medfreq(Pzz,f);

Сгенерируйте сигнал, PSD которого напоминает частотную характеристику КИХ-фильтра полосы пропускания 88-го порядка с нормированными частотами среза $$0.25\pi$$ рад/отсчет и $$0.45\pi$$ рад/отсчет.

d = fir1(88,[0.25 0.45]);

Вычислите медианную частоту сигнала между $$0.3\pi$$ рад/отсчет и $$0.6\pi$$ рад/отсчет. Постройте PSD и аннотируйте интервал измерения и медианная частота.

medfreq(d,[],[0.3 0.6]*pi);

Выведите медианную частоту и мощность полосы интервала измерения. Определение частоты дискретизации $$2\pi$$ эквивалентно отъезду сброса уровня.

[mdf,power] = medfreq(d,2*pi,[0.3 0.6]*pi);

fprintf('Mean = %.3f*pi, power = %.1f%% of total \n', ...
    mdf/pi,power/bandpower(d)*100)

Mean = 0.371*pi, power = 77.4% of total 

Добавьте второй канал с нормированными частотами среза $$0.5\pi$$ рад/отсчет и $$0.8\pi$$ рад/отсчет и амплитуда, которая является одной десятой тот из первого канала.

d = [d;fir1(88,[0.5 0.8])/10]';

Вычислите медианную частоту сигнала между $$0.3\pi$$ рад/отсчет и $$0.9\pi$$ рад/отсчет. Постройте PSD и аннотируйте медианную частоту каждого канала и интервала измерения.

medfreq(d,[],[0.3 0.9]*pi);

Выведите медианную частоту каждого канала. Разделитесь на $$\pi$$.

mdf = medfreq(d,[],[0.3 0.9]*pi)/pi

mdf = 1×2

    0.3706    0.6500