Отправлено ВН 20 февраля 2003 г. 13:12 В ответ на: Блин совсем запутался. отправлено Newka 20 февраля 2003 г. 11:46

Не видно команды SSBX FRCT, но это не главное.
Похоже Вам нужно сменить знаки у y_c1 и y_c2. Для y_c1, точнее, использовать mas вместо mac.
Я не пользуюсь матлабом, не знаю, что он генерирует при расчете к-тов, но похоже y_c1, y_c2 даны с теми знаками, с которыми они входят в H(z), а не в разностное уравнение.
И вот проверка.
Для упрощения из x_c0-x_c2 оставлю только x_c0=1.
К-ты y_c0,y_c1,y_c2 в формате Q15 равны 1,-1,0.69
Фильтр описывается разностным уравнением:
y(n)=x(n)+y_c1*y(n-1)+y_c2*y(n-2). И, соответстенно, передаточной функцией:
H(Z)=1/(1-y_c1*Z^(-1)-y_c2*Z^(-2)). И при данных y_c1,y_c2:
H(Z)=1/(1+Z^(-1)-0.69*Z^(-2)).
Наиболее распространенный случай ФНЧ второго порядка - 2 комплексно-сопряженных полюса.
В общем виде знаменатель H(z), для такого ФНЧ, имеет вид:
(1-b*exp(j*W0)*Z^(-1))*(1-b*exp(-j*W0)*Z^(-1)). Т.е. 2 полюса, комплексно-сопряженных: b*exp(j*W0) и b*exp(-j*W0). b - "модуль" полюса, положительно и меньше 1. Определяет добротность полюса.
W0 - цифровая частота. Т.е. знаменатель H(z):
1-2*b*cos(W0)*Z^(-1)+b^(2)*Z^(-2). Пользуясь Вашими обозначениями:
y_c2=-b^(2), y_c(1)=2*b*cos(W0). Поскольку b^(2) больше 0, то y_c2 должно быть меньше 0. У Вас наоборот.
Про ssbx frct только не забудьте.

• Спасибо за ответ — Newka   (20.02.2003 14:36, 2264 байт)
  • Ответ: — ВН   (20.02.2003 16:45, 4791 байт)
    • Спасибо огромное Вам BH! — Newka   (20.02.2003 19:16, 381 байт)
      • А голова зачем тогда, если матлаб все сам делает. BH у нас профессор, только он оценки не ставит :-) — AntZ   (20.02.2003 19:27, пустое)