Огибающая сигнала
[ возвращает верхние и более низкие конверты входной последовательности, yupper,ylower] = envelope(x)x, как величина его аналитического сигнала. Аналитический сигнал x найден с помощью дискретного преобразования Фурье, как реализовано в hilbert. Функция первоначально удаляет среднее значение x и добавляет его назад после вычисления конвертов. Если x матрица, затем envelope действует независимо по каждому столбцу x.
[ возвращает конверты yupper,ylower] = envelope(x,fl,'analytic')x решительное использование величины его аналитического сигнала. Аналитический сигнал вычисляется путем фильтрации x с Гильбертовым КИХ-фильтром длины fl. Этот синтаксис используется, если вы задаете только два аргумента.
[ возвращает верхние и более низкие среднеквадратичные конверты yupper,ylower] = envelope(x,wl,'rms')x. Конверты определяются с помощью раздвижного окна длины wl выборки.
[ возвращает верхние и более низкие пиковые конверты yupper,ylower] = envelope(x,np,'peak')x. Конверты определяются использование интерполяции сплайна по локальным максимумам, разделенным, по крайней мере, np выборки.
envelope(___) без выходных аргументов строит сигнал и его верхние и более низкие конверты. Этот синтаксис принимает любой из входных параметров от предыдущих синтаксисов.
Сгенерируйте квадратичный щебет, модулируемый Гауссовым. Задайте частоту дискретизации 2 кГц и длительность сигнала 2 секунд.
t = 0:1/2000:2-1/2000; q = chirp(t-2,4,1/2,6,'quadratic',100,'convex').*exp(-4*(t-1).^2); plot(t,q)
Вычислите верхние и более низкие конверты щебета с помощью аналитического сигнала.
[up,lo] = envelope(q); hold on plot(t,up,t,lo,'linewidth',1.5) legend('q','up','lo') hold off
Сигнал асимметричен из-за ненулевого среднего значения.
Используйте envelope без выходных аргументов, чтобы построить сигнал и конверты в зависимости от демонстрационного номера.
envelope(q)
Создайте двухканальный сигнал, произведенный на уровне 1 кГц в течение 3 секунд:
Один канал является экспоненциально затухающей синусоидой. Задайте частоту 7 Гц и постоянную времени 2 секунд.
Другой канал является перемещенным временем Модулируемым гауссовым образом щебетом со значением DC 2. Задайте начальную частоту щебета 30 Гц, которая затухает к 5 Гц после 2 секунд.
Постройте сигнал.
t = 0:1/1000:3; q1 = sin(2*pi*7*t).*exp(-t/2); q2 = chirp(t,30,2,5).*exp(-(2*t-3).^2)+2; q = [q1;q2]'; plot(t,q)
Вычислите верхние и более низкие конверты сигнала. Используйте Гильбертов фильтр с длиной 100. Постройте каналы и конверты. Используйте сплошные линии для верхних конвертов и пунктирные линии для более низких конвертов.
[up,lo] = envelope(q,100,'analytic'); hold on plot(t,up,'-',t,lo,'--') hold off
Вызовите envelope без выходных аргументов, чтобы произвести график сигнала и его конвертов в зависимости от демонстрационного номера. Увеличьте длину фильтра до 300, чтобы получить более сглаженную форму. 'analytic' флаг является значением по умолчанию, когда вы задаете два входных параметра.
envelope(q,300)
Вычислите и постройте движущиеся конверты RMS записи гудка поезда. Используйте окно с продолжительностью 150 выборок.
load('train') envelope(y,150,'rms')
График верхние и более низкие пиковые конверты речевого сигнала смягчил интервалы с 30 выборками.
load('mtlb') envelope(mtlb,30,'peak')
Создайте и постройте сигнал, который напоминает начальное обнаружение распространения светового импульса через дисперсионный носитель.
t = 0.5:-1/100:-2.49; z = airy(t*10).*exp(-t.^2); plot(z)
Определите конверты последовательности с помощью величины ее аналитического сигнала. Постройте конверты.
envelope(z)
Вычислите аналитический конверт сигнала с помощью Гильбертова фильтра с 50 касаниями.
envelope(z,50,'analytic')
Вычислите конверт RMS сигнала с помощью движущегося окна с 40 выборками. Постройте результат.
envelope(z,40,'rms')
Определите пиковые конверты. Используйте интерполяцию сплайна с условиями не-узла по локальным максимумам, разделенным по крайней мере 10 выборками.
envelope(z,10,'peak')