Экспоненциальный развертываемый синус
excitation = sweeptone()
excitation = sweeptone(swDur)
excitation = sweeptone(swDur,silDur)
excitation = sweeptone(swDur,silDur,fs)fs Гц.
excitation = sweeptone(___,Name,Value)Name,Value пары аргументов, в дополнение к входным аргументам в предыдущих синтаксисах.
Создайте сигнал возбуждения тона свип-сигнала с помощью sweeptone функция.
excitation = sweeptone(2,1,44100);
plot(excitation)
title('Excitation')
Пропускайте сигнал возбуждения через фильтр с бесконечной импульсной характеристикой (БИХ) и добавляйте шум для моделирования реальной записи (системной характеристики).
[B,A] = butter(10,[.1 .7]);
rec = filter(B,A,excitation);
nrec = rec + 0.12*randn(size(rec));
plot(nrec)
title('System Response')
Передать сигнал возбуждения и отклик системы на impzest функция для оценки импульсной характеристики. Усечь оценку до 100 пунктов. Использовать impz для определения истинной импульсной характеристики системы. Постройте график истинной импульсной характеристики и расчетной импульсной характеристики для сравнения.
irEstimate = impzest(excitation,nrec); irEstimate = irEstimate(1:101); irTrue = impz(B,A,101); plot(0:100,irEstimate, ... 0:100,irTrue,'ro') legend('True impulse response','Estimated impulse response')
Создайте экспоненциальный сигнал синусоидального сдвига (ESS) с 3-секундным сдвигом, который продолжается от 20 Гц до 20 кГц и заканчивается 2-секундным молчанием. Укажите частоту дискретизации 48 кГц.
fs = 48e3;
excitation = sweeptone(3,2,fs,'SweepFrequencyRange',[20 20e3]);Визуализируйте возбуждение во времени и частоте.
t = (0:numel(excitation)-1)/fs;
plot(t,excitation)
xlabel('Time (s)')
spectrogram(excitation,512,0,1024,fs,'yaxis')
[1] Фарино, Анджело. «Достижения в измерениях импульсной характеристики с помощью синусоидальных сдвигов». Представлен на 122-й конференции Общества аудиотехники, Вена, Австрия, 2007 год.