|
|
 |
Ответов
|
Jan 26 2006, 07:52
|
Частый гость
 
Группа: Свой
Сообщений: 140
Регистрация: 18-10-05
Пользователь №: 9 792

|
Еще раз о коэффициентах фильтра.
Если рассчитанные коэффициенты: A1=1.584584 A2=-0.987981 B0=-B2=0.006009
То для рассчетных формул вида:
1. D0=x[t]+a1*D1[t]+a2*D2[t] 2. y[t]=b0*(D0-D2)
где D1[t]=D0[t-1], D2[t]=D1[t-1], x[t] - входной сигнал, y[t] - выходной.
Коэффициенты в целочисленных процессорах коэффициенты обычно сдвигают так, чтобы максимальный из них занимал все значащие разряды.
В нашем случае, для 16-разрядной арифметики сдвиг на 14 разрядов,
a1=25961 (6569h) a2=-16187 (С0С5h) b0=98 (62h) конечный результат надо сдвинуть на 14 разрядов вправо
Надо еще внимательно посмотреть на разрядность D0,D1,D2, чтобы не было переполнений. В любом cлучае, результат умножения a1*D1 и a2*D2 - как минимум 32-разрядный, для дальнейших вычислений младшие 16 бит отбрасываются.
Возможно, для оптимизации разрядности промежуточных результатов, сдвигать надо не на 14 бит, а на 13 или 12 (соответственно рассчитанные целочисленные коэффициенты сдвинуть на 1 или 2 бита вправо).
Сообщение отредактировал Johny - Jan 26 2006, 08:18
|
|
|
|
|
Jan 26 2006, 08:34
|
Частый гость
 
Группа: Свой
Сообщений: 140
Регистрация: 18-10-05
Пользователь №: 9 792

|
Цитата(Johny @ Jan 26 2006, 10:52)  Надо еще внимательно посмотреть на разрядность D0,D1,D2, чтобы не было переполнений. В любом cлучае, результат умножения a1*D1 и a2*D2 - как минимум 32-разрядный, для дальнейших вычислений младшие 16 бит отбрасываются. Почему-то не смог отредактировать. У результата умножения a1*D1 и a2*D2 не отбрасываются младшие 16 бит, он сдвигается вправо на 14 бит. Кстати, надо следить за разрядностью входного сигнала. Если хотите, чтобы сумматоры были 16-разрядными, входной сигнал навскидку должен быть не более 14 разрядов. Иначе надо переходить на более высокоразрядную арифметику (например на 32-разрядную)
|
|
|
|
Сообщений в этой теме
ASV Цифровой фильтр Dec 8 2005, 20:08 mse Не всё так просто. ЦФ бывают, как минимум, FIR и I... Dec 8 2005, 21:47 Tiro ASV
По английски читаешь? Посмотри здесь, наприме... Dec 8 2005, 22:12 Dr.NoA Цитата(ASV @ Dec 8 2005, 23:08) Вопрос: и... Dec 9 2005, 09:22 mse ЦитатаКак раз все 10 отсчетов и надо пропускать че... Dec 9 2005, 19:58 GrayCat Цитата(ASV @ Dec 8 2005, 22:08) На скольк... Dec 12 2005, 08:58 aal В http://vlab.netsys.ru/forum/showthread.php?t=388... Dec 15 2005, 12:28 nano Крик о помощи
Начальные условия : топик прочел... Dec 15 2005, 16:41 ASV У меня задача следующая. Мне нужно померить напряж... Dec 20 2005, 21:19 mse Цитатаи умеет быстро умножать только 8bit
Ну не ст... Dec 20 2005, 21:43 ASV Цитата(mse @ Dec 21 2005, 01:43) Цитата и... Dec 20 2005, 21:48 mse Цитатаа не в строчку как это делают DSPIC-и.
Хех. ... Dec 20 2005, 21:54 ASV Цитата(mse @ Dec 21 2005, 01:54) Цитатаа ... Dec 21 2005, 22:46 urasinov nano
В MATLAB в DSP Blockset имеется Filter Realiz... Dec 22 2005, 07:17 Sokol Цитата(nano @ Dec 15 2005, 19:41) Крик о ... Dec 22 2005, 12:34 skopus Уважаемые, скажите пожалуйста насколько сильно мож... Jan 23 2006, 16:52 BVU Цитата(skopus @ Jan 23 2006, 19:52) Уважа... Jan 25 2006, 08:06 Dr.NoA Я правильно понял проблему? Вы подаете синус с амл... Jan 23 2006, 20:54 skopus нет. Вы не совсем правильно поняли. Фильтры не сое... Jan 24 2006, 05:02 skopus поисследовал ситуацию поподробнее и обнаружил, что... Jan 24 2006, 08:59 Dr.NoA Какие-то странные графики. Почему, например, для ф... Jan 24 2006, 10:54 skopus фильтр на 316гц по идее должен пропускать частоту ... Jan 24 2006, 15:08 Johny Цитата(skopus @ Jan 24 2006, 18:08) для ф... Jan 25 2006, 16:25 Dr.NoA Я бы все таки проверил характеристики фильтра в Ma... Jan 24 2006, 19:54 skopus я уже неоднократно проверял. И в матлабе, и в про... Jan 25 2006, 07:06 skopus я не говорил что у меня НЧ фильтры. У меня узкопол... Jan 25 2006, 08:34 BVU Для выяснения причин данной проблемы нужно подать ... Jan 25 2006, 10:20 skopus Цитата(BVU @ Jan 25 2006, 13:20) Для выяс... Jan 25 2006, 14:54  BVU Цитата(skopus @ Jan 25 2006, 17:54) Цитат... Jan 25 2006, 15:46   Stanislav Цитата(BVU @ Jan 25 2006, 18:46) ...Судя ... Jan 25 2006, 16:05  Johny Цитата(skopus @ Jan 25 2006, 17:54) http:... Jan 25 2006, 16:56 michael34 Вряд ли ты прав... Хотя точно сказать не могу... Jan 25 2006, 11:11 BVU Цитата(michael34 @ Jan 25 2006, 14:11) Вр... Jan 25 2006, 13:34 Dr.NoA Как я понял, при расчете в таком виде
Цитата// 1. ... Jan 26 2006, 08:44 skopus Спасибо за совет. Но удивительное в том, что перво... Jan 26 2006, 11:17 Dr.NoA Наоборот, при такой записи уравнения
Цитатаx[kT]+a... Jan 26 2006, 13:12 skopus так в том-то и дело что этих b0 и b2 в уравнениях ... Jan 26 2006, 14:26 Dr.NoA Цитата(skopus @ Jan 26 2006, 17:26) так в... Jan 26 2006, 16:14
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|