envelope
Огибающая сигнала
Синтаксис
Описание
[ возвращает верхние и более низкие конверты входной последовательности, 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 аудиозаписи
Вычислите и постройте движущиеся конверты 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')