f
package inter; import ca.uol.aig.fftpack.*; public class fourier{ int i,n; double dx,x; double a; double alpha=new double[1024]; double beta=new double[1024]; double data=new double[1024]; void makedata(){ n = 1024; a = new double [n]; RealDoubleFFT my = new RealDoubleFFT(n); dx = 2 * Math.PI / n; for (i = 0; i < n; i++) { a[i]=data[i]; } my.ft(a); alpha[0]=a[0]/n; int s; for(s=1;s<512;s++){ alpha[s]=2*a[2*s-1]/n; beta[s]=-2*a[2*s]/n; } } }
</pre>