Как построить спектр сигнала, график которого в TChart?

[SantilaS]

Ребята, помогите пожалуйста. Давно не программировал((( Сколько не читаю информации, не получается. Задача такая: есть график сигнала в TChart, нужно построить его спектр, можно через преобразование Фурье, можно через вейвлет анализ. Подскажите какие-нибудь алгоритмы для выполнения этой задачи или киньте исходный код примера решения похожей задачи на avfelichev@mail.ru Заранее благодарен откликнувшимся!

[raxp]

http://pblog.ru/?p=658

[саша40]

Скажите только каких сигналов? Если как в плеере, то надо определять FFT-параметры канала. Функция Bass_ChannelGetData. Если же других, то в массивe Single, затем по цыклу, рисуете на канвасе PaintBox.

[SantilaS]

В приложении пример данных, с которыми я работаю... саша40, если вы можете помочь, то напишите мне совет в письме на avfelichev@mail.ru Могу прислать исходник программы открывающей эти файлы...

[raxp]

...выкладывайте в теме исходник и структуру ваших DAT-файлов. А приложенный exe-ник выкиньте, никто его запускать не будет.

[SantilaS]

Исправляюсь) Исходник...

[саша40]

У меня есть функция, которая работает через FFT-параметры. Вот код:

TFFTData = array [0..512] of Single;
procedure Draw(HWND : THandle; FFTData : TFFTData; X, Y : Integer);
FFTPeacks  : array [0..128] of Integer;
 FFTFallOff : array [0..128] of Integer;
implementation

//Timer
var
FFTFata: TFFTData;
begin
if BASS_ChannelIsActive(Stream)<>BASS_ACTIVE_PLAYING then exit;
Bass_ChannelGetData(stream,@FFTFata,BASS_DATA_FFT1024);
Draw(pb1.Canvas.Handle,FFTFata,0,-5);

//Draw
var i, YPos : LongInt; YVal : Single;
begin
//очистка области визуализации
          pb1.Canvas.Pen.Color := clWhite;
          pb1.Canvas.Brush.Color := clWhite;
          pb1.Canvas.Rectangle(0, 0, pb1.Width, pb1.Height);
       //визуализация
       // pb1.Canvas.Pen.Color := clRed;
         for i := 0 to 127 do begin // количество столбиков
           YVal := Abs(FFTData[i]);
           YPos := trunc((YVal) * 500);
           if YPos > pb1.Height then YPos := pb1.Height;

           if YPos >= FFTPeacks[i] then FFTPeacks[i] := YPos
             else FFTPeacks[i] := FFTPeacks[i] - 1;

           if YPos >= FFTFallOff[i] then FFTFallOff[i] := YPos
              else FFTFallOff[i] := FFTFallOff[i] - 3;
                     //рисование пика
                     pb1.Canvas.Pen.Color := clRed;
                     pb1.Canvas.MoveTo(X + i*(3+1) , Y + pb1.Height - FFTPeacks[i]);
                     pb1.Canvas.LineTo(X + i*(3+1) + 3, Y + pb1.Height - FFTPeacks[i]);
                     //рисование столбика
                     pb1.Canvas.Pen.Color := clBlue;
                     pb1.Canvas.Brush.Color := clBlue;
                     pb1.Canvas.Rectangle(X + i*(3+1) , Y + pb1.Height - FFTFallOff[i], X + i*(3+1) + 3, Y + pb1.Height);

              end;

//paintbox.paint
 pb1.Canvas.Pen.Color := clWhite;
          pb1.Canvas.Brush.Color := clWhite;
          pb1.Canvas.Rectangle(0, 0, pb1.Width, pb1.Height);

Преобразование Фурье используется в визуализации. Это звук. Но у вас же другое. Попробуем использовать другие данные.

[raxp]

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

По данным топикстартера видно, что чипы сигнала (укладывается по 6 периодов в чипе) в архивах имеют период 40 мс, что соответствует частоте канального сигнала 25 Гц (вероятно медицинская какая аппаратура, а может еще что). Тут важнее - с какой частотой оцифровывалось это все, ибо это нужно для правильной визуализации (восстановления частоты гармоник на графике). По моим прикидкам частота дискретизации записанных сигналов в архивах в районе около 2700 Гц, что подтверждается переделкой.

Вот полученные спектры для двух каналов из архивов (для облегчения работы конечного пользователя мною была разработана 'FFTR.DLL'):










Вот видео работы:


p.s.: 'Описание API вызова библиотеки FFTR.DLL' приложил к архиву. Кстати, ТС, что за прибор у вас все-таки?

[саша40]

Говорил же сказать что за сигнал. Может автор имеет разные виды сигналов(звук, инфазвук, ультразвук, излучение)?

[raxp]

Саша, топик-стартер приложил сам сигнал в посту #4. Да и я уже визуально его показал и рассказал за его характеристики.

...разные виды сигналов(звук, инфазвук, ультразвук, излучение)

не разачаровуйте меня, вас иногда заносит в сторону. Первые три всего лишь по частотному диапазону условно разделены, а "излучение" - вообще довольно расплывчатая формулировка, ничего не говорящая о характере и сигнатуре самого сигнала.