Не вдаваясь в анализ Ваших результатов, могу только еще раз повторить, что к системам типа дециматор-интерполятор совершенно не применимы такие понятия, как АЧХ и передаточная ф-я. Так что вопрос "что правильно 6 или 30?" не имеет смысла, ибо не правильно ни то ни другое.

В дополнение могу порекомендовать ознакомиться с работами Vaidyanathan P. P., например Quadrature mirror filter banks, M-band extensions and perfect-reconstruction techniques, цитата:

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

И посмотреть (хотя бы на Wikipedia) что такое Time-variant system:

"There are many well developed techniques for dealing with the response of linear time invariant systems, such as Laplace and Fourier transforms. However, these techniques are not strictly valid for time varying systems."

"It does not have an impulse response in the normal sense"

Ну уж совсем ничем характеризоваться они не могут наверное. Что-то должно быть, что?

- а какое подавление должно быть на самом деле если на вход системы будет подана частота 0.25?

Далась Вам эта частота)). Полистайте форум, пару-тройку месяцев назад тут были серьезные диспуты по поводу теоремы Котельникова (в который раз уже) и частоты 0.5*Fd (ну или 0.25*Fd до децимации). С практической точки зрения говорить про эту частоту не имеет смыла, в этом легко убедиться если посмотреть на результат фильтрации и децимации двух синусоидальных сигналов с одной и той же частотой 0.25*Fd но с разными фазами.

Если же речь идет о частотах меньше 0.5*Fd (0.25*Fd), то насколько я помню, эта тема уже обсуждалась: Halfband фильтр...

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

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

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

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

Ясно же сказали, что к системам с дециматором неприменимо (в полной мере) понятие ЧХ.
Исходные ЧХ (до децимации, синие) совпадают.
После децимации (малиновые) различаются из-за того, что выброс входного сигнала попадает под децимацию. Если конкретно для децимации в 2 раза иметь входной сигнал с ДВУМЯ единицами подряд, малиновые графики должны совпасть более точно (не проверял, но так думаю). Но это - уже не частотная характеристика.

Еще раз:

"...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

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

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


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

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

s = [zeros(1,127) 1 1 zeros(1,127)];

s = [zeros(1,127) 1 1 zeros(1,128)];

Как видите, не устойчив... А, нет, одинаковый.