import?sys
import?os
from?math?import?*

import?numarray
import?numarray.fft

import?rsoft.rspytools?as?rspy

if?len（sys.argv）<2:
??print?‘usage:?filter?‘
??sys.exit（1）

fname=sys.argv[1]
if?not?os.path.exists（fname）:
??print?“Can‘t?find?%s“%fname
??sys.exit（1）

prefix=fname[:-4]

#load?.tmn?file
tmn_data=rspy.file_to_numarray（fname）

v_t=tmn_data[:0]
v_ey=tmn_data[:1]
clight=2.9979e14

#trim?off?start?of?.tmn?file?containing?CW?ramp
npts=len（v_t）
highpow2=int（floor（log（npts*1.0）/log（2.0）））
goodpts=int（pow（2?highpow2））

#apply?raised?cosine?window?function
dt=（v_t[1]-v_t[0]）/clight
v_ey=v_ey[-goodpts:]
cos_domain=numarray.arange（goodpts?typecode=‘d‘）/（goodpts-1）*pi-pi/2
filt_ey=v_ey*numarray.cos（cos_domain）
??
#make?fft
v_ey_fft=numarray.fft.fft（filt_ey）
v_ey_fft_pos=numarray.abs（v_ey_fft[0:goodpts/2]）/goodpts
v_nu=numarray.arange（goodpts/2?typecode=‘d‘）/（goodpts*dt）/1e12??#?freq?in?THz

#output?windowed?fft?data
rspy.numarray_to_file（‘%s.fft2‘%prefix?v_ey_fft_pos?x=v_nu）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????9298928??2006-08-03?11:00??RSOFT7\docs\bandsolve.pdf

?????文件????3543151??2009-11-10?23:56??RSOFT7\docs\beamprop.pdf

?????文件????1272279??2006-08-02?11:00??RSOFT7\docs\diffractmod.pdf

?????文件?????896311??2006-08-03?11:00??RSOFT7\docs\femsim.pdf

?????文件????2557856??2010-03-11?11:39??RSOFT7\docs\fullwave.pdf

?????文件????1299668??2006-08-21?11:00??RSOFT7\docs\gratingmod.pdf

?????文件?????174591??2006-08-21?11:00??RSOFT7\docs\install.pdf

?????文件????3921843??2006-04-05?13:48??RSOFT7\docs\lasermod.pdf

?????文件????4198886??2006-08-03?11:00??RSOFT7\docs\most.pdf

?????文件?????540510??2002-04-03?11:00??RSOFT7\docs\rplot.pdf

?????文件?????889326??2004-10-20?11:00??RSOFT7\docs\rsdatabrowser.pdf

?????文件????1965925??2006-08-03?11:00??RSOFT7\docs\rsoftcad.pdf

?????文件????????288??2006-03-30?11:00??RSOFT7\etc\bcadw32.reg

?????文件????????766??2003-04-15?11:00??RSOFT7\etc\setup.ico

?????文件???????2252??2006-03-30?11:00??RSOFT7\etc\winplot.reg

?????文件????????304??2006-03-30?10:00??RSOFT7\etc\xp\bcadw32.reg

?????文件???????2268??2006-03-30?10:00??RSOFT7\etc\xp\winplot.reg

?????文件????????133??2005-08-05?17:30??RSOFT7\example\BandSOLVE\ApplicationNotes\makefile

?????文件????3927144??2008-09-17?16:19??RSOFT7\example\BandSOLVE\ApplicationNotes\pcslabs\bstmp_domain.vtk

?????文件??????15318??2005-06-15?17:51??RSOFT7\example\BandSOLVE\ApplicationNotes\pcslabs\holes_solid_cladding.ind

?????文件????????573??2005-08-05?17:28??RSOFT7\example\BandSOLVE\ApplicationNotes\pcslabs\makefile

?????文件??????15868??2003-12-10?21:56??RSOFT7\example\BandSOLVE\ApplicationNotes\pcslabs\nosymmetry.ind

?????文件???????1304??2006-04-03?14:23??RSOFT7\example\BandSOLVE\ApplicationNotes\pcslabs\nosymmetry.kpts

?????文件??????16012??2005-06-15?17:48??RSOFT7\example\BandSOLVE\ApplicationNotes\pcslabs\rods_periodic_cladding.ind

?????文件??????10990??2005-05-12?13:38??RSOFT7\example\BandSOLVE\ApplicationNotes\pcswaveguides\cubic\circref_clad.ind

?????文件???????9358??2005-05-12?13:38??RSOFT7\example\BandSOLVE\ApplicationNotes\pcswaveguides\cubic\clad_lightline.ind

?????文件???????9414??2005-08-15?11:18??RSOFT7\example\BandSOLVE\ApplicationNotes\pcswaveguides\cubic\clad_lightline_full.ind

?????文件????????343??2005-08-12?04:28??RSOFT7\example\BandSOLVE\ApplicationNotes\pcswaveguides\cubic\combined.pcs

?????文件???????1334??2005-08-15?12:29??RSOFT7\example\BandSOLVE\ApplicationNotes\pcswaveguides\cubic\combined_scan.pcs

?????文件???????9406??2005-05-12?13:38??RSOFT7\example\BandSOLVE\ApplicationNotes\pcswaveguides\cubic\cubicslab.ind

............此處省略1121個文件信息