フーリエ

package inter;

public class pro{

double  alpha=new double[1024];
double  beta=new double[1024];
double  data=new double[1024];
int s;

    public static void main(String args) {
 pro test=new pro();
    }

    
    
    pro(){


        for(s=0;s<1024;s++){
        data[s]=2*s+1;
        }

        fourier sub=new fourier();
        sub.data=data;
        sub.makedata();
        alpha=sub.alpha;
        beta=sub.beta;

        System.out.println(fou(150));
        
    }

double fou(double x){
    double f1;
    double dx;
    int n,nx;
    n=1024;
    dx=2*Math.PI/n;
    
    f1=alpha[0];
    for(nx=1;nx<512;nx++){
f1=f1+alpha[nx]*Math.cos(nx*x*dx)+beta[nx]*Math.sin(nx*x*dx);
    }
return f1;
}

}