24小时热门版块排行榜

返回列表

weivichen

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 249.5
红花: 1
帖子: 20
在线: 8.9小时
虫号: 1556819
注册: 2011-12-30
性别: MM
专业: 计算数学与科学工程计算

[求助] sin函数的傅里叶变换

我想做sin函数的傅里叶变换可是程序没有error,但是运行就说那个exe文件停止工作，程序如下，好像很多问题的样子，新手求指教！
program main
use textgraphlib
implicit none
real,parameter::endx=3.1,startx=0.0,xinc=0.1
integer,parameter::maxx=32,maxy=20
real ::x,px=0.0,py=0.0
integer i,num,fx,fy,xx
real,allocatable::datax(

,datay(

x=startx

do i=1,maxx
x=px
py=sin(x)
datax(i)=px
datay(i)=py
px=px+xinc
end do

call four1(datay,16,1)
call setscreen(maxx,maxy)

  do xx=1,maxx
fy=datay(xx)
call putchar(xx,fy)
  end do

  call updatescreen()
stop
end program

subroutine four1(data1,nn,isign)
integer isign,nn
real data1(2*nn)
integer i,istep,j,m,mmax,n
real tempi,tempr
double precision theta,wi,wpi,wpr,wr,wtemp
n=2*nn
j=1
do  i=1,n,2
if(j.gt.i)then
tempr=data1(j)
tempi=data1(j+1)
data1(j)=data1(i)
data1(j+1)=data1(i+1)
data1(i)=tempr
data1(i+1)=tempi
end if
m=n/2
1 if((m.ge.2).and.(j.gt.m)) then
   j=j-m
   m=m/2
   goto 1
   end if
   j=j+m
   end do
   mmax=2
2 if(n.gt.mmax) then
   istep=2*mmax
   theta=6.28318530717959d0/(isign*mmax)
   wpr=-2.d0*sin(0.5d0*theta)**2
   wpi=sin(theta)
   wr=1.d0
   wi=0.d0
   do  m=1,mmax,2
   do  i=m,n,istep
      j=i+mmax
   tempr=sngl(wr)*data1(j)-sngl(wi)*data1(j+1)
   tempi=sngl(wr)*data1(j+1)+sngl(wi)*data1(j)
   data1(j)=data1(i)-tempr
   data1(j+1)=data1(i+1)-tempi
   data1(i)=data1(i)+tempr
   data1(i+1)=data1(i+1)+tempi
   end do
   wtemp=wr
   wr=wr*wpr-wi*wpi+wr
   wi=wi*wpr+wtemp*wpi+wi
   end do
   mmax=istep
  goto 2
  end if
  return
  end

回复此楼