Быстрое преобразование Фурье

@sinsin · Регистрация: 21.01.2010

Студворк — интернет-сервис помощи студентам

Доброго времени суток.
Пытаюсь реализовать на С++ создание цифрового фильтра методом свертки, при этом использую БПФ, реализованную на http://algolist.manual.ru/maths/fft.php, но получаю странный результат:
после прогона заданной частотной характеристики фильтра через обратное БПФ(пока использую FFT_T) в выходном массиве появляеются ненулевые комплексные составляющие.

C++

  int N = FDiscr*2;
    Complex *FHar = new Complex[N];
    double *tmp_arr = new double[N];
    Complex tmp; 
    // частотная характеристика
    for (int i = 0;i<FDiscr;i++)
    {
        FHar[i].i = 0;
        if ((i<=f2)&&(i>=f1))
            FHar[i].r = 1;
        else
            FHar[i].r = 0;
    }
    FFTReOrder(FHar,FDiscr);
    // обратное бпф
    FFT_T(FHar,FDiscr,-1);

Как это можно побороть?
Чтобы насчет правильности действий не возникалли вопросы, синтез фильтра в mathcad прикреплен в Filter_exper.rar.

@vital792 · 27.02.2011, 23:11

Сообщение от sinsin

после прогона заданной частотной характеристики фильтра через обратное БПФ(пока использую FFT_T) в выходном массиве появляеются ненулевые комплексные составляющие.

ниче не понял, какие еще ненулевые комплексные составляющие? После ifft вы должны были получить импульсную характеристику фильтра и не в комплексном виде(мнимую часть просто можно отбросить).
И естественно она не будет нулевой. Судя по коду вы пытались сделать полосовой фильтр с полосой пропускания от f1 до f2. Посмотрите в учебнике импульсную характеристику полосового фильтра и сравните со своей. Я думаю не будет сильно отличаться

@sinsin · 28.02.2011, 07:10 **[ТС]**

Проблема решена, см. [ссылка удалена]

Комментарий модератора

Ссылки на другие форумы запрещены.
Перечитайте правила форума.
Скопируйте решение сюда

@vital792 2013 / 1285 / 61 Регистрация: 05.06.2010 Сообщений: 2,213
	27.02.2011, 23:11	2
	Сообщение от sinsin после прогона заданной частотной характеристики фильтра через обратное БПФ(пока использую FFT_T) в выходном массиве появляеются ненулевые комплексные составляющие. ниче не понял, какие еще ненулевые комплексные составляющие? После ifft вы должны были получить импульсную характеристику фильтра и не в комплексном виде(мнимую часть просто можно отбросить). И естественно она не будет нулевой. Судя по коду вы пытались сделать полосовой фильтр с полосой пропускания от f1 до f2. Посмотрите в учебнике импульсную характеристику полосового фильтра и сравните со своей. Я думаю не будет сильно отличаться 0

Опции темы