Сгенерируйте 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) и аннотируйте среднюю частоту.

meanfreq(x,Fs)

ans = 7.5032e+04

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

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

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

y = meanfreq([x x2],Fs)

y = 1×2
105 ×

    0.7503    2.4999

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

meanfreq([x x2],Fs);

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

meanfreq(x+x2,Fs)

ans = 2.1496e+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);

meanfreq(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 = meanfreq(Pyy,f)

y = 1×2
105 ×

    1.0013    2.5732

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

meanfreq(Pyy,f);

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

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

meanfreq(Pzz,f);

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

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

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

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

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

[mnf,power] = meanfreq(d,2*pi,[0.3 0.6]*pi);

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

Mean = 0.373*pi, power = 75.6% of total 

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

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

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

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

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

mnf = meanfreq(d,[],[0.3 0.9]*pi)/pi

mnf = 1×2

    0.3730    0.6500