Простейший БИХ-фильтр, Не пойму как коэффиценты сделать целочисельными Опции

[condor]

Надо на 16-разрядном контроллере сделать ФНЧ.
Пытаюсь в МАТЛАБе смоделировать хотя бы фильтр второго порядка в виде канонического звена.
Коэффициенты получил с помощью fdatool.
С дробными коэффициентами все правильно работает, а как сделать их челочисельными?
Везде в литературе об этом вскользь упоминается как о само собой разумеющемся.
Если их просто умножить на 2^14, то из-за рекурсии при вычислении
w(n) = x(n) - a1*w(n-1) - a2*w(n-2)
w(n) быстро уходит в бесконечность. (На всякий случай: w - это промежуточное значение в каноническом звена; x - входной сигнал; ai - коэф.).
Подскажите, а то я сегодня туго соображаю

[condor]

Я все это делал.
Квантовал в fdatool.
Написано Stable: Yes

Quantized Numerator:
1
1.99591064453125
1
Quantized Denominator:
1
0.06439208984375
0.3143310546875
Quantized Gain:
0.3074951171875

С этими квантованными коэффициентами все работает отлично.
Но если их умножить на 2^14, то получим.
a1 = 0,0644*2^14 = 1055;
a2 = 0,3143*2^14 = 5150;
b1 = 1.99591*2^14 = 32701;
b2 = 1*2^14 = 16384;
sk = 0,3075*2^14 = 5038
Такие же коэффициенты генерятся в Сишном хедере.

w(n) = sk*x(n) - a1*w(n-1) - a2*w(n-2) ;
w(n) = 5038*x(n) - 1055*w(n-1) - 5150*w(n-2);
Начальные w я принял равными 0, а x(n) - какое-нибудь постоянное значение.
Теперь даже проделав вручную нескольно итераций видно, что w(n) нестабильно.
Думал что в формуле опечатка, но ведь с дробными кавантованными все работает, и даже АЧХ примерно совпадает.

Вот МАТЛАБовский код. Может я где-то ошибся, но кажется тут особо не где.

Nsamp = 2000;
wk = 0;
w1 = 1;
w2 = 1;
sk = 0.3075 * 2^14;
a1 = 0.0644 * 2^14;
a2 = 0.3143 * 2^14;
b0 = 1;
b1 = 1.996;
b2 = 1;

freq = 10;
t = 0:0.0005:1;
t (end) = [];

x = sin (2*pi*freq*t);

for n = 1:Nsamp
xn = x(n);
wk = sk*x(n) - a1*w1 - a2*w2;
y(n) = b0*wk + b1*w1 + b2*w2;
w2 = w1;
w1 = wk;
end