Форум программистов [Powered by Invision Power Board]

Код

int nn, ii, jj, n, mmax, m, j, istep, i, isign;
double wtemp, wr, wpr, wpi, wi, theta, tempr, tempi;
nn = Memo1->Lines->Count;
double* a = new double[nn];
for (int i=0; i<Memo1->Lines->Count; i++)
        {
        a[i] = StrToFloat(Memo1->Lines->Strings[i]);
        }


n = 2*nn;
j = 1;
ii=1;
while (ii<=nn)
        {
        i = 2*ii-1;
        if (j>i)
                {
                tempr = a[j-1];
                tempi = a[j];
                a[j-1] = a[i-1];
                a[j] = a[i];
                a[i-1] = tempr;
                a[i] = tempi;
                }
        m = n/2;
        while (m>=2 && j>m)
                {
                j = j-m;
                m = m/2;
                }
        j = j+m;
        ii=ii+1;
    }
    mmax = 2;
    while (n>mmax)
    {
        istep = 2*mmax;
        theta = 2*M_PI/(isign*mmax);
        wpr = -2.0*(sin(0.5*theta))*(sin(0.5*theta));
        wpi = sin(theta);
        wr = 1.0;
        wi = 0.0;
        ii=1;
        while (ii<=mmax/2)
        {
            m = 2*ii-1;
            jj=0;
            while (jj<=(n-m)/istep)
            {
                i = m+jj*istep;
                j = i+mmax;
                tempr = wr*a[j-1]-wi*a[j];
                tempi = wr*a[j]+wi*a[j-1];
                a[j-1] = a[i-1]-tempr;
                a[j] = a[i]-tempi;
                a[i-1] = a[i-1]+tempr;
                a[i] = a[i]+tempi;
                jj=jj+1;
            }
            wtemp = wr;
            wr = wr*wpr-wi*wpi+wr;
            wi = wi*wpr+wtemp*wpi+wi;
            ii=ii+1;
        }
        mmax = istep;
    }
  for (int i=0; i<nn; i++)
        {
        Memo2->Lines->Add(FloatToStr(a[i]));
        Form2->Series1->Add(a[i],2,Color1);
        }
delete(a);