|
|
  |
быстрый FIR для комплексного сигнала |
|
|
|
Apr 14 2008, 08:04
|
Участник

Группа: Новичок
Сообщений: 37
Регистрация: 30-07-07
Пользователь №: 29 445

|
Цитата(Михаил_K @ Apr 14 2008, 09:31)  Поясните плиз. Эффект Гиббса возникает при наглом усечении ИХ линейной системы. Какое отношение это имеет при фильтрации в частотной области, где просто зануляются ненужные частоты. У нас есть средства программной обработки сигналов, где выполняется фильтрация в частотной области, и значения АЧХ фильтров +1 и 0. дело в том что при наглом занулении фурье-коэф в частотной области и обратном преобразовании фурье получаются артефакты в сигнале. а т.к. этот сигнал в конечном счете переносится в звуковой спектр, то это все приходится слушать ушами. поэтому и такое сложное формирование ачх, а так же применение overlap&add. странная дискуссия получается - я задал вопрос, по существу никто ничего не сказал, зато все спрашивают.
|
|
|
|
|
Apr 14 2008, 08:14
|
Гуру
     
Группа: Свой
Сообщений: 3 106
Регистрация: 18-04-05
Пользователь №: 4 261

|
Цитата(Relayer @ Apr 14 2008, 11:04)  странная дискуссия получается - я задал вопрос, по существу никто ничего не сказал, зато все спрашивают. А что свойтво линейности фильтров уже отменили? Т.е., буквально, фильтруем отдельно действительную часть сигнала, затем комплексную, складываем и получаем комплексный результат. Фильтр FIR, ессно, один и тот же, что для Re, что для Im.
|
|
|
|
|
Apr 14 2008, 08:24
|
Участник

Группа: Новичок
Сообщений: 37
Регистрация: 30-07-07
Пользователь №: 29 445

|
Stanislav, я не намерен отвечать на ваше пустословие в этой теме. вы не написали ничего по существу вопроса. общие фразы и только. вдобавок вы позволяете себе хамское поведение. если вам не понятно - я не собираюсь вступать с вами тут в бессмысленную перепалку. хотите говорить - говорите. до тех пор пока модератор не сочтет ваши разговоры раздуванием флейма Цитата(blackfin @ Apr 14 2008, 11:14)  Т.е., буквально, фильтруем отдельно действительную часть сигнала, затем комплексную, складываем и получаем комплексный результат. Фильтр FIR, ессно, один и тот же, что для Re, что для Im.  это все применимо к симметричным относительно 0Hz фильтрам. и еще один момент - мне фильтровать надо именно в частотной области, без перехода во временную. дело в том что перед и после фильтрации присутствует еще некоторая обработка в частотной области, поэтому лишние "телодвижения" туда-сюда из частотной во временную и обратно скажутся на производительности системы.
|
|
|
|
|
Apr 14 2008, 09:12
|
Гуру
     
Группа: Свой
Сообщений: 3 106
Регистрация: 18-04-05
Пользователь №: 4 261

|
Цитата(Relayer @ Apr 14 2008, 11:24)  это все применимо к симметричным относительно 0Hz фильтрам. Здесь Вы, мягко говоря, ошибаетесь.. Сможете доказать, что у фильтров с несимметричной АЧХ отсутствуют линейность и стационарность.. ??? Цитата(Relayer @ Apr 14 2008, 11:24)  и еще один момент - мне фильтровать надо именно в частотной области, без перехода во временную. дело в том что перед и после фильтрации присутствует еще некоторая обработка в частотной области, поэтому лишние "телодвижения" туда-сюда из частотной во временную и обратно скажутся на производительности системы. Да я не навязывал Вам никакой конкретной реализации фильтра. Просто хотел лишний раз напомнить о линейности операции фильтрации в её математическом смысле..
|
|
|
|
|
Apr 14 2008, 09:44
|
.
     
Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753

|
Цитата(blackfin) Т.е., буквально, фильтруем отдельно действительную часть сигнала, затем комплексную, складываем и получаем комплексный результат. По одной части сигнала кажется невозможно отличить частоты +F и -F (симметрия от 0) т.к. они совпадают с противоположной, сдвинутой на 180 град. Так что фильтр должен работать сразу с комплексным представлением. Либо предворительно сдвинуть спектр в части, подлежащей фильтрации. Цитата(Relayer) дело в том что при наглом занулении фурье-коэф в частотной области и обратном преобразовании фурье получаются артефакты в сигнале. а т.к. этот сигнал в конечном счете переносится в звуковой спектр, то это все приходится слушать ушами. поэтому и такое сложное формирование ачх, а так же применение overlap&add. Из чистого любопытства, для чего понадобилась такая необычная (ИМХО) фильтрация?
Сообщение отредактировал GetSmart - Apr 14 2008, 09:45
--------------------
Заблуждаться - Ваше законное право :-)
|
|
|
|
|
Apr 14 2008, 10:24
|
Знающий
   
Группа: Свой
Сообщений: 552
Регистрация: 29-02-08
Пользователь №: 35 481

|
Не знаю, поможет или нет. Но мы для расчета комплексного фильтра (чтобы не было симметрии относительно 0) делали так. 1. Вычисляли обычный полосовой фильтр необходимой конфигурации в матлабе. При этом получался фильтр с действительными коэффициентами, который симметричен относительно 0. 2. Делали преобразование Гильберта над фильтром. При этом одна половинка исчезает.  .
|
|
|
|
|
Apr 14 2008, 12:51
|
Участник

Группа: Новичок
Сообщений: 37
Регистрация: 30-07-07
Пользователь №: 29 445

|
Цитата(GetSmart @ Apr 14 2008, 12:44)  По одной части сигнала кажется невозможно отличить частоты +F и -F (симметрия от 0) т.к. они совпадают с противоположной, сдвинутой на 180 град. Так что фильтр должен работать сразу с комплексным представлением. Либо предворительно сдвинуть спектр в части, подлежащей фильтрации. о чем я и говорил Цитата(GetSmart @ Apr 14 2008, 12:44)  Из чистого любопытства, для чего понадобилась такая необычная (ИМХО) фильтрация? вот для этого: http://code.google.com/p/sdrlabЦитата(Михаил_K @ Apr 14 2008, 13:24)  Не знаю, поможет или нет. Но мы для расчета комплексного фильтра (чтобы не было симметрии относительно 0) делали так. 1. Вычисляли обычный полосовой фильтр необходимой конфигурации в матлабе. При этом получался фильтр с действительными коэффициентами, который симметричен относительно 0. 2. Делали преобразование Гильберта над фильтром. При этом одна половинка исчезает.  . Михаил, спасибо. буду обязательно пробовать! за две страницы первый солюшн. насколько я понимаю такой методикой можно фильтр с вещественными коэфф отобразить только либо на положительные либо на отрицательные частоты. правильно? PS 2 fontp и blackfin - господа, вы наверное не совсем поняли - я тут никому ничего доказывать не собираюсь. высшее математическое есть, но по цос у меня самообразование. и по работе я специализируюсь в областях далеких от электроники и цос. а ветку открыл именно с целью чтобы профессионалы не зубоскалили и выпендривались друг перед другом, а помогли реальным советом
|
|
|
|
|
Apr 14 2008, 13:21
|
Знающий
   
Группа: Свой
Сообщений: 552
Регистрация: 29-02-08
Пользователь №: 35 481

|
Цитата(Relayer @ Apr 14 2008, 16:51)  насколько я понимаю такой методикой можно фильтр с вещественными коэфф отобразить только либо на положительные либо на отрицательные частоты. правильно? Ну вроде как да. В принципе можно делать и в лоб. Т.е. задаете требуемую АЧХ, делаете преобразование фурье комплексное, потом усекаете окном. Тут главное правильно задать наклон ФЧХ. Но это нужно делать ручками, матлаб к сожалению расчет комплексных фильтров не поддерживает. У способа с преобразованием Гильберта есть недостаток - увеличения порядка фильтра на порядок преобразователя гильберта  . Мы синтезировали большие фильры с числом коэффициентов более 200, поэтому +- 14 коэффициентов было не так страшно. В вашем случае не знаю как это будет.
|
|
|
|
|
Apr 14 2008, 15:03
|
Участник

Группа: Новичок
Сообщений: 37
Регистрация: 30-07-07
Пользователь №: 29 445

|
Цитата(Михаил_K @ Apr 14 2008, 16:21)  Мы синтезировали большие фильры с числом коэффициентов более 200, поэтому +- 14 коэффициентов было не так страшно. В вашем случае не знаю как это будет. у меня порядки 500 и выше. так что некритично
|
|
|
|
|
Apr 15 2008, 22:37
|

Гуру
     
Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987

|
Цитата(Relayer @ Apr 14 2008, 12:24)  Stanislav, я не намерен отвечать на ваше пустословие в этой теме. вы не написали ничего по существу вопроса. общие фразы и только. вдобавок вы позволяете себе хамское поведение. если вам не понятно - я не собираюсь вступать с вами тут в бессмысленную перепалку. хотите говорить - говорите. до тех пор пока модератор не сочтет ваши разговоры раздуванием флейма Ошибаетесь, уважаемый. Я дал здесь Вам целых два полезных совета по существу вопроса. Ваша нежелание/неспособность принять их к сведению нисколько не умаляет их сути, а Ваше реакция на добрые намерения напрочь отбивает всякую охоту помогать Вам хоть в чём-то впредь.  Цитата(Relayer @ Apr 14 2008, 16:51)  ...и по работе я специализируюсь в областях далеких от электроники и цос... Это заметно. Особенно в темах о компрессии сигнала.
--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
|
|
|
|
|
Apr 16 2008, 05:51
|
.
     
Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753

|
Цитата(Stanislav) Я дал здесь Вам целых два полезных совета по существу вопроса. Перечитал снова. Не нашёл ни первого ни второго  Может носом ткнёте? ЗЫ. Я не автор, но мне любопытно.
--------------------
Заблуждаться - Ваше законное право :-)
|
|
|
|
|
  |
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|