ふぉうｃｘ - プログラム倉庫

class fou{

double delta;
double[] data=new double[200];
double[] alpha=new double[200];
int datanumber;
int qx,m,n;  

double makedata(double nx){

double fx;

seekalpha(datanumber,data);

fx=fou(datanumber,nx,alpha);

return fx;
}

double fou(int datanumber,double nx,double[] alpha ){
double f1;
double delta;
int n,qx;
delta=2*Math.PI/datanumber;
f1=alpha[0];
for (qx=1;qx<datanumber;qx++){
f1=f1+alpha[qx]*Math.cos(qx*delta*nx);
}
return f1;
}

void seekalpha(int datanumber,double[] data){
double delta,z;
double[][] a=new double[datanumber][datanumber];

delta=2*Math.PI/datanumber;

for(m=0;m<datanumber;m++){
for(n=0;n<datanumber;n++){
a[m][n]=Math.cos(-m*n*delta);
}
}

for(m=0;m<datanumber;m++){
z=0;
for(n=0;n<datanumber;n++){
z=z+a[m][n]*data[n];
}
alpha[m]=z/datanumber;
}

}






}