Чо там читать? Берешь и ебашишь!

function wave(freq, phas, amp, t)
  phas2 = (phas*2*pi) / freq;
  w = amp .* sin(freq .* (t .+ phas2));
  return(w);
  end;

(Reply to this) (Thread)

отсюда:

void FFT(int dir, long m, complex  x[])
{
   long i, i1, i2,j, k, l, l1, l2, n;
   complex  tx, t1, u, c;

   /*Calculate the number of points */
   n = 1;
   for(i = 0; i < m; i++) 
      n <<= 1;   

   /* Do the bit reversal */
   i2 = n >> 1;
   j = 0;

   for (i = 0; i < n-1 ; i++)
   {
      if (i < j)
         swap(x[i], x[j]);

      k = i2;

      while (k <= j) 
	  {
         j -= k;
         k >>= 1;
      }

      j += k;
   }

   /* Compute the FFT */
   c.real(-1.0);
   c.imag(0.0);
   l2 = 1;
   for (l = 0; l < m; l++) 
   {
      l1 = l2;
      l2 <<= 1;
      u.real(1.0);
      u.imag(0.0);

      for (j = 0; j < l1; j++) 
	  {
         for (i = j; i < n; i += l2) 
		 {
            i1 = i + l1;
            t1 = u * x[i1];
            x[i1] = x[i] - t1; 
            x[i] += t1;
         }

         u = u * c;
      }

      c.imag(sqrt((1.0 - c.real()) / 2.0));
      if (dir == 1)
         c.imag(-c.imag());
      c.real(sqrt((1.0 + c.real()) / 2.0));
   }

   /* Scaling for forward transform */
   if (dir == 1) 
   {
      for (i = 0; i < n; i++)
         x[i] /= n;      
   }   
   return;
}

и фсё

(Reply to this) (Parent)

А первые два коммента красивые, правда?

(Reply to this) (Thread)

не по делу же

(Reply to this) (Parent) (Thread)

Именно этим.

(Reply to this) (Parent)

А по делу -- AD DSP Guide, http://www.analog.com/en/education/education-library/scientist_engineers_guide.html , chapters 8-12,31.

Функан преобразования Фурье как таковой (ортогональные системы функцый и вот это всё) там не даётся, но по дискретному -- в общем, достаточно дажэ теории чтобы это самому вывести, а практических советов так прямо норм.

(Reply to this) (Thread)

Ну, я посмотрю конечно. Но вот к книжкам со словами "signal processing" в названии я после пары попыток отношусь с некоторым предубеждением. Оно-то конечно про преобразование Фурье, но под таким специфическим углом зрения, который мне мало интересен.

Мне, в смысле, преобразование Фурье интересно как средство анализа данных, а не как средство преобразования сигнала.

(Reply to this) (Parent) (Thread)

Я вообще не понял о чём ты, но подозреваю, что этот гайд тебе нужно посмотреть с первой части. А потом ещё часть Applications.

(Reply to this) (Parent) (Thread)

книжка очень суперузкоспециальная же. аудиосигналы и обработка растровых изображений его явно не интересуют. его интересует типа обработка геопространственных данных с примерами, других вещей, все такое. меня в общем тоже.

(Reply to this) (Parent)

Анализ данных это и есть разновидность обработки сигналов.

(Reply to this) (Parent)