[dennou-ruby:000897] Re: numfft

[Takeshi Horinouchi <horinout@xxxxxx>]

��Ƿ��Ǥ���

�ⶶ(����)�͡�

�ɤ��ɤ�Ǥʤ��ä��Τǡ����Τ������Ǻ��Ǥ�ȤϻפäƤޤ���Ǥ�
����

> ����ܤμ���Ǥ���
> �����Ĥμ������Ф��ơ�fft ��ݤ�����ϡ�
> ¿��������� shape ���Ѥ��ʤ��Ǥ�����
> ʣ���μ����� fft �򤫤ޤ������Ǥ���

���ξ����Ѥ�ޤ���

> �㤨�С�
>         na = NArray.new(Float,5,4,3)
> �Ȥ��ơ�
>         na.rfft(0)
> �Ȥ�������Ĺ��5�ΰ켡�������4*3=12��fft ����櫓�Ǥ���͡�
> ���λ��ϡ�Ĺ��5�� fft ���줿����12�ܡ�5*4*3�ζ��֤˼��ޤ�ޤ�����
> ʣ���μ����� fft �������Ȥ��ϡ��⤦����ʾ她�ڡ�����̵���櫓�Ǥ����顢
> ����η��򿷤������ꤷ�ʤ���Фʤ�ʤ���

�Ĥ����裲�����ܤˤ������Ȥ��ϡ�ñ��Ĺ������fft��5*3����ޤ���
��̤Ϥ�Ȥ�Ʊ�� 5*4*3 ������Ǥ������᡼���Ȥ��ƤϤ���ʴ�����
             
    f(x,y) ==> F(k,l)

  where

    F(k,l) �� \int \int f(x,y) exp(-ikx-ily) dx dy

(����ܤ���ʬ��̵�뤷�Ƥޤ���Ϣ³�ؿ����Ф��뼰��񤤤Ƥޤ���
�ǼºݤȰ㤤�ޤ������줬������Ȼפ�ʤ��Ǥ�����������
����ϰʲ��Τ褦�ˡ�x �ˤĤ���1����Ū�˷׻����ơ����η�̤� y ��
�Ĥ��Ʒ׻������ OK������ y ����ˤ�äƤ�Ʊ���ʤΤ�ʬ����ޤ��͡�

    f~(k,y) = \int f(x,y) exp(-ikx) dx
    F(k,l) = \int f~(k,y) exp(-ily) dy

    f'(x,l) = \int f(x,y) exp(-ily) dy
    F(k,l) = \int f'(x,l) exp(-ikx) dx

Ǽ���Ǥ��礦����