Экспоненциал развернул синус
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')
Передайте сигнал возбуждения через бесконечную импульсную характеристику (IIR), фильтруют и добавляют шум, чтобы смоделировать реальную запись (отклик системы).
[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.