Re: общие вопросы по ЦФ
Добавлено: 23 ноя 2011, 16:55
удалить
Теория и практика цифровой обработки сигналов
http://ru.dsplib.org/forum/
Код: Выделить всё
f=0
q=0
while q<=n-1:
f+=sn[q]*np.cos(xa*(q))-sn[q]*np.sin(xa*q)*1j
q+=1
Если зеленая это мнимая часть то реальная часть должна совпадать с модулем (красный с синим) ищите ошибку. Физически фильтр с нулевой фазочастотной характеристикой (нулевая мнимая часть H(w)) реализовать нельзяeiyawii писал(а):
на которой чётко видна жирная зеленая линия, которая на всей оси частот равна нулю - это мнимая ЧХ(!), то есть мнимая часть равна нулю. Синяя кривая - это действительная ЧХ, точнее её модуль. Красная кривая - АЧХ. Теперь всё-таки если мнимая часть равна нулю, то она не должна вносить вклад в АЧХ, но если строить без неё, то получается не совсем то. Из этого следует, что где-то, возможно, ошибка.
Бахурин Сергей, знаете ли КИХ фильтры, у которых мнимая часть ЧХ равна нулю? И не может ли в си
ну я найду этому применение есть задачи довольно прагматичные, которым не помешал бы такой фильтр.Бахурин Сергей писал(а):вы представляете себе что такое 500 дБ? Зачем вам такое подавление? Я думаю что 500 дБ на сегодняшний день это не реально
300 дБ я получил, в результате не сомневаюсь. Думаю с этим и связано осыпание АЧХ - если бы не ошибки округления и подобные вещи, то можно и в 400 залезть. Если я дам вам коэффициенты фильтра, можете проверить? Или скажите как проверить в каком-нибудь матлабе, зная коэффициенты фильтра(это даже лучше будет)?ivan219 писал(а):500 дБ это из области фантастики вы что бозон Хиггса собрались искать Тут даже 300 дБ получить с числами двойной точности а это 64 бит проблематично из за ошибки округления.
только вчера об этом думал, попробую.ivan219 писал(а):Но если нужно то поступите так. Создайте фильтр с желаемыми вами характеристиками и сверните его самим собой несколько раз и вы получите желаемое подавление в 500 дБ
ну это был сарказм над самим собой, по поводу мнимой характеристики) Ошибка состояла в том, что на вещественных осях нельзя вывести комплексное число, и поэтому мнимая ЧХ совпадала с осью абцисс, типа нулю была равна. Поэтому взял абсолютное значение мнимой ЧХ.Бахурин Сергей писал(а):Если зеленая это мнимая часть то реальная часть должна совпадать с модулем (красный с синим) ищите ошибку. Физически фильтр с нулевой фазочастотной характеристикой (нулевая мнимая часть H(w)) реализовать нельзя
Вот это я не понял.Бахурин Сергей писал(а):Физически фильтр с нулевой фазочастотной характеристикой (нулевая мнимая часть H(w)) реализовать нельзя
Не знаю на счёт матлаба не пользуюсь им.eiyawii писал(а): 300 дБ я получил, в результате не сомневаюсь. Думаю с этим и связано осыпание АЧХ - если бы не ошибки округления и подобные вещи, то можно и в 400 залезть. Если я дам вам коэффициенты фильтра, можете проверить? Или скажите как проверить в каком-нибудь матлабе, зная коэффициенты фильтра(это даже лучше будет)?