Запуск библиотеки FFT C#

[Krasi]

Решил я воспользоваться этой библиотекой, хочу сначала просто сделать БПФ на любое число выборок любого сигнала.

http://www.fftw.org/install/windows.html
http://www.fftw.org/fftw2_doc/fftw_2.html

Я знаю, как подключить к проекту dll и добавить h-файл. Теперь мне надо понять, как будет выглядеть основной код для этой задачи.

У меня к Вам 2 вопроса.
1. Я вообще программирую с помощью Windows Forms, поэтому в идеале я жду простейший вариант кода, либо подскажите хоть как-то, чтобы разобраться, что писать.

2. Я посмотрел примеры, ориентируюсь на ссылки выше, но нашел только код

#include <fftw.h>
...
{
     fftw_complex in[N], out[N];
     fftw_plan p;
     ...
     p = fftw_create_plan(N, FFTW_FORWARD, FFTW_ESTIMATE);
     ...
     fftw_one(p, in, out);
     ...
     fftw_destroy_plan(p);  
}

Может быть, Вы знаете, как сделать хотя бы консольное приложение, возможно на основе этого кода, которое бы могло из одного массива данных с помощью БПФ получать другой массив?

[Krasi]

Нашел, как это уже делали. http://ru.stackoverflow.com/question...80%D1%8C%D0%B5

Проблема в том, что я запустил этот проект, понял, что в fin входные данные, в fout - выходные. Я решил перенести это в свой проект. Все заступорилось с подключения библиотеки. Смотрите, я добавил библиотеку fftwlib.dll в проект, написал "embedded resource", затем ищу в Reference fftwlib и не могу найти. В проекте пишу строчку

using fftwlib;

, так компилятор не может понять, откуда взять fftwlib. Подскажите, как я могу хотя бы успешно подключить эту библиотеку к себе в проект? Почему может не получаться?

[Stilet]

Насколько я помню, сборки подключаются так:
Безымянный.jpg

[p51x]

1. Лично я вам несколько раз, что не надо сборки добавлять как файлы в проект или как ресурсы. Используйте меню для референсов.
2. Либа с http://www.fftw.org - это действительно просто либа, не сборка
3. Вы в первом посте ссылаетесь на нативную либу, а во втором сетуете, что взяли код из примера использования обертки либы, а обертку не взяли

Может вам начать с чего-нибудь попроще?..

[Krasi]

Я долго мучался, и в итоге библиотеку он видит. Через референсы и добавляю, но проблема в том, что когда я переносил dll в корень проекта, среда Visual его не видела в референсах. Проблема была в том, что я хотел добавить через Reference, но не мог найти библиотеки, которую надо добавить. Ну ок.
Посмотрите, пожалуйста, пример, я руководствовался ссылкой, которую указал выше, но я не знаю, как из консольного приложения переписать windows forms. Что мне надо изменить в исходном коде, чтобы у себя не иметь ошибок и тп. Когда я, например, скопипастил класс public class FFTWtest перед стандартным классом windows forms, начало писать ошибки по поводу дальнейшей инициализации windows forms.

Я просто еще не соображаю, где брать классы и примеры, чтобы использовать функции из dll. На математическом уровне-то все ясно - мне сейчас надо сделать просто входной и выходной массивы, на входе пилу, например, подать входной сигнал на обработку библиотеке FFTW, и результат записать в выходной массив. А вот подробнее - я вообще не знаю, с чего начать, с каких функций! Даже этот пример меня не очень пока вдохновил, не могу продвинуться http://forum.sources.ru/index.php?showtopic=368270

Следующий код я переосмыслил и записал к себе в приложение Windows Forms

int n = 16348; float [] fin = new float[n * 2]; float [] fout = new float[n * 2];
                for (i = 0; i < n * 2; i++) fin[i] = i % 50;
                for (i = 0; i < n * 2; i++) fout[i] = i % 50;
                IntPtr pin = fftwf.malloc(n * 8); IntPtr pout = fftwf.malloc(n * 8);
                
                IntPtr fplan1 = fftwf.dft_1d(n, pin, pout, fftw_direction.Forward, fftw_flags.Estimate);
                
                fftwf.execute(fplan1);
                double a = 0, b = 0, c = 0;
                fftwf.flops(fplan1, ref a, ref b, ref c);
                Marshal.Copy(fin, 0, pin, n * 2); Marshal.Copy(fout, 0, pout, n * 2);

Проблема в том, что возникает Exception для строчки

IntPtr pin = fftwf.malloc(n * 8); IntPtr pout = fftwf.malloc(n * 8);

Почему это может происходить?

[p51x]

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

Я долго мучался, и в итоге библиотеку он видит. Через референсы и добавляю, но проблема в том, что когда я переносил dll в корень проекта, среда Visual его не видела в референсах. Проблема была в том, что я хотел добавить через Reference, но не мог найти библиотеки, которую надо добавить. Ну ок.

Вы читаете, что вам пишут? Это не СБОРКА, это НАТИВНАЯ ЛИБА.

Посмотрите на https://github.com/tszalay/FFTWSharp...arp/Interop.cs
Обратите внимание на всякие там DllImport...

[Krasi]

ООО, я просто 2 dll добавил, которые в вашей сборке были указаны, и у меня перестало выдавать ошибку!! Серьезно что ли? А как правильно тогда вообще эту сборку в проект добавлять? Получается, Reference здесь нужна не на dll, а на проект FFTWSharp, который я должен добавить в свой проект? Кстати, здесь попроще код http://ru.stackoverflow.com/question...81-idisposable

[Krasi]

Уже ближе, тогда что у меня в коде может быть неверно сейчас? Это от начала до конца операция БПФ, и исходные данные, и выходные можете видеть.

IntPtr pin, pout;
                int n = 16348; float [] fin = new float[n * 2]; float [] fout = new float[n * 2];
                
                pin = fftwf.malloc(n * 8); pout = fftwf.malloc(n * 8);
                fin = new float[n * 2]; fout = new float[n * 2];

                for (i = 0; i < n * 2; i++)
                    fin[i] = i % 50;

                
                //copy managed arrays to unmanaged arrays
                Marshal.Copy(fin, 0, pin, n * 2);

 
                IntPtr fplan1 = fftwf.dft_1d(n, pin, pout, fftw_direction.Forward, fftw_flags.Estimate);
                
                fftwf.execute(fplan1);
                double a = 0, b = 0, c = 0;
                fftwf.flops(fplan1, ref a, ref b, ref c);

                this.Controls["richTextBox2"].Text = fin[80].ToString();

                fftwf.free(pin);
                fftwf.free(pout);
                fftwf.destroy_plan(fplan1);

fin, вроде, пила, но вот хочу я вместо fin вывести fout - выдает 0, а по-логике у меня результат должен быть в pout. Может что-то не дописал? Массив fout просто заполнен нулями.

[Krasi]

Во, вот с этой строчкой заработало

Marshal.Copy(pout, fout, 0, n * 2);

[Krasi]

Теперь вопрос попроще. Чтобы протестировать работу этого кода, я хочу использовать гармонический сигнал. И тут я забыл, что синус выдает double, а мне надо float.

fin[i * 2] = 5 * Math.Sin(2 * Math.PI * 10 * i / 100);

Что я могу сделать?