Никто не запрещает сделать CIC 12го порядка но
1) СIC дециматор в 10 раз 12го порядка будет уметь усиление 10^12 это потребует максимальную разрядность
в фильтре = число_разрядов_АЦП + 40
2) Полоса пропускания ну очень кривая будет (представляете sinc(x)^12 )
3) Почитайте описание микросхем DDC , врядли там найдете СIC больше 5го порядка.

p.s. Для Вашей задачи лучше всего FIR будет или CIC + FIR

интеграторы неправильно сделали, нарисуйте реализацию сic фильтра в симулинке и наложите ее на свое железо. увидите разницу %)

и разрядность у вас не правильно выбрана, работать не будет. учите матчасть.


+1. Тем более что полоса фиксирована и всё равно будете ставить компенсационный FIR %)

Вам подойдет и 11-го порядка CIC - неравномерность на краю полосы (на 8 МГц) -8.21дБ, максимальная помеха от децимации на краю полосы -108дБ.

Прошу прощения - почему-то сначала увидел 100 дБ максимальной помехи. Действительно, только 12-ый порядок подойдет.


Сообщение отредактировал 729 - Aug 13 2010, 08:39

Обратите пристальное внимание на то, как вы делаете округление. Меня всегда ошибки преследовали именно в этом.

Мы не делали округление, а просто отбрасывали, что вело к появлению некоторой постоянки. Но её легко замерить, а потом вычесть - работает вполне прилично.

C замечанием согласен. Исправлю, получившийся код выложу! Про разрядность - взял наобум. Потом планировал отрезать совсем другие разряды. Спасибо!

Вы глубоко заблуждаетесь, CIC фильтр эквивалентен N каскадам КИХ фильтров с прямоугольной импульсной характеристикой, соединенных последовательно (разница только в физической реализации). А суммарная импульсная характеристика фильтра определяется как свертка всех импульсных характеристик последовательно соединенных каскадов. Для примера, импульсная характеристика 2-х каскадов - равнобедренный треугольник, 3-х каскадов - параболический вид (сложно в уме прикинуть) и т.д.

Хотелось бы прояснить для себя следующий момент касаемо в части реализации CIC-фильтров. К примеру есть 3 каскада интеграторов + гребенчатых фильтров, т.е. 3-х каскадный CIC-фильтр. При этом разрядность входных данных больше чем выходных (например in_data = 20 бит, out_data = 12 бит). Допускается ли в этом случае отбрасывать младшие биты в интеграторах согласно методики Хогенауера (Hogenauer, Eugene. An Economical Class of Digital Filters For Decimation and Interpolation, IEEE Transactions on Acoustics, Speech and Signal Processing, Vol. ASSP-29, pp. 155-162, April 1981.). Согласно алгоритму в данной статье, фильтр работает при разрядности входных данных меньшей чем разрядность выходных. В противном случае моделирование показывает что фильтр работает неправильно! Есть соображение, что когда разрядность входных данных больше чем выходных, отбрасывать младшие биты в интеграторах вообще нельзя. Так ли это, или я ошибаюсь?
Заранее благодарен за ответы.

нет, простое "обрезание" разрядности не подойдет

Возьмите промоделируйте - смотрите на выход СIC фильтра (полный выход со всеми разрядами) - уменьшите масштаб при моделировании (чтобы много выходных данных на экране видели) - смотрите какие биты меняются, а какие нет или практически не меняются.
В итоге берите только те которые меняются - значащие биты в необходимом количестве (чаще делают разрядность входа и выхода равными).
Если СIC фильтр перестраиваемый, то такую операцию нужно делать для каждого коэффициента. Можно взять несколько точек, а потом аппроксимировать...

Посмотрите как меняются входные данные и какие значащие биты там, постарайтесь чтобы выходные данные с СIC фильтра были аналогичными...

PS СIC фильтр делайте полный без "обрезания" на каскадах - совет
PS PS прошу прощения - ошибся - говорил про децимирующий СIC фильтр, но думаю в Вашем случае (интерполяции) все аналогично с точность до наоборот...

Ну делать подбор усечения в каждом каскаде с последующим моделированием это слишком расточительно по времени. Я спрашивал про методику усечения описываемую Хогенауером. Согласно этой методики, можно заранее вычислить сколько бит и в каком каскаде можно отбрасывать с минимальным ущербом для конечного результата.
Другое дело что фильтр перестает правильно работать, когда разрядность входных данных больше чем выходных.
К примеру: есть два трехкаскадных фильтра: 1)N=3 кол-во секций, R=32 коэфФ. прореживания, М=2 разностная задержка. Bin= 8 бит, Bout =10 бит.
2)N=3 кол-во секций, R=32 коэфФ. прореживания, М=2 разностная задержка. Bin= 18 бит, Bout =10 бит.
Согласно приведенной в статье методики: для первого фильтра имеем следующие разрядности в промежуточных каскадах: 1 каскад = 26 бит, 2-ой =21 бит,
3-й = 16 бит, 4-й = 14 бит, 5-й = 13, 6-й = 12 бит.
для второго фильтра имеем такие же результаты, но первый фильтр работает правильно, второй нет.
Возникает вопрос ПОЧЕМУ? Для первого фильтра максимальная разрядность 1-го интегратора вычисляется по формуле Bint1=Bin + N*log2(R*M) = 8 + 3*log2(64)= 26 бит.
Для второго 20 + 3*log2(64) = 38 бит из которых мы учитываем только 26 старших (согласно Хогенауеру). Так вот я хотел спросить у знающих
людей: можно ли пользоваться данной методикой когда входная разрядность больше выходной. У Altera ядро работает корректно в обоих случаях. У меня нет.

нет Вы меня не поняли...



а потом уже применяете только ко входу СIC фильтра (потому что у Вас интерполирующий) то, что я советовал ранее...

Выбирать значащие биты Вы будете делать столько раз - сколько каскадов фильтров... Никуда не денитесь...

Ваш метод понятен. Но во первых он основан на эмпирических подходах, что далеко не всегда удобно, а во вторых я все же спрашивал про метод описанный в статье.
У меня есть подозрение что я все таки что то не так: реализую, понимаю или может что-то упустил. Но все равно спасибо!
Фильтр децимирующий.

Сообщение отредактировал PavPro - Jan 15 2015, 20:34

тут смотрели? я по Хогенауэру делал, но не помню, в каких условиях проверял