Характеристика линейной системы Опции

[Leshii]

Предположим, что у нас имеется некий черный ящик на вход которого можно подавать какой-либо сигнал.
Задача оценить его (черного ящика) частотно-фазовую характеристику.
Если я подам на вход системы импульс, то вроде бы на выходе ЛС будет отклик, соответствующих ИХ данной ЛС и если взять Фурье от полученной ИХ то и получим частотно-фазовую характеристику.
Вопрос - насколько правомочно применение данного подхода, если одним из компонентов в черном ящике является устройство децимации?

Сообщение отредактировал Leshii - Apr 17 2009, 20:23

[ViKo]

Перескочил сюда из недавней темы. Попробовал то, что здесь предложено. Но изменил длину входной последовательности, чтобы было кратно степени 2.

Код

clear;
clf;

% Единичный импульс
s = [zeros(1,127) 1 zeros(1,128)];

% Импульсная характеристика фильтра
h = [0.042041 0 -0.087017 0 0.31156 0.5 0.31156 0 -0.087017 0 0.042041];

% Фильтр
s_f = filter(h,1,s);

% Дециматор в два раза
s_d = s_f(1:2:end);

% Преобразование Фурье от ИХ
grid on
hold on
plot(20 * log10(abs(fft(s_f)) + eps), 'b');
plot(20 * log10(abs(fft(s_d)) + eps), 'm');
hold off

Имею картинку, отличную от показанной топикстартером.

Если же последовательность другая, то и картинка другая.
s = [zeros(1,128) 1 zeros(1,127)];

Как это понимать?

[Самурай]

...
Как это понимать?

Еще раз:

"...decimator is a time-varying device, even thought it is linear. Accordingly, it cannot be represented by a transfer function..."

И для закрепления материала вот так попробуйте:

% Дециматор в два раза
s_d = s_f(2:2:end);

Сообщение отредактировал Самурай - Oct 20 2014, 09:30

[ViKo]

И для закрепления материала вот так попробуйте:
% Дециматор в два раза
s_d = s_f(2:2:end);

И что? Первый вариант поменялся со вторым.

С двумя единицами уже что-то бессмысленное вылазит...

После децимации (малиновые) различаются из-за того, что выброс входного сигнала попадает под децимацию.

Очень уж удачно попадает. Так, что плоская АЧХ получилась. Ну, пусть не АЧХ...

[V_G]

С двумя единицами уже что-то бессмысленное вылазит...

Вот с этим сигналом:

Код

% Единичный импульс
s = [zeros(1,127) 1 1 zeros(1,128)];

?
Ничего не бессмысленное, а спектр отклика на сигнал из двух единиц подряд. Он уже устойчив к сдвигу при децимации в 2 раза.