Экспоненциальный стреловидный синус
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 год.