フーリエ変換サブルーチン
V2 Fourier Transformation Subroutine
Closed; even; 50 storage locations + 34 storage locations for S2
and 16 storage locations for working spaces.
Working spaces: 0w 〜 15w.
Puts Cm in 12wL and Sm in 14wL, where 2m = C(3)
and
Cm = 2-15Σr=0n-1arcos(2mrπ/n),
Sm = 2-15Σr=0n-1arsin(2mrπ/n),
where ar=C(A+r).
Presets parameters:
0n = number of points,
0a = A location of a0
0h = 15, ┐
0s = 50f, │ set by the subroutine tape.
0w = 34s = 84f,┘
Notes:
1. Scale factor is set to 2-15, but it can be changes to 2-h
by re-setting the h-parameter.
2. The subroutine places -2m/n ,2-3cos(2mπ/n) and
2-3sin(2mπ/n) in 6wL, 8wL, 10wL respectively.
0f: a 40 ┐plant link
1 x 46r ┘
2 q 2 ┐
3 n 3 │-2m/n to 6wL
4 d 62 │
5 tl 6w ┘
6 tlo ┐
7 it │2-3cos(2mπ/n)
8 jls │to 8wL
9 tl 8w ┘
10 p 2 ┐
11 tlw │clear working spaces
12 tl 4w ┘
13 p 47r ┐
14 x 25r │
15 s 48r │
16 z 32r │
17 q 49r │
18 pl 4w │
19 wl 8w │
20 ll 4 │
21 slw │
22 tl 2w │iterative cycle
23 pl 4w │
24 tlw │to compute
25 p ( ) │Si+1=
26 rh │2Sicos(2mπ/n)
27 al 2w │-Si-1+2-hai
28 tl 4w │
29 p 25r │
30 a 49 │
31 jl 14r ┘
32 pl 4w ┐
33 vl 8w │Cm to 12wL
34 ll 3 │
35 slw │
36 tl 12w ┘
37 nl 6w ┐
38 a 21s │
39 tlo │2-3sin(2mπ/n)
40 it │to 10wL
41 jls │
42 tl 10w ┘
43 vl 4w ┐
44 ll 3 │Sm to 14wL
45 tl14w ┘
46 jl ( ) link
47 | 0a
48 | 0na
49 | t round off constant
Date 60-01-11
Author H. Takahasi