фильтр Баттерворта

[addddd]

Рассчитать реакции ФНЧ Баттерворта n-го порядка на белый шум.

Помоги те разобраться!!!

Я расчитал фильтр 8 порядка, это раз, нашел его ачх.
Как я знаю для нахождения спектра выходного сигнала нужно ЧХ поделить на спектр входного сигнала, но как задать белый шум функцией незнаю, и незнаю как вычислить спектры в дельфи.

[addddd]

repeat

h1:=(w0/sqrt(sqr(w0)+sqr(w*w*w0*0.195)));
h2:=(sqr(w0)/sqrt(sqr(sqr(w0)-sqr(w))+sqr(2*w*w0*0.556)));
h3:=(sqr(w0)/sqrt(sqr(sqr(w0)-sqr(w))+sqr(2*w*w0*0.831)));
h4:=(sqr(w0)/sqrt(sqr(sqr(w0)-sqr(w))+sqr(2*w*w0*0.98)));
h:=h1*h2*h3*h4;

Series1.AddXY(w,h1,'',clTeeColor);
Series2.AddXY(w,h2,'',clTeeColor);
Series3.AddXY(w,h3,'',clTeeColor);
Series4.AddXY(w,h4,'',clTeeColor);
Series5.AddXY(w,h,'',clTeeColor);

w:=w+0.001;
until w>ww;

[KobolD]

Белый шум задается случайными величинами (random), ты же входной сигнал всеравно квантовать будешь. А для нахождения спектра сигнала можно воспользоваться дискретным преобразованием фурье. Кстати гдето на форуме была тема про него.
Могу пример ДПФ на C# кинуть.

[addddd]

Дайте на код взглянуть пожалусто

[addddd]

тоесть нужно прямое преобразование фурье над самим сигналом???? на бумаге легко а в коде не получатся...

[KobolD]

        const int N = 50; //количество разбиений
        double[] InnerArray = new double[N];//Массив входных данных
        Complex[] Furie = new Complex[N];//Массив значений ДПФ

        public class Complex //Класс комплексных чисел
        {
            public double Re;//Реальная часть
            public double Im;//Мнимая
        }

        private void main(object sender, EventArgs e)
        {
            for (int i = 0; i < N; i++)
                InnerArray[i] =Math.Cos(2 * Math.PI * i / N);//задаем форму сигнала
            DPF(InnerArray);
        }

        void DPF(double[] Inner)//Прямое преобразование фурье
        {
            int N = Inner.Length;
            double Arg;

            for (int k = 0; k < N; k++)
            {
                Furie[k] = new Complex();
                for (int n = 0; n < Inner.Length; n++)
                {
                    Arg = 2 * Math.PI * k * n / N;
                    Furie[k].Re += Inner[n] * Math.Cos(Arg);
                    Furie[k].Im -= Inner[n] * Math.Sin(Arg);
                }

            }
        }

Вот тут немного теории плюс формулы как посчитать амплитуду и фазу для АЧХ

[MypkAnn]

А как тогда будет выглядеть функция обратного преобразования Фурье что бы можно было обратно в массив double поместить значения и вывести полученный график?