program fourierparameter(pi=3.141592653,max=100)complex:i=(0,1) real:xm=0.0,x,n,m,dx,ddx=2*pi/maxdo m=0,99do n=0,99d=5.0/nfm=m/(max*d)x=cos(5*dx*n)xm=xm+x*exp(-i*n*m*2*pi/max)end doopen(1,file=1.txt)write(1,*)fm,xmxm=0.0end doend0 4 8 12 16 20-200204060有限离散傅氏正变换
