Как решить разностное уравнение с использованием Z-преобразования в Matlab
Я пытаюсь найти импульсную характеристику фильтра, выраженную уравнением разности, как показано ниже:
y [n] - y [n-1] + 0,9*y[n-2] = x[n] - x[n-2]
Я хочу использовать Z-преобразование, чтобы решить уравнение. Я использую этот код в Matlab:
syms y(n) z
eq = y(n) - y(n - 1) - 0.9*y(n - 2) - impulse_me(n) + impulse_me(n - 2);
Zeq = ztrans(eq, n, z);
и я объявил функцию impulse_me, как показано ниже:
function [ y ] = impulse_me( n )
y = (n==0);
end
но когда я запускаю код, результат:
Zeq =
ztrans(y(n) - (9*y(n - 2))/10 - y(n - 1) - 2 == y(n) - (9*y(n - 2))/10 - y(n - 1), n, z)
Кто-нибудь может мне помочь с проблемой?