Функция freqz.

[Ivan52]

Добрый вечер.

Мной строится ФНЧ Чебышева 1-го рода. При реализации на заключительной стадии необходимо применить команду:

[h, w] = freqz(a, b);

Подскажите, пожалуйста, как реализовать эту функцию? То есть не использовать готовую freqz, а реализовать свою собственную.
Мной был получен порядок ФНЧ:

[N, W_p] = cheb2ord(W_p, W_s, R_p, R_s);

Затем был рассчитан ФНЧ:

[a, b] = cheby1(N, R_p, W_p);

Теперь необходимо найти АЧХ рассчитанного фильтра. Это можно было бы сделать с помощью встроенной функции freqz, но по заданию требуется ее реализовать. То есть не использовать готовую функцию freqz, а написать собственную. Помогите, пожалуйста, это реализовать.
Спасибо.

[Ivan52]

Обнаружил в сети реализацию функции freqz, но, либо я наверно что-то не так делаю, так как она выдает иной результат, нежели функция freqz. Подскажите, как оптимизировать этот код, чтобы он работал корректно?
Код функции freqz:

N = 1024; % Number of points to evaluate at
upp = pi; % Evaluate only up to fs/2
% Create the vector of angular frequencies at one more point.
% After that remove the last element (Nyquist frequency)
w = linspace(0, pi, N+1); 
w(end) = [];
ze = exp(-1j*w); % Pre-compute exponent
H = polyval(b, ze)./polyval(a, ze); % Evaluate transfer function and take the amplitude
Ha = abs(H);
Hdb  = 20*log10(Ha); % Convert to dB scale
wn   = w/pi;
% Plot and set axis limits
xlim = ([0 1]);
plot(wn, Hdb)
grid on

[Ivan52]

Тема более не актуальна.