IIR Matlab Опции

[Mad_max]

Доброго времени суток!

Сообразил в матлабе фильтр

Код

Fs=8000;        % Sampling frequency 8 kHz
fc=800;        % Passband cutoff frequency 800 Hz
fs=1600;        % Stopband frequency 1.6 kHz
Rp=1;            % Passband ripple in dB
Rs=20;        % Stopband attenuation in dB
Wp=2*fc/Fs;    % Normalized passband edge frequency
Ws=2*fs/Fs;    % Normalized stopband edge frequency

[N,Wn]=buttord(Wp,Ws,Rp,Rs);    % Filter order selection
[b,a]=butter(N,Wn);        % Butterworth filter design
[Z,P,K]=tf2zp(b,a);        % Transfer function to zero-pole
[sos,G]=zp2sos(Z,P,K);        % Zero-pole to 2nd-order section

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

Код

sos =

    1.0000    2.0000    1.0000    1.0000   -0.8659    0.2139
    1.0000    2.0000    1.0000    1.0000   -1.1249    0.5770

Коэффициенты нерекурсивной части фильтр почему-то округлились до целых чисел.
С матлабом работаю недавно и ни как не соображу, что делаю не верно.

[Чиповод]

На самом деле они не совсем округлились. Это их матлаб округлил при выводе в формате по умолчанию.

А если вы наберете в командной строке

Код

format long

То результат будет такой.

Код

sos =

  Columns 1 through 4

   1.000000000000000   2.000000057757208   1.000000007739223   1.000000000000000
   1.000000000000000   1.999999942242791   0.999999992260774   1.000000000000000

  Columns 5 through 6

  -0.865894897828775   0.213946368900228
  -1.124874505762796   0.577024330739516

Он правильный?

[Mad_max]

Не могу прямо сейчас проверить, но как-то результат больше чем 2.0 настораживает.
До Q13 масштабировать что ли?