Как fdatool преобразует нули и полюса в numerator/denominator Опции

[coding4dsp]

Добрый день!

Пример.

У фильтра только один ноль, равный +1j*0.5. Если создадим эти условия в fdatool и перейдем на вкладку FilterCoefficients, то увидим FIR фильтр с коэффициентами: [1, -0.5j]. Как FDATool это делает?
Проверим, получим нули и полюса от numerator: tf2zpk([1, -0.5j],1) выдаст ноль +1j*0.5. Правильно!

Теперь я хочу в m-коде, используя zp2tf() получить numerator, как в fdatool.
[a,b] = zp2tf((1j*0.5),0,1); не выдает комплексный numerator.
zp2tf использую не по правилам? Других функций, подходящих по смыслу, не наше.
Не понимаю, что творю, заблудился?!

Update: реализовал tf2zpk через функцию poly. Все равно интересно, почему tf2zpk не выдает комплексных чисел.

Сообщение отредактировал coding4dsp - Dec 21 2017, 16:20

[andyp]

Не понимаю, что творю, заблудился?!

Да, имхо. Коэффициенты передаточной функции предполагаются вещественными если охота остаться в рамках классических LTI систем. Отсюда и real() в коде zp2tf (см. type zp2tf).

Если коэффициенты передаточной функции - комплексные, то придется столкнуться с рядом неожиданных эффектов - например, попробуйте freqz(1,[1 (1+8j)], 'whole'). АЧХ перестала быть чётной. упс.

Можно почитать например
http://portal.research.lu.se/portal/files/...ACC_0939_FI.pdf

если тянет к странному.

[coding4dsp]

Коэффициенты передаточной функции предполагаются вещественными если охота остаться в рамках классических LTI систем.

Про это. Я не понял смысл. FIR фильтр с комплексными коэффициентами не соответствует LTI системе?

Сообщение отредактировал coding4dsp - Dec 22 2017, 13:38