Professional Documents
Culture Documents
Research Online
University of Wollongong Thesis Collection University of Wollongong Thesis Collections
1977
Recommended Citation
Montagner, G. John, Identification and optimisation of hopper-discharge chute systems for bulk granular materials, Doctor of
Philosophy thesis, Department of Mechanical Engineering, University of Wollongong, 1977. http://ro.uow.edu.au/theses/1583
Research Online is the open access institutional repository for the University of Wollongong. For further information contact the UOW Library:
research-pubs@uow.edu.au
IDENTIFICATION AND OPTIMISATION OF HOPPER-DISCHARGE CHUTE
by
of Doctor of Philosophy
May, 1977
808157
TABLE OF CONTENTS
PAGE
ABSTRACT 6
ACKNOWLEDGEMENTS 8
LIST OF FIGURES 9
LIST OF PHOTOGRAPHS 14
LIST OF TABLES 15
CHAPTER 1: INTRODUCTION 18
1.1 Genera I 18
2.1 General 25
2.3 Model 33
2.3.1 Introduction 33
2.4.1 Introduction 43
2.4.2 Apparatus 43
2.4.3 Procedure 45
' PAGE
CHAPTER 3: OPTIMISATION 52
3.1 General 52
3.3.1 Introduction 59
3.4.1 Description 61
Calculations
4.7 Conclusions 95
PAGE
(Assumptions, Limitations)
PAGE
Perturbations 155
PAGE
ABSTRACT
dynamic characteristics.
The design procedure for minimum transit time is presented in the form
chute system under gravity flow are examined. Using the P.R.B.S. and
for the system under varying conditions of initial flow and chute
model were determined, enabling both the system bandwidth and the
ACKNOWLEDGEMENTS
John Montaqner.
9
LIST OF FIGURES
chute.
chute.
metre chute.
10
metre chute.
chute.
chute.
chute.
metre chute.
metre chute.
chute.
valve perturbations.
6 = 35°.
6 = 40°.
6 = 45°.
and Watanabe.
LIST OF PHOTOGRAPHS
3 2.2 FIowmeter.
perturbat ions.
Transducer.
15
LIST OF TABLES
alphabetical order.
a coefficient
a acceleration
B chute width
f frequency
j /-T
at zero velocity
m mass
n normal co-ordinate
n number of samples
R correlation function
s d isp 1acement
t t ime
t tangential co-ordinate
x» y Cartesian co-ordinates
X (t ) input signal
P radius of curvature
U) computation variable
CHAPTER 1: INTRODUCTION
1.1 Genera I
usually.
vi ab i Iity .
per year. For 1973 this amounted to 20,135,000 tonnes with a dollar
optimising the design for gravity flow under steady state conditions
as blending.
by Roberts et a I (References 4 to 6 ).
problem.
For the controlled flow case where flow transients are caused
methods are now well known, being more reliabile than conventional
stochastic analysis.
insensitivity.
(Reference 10) efficient only when Coulomb friction was ignored. The
the low exit angle, instabilities in flow occurred and the method
(Reference 15).
22
equipment cost for on-line analysis, and large data reductions for
off-line work. The former was not a problem in this case as the
of the ’fast’ flow phase. In the ’fast’ flow mode the material makes
contact with the chute bottom and side walls but does not make
contact with the top of the chute. In the ’slow’ flow mode the chute
the complete perimeter of the chute. ’Fast’ flow is more efficient with
the flow rate out of the hopper being entirely governed by the
retardation effects are encountered which impair the flow rate and
onto the plant under test whi 1st the plant is in normal operation.
2.1 Genera 1
mode of flow is the ’fast’ flow mode in which the grain stream
receiving hopper, which is the feed hopper for the grain conveyor.
The recycling ability is clearly seen (in Figure 1 and Plate 1).
P.R.B.S,
NOISE MILLET RETURN LINE
GENERATOR
POWE R
AMPL IFI ER
.
k g /m i n
Fo w r a t e
1 HOPPER
P.R.B.S.
TRANSDUCER
HOPPER
GRAIN
ELEVATOR
CHUTE
FLOW
METER
GRAIN ELEVATOR
RECEIVING HOPPER
POWER
AMPLIFIER
ControI-Perturbinq Valve
hopper orifice.
Chutes
were cut using a bandsaw from full sized drawings glued to the
IncIi nometer
FIowmeter
2.3 Model
Pn
F igure 3 Pr e s s u r e distribution on
CHUTE WALLS ( ROBERTS )
the bottom and side wall pressure distribution and found that the
depth of gra in .
F igure 4 Pr e s s u r e distribution on
CHUTE WALLS ( PARLOUR )
Te = T(1 + k jj)
flowing stream was undertaken and the ratio of work done due to
following proportions.
This showed that the major energy loss was against the chute
lumped model takes into account the friction drag force based on Us
found that on this basis the ratio of work done by the various
components was
= 0.306
xe = t [ 1 + keo | (1 + c v 2)] (2 .1 )
r
S•inee vO = 7—
H
v Ho
where
_ keo v q Hq
1 B
= keo vp Hp c
C2
keo = 0.346
c = 0.0215 (sec/ft2)
= 0.001997 (sec/m2)
Ci = 0.346 Vo "Tp*
b
c 2 = 6.9096 x ( 10- “) x
b
40
T = 0.306
0f = 90°-tan ^ T q ) (2 .4 )
with the vertical was increased whilst ’fast’ flow was in progress.
stream flow within the chute and upon leaving the chute were made.
measure the transit time for a single grain particle moving down the
with the grains sliding on the bottom surface and striking the
TRACE 1
FLOWMETER RESPONSE
TRACE 2
valve within the empty hopper and the valve was then stepped open
parameter model, tests were carried out as discussed above. The mean
descent time was found to be 0.704 seconds* for the model chute with
shown (Section 5.2) that the cross-correI at ion determined under these
integrated, yields the step response. The transit time was taken to be
the delay time to the midpoint of the step rise, as shown in Figure 6.
As the transit time so calculated for the bulk flow was found to
assumpt io n .
2.4.1 Introducti on
method.
2.4.2 Apparatus
All the uniform flow tests were carried out on the model
Figure 7 and Plate 5 were, in turn connected to the testing rig such
that all of them received from and discharged to the same end points.
That is, the straight chute at 40° to the vertical, the parabolic
Xo
where c =
(Yo)
and both the unconstrained and the constrained (as detailed later)
RX IS (METRES)
F igure 7 Ch u t e p r o f i l e s f o r m o d e l c h u t e
£*Il '’
¿g& - t< »»v *
Jt.î¿¿ *'"T ‘'‘■v’
";.’
;j''i
*J»Í
&5$
OPTIMUM CHUTE
(t = O.O)
PARABOLIC
CHUTE STRAIGHT CHUTE
Pl a t e 5 Co m p a r i s o n of c h u t e p r o f i l e s
46
form a closed loop enabling the static capacity of the model hopper
2.4.3 Procedure
In each case millet seed was set cycling round the testing loop
to condition the chute surface, and to filter out husks which appeared
the transport delay, of the time taken for the grain stream to pass
from the perturbing gate to the flowmeter. The main ’flow’ gate was
4 kg/min. whilst the perturbing gate stroke was tested at both ±2mm
and ±4mm. The capability of the correlation method to cope with this
yx‘
Vp = (2.5)
e |_2cos 2 « (y-x tan °0 _
knowing the film scaling factor and the framing rate. In all
interest, and the camera enabled timing pips, derived from a timing
chute runout.
section the average velocity of the grain stream is v and the mean
stream thickness is H.
F igure 8 Ac c e l e r a t i o n m o d e l - f o r c e b a l a n c e
i.e., ? m v 2 2 - m g y - W . D . 1^ 2 + i m v i 2 (2.6)
D e f i n e a dr ag fo rc e F q _ |n |+ yv (2.7)
w h e r e N = dmg SinG + dm a^ ( 2. 8)
49
Now
a = /a)'
n
= /(— ) 2
f dt
= )(— )
J at; at;
1 ,dsN ,d0 .
= / - nf dt (For this orientation ds = -A10)
ds
2 d8 Q (2.9)
= -v — COS0
dy
dfl
N = dm g Sin0 + dm (-v2 — cos0) (2. 10)
a dy
: _ g Sin0 - v 2
N11 cos0 (2 . 11 )
dy
where
Now
ds = / l + ( ~ ) Z dy
dy '
= J 1 + tan 20 dy
= sec 0 dy (2.14)
50
W .D .i->2 = /2 Fd dy
i cosG
dy + I (2.15)
vz 2 = gy - /v fd Vl
o COS0
re-arranging
FP
v 22 = 2gy - 2 f dy t vi2 (2.16)
cos0
and final Iy
v 2 = *4gy - 2 /Y dy + Vj 2 (2.17)
37 O COS0 7
friction components.
yields
and as before
Now
d (v)
A+ =
dt
_d_
dt V « ’]
pd20 6f_ d0
( 2 . 21 )
Jdt 2 dt dt
and An = /(di ) 2 ( 2 . 22 )
J dt;
yields:—
“ 2 ^ ^ od20 d/ d0 Q n
Te g Sin© + f(^rr) + yv + + dt * dt " 9 cos6 = 0
a dt
(2.23)
re-arrangi ng
d 20 1 ,df . d0
(^ 0 + ^ (Te Si n0 - cos0 ) + = 0
dt 2 f dt dt ) + Te
(2.24)
friction forces.
CHAPTER 5 : OPTIMISATION
3.1 Genera I
alluded to in Section 1.1 and once the material parameters have been
discharging material from a hopper into a rail car, one may wish to
the most favourable flow pattern, i.e., ’fast’ flow, which helps to
speed and direction of material falling onto a conveyor belt has been
co-ordinates 0 and y are the most satisfactory. 0 (y) is the chute slope
(2.17)
repeated
The descent time between two successive cross-sections A and B is
given by:
+AB = /B dt
A
jB d£
A v (3.1)
54
, rB sec0 ,
+AB = f Z dY (3.2)
A
Since the end points (X, Y) are fixed, the function 0(y ) must
sati sfy
,Y
tan0 dy = X (3.3)
,Y d
w(y) = -/ yv + Te |N| (3.4)
cos
0 L
Note that u)(y) satisfies the differential equation
dio = -[ yv + xe IN11
dy cos0 (3.5)
with w(o) = 0 .
.Y dy (3.7)
^ /vi 2 + 2 (gy + w) cos0
Equation 3.3.
55
framework one may aim to optimise either the exit velocity or some
losses;
Now ds = Sec0 dy
N = dmg Sin6 + dm An
56
d0
(V2 2 - v i2) = 2gy - 2^S {ie (g Si n0 - V2 2 C O S0
dy
(3.12)
and exit angles, entry and exit co-ordinates one could complete
the exit velocity for the exit angle specified. In the present study
angle that the chute bottom makes with the vertical, 0 , needs to be
and can exceed the limiting angle, 0f. For ’fast’ flow this angle
determines the stability of the chute system. Should the stream flow
would stall and choke or at the very least flow would change from
’fast’ flow to ’slow’ flow. Once stalled, the chute would not self
was detected whilst the chute was discharging under a light flow
the material being handled and to ensure that the combination of inlet
ppa*” *,*»
,;^r f '
’" miijr*— •-
tiËÊÊÊm ''*.-- :
FRAME 1 FRAME 2
FULLY
í*** -.'\;'.'••/ ■/- í, v •’\ W CHOKED
.- • . .r >'• r .
CHUTE
INITIAL FORMATION OF
GRAIN FLOW INSTABILITY |$P§pp
FRAME 3 FRAME 4
Pl a t e 6 C h u t e s t a l l i n g (u n c o n s t r a i n e d )
3.3 Methods of Solution
3.3.1 Introduction
o m
i.e., tan0(y) = Ao + Aiy + A 2y + .... + Am y
60
d 0_
An expression for , required in Equation 2.18 is obtained
dy
by differentiating Equation 3.13
M
d9 m-i
i.e., sec 20 £ m Am y (3.14)
dy
m=i
re-arranging
M _
d0
cos 2e £ m Am ym 1 (3.15)
dy
m= i
Now substituting Equation 3.13 into the end condition Equation 3.3
M ,
y
L A v m+1
Mm,.y____ X = 0 (3.16)
m=o m +1
for which
M-1
M+1
X £ Am Ym+1 (3.17)
Am y M+1 m+i
m=o
This means that only the coefficients {A0 , Ai, A2, ... A^_^}
M coefficients.
Equation 3.3. .
0O = tan 1 A0
61
3.4.1 Description
where K = 1, 2, 3, . ..., M.
62
proceeds as follows:
M M
$ = F - E AK G k + B E G K2 (3.21)
K=1 K=1
is set equal to the identity matrix for the first iteration. The
N M N
3G j 3G k 3Gk 8F . ( 3 . 22 )
E E (A E
3xj 3x{ 3xj 3xj
¡=i j=i i= i
where K = 1, 2, ..., M.
constraining equations.
Initial values for the chute parameters are read in, including
data.
During the course of the research programme the need for improved
which was carried out in parallel with the main experimental programme
and tested towards the end of the research programme. This ultrasonic
say, inlet to and outlet from the chute. Having established that the
exit velocity from the model hopper into the chute entry was
approximately 0.5 metre per second and using the previously derived
this model scale (see Section 4.5 for a detailed sensitivity analysis)
this variation does affect both the descent time and the exit
Coefficient of Viscous
Drag y - 0 1 .0 1.5 2.0
Aiy2+A 2y 3+A3y‘t+Aity 5
Equation of Profile x -
2 3 4 5
66
for this chute experimental results were available for inlet velocity,
exit velocity and descent time, as well as the length of the chute,
open chute subject to friction on the chute bottom and waI Is. It is
forces.
F igure 9 Fo r c e balance
67
dv
= g cos0 - T0 N - yv (4.1)
dt
= g cos0 - Te g sin0 - yv
dv
Note that ppp - 0 for terminal velocity.
f = A - B - f - Dv - uv
where
A = g cos0
B = tg sin0
C = xg sin© ci
D = t g sin0 C2
E = A- B
F = D+ y
which yields
± L = e - £. - F (4.3)
dt v '
| 2 dv = /(E ~ ~ ~ Fv )dt
After re-arranging
ultimately results in
E J_
f 2Fv i v 2 - E(vi + y2 ) - 2Fa(v2 - V i ) + 2c ] ^ 3 ["Fvi2 - Evi + c l ^
L 2F v i v 2 - E(vi + v2) + 2Fa(v2 - Vi) + 2c J L F v 22 - Ev2 + cj
v2 = 2.7 m/s
t = 0.735 sec.
This result also needs to be considered in the context that all the
to yield the ’best case’ and ’worst case’ value of y. This result
would seem to support the assumption that for the test environment the
69
Using the same right hand side of Equation 4.3 we rewrite the
acceleration expression as
dv c c _
ds v
Ev - c - Fv2
• 2 dv r •
i .e., v T - = Ev - c - Fv (4.8)
ds
Re-arrang ing
v2dv
= ds
Ev - c - Fv2
v2dv
= ds
C 2\
F(— v - -p - vz )
v2dv
= ds
Ev.
-F(v2 + ^ -
vzdv (4.9)
= -Fds
in _E_
2F
(V 2-V 1 ) ifl£i + c - Ev 2 2
I f v1 2 + c - E 1
2a \2¥
r 2Fy1v2 ~ E (v i+ V 2 ) - 2 F a (V 1 - V 2 ) + 2c 1
[2FV1V2 ” E(v i + V 2 )) +
+ 22Fa(V
F a ( V 11--V
V22 )) +
+ 2c
2c J
= e ~ F ( s 2- S i ) (4.15)
coefficient such that Equation 4.15 is satisfied. Use was again made of
vi = 0.5 m/s
v 2 = 2.7 m/s
s 2 - Si = 1.145 metres
±3$ in just these three parameters revealed the spread in the value
parameter variations to yield the 'best case' and 'worst case’ value
y 1.83* 2.28**
CENTROID OF IMPULSE
above interprefation for chutes of known form. The descent time for
It can be seen from the above that even at this small scale level
designed chutes when compared to the more commonly used chute forms.
Parabo1ic
Chute .712 .707 .721 .72 .716 2.7 2.4
Constra ined
Chute .704 .700 .707 .70 .702 2.7 2.3
varied over a range that was considered practicable for the model.
this angle precisely. For this analysis and that of Section 4.5 the
and hence performance,and the results from the 60° case only are
reported here.
73
For each analysis one variable was held constant whilst each
purposes.
drag (y) between 0.5 and 1.5. However, this variation in y does
74
- 0 . 1
- 0 2 -
-0- 3
AXIS (METRES)
-0. 4
-0. 5
- 0 6 .
- 0 . 7
- 0. 8
-0. 9
D-Ü 0-1 0-2 0.3 0.4 0.5 0.6 0.7
X AXIS (METRES)
AXIS (METRES)
Q . Q -4____ i____ 1____ i____ I____ i____ I____ i____ I____ i____ I____ i____ I____ i____ L
RXIS (METRES)
time. The maximum spread in exit velocity was 34$ while the
limiting angle (0f) between 55° and 65°, together with only minor
velocity was 1$ and the spread in descent time was only 0.3%.
velocity and descent time. The spread in exit velocity was 6.5$
height to breadth ratio (HB) but are associated with only a slight
exit velocity was 2.8$ and the spread in descent time only 0.9$.
if they were seen to hold for larger scale chutes. While it was not
essentia I.
80
• — \—
1 X = 0,75 m
Y = 0.9 M
SCALE : 20 MM TO 1.0 M
X = 3.0 M
, Y = 4.0 M
X = 5.0 m
! Y = 5.0 m
f Y
F i g u r e 16 Sc a l e u p o f c h u t e g e o m e t r i e s
and a parabolic chute, and profiles for both the unconstrained and the
(METRES)
T RX1S
purposes.
fi % % % % %
Viscous Drag
Coefficient 17.6 33.4 46.1 34.2 60 15 A
Limiting
Angle .3 .18 .54 1 1 .6 .3
Inlet
VeIoci ty 16 8.8 7.4 1 .8 .28 .6
Coeff icient
of Coulomb
Friction 2.3 .72 1.6 6.5 .64 3.5
Height to
Breadth Ratio .92 .08 .24 2.8 .49 .69
with y > 1.0. The descent time variation increases from a spread
the same time the spread of exit velocity increases from 34.2%
to 75.4%.
84
(METRES)
AXIS
F i g u r e 19 Va r i a t i o n in limiting angle 0f
FOR (3x4) METRE CHUTE
(METRES)
R XI 5
-4. 0
0. ü 1 .ü 2. D 3. ü
X RX I S (METRES)
Figure 20 Va r i a t i o n in i n i t i a l v e l o c i t y
FOR (3x4) METRE CHUTE
86
(METRES)
AXIS
AXIS (METRES)
2. For changes in the limiting angle (0f) between 55° and 65° there
0.5^, while the exit velocity spread has a maximum across the
range of 1.6$.
3. Changes in inlet velocity (VQ ) from 0.5 to 1.5 m/s across the
The descent time spread drops from 16$ to 7.4$ and the exit
( 3 x 4 ) metre chute.
0.306 across the size range of the chutes have a decreasing effect
across the whole range being always less than 1$. The exit
( 3 x 4 ) metre chute.
for chutes operating under uniform flow conditions in the TfastT flow
95
parameters are read into and depending on the computer used, the
execution time and cost are relatively small. In the case reported
comp i Ier.
constra ints.
savings.
of the chute.
minimise the waste. Pressure for better quality control and increased
test input characteristics and wide band and output drift. His
2. Set the digit interval (AT) to one quarter of the smallest time
constant of interest.
presentat io n .
input signals x(t) and input noise n(t) and with an output signal or
The input signal x(t) is the test signal (P.R.B.S.) which is injected
X (t) --
SYSTEM
INPUT
operating signal and any random noise signal originating from some
„00
y(t) = f x (t - A) h(A) dA (5.1)
o
correlation function Rxy(t ) of the input and output signals x(t) and
y (t).
expressed by
AIternateIy
Rxy (T ) = E {x (t ) y (t )} (5.4)
103
00
Rxy(x) = f Rxx^T — X ) h (X ) dA (5.5)
k is norma I Iy <
10
where
For any arbitrary random signal input x(t) one can identify the
Rxv (5.8)
h = [Rxx] AX
mean that the input signal bandwidth is very much greater than the
Rxx = SQ 6 (t ) (5.9)
leads to
RXy(T) = SQ h (t ) (5.10)
re-arrangi ng
n(t) and the respective 'weighting function' of the signal and the
,0 0 CO
Rxy(T) = £ Rxx(x - A) h (A) dA + £ Rxn (x - A) hn (A) dA
(5.12)
signals x(t) and the noise input n(t). Now provided these two signals
Rxn<T - A) = x n (5.13)
Rxy(T ) - So h (x ) + c (5.14)
105
Thus by ensuring that the input signal has a zero mean the
of a clock period AT. In this case NAT the sequence time is 15AT
and the signal pattern repeats itself after this time interval.
of this P.R.B.S. are also shown in Figure 30. Leary (References 56,
X(t)
P •R , B .S • ( N = 1 5 )
+ a
- a
Aï
1 SEQUENCE = SEQUENCE
R x x (t) N. ¿T REPEATS
Autocorrelation
P o w e r (log)
k
power spectrum
(sin x / x ) 2
3 dB
Shaped spectrum
1 >
1
1
1
1 log f
0.45 f f 2f 3f
c C c c
F igure 30 Te s t s i g n a l c h a r a c t e r i s t i c s
i CORRELATOR
SPECTRUM DISPLAY
(fissa P ■R •B IS I
fSpI NOISE
.u,J GENERATOR
o fa fo ’o l [ -'—
^ J>
STRAIN GAUGE
BRIDGE AMPLIFIER
m w
X -Y RECORDER
Pl a t e 7 Id e n t i f i c a t i o n hardware
108
experimental test run that was faulty, particularly when carrying out
system input and output signals are sampled on-line via an analogue
with that obtained by convoluting the system’s input with its impulse
H.P. 7DTDC, and using this as the P.R.B.S. input to the system in
P iR •B iS ì
NOISE MILLET RETURN LINE
GENERATOR
E
::rC - r l v HOPPER
*• V
* v'MILLET--',
, ^ f • ,
POWE R >'^ .‘ ‘ / --
AMPL IFIER
^ s'
. - i - j 1 7 < , v *• - FLOW CONTROL VALVE
zn
PERTURBING / X
TRANSDUCER
STRAIN GAUGE
FLOW METER
TO NOVA
A/D CONVERTER
OUTPUT
SIGNAL
that a sequence length of about five times the longest time constant
that more reliable results are achieved when (AX) is set equal to
(AT).
particularly where the signal to noise ratio is very poor. This may
normally.
5.3.1 Introduction
soon complemented by a small mini, i.e., NOVA 1200 with 8K-16 bit
was used. Just prior to completing this project a large Uni vac 1106
functions to be computed:
size, so that the maximum amount of memory was available for data
digital converter, the teletype reader, the high speed paper tape
channels, representing the input to the system and the response from
only limited by the core size, while the latter is limited by the
amplitude to limit the ordinate to six and a half inches across the
roll form.
Firstly, the data is loaded into arrays, counted, and the number of
Section 5.2.1.
Figure 33 and have been programmed from the digitised versions of the
10.5.1.
the chute system had non-linear effects due to the motion of the flow
out.
F igure
33
Op e r a t i n g
system
subroutine
hierarchy
ir
T
jk
6.1 Genera I
cross-correI ation analysis the impulse flow responses for the system
linearities were associated with the motion of the flow control valve.
run had been completed and the data collected, the package also
digital filtering, the removal of any bias; and made the data itself
or paper tape.
6.2 Apparatus
times for the uniform flow case. The flow control gate was used to
set the required datum flow of material, while the perturbing gate
was actuated by the P.R.B.S. generator via the power amplifier and
meter located near the end of the chute detected momentum changes of
Both the cross-correI at ion signal RXy(T) from the correlator and
(N), clock period (AT) and amplitude (A) could be varied over a wide
range. This low energy signal was amplified in the Philips power
amplifier (GM 5535) before passing to the Philips solenoid (PR 9270).
This was capable of exerting a maximum force of 3.5 newtons per amp
was used.
delay points being displayed at any one time and with a total delay
was directly tied to the selected correlation delay time and the
total range of frequency displayed was always two decades. Both the
correlator and the spectrum display unit were fitted with direct
faci Iity.
120
6.3 Procedure
datum flows and straight chute s Iope angIes were incorporated in the
test programme.
main flow control gate was preset to a fixed datum flow rate and the
satisfied that the result was stable and repeatable. At this point
the discrete correlation function in dot form and the spectral function
the tests, the final experimental run using the correlator package was
package. This duplication made the results not only useful for cross
correlator package.
insight into the flow characteristics associated with this valve motion
chute system and the results compared to those obtained from the
stepped, open and shut, while the main flow control gate was open
frequency square wave generator, while the flow response was recorded
no main flow and with a main flow of 4 kg/min were tested. In order
D.C. level of the main flow for that test. To improve the keeping
solution.
from the strain gauge bridge amplifier to the U.V. recorder. This
frequency.
flow through the bin orifice during the motion of the perturbing gate.
was used. The camera incorporates electronic speed control and has
projected onto the glass top table. Since a reference grid was
film was projected one frame at a time onto the glass top table a
selected grains for each time increment. The framing rate was
precisely known since the film edge had timing pips as detailed below
neon discharge tubes that produce timing marks on either side of the
film edge. A timing light generator to provide one lamp with 50Hz
pulses and the other with 1000Hz pulses was designed and built
Weston Mark IV light meter of the reflective reading type. The more
appropriate spot type light meter was not available necessitating some
trial and error to obtain reasonable exposures using the Weston meter.
observe the grain stream at the orifice exit, at the chute exit and,
depending on the framing rate required and the observation time sought.
trigger one of the timing lights (described above) in the high speed
The high speed cine photographic analysis of the flow through the
Figure 36 and its stroke was halved. The effect persisted and the
high speed film revealed that the grain profile flowing out of the
PERTURBING GATE
AFTER SHARPENING
F igure 36 Pe r t u r b i n g gate e n d -s e c t i o n
MODIFICATION
OPENING CLOSING
PERTURBATION PERTURBATION
F igure 37 Gr a i n f l o w s h a d o w p r o f i l e s s h o w i n g
EFFECT OF VALVE PERTURBATIONS
extracted from the grain profiles shown in Figure 37. The overshoot
co
co
LU
LU
CC
I—
CO
Figure 38 St r e a m t h i c k n e s s t r a n s i e n t r e s p o n s e c h a r a c t e r i s t i c s
AT BIN OUTLET RESULTING FROM VALVE PERTURBATIONS
transparent hopper.
In order to observe the flow pattern above and around the bin
blackened and plain millet seed was used to aid in the photographic
shutting of the perturbing gate on the '2D' hopper as was used for the
129
Pl a t e 8 "Tw o -d i m e n s i o n a l " h o p p e r
i.e., only the perturbation, and with main flow present were tested.
Using a reference grid behind the falling grain enabled the grain
closing.
130
analyse the grain flow from the main hopper orifice. in addition,
in the particular grain that was being plotted on the glass top table
disappearing from view into the body of the hopper, requiring the
the flow pattern within the funnel flow bin but this was considered
the flow behaviour through the chute is quasi- 1inear. That is, for
uncorrelated with the test signal so that the simpler Equation 5.14
of chute inclination angles for the straight chute each with eight
initial flow settings. This linear modelling in the time domain was
function RXy(t ).
inclined chute under datum flows ranging from 12mm to 57mm valve
for inclinations of 35°, 40° and 45°. Some general observations can be
made.
2. For all flow rates and inclinations tested the flow experiences
during the transient flow through the orifice of the bin and is
4. For the 35° inclination with its relatively high flow velocity
the response shows more oscillatory behaviour for all flow rates.
variations in the preset flow. This last fact supports the quasi-
rp . -
in
8.9 kg/min.
•l#..**#**««*,
in
'•‘• / . i *• * •*, .«.* • t - • ••**»'*' *‘ *,#**#l**^l «Nil •**«•••*
6.4 kg/min.
Va l v e o p e n i n g m m
CO
A m p l i t u d e v /d i v
• **«• •••«*• « CO
CM
................................................................. ^ •• ••h ,
*• • ....... .
4.0 kg/min.
-«-•«• 4##. » -
m
h- CM
...... .. *•...+
or
^ •••• #
1.6 kg/min,
•'♦,**•* *»•*»
nj
• U«#*•<»,< /
r 1 I
0.0 0.25 0.5 1.0 1.25 1. 5
T ime ( t ) $e c ,
W + *««r1 ....... ..
r
in
8.9 kg/min.
. 6.4 kg/min.
/d i v
mm
.... ••• .... 00
•»« ....
-iI fV)
opening
v
Am p l i t u d e
Va l v e
*'**%#*^##«%\**#»*#<#*4
CM
H1■«*<‘n ^ *nr‘,*t>»<'t‘><**» »**'*••' ro
f
• ** •• • « • • • * *
4.0 k g /min.
LD
V S ...* * ...... 04
CT>
* *./* ***.... .....«•»••».» ,...**•*** **..«.**.
1.6 kg/min.
•— ............#
r T
------- 1
0.0 0-25 0.5 0.75 1.0 1.25 1.5
T ime ( t ) $e c ,
..... r
4 •« in
8.9 k g / m i n .
* . . . ...»
in
** Bfc . . ^ • •• •• •• • •• /•' 4 in
6.4 kg/min
/d i v
mm
........ •4 •• . _ •. 4- CO
oo
opening
v
Am p l i t u d e
Va l v e
.*....... Osi
w. ••• 4*« H CO
4.0 kg/min
in
CN
1.6 k g / m i n .
%* ....s*•
. CM
r T
0.0 0.25 0 .5 0.75 1.0 1.25 1.5
T i m e (t ) S e c ,
travel through the chutes. The approximate transit times for the
three inclinations 35°, 40° and 45° are 0.70, 0.77 and 0.89 seconds
this fact for an open channel bulk granular material flow meter is
ind icated.
7. Since the flow through the bin orifice under opening and closing
shows the h(x) records of the straight chute for a flow rate of
4.0 kg/min for inclinations of 35°, 40° and 45°. Clearly the pattern
••• 0 = 45°
kg/min . flowrate
A m p l i t u d e v /d i v
• • * • ,
0 =40°
»•V
.............. I,
» • ..............Z’““
••
4.0
0 =35°
*,••••••••••• .•••j
r ---------- 1----------
1 ------------------------------------------------ j
the hardware correlator to obtain the Bode magnitude and phase versus
phase plots.
1. The frequency response based on the -3dB point for each inclination
as the zero asymptote the spread in the -3dB points was from
This means that step flow adjustments with a period of less than
response of the form e^"*". Consequently this type of plot does not
4.0 k g /m i n .
FLOWRATE
-18Q
-360
LU
CO
<
X
Q_
-540
-720
-900
•5 5 50
LOG FREQUENCY Hz
fourth order.
provide a useful measure of the transit time and hence average velocity
enabling both the system bandwidth and the behaviour within that
bandwidth to be analysed.
its inclusion in this comparison. Each of the three chutes were tested
of known curved form each operating at a datum flow rate of 4.0 kg/min.
4.0 kg /m i n . flo w r a t e
***t»«• * •*
•— « . ,1
,* / **•«•** ** , ,*• *' 1
/d i v
v
PARBOLIC CHUTE
Am p l i t u d e
•**•» .. „ i*•
ai••
.• § ,*
.... * •*
«• • ♦«•**•••**«*
OPTIMUM CHUTE
( t = 0.0)
• ><••• ***.,.,....•••.. V.../**
••• •
r ---------- 1---------- T 1
0, 0 0.25 0,5 0.75 1,0 1.25 1.5
T ime ( T) Sec .
again present.
that for stable flow the designed flow rate would need
Va l v e opening mm
9V l
147
"7-2.4 Bode Diagrams for the Optimum Chute with Varying Flow Rate
1.6, 4.0, 6.4 and 8.4 kg/min the respective bandwidths are
a number of respects. The optimum chute profile has the fastest transit
chute has the highest bandwidth and the most oscillatory response
transient behaviour not obvious from the h(x) curves, can be gained by
MAGNITUDE dB
ratio of one, with a 25mm datum flow setting, the step response
trends observed above again apply; the ’dip’ is much more pronounced
and further to the right indicating the longer transit time. The
significant difference being that the 38mm datum flow rate now
flow patterns. Again the 25mm datum flow rate corresponds to the
12mm datum flow has the most attractive step characteristic. The
25mm datum flow corresponds to minimum descent time but the levelling
NORMALISED AMPLITUDE
F i g u r e 49
I I I 1 > I I ¡ I i l l 1 I » I 1 I 1 I i 1 f I « 1 I «- » J I ^ I 1 1 1 I f
= 38 M M
r~ n ~ r- T - r » r i i i i m t i -t t t t t m i i i T T T 't i r r n i t i n —
0.8 0.9 1.0 1.1
TIME ID SEC.
U1
N3
153
F igure 51 Step r e s p o n s e s t r a i g h t c h u t e 12 mm fl o w
NORMALISED RMPLITUDE
154
Figure 52 Step response straight chute 25 mm flow
NORMALISED AMPLITUDE
Oí
influences in the impulse response. The 38mm datum flow rate was
reduce the very large oscillatory behaviour and should not strictly
than one. It should be noted that the curve for 45° inclination is the
same one as was mentioned earlier and because of the necessary pre
chute is half that for the straight chute inclined such that the
setting or 4.0 kg/min. The designed flow rate exhibits the minimum
NORMRLISED AMPLITUDE
F igure 54 St ep r e s p o n s e c h u t e s of know n fo r m
15^
NORMALISED AMPLITUDE
again implies that for 'fast* stable flow the design point ought to
transient mode.
Using Equation A5.14 the computed step response for the straight 40°
also shown in the same figure. The ’dip’ and the rise-time
the transport lag, before the flow finally increases to its new steady
state value.
the response to any input. The corollary, for predicting the required
The numerical difficulties are associated with the use of Equation A5.15
repeated here
r J ’ 1 -
1
x(j) 2 H(j-k) x(k) (A5.15)
H(0)
AA k=0 repeated
j = 0, 1, ... m-1 .
NORMALISED AMPLITUDE
162
F igure 56 Co m p a r i s o n of c o m p u t e d a n d m e a s u r e d step
RESPONSES FOR STRAIGHT INCLINED CHUTE
163
and the impulse function model for the optimum chute for the
designed flow rate of 4.0 kg/min, the required control valve motion
step increase in flow above a given datum flow rate. The curve
approximately three times this required step and then be closed down
•
1 4 1 1
.1 ’.2 .3 .¿I
• ELAPSED TIME (SEC)
•
•
ft
•
•
'**••••• ft------------------------------------- ------
•ft •
HOPPER
ORIFICE
INITIAL
DATUM PREDICTED
INITIAL
OPENING
TARGET
DATUM
the difference between the desired and the actual hopper chute
grain flow from the bin orifice. Figure 37 and Plate 9 depict the
time the stream thickness in reducing to the new steady state value
behaviour the perturbing valve end section was ground down to a fine
grain forces within the hopper and discharge orifice resulting from
166
discussed below.
the bin have been plotted as described in Sections 6.5 for typical
contours at the bin orifice are plotted for the grain stream flow
at the bin orifice have been computed and four sets of velocity
profiles are shown in Figure 59. The four flow regimes depicted
position.
pos it io n .
For each flow regime two profiles are shown, the short dashed
5 milliseconds after exit from the bin outlet and the long dashed line
opening CLOSING
PERTURBATION PERTURBATION
VELOCITY METRES/SEC
VALVE MOTION
for all cases after the grains have fallen for 20 milliseoncds
grain flow from the bin orifice to be determined before, during and
deI ineated.
felt that any additional insight into the flow patterns immediately
valve motion influenced the grain flow well into the body of the
illustrating the shift in the right hand shear plane when the
slip occurs between the two pairs of lines associated with the same
datum point at the orifice. As for the case discussed in Section 7.5.1,
the path traced by the individual grains was plotted for successive
frames and from these the velocity profiles across nominated horizontal
pos it io n .
above the bin orifice. The approximate uniform velocity profile before
o
in
m m /sec
o
o
o
m
□ BEFORE CLOSING
O AT CLOSING
> AFTER CLOSING
O AT OPENING
VALVE OPEN
VALVE CLOSED
above the bin orifice and the effect of the valve motion is becoming
steady state velocity profile with the perturbing valve closed now
the grains are leaving the influence of the hopper and are becoming
of 18$ while the far side grain stream velocity drops only 7$,
the average grain stream velocity. At the centre line this increase
in velocity from the 100mm point to the bin orifice within the shearing
section of the body of the hopper. The average centre line velocities
o
O
m m /se c
O
O
CN
O
O
co
2 X FULL SIZE lu
UJ
(X
LU
O
VALVE OPEN
VALVE CLOSED
□ BEFORE CLOSING
O AT CLOSING
> AFTER CLOSING
O AT OPENING
LU
FULL SIZE LU
cr
LU
U
O
O
m
o
o
LD
m m / sec
C] O
O
r-
o
o
oo
□ BEFORE CLOSING
o AT CLOSING
A AFTER CLOSING
0 AT OPENING
The !2-DT hopper results have shown that the flow behaviour
above the bin orifice during perturbing valve movement. High speed
Un iform FIow
6. Even at the small scale tested the optimum chute geometry exhibited
Transient Flow
foI lowing:-
enabling both the system bandwidth and the behaviour within that
bandwidth to be determined.
179
product.
the other hand, the results clearly show that the optimum chute
stable flow.
was relatively minor and the use of the linear systems theory
not due to the perturbing valve end section but rather to the
hopper.
10. The overall performance of the optimum chute was clearly superior
to this study and would broaden the scope for application of the
present findings:
presented is supported.
cohesive materials.
182
Curved Chutes".
Mech. & Chem. Eng. Transactions, I.E.Aust., Vol. MC3, No. 2 (1967).
Publ.73/4 pp.156-161.
183
1972 pp.149-153.
p p .39-45.
13. CHARLTON, W.H. and ROBERTS, A.W. "Chute Profile for Maximum
1970 pp.292-294.
AIgori thms".
System Design".
Fortran".
AIgori thms".
Process Identification".
34. PHUDKE, M.S. and WU, S.M. "Function and Noise Model of a
Private Communication.
Industry".
Private Communication.
39, BRYANT, G.F. et al. "The Design and Tuning of Steel Mill
Control Systems".
Private Communication.
40. BENDAT, J.S. and PIERSOL, B.A.G. "Random Data: Analysis and
Measurement Procedures".
Wiley-Interscience 1971.
187
Wi ley-Interscience 1970.
Seguenc ing".
Control System".
pp.205-211.
52. GOODWIN, G.C. et al. "Optimal Test Signal Design for Linear
53. GOODWIN, G.C. et al. "Optimal Test Signal Design for Linear
Conf. on Comp. Aided Cont. Syst. Design Cambr. April 2-4, 1973.
Sequences".
Techn iques".
Dev ices".
of Particulate Solids".
pp.224-228.
; 90
74. OKI, K. et al. "A New Method for Evaluating the Size of
10.1.1 Introduction
measuring grain stream velocity under the same design criteria was
■
i .e ., , mv.)
(m = —
v
turn, gives an output signal that for small angular displacements would
COUNTE RWEIGHT
__ >iti>_____ a n I j L __
l---- - ------------- 1------- .- 1 -7^------ 1
RESERVOIR
on the plate and more fundamentally by the accuracy with which the
Nolte (References 64, 65) has been associated with the design
His refinement in ensuring that the sensor moved only vertically and
vertical velocity. The material then falls a set height within the
RIGID FRAME
194
SPREADER
BELLOWS
SEAL
SPRING
HOUSING
L PNEUMATIC
TRANSDUCER
BOTTOM FLANGE
through the meter, causes a torque in the supporting shaft which can
Material enters the meter and strikes the cone at the top of the
meter and is then directed into the top of the rotating vanes near
the axis of rotation. The angular rotation of the vanes causes the
torque only and ignore any bending moments which may be present.
As with most such systems measuring small toques from rotating shafts
with the usual high noise present due to the slip-ring subsystem,
time delay.
measuring range is limited and the in-meter transit time would make
this investigation.
ROTOR
IMPELLER
SYNCHRONOUS
MOTOR
considered and some of the more promising devices were tested. The
FIowmeters
drift, requiring frequent recaI ibration and was thus not considered
special precautions for storage and the need for more than the usual
PIEZOELECTRI
TRANSDUCER
IMPACT
PLATE
Velocity Detectors
decade counter with each digit driving a separate light via switching
the reflection of light from the perspex chute, the randomness of the
201
the grain head provided the mass flow rate was constant
i .e.
_H_ = Vo
Ho V
the head and, knowing the inlet conditions and assuming uniform flow,
Problems with light reflections off the perspex chute and with
Figure A 1.10. This device would enable the surface and subsurface
friction bearings and low rotational inertia, together with the need
overshadowed this unit and further work on the wire paddle wheel' unit
was deferred.
experimental investigation.
deflection away from the grain stream for a high datum flow rate
n = 1, 2, 3 respectively.
devices have been described in the literature (References 69, 73, 74)
which rely for their operation on the Doppler shift principle. The
high equipment cost associated with both laser and microwave based
target is given by
.. = V(d fd (A 1 -1 )
i" 2foj cos0
fd = Doppler frequency
206
Figure Al.ll Calibration curves for the strain gauge flowmeter
0 = a n g l e b e tw een t h e c e n t r e line of th e e m it t e r ve
velocity di r a t i o n .
vt = Kfd (A i.2)
wave Iength A =
fu)
requ ir e d .
The prototype was built at a cost of less than S100 with the
1. the range
Results of these tests are too lengthy to be presented here and will
gear used previously. Frame 3 shows the fitting of the plastic tubes
present. Its portabi Iity and digital readout are clearly evident
wheat covered rubber belt. Frame 2 shows the layout of the printed
as part of the testing programme for the uniform flow of grain down
chutes.
required between the meter and the ultraviolet recorder. As none was
• ©
CALIBRATING
MOTOR SPEED TEST RIG
CONTROLLER
Fr a m e 1 Ex p e r i m e n t a l rig Fr a m e 2 Fr e q u e n c y monitoring equipment
LED READOUT
- .^>■.^ ■ , . , to.
¡ULTRASONIC TRANSDUCERS!
¡SWIVEL MOUNTED ssv !
?7^Ns^^wweTWi'* w * œ s
‘fe
■S I
••îsittWS» PC BOARDS
IN POSITION
ULTRASONIC
TRANSDUCERS
Í WHEAT ON"BELT
Figure A 1.12 shows the circuit of the filter with Ri, R2, Ci and C2
were achieved. These were considered adeguate for the range of tests
degree of filtering.
r2
FREQ R2 R2j r3 Cl c2
5 750 k 150 k ,27 .056
10 768 k 159 k .1 .0208
25 768 k 159 k ,047 .000828
50 768 k 159 k .022 .000414
F i g u r e A l , 12 F lo w meter f i l t e r
212
The two neon timing lights fitted on either side of the film
expose, the edges of the film during a filming run and thus provide a
useful to have the motion of the flow control gate, in this case the
perturbing gate, actually switch one of the timing lights on and off
square wave applied to the other timing light the film is marked into
one millisecond sections. From these timing marks the portion of the
film correctly exposed after the camera has reached operating or set
unit provides the timing reference for all subsequent film analysis
film spools up to 120 metres in length. The optical head used was
suitable for exposing 16mm full frame positive or negative movie film
and would accept any standard !C T mount lens. The quoted resolution
edges.
an f stop lower than this. The exposure time for this camera is
calculated as follows
For a given light level and film speed this exposure time
relatively high f stop setting which can only be obtained, all other
parameters held constantly increasing the film speed. The Ilford Mark V
negative film rated at 400 ASA was forced developed to double that ASA,
analysis of the projected image when the film was forced beyond this
value. It should be noted that for framing rates above about 1000
214
spot type exposure light meter becomes essential when the field to be
CorreI at ion
caIcuI at ions
N-j + 1
R (jAT) = . E x (i)x (i+j ) (A5.1)
J N-j ¡=0
signals under test, use is made of the fact that the power spectral
f = kfc
(A5.2)
The final digitised formula for the power spectral density becomes
m-1
G (k ) 2AT R(0) + 2 Z D(r)R(r) cos(-^-^-) (A5.4)
m
r= 1
of any time delays and any phase shift. In a similar manner to the
kfc
f
m
Q(k)
0xy(f) - tan 1 (A5.6)
C(k)
and
m- 1 ~]
4AT trr^
Q(k) Z (r) s in (A5.3)
m
r-1
The odd and even parts of the cross-corre Iat ion function M r ) anc
and
B(r) = i [RXy (r ) - Ry x (r )] (A5.10)
for r = 1 , 2, ... m - 1
Syster Computation
for the input and output signals may oe expressed in terms c-f the
HO)
Rxy(O)
Rx (0)AX
(A5. H )
J-l
h (;N 1 Rx v (j ) I Rx (j-k)H(K)
Rx (0) ' AX
L k=0
for j = 1, 2, ... m - 1 .
H = [Rx]_ 1 [Rxy]
AX
y (c) (x(O)H(O))AX
1 j-1
x(j) Z H(j-k) x(k) (A5.15)
H (0 ) AA
k=0
N-1 N-1
G(ja)) = AT E H( iAT)cos(iwAT) jAT E H (iA T )s in (iwAT)
i=0 ¡=0
(A5.16)
data processing tool. The Data General User Manual (Reference 77)
was found useful here. The intimate interaction between operator and
this nature on a newly installed system having only slow speed input
consumi ng.
logging and data reduction for student projects within the Department
of Mechanical Engineering.
219
etc.
hierarchy
223
-00 1 SI GAM
; DATA A C C E P T E D F R O M T E L E T Y P E , H I G H S P E E D P A P E R
* T A P E R E A DER, C A R T R I F I L E , A TO D C O N V E R T E R .
5 N O R M A L I S E D DATA M A Y BE P L O T T E D OR LISTED.
; S I G N A L C O M P U T A T I O N S I N C L U D E A U T O C O R R E L A T I O N OF
I I N P U T OR OUTPUT, C R O S S C O R R E L A T I O N , P O W E R S P E C T R A L
; FUNCTIONS.
; SYS T E M C O M P U T A T I O N S I N C L U D E I M P U L S E F U N C T I O N
; C O N V O L U T I O N (TO P R E D I C T T H E OUTPUT)
; I N V E R S E C O N V O L U T I O N CTO P R E D I C T TH E INPUT)
; F R E Q U E N C Y R E S P O N S E . (B O D E D I A G R A M )
5 PRBS PA R A M E T E R S : -
I DT = T I M E D I S P L A C E M E N T ( D I G I T INTERVAL)
; S E L E C T N O I S E B A N D W I D T H TO S UIT S Y S T E M
I U S U A L L Y TEN T I M E S T H E S Y S T E M ' S B A N D W I D T H .
5 P R B S P E R I O D TO B E A P P R O X . 1 . 5 * S E T T L I N G TIME.
5 T C P R B S )=1.5*TCSETTLING)
; TC S E T T L I N G ) = T ( D E A D T I M E ) + 4 * T ( L O N G E S T
i T I M E CONST. OF I N T E REST)
; S E L E C T DT < 0 . 5 * T C S H O R T E S T T I M E CONST.
J OF I N T E R E S T )
I CORRELATION PARAMETERS:-
5 N = NO OF D ATA P O I N T S (SE L F COU N T I N G )
; M = NO OF C O R R E L A T I O N L A G S (SET TO N/10)
I DELAM = T I M E I N T E R V A L B E T W E E N L A G S
5 ( SET E Q U A L TO DT)
. TITL S I G A N ; G E N E R A T E O P E R A T I N G SY S T E M
. ENT W S A , G E T C , P U T C , . C R L F , . T X T T , .N E X T P , . R E S L T
. ENT T E M P O , T E M P 1 , C 1 2 , C I 5,N U M , N U M P 1 , J , M 1 , N
. ENT S U M , S U M P 1,R X X O , D E L A M , O N E , T W O , S T A R T
. ENT C 5 P 1,C 9 P 1, M S K , N L P 1» N 1 , I C H N , S R A T E , S C N T R
. ENT • PLOT, N L , C 9 , C 5 , W, WW, D, DD, K,M, T W NTY
. f:n t a m e s o , a m e s i , . g t c t , . g t c p , . g t c c ,p i
• ENT DEG, LNLOG, TEN
. E X TN USAR,TTIN,TTWR,RR,H,G,PLOT,GXY,,CON,FRESP
. EX TN RAY 1, RAY2, RAY 3, RAYA, RAY 5, RAY 6, RAY 7 , RAY 8
. F^XTN DEL, DATAIN, INPUT, OUTPUT, AUTOI , AUTOU, C R O S S
. EXTN MESO,MES1,MES2,MES3,MESA,MES5,MES6
. FXTN MES7,MES8,MES9,MEAN,A2DC,LIST1,LIST2,QUES
•EXTN TTYR,HSPTR,CARTF,A2DI,CRLF,TXTT
. EXTN N E X T P , R E S L T , N P A G E , H E A D , F I N T , F ENT
. p:x t n G T C T , G T C P , G T C C , INPOW, O U T P O W , X P O W
. EXTN C O N V,INCON,ICON,FREQ,STRTM
224
000 2 SI GAM
• ZREL
00000-177777 WSAS USAR % P O I N T E R TO U S E R A R E A
0 0 0 0 1 - 1 7 7 7 7 7 GETC: TTIN 0
9 P O I N T E R TO G E T A CHAR« R O U T I N E
0 0 0 0 2 - 1 7 7 7 7 7 PUTC: TTWR 0
9 P O I N T E R TO W R I T E A CHAR. R O U T I N E
0 0 0 0 3 - 0 0 0 0 0 l-.TTIN: TTIN 0
9 P O I N T E R S IN P A G E Z E R O TO R O U T I N E S
0 0 0 0 4 - 1 7 7 7 7 7 • RR: RR
0 0 0 0 5 - 1 7 7 7 7 7 .Hi H
0 0 0 0 6 - 1 7 7 7 7 7 • Gl G
0 0 0 0 7 - 1 7 7 7 7 7 . GXY: GXY
0 0 0 1 0 - 1 7 7 7 7 7 • QUESJ Q UES
0 0 0 1 1 - 1 7 7 7 7 7 •GTCTl GTCT
0 0 0 1 2 - 1 7 7 7 7 7 •GTCPI GTCP
0 0 0 1 3 - 1 7 7 7 7 7 .G T C C : G T C C
0 0 0 1 4 - 1 7 7 7 7 7 • TTYR: TTYR
0 0 0 1 5 - 1 7 7 7 7 7 • H SPTRi H S P T R
0 0 0 1 6 - 1 7 7 7 7 7 • CARTFl C A R T F
0 0 0 1 7 - 1 7 7 7 7 7 •A 2 D I 2 A2DI
0 0 0 2 0 - 1 7 7 7 7 7 . A2DC2 A2DC
0 0 0 2 1 - 1 7 7 7 7 7 •PLOT: PLOT
0 0 0 2 2 - 1 7 7 7 7 7 • CON: CON
0 0 0 2 3 - 1 7 7 7 7 7 . ICON; ICON
0 0 0 2 4 - 1 7 7 7 7 7 . FRESP: F R E S P
0 0 0 2 5 - 1 7 7 7 7 7 . D ATAIN 2 D ATAIN
0 0 0 2 6 - 1 7 7 7 7 7 .L S T 12 LISTI
00027-177777 •LST2; LIST2
0 0 0 3 0 - 1 7 7 7 7 7 •CRLF; CRLF
0 0 0 3 1 - 1 7 7 7 7 7 . TXTTl TXTT
00032-177777 •NEXTP2 NEXTP
0 0 0 3 3 - 1 7 7 7 7 7 • RESLT2 R E S L T
00034-177777 AMESO2 - 1 «
9 DATA SOURCE
00 0 3 5 - 177777 AM ESI 2 - 1 0
9 DATA L I S T I N G
0 0 0 3 6 - 1 7 7 7 7 7 AMES2I - 1 0
9 FREQUENCY RESPONSE
0 0 0 3 7 - 1 7 7 7 7 7 AMES3: - 1 0
9 C O N V O L U T I ON
0 0 0 4 0 - 1 7 7 7 7 7 AMES4: - 1
•
9 A U T O C O R . OF INPUT
0 0 0 4 1 - 1 7 7 7 7 7 AMES5: - 1 A U T O C O R . OF O U T P U T
0 0 0 4 2 - 1 7 7 7 7 7 AM E S 62 - 1 0
9 CROSSCORRELATION
0 0 0 4 3 - 1 7 7 7 7 7 A M E S 72 - 1 •
9 IMPULSE FUNCTION
0 0 0 4 4 - 1 7 7 7 7 7 AMES82 - 1 •
9 POWER SPECTRAL ROUTINES
•
0 0 0 4 5 - 1 7 7 7 7 0 NL: -10 9 NL
0 0 0 4 6 - 1 7 7 7 0 4 N L P 1: - 7 4
0 0 0 4 7 - 1 7 7 7 6 6 C9 2 - 12 p C9
0 0 0 5 0 - 1 7 7 7 6 6 C 9 P 1: - 12
•
0 0 0 5 1 - 1 7 7 7 7 3 C5: -5 9 C5
0 0 0 5 2 - 1 7 7 7 7 3 C 5 P 1: -5
l 'u u J O y ( . .1u \w*
m T E M P O 2 V-<
A 0
9 TEMPO
A A A C / i_ A A A A A A
V-' V ' v.» O H TEMPI:
A
V
0
9 TEMPI
A A A C C _ A A A A
\ Zj Zj V_< V.* 1
f Q
C 12: 12 9 C 12
v'uv.» J C v.»v-' v.' v.» 1 C 1 5: 15 0
9 C 15
1 NUMI MU M
A A C 7 ^ rs A A A A
V.I J f “ v.» V_' V.» V
_*V \w* 9
0 0 1 0 3 - 0 0 0 0 0 0 DELAM: A
KJ ; DELAM
V.» U iU H
A
A A 1 A A A A A A A
U U L » V.* V.» i ARITHMETIC CONSTANTS
0 0 1 0 5 - 0 4 1 0 7 1 DEG: 57* 2 9 5 7 7 9 5 •
9 RADS. TO DEGS.
00106-045670
0 0 1 0 7 - 0 4 0 4 6 2 PI: 3*14159265
00110-041766
00 1 1 1 - 0 4 0 4 4 4 L N L O G S 2.30258509 •
LN TO LOG
00112- 1 53543
00113- 0 4 0 4 2 0 OME: 1.0 ; OME
uu 1 1 *4
A A 1 | /|_ A A A A A A
00 1 1 7 - 0 4 0 6 4 0 TEM:
A A $ A A _ A A A A A A
0 0 1 2 1 - 0 4 1 0 2 4 TWMTYS cvy • u
A O AN
U U 1C C
A A | O O — A A A A A A
V.»V^L*UV/V*
•N R F L
0 0 0 0 0 * 0 6 2 6 7 7 START: I0R S T •
0 RESET START
A A A A 1 f A / ; O A 7 7
L » U V-Mw» 1 Y D O u /f HALT •
0 WAIT # PRESS CONTINUE
A A A A A t
V. » VJ\. » V ' <o
A A
U U U U J C
O-*
JSR § • NEXTP
A A A A Q 9 A A ^ A O | _
u u u u O u u O u O 1 JSR § • TXTT •
0 STARTUP MESSAGE
0 0 0 0 A * 177777 STRTM
A A A A C
V.' V» S.» J
• A A A /i A O
U U M U J JMP . +3
A A A A / :
u u v ^ u C
t A A / ;A A
u u u u O
1 _
1 JSR § . TXTT
«
0 RESTART HEADING
00007*177777 HEAD
V.IUV.I
i r> #
1l'
a *1/ i « ^ ^ .
L 'O ^ U U V . ' LD A 3# WSA *
0 SET F O R M A T F O R
00011*020066- LDA c*ww : FLOATING POINT
00012*041521 STA 0# 121# 3
r \ r \ n <o 1
V. < U U IJ
Q A A 7
\J t\J C .\ J \J LDA 0#DD
00014*041522 STA 0#122# 3
000 15* 177777 FI NT
•
0 INI T A L I S E FLT* PT.
A 1 t 1 A/> 1 A _
JSR §.QUES 0 SET THE RUN S T R E A M BY
00 0 1 7 * 1 7 7 7 7 7 MES0 0 A S K I N G AL L R E L E V E N T
v.*\~i
A A A 9VA O
Q y.% A A , , ,
*4 1 sj \J LDA 3# PFT 0 QUESTIONS.
v_* ¿
L1 i I iO U U
A A A A 1 9 1 | T A A A
W ADD 0# 3
A A A
v.1u u
A A
c
9 A / l O A O / j .
c* H u u O * * STA 0 # AME S 0 0 DATA S O U R C E
V.»V.* c J t uu
A A A A Q
V.*
A A 7 / 1A A
f JSR §0.3
A A A O/ l• JSR § . TXTT • AMS. F0LL. QUES.
V. <v . *U u u O u J 11 » 0
00025*177777 MES2
i.>(.'v.'CO•0
A A r t A A £ A 1O _
uCu K 1 JSR § . Q UES
00027*177777 MES1
• DATA L I S T I N G
W
A A A
0 J uA •u^luyOJ
A / i A A O C « .
STA C#AM E S I 0
A A A A 1 1A A £ A (A
v_iv.<v.1O l y y D y I k JSR §• QUES
00032*177777 MES9
00033*040036- STA 0 . A M E S 2 0• F R E Q U E N C Y R E S P O N S E
00034*100015 COM# 0# 0# SNR
00035*000404 JMP • +4
00036*040040- STA 0# A M E S 4 0• AUT0C0R. OF INPUT
0 0 0 3 7 *0 4 0 0 4 2 - STA 0# A M E S 6 0• C R O S S C O R R E L A T I O N
0 0 0 4 0 *0 4 0 0 4 3 - STA 0# A M E S 7 0• I M P U L S E F U N C T I O N
226
0004 SI GAN
00041•006010- JSR 0.QUES
00042*177777 MES3
STA 0 * AM ES 3 « CONVOLUTION
»
00043*040037-
00044*024036- LDA 1# A M E S 2
00045*124015 COM# 1 » 1#SNR
O L ' u H D l u u u l j COM# 0*0# SNR
JMP • +5
A r t A C A • A A A / | A _
STA 0# A M E S 4 I A U T 0 C 0 R « OF INPUT
U O V.' 1 U H U U H C
A A A C
□
1 I A / ,A A / . A .
STA C*AMES6 #• C R O S S C O R R E L A T I O N
00052*040043- ST A 0* A M E S 7 •
# IMPULSE PONCTION
U O U D O U U U*J 1 d
A A A C ' l f A A A A I A
JMP PS
00054*124014 COM# 1#1# S Z R
U O V H 1 V» JMP
A A A C C • A A A / . « A
J O PS
A A / -A I A , .
rs r% c
o c » v.» J O
/L 1
u v D v i o JSR e.QUES
00057*177777 MES7
u u u D o W ^ U W 'i O
A A A £ A « A / . A A /. Q _
STA 0#AMES7 ; IMPULSE FUNCTION
A A A £ I « 1A A A 1 C
u c * v »0 1 i u u u i D C OM# 0 * 0 * SNR
u v . ' o o c u v - ' i ' ^ u J JMP PS
A A A £ A . A
v . ' u v .'00
/. A A / , A _
STA 0# A M E 5 4 #« AUT0 C 0 R . OF INPUT
00064*040042- STA 0*AMES6 #* C R O S S C O R R E L A T I O N
r » A r y / ;c
c» c» o O J
• a a /• a
u u O u
i
1v
A ^ n c
r j
•
« JSR @.QUES
00066*177777 M ES8
00067*040044- STA 0 # AMES8 • POWER SPECTRAL ROUTINES
#
A A A T A
0 C »0 i O
• J A A A
I V U V
1 C
1J COM# 0 * 0 # SNR
A A A T
0 0 0 f
1 f A A A /. | A
1 O O O H I C JMP PS1
A A A T A • A /i A A / . A _
v.# ^ K.» f r z L » *4 %.» v *4 L» STA 0* A M E S 4 #
•
A U T 0 C 0 R R . OF INPUT
00073*040041- STA 0 # AM ES 5 # A U T 0 C 0 R R . OF O U T P U T
00074*040042- STA 0 # AMES6 #• C R O S S C O R R E L A T I O N
00075*024043- LDA 1# A M E S 7
00076*124014 COM# 1* 1 # S Z R
A A A A A
V.« V.* I I
• A A A / i
O U V -'H C v.' JMP PS2
A A
O W
« A A f A /, A A
1 V.HJ U ^ U U
/| A _
t U STA 0# A M E S 4 « AUT0 C 0 R . OF INPUT
#
00101 *040042- STA 0# A M E S 6 5 CROSSCORRELATION
A A
V.* V.I
I A A *
1 V_i C .
A A A /, 1 C
t ' O O H 1 J JM P PS2
A A
i .1u
J A A
1 V.1O
I A A £ A | A
o y u w K *
. n C
" O
1
I
•
• JSR @ •Q U E S
00104*177777 MES5
00105*040041- STA 0* A M E S 5 ; AUT0 C 0 R . OF O U T P U T
00106*024043- LDA 1# A M E S 7 * IMPULSE FUNCTION
00107*124014 COM# 1*1# S Z R
A A
v_< O
t
1 1U
I <A t A A A /| A 7
LM -M -'H o I JMP PS2
«_<\ J I 1 1» t UA LA 'D
A A I I
o iu
£ A 1 A _
JSR @.QUES
00112*177777 MES4
00113*040040- STA C# A M E S 4 A U T0C0R. 0F INPUT
00114*006010- JSR 0 . Q UES
00115*177777 MES6
00116*040042- STA 0* A M E S 6 CROSSCORRELATION
001 17 *0 2 0 0 37-PS2J LDA 0# A M E S 3 CONVOLUTION
00120'024044- LDA 1# AMES8
00121*100015 COM# 0# C# SNR
00122*124014 COM# 1»1# SZ R
r \ r\ « O Q t r \ r \ r \ J \ r \ c\
v v 1 c J c*u JMP •+ 2
00124*000413 JMP . + 13
00125*006031- JSR @ • TXTT ; TIME DISPLACMENI
00126*177777 DEL
A A
\
j \
ji
I O T • ^ O n ,>
1 c / u c u u o
1 «
I LDA 0# G ETC
00 130* 0 2 4 0 0 3 LDA 1# • T T I N
0 0 131 * 0 4 4 0 0 1 STA 1»GETC
0 0 132* 177777 PENT
00133*120000 F DF C GET DELAM
227
0005 SIG AN
00134* 040103- FSTA 0# DEL AM
00135*100000 F EXT
001 3 6 * 0 4 0 0 0 1 - ST A 0# G E T C
00137*126120 ADCZL 1# 1
0 0 1 4 0 *0 2 0 0 3 4 - LDA 0# AM ESO ; DATA SOURCE
00141* 107014 ADD# 0#1 # SZR ; A TO D ?
00142*000403 JMP •+ 3 ; NO
00143*006020- JSR 0. A 2 D C ; YES» J U M P TO SUB.
00144*000413 JMP . + 13
00145*024035- LDA 19 A M E S 1
00146*020040- LDA 0# AME S 4
00147*030042- L DA 2» AME S 6
00150*150015 COM# 2 # 2 # SNR
00151*124014 COM# 1#1 * S Z R
UUIJC 1 C O I A A A / | / > Q
JMP • +3
00153*100015 COM# 0 * 0 » SNR ; AUTOCOR. OF INPUT ?
00154*000415 JMP NO ; NO
1J J
V. H. ’ Ou CÍD V.I JSR 0 . D A T A I N; INPUT D A T A C A L L
0 0156* 177777 RAY 1
00157*006032- JSR e.NEXTP
00160*006031- JSR 0.TXTT ; INPUT H E A D I N G
00 161 * 177777 INPUT
00162*006026- JSR a.LSTl
00163*000156* RAY 1
00164*020035- LDA 0# AMESI
00165*100015 COM# 0# 0# SN R ; LIST. A N D OR PLOT. ?
00166*000403 JMP NO
00167*006027- JSR «.LST2
Ars 1 7 A f « r 7 |
I / v.'
a a a
c ' v_' v j 1 O O RAY 1
00171 * 126120 M 0 : ADCZL 1# i ; S T A R T TO E X E C U T E THE
00172*020034- LDA 0 AM ESO ;
• RUNSTREAM
00173*107015 ADD# 0 # 1 SNR ; # DATA S O U R C E ?
00174*000413 JMP . + 13
00175*024035- LDA 1 AMESl »
A A A 1 A
V*Vȣ-I^*
f A A £ A 7
V. *V. » UU J
1 _
I JSR P.TXTT ; OUTPUT HEADING
0021 1 * 177777 OUTPUT
00212*006026- JSR 0.LST1
00213*000206* RAY 2
00214*020035- LD A 0 » AMES1
A A A
S»'V.* ¿
1 C
1 J
1 1 A A A
1 O V.» U
| C
I J COM# O p SNR ; LIST.
O p A N D OR PLOT. ?
A A A
VJ VJ
1 £
1 O
t A A A
s.» u
/i A 7
J JMP N1
00217*006027- JSR 0.LST2
f > n O O ' 1 •
L.» l < ez. C i L< V.I V.I V.I &
1 O
1 O
*
R AY2 *
00221 *0 2 0 0 6 3 - N 1 : LDA o *m ; SET P L O T T I N G C O U N T E R
00222*040062- STA 0» M 1
0 0 2 2 3 *0 2 0 0 4 0 - LDA O A M E S 4 ; A U T O C O R R E L A T I O N OF
p
0006 SI GAM
00227*006031- JSR » . TXTT
00230*177777 AÜT0I
00231*006004- JSR G.RR »«
J U M P TO C O R R E L A T I O N
00232*000170* RAY 1 »• R O U T I N E
00233*000232* RAY 1
00234*177777 RAY3
A A
UUCJ J
tO
A O A A / iA _
G tK 'n U LDA 0# A M E S 4
AAOOA •I
UUCÜÜ &A A / . A A
u u h u o MEG 0 #0
A A O O * 7 t | A A A A /.
COM 0 #0 # SZ R »•
PLOTTING ?
^ ^ O /I A • A A ‘A /»A A JM P N2
0 0 2 A 1 *0 0 6 0 2 1 - JSR #.PL0T
00242* 0 0 0 2 3 4 * RAY3
00243* 0 2 0 0 4 1 - M 2 J LDA 0» AM ES 5 # A U T O C O R R E L A T I O N OF
00244*100005 COM 0 * 0 # SNR » O U T P U T D A T A C A L L
00245*000416 JM P N3
««
ü u co/i A •Art
*40 u w O£rt»jn wJd JSR ».NEXTP
00 247 *0 0 6 0 3 1 - JSR ».TXTT
00250*177777 AUT0U
A r t oC 1 t r t r t / r r t f t
u u c J 1 uuOuu*i* JSR ». RR #■ J U M P TO C O R R E L A T I O N
rtrtrtCrtl
UUC J c
rt rt rtrtrtrt •
U L 'U C G U RAY2 #• R O U T I N E
rtrtrtC «
C 'U C J O
J« rtrtrtrtÇ O *
UV.H.ICJC RAY2
00254*177777 R AY4
00255*020041- LDA 0* AM ES 5
rt />
r t r t r t C Z 1 1 rtrt /irtrt
(j u c j D I u u *<u u MEG u * u
r t r t r t C 7 1 | rtrtrtrt
u u c j i Iw y u u *» COM 0 * 0 # SZR •
» PLOTTING ?
r t r t O i r t • r t rtrt /.rtrt
ui.ioOv.1 J JM P N3
r t r t r t i 1 * rtrt £ rt o 1 _
WV.ICU 1 v .'ü D w c l JSR ».PLOT
r t r t O A O * rtrtrtrtc/. •
c l ' c UC u u w g j *J RAY4
0 0263* 0 2 0 0 4 2 - M 3 : LDA 0 # AMES6 •
» CR0SSC0RREL. CALL
1A A A A C
COM 0# 0# S M R
00265*000416 JMP M4
rtrt O i £ • r t r t i r t O O «
l/ U tO O Du OG JSR ».NEXTP
rtrt O A 7 • o o A* * ^ t •
o u g D I v^yO uJI JSR ».TXTT
00270*177777 CROSS
00271*006004- JSR » . R R
r t r t r t T O 1 rtrtrt Q t 7 »
VJUG 1 G UUVJ G J O RAY2
r\ri 0 7 *1 • r t r t r t r t l l *
U C C I O 1/yv.iGJJ RAY 1
00274*177777 RAY 5
00275*020042- LDA 0# A M E S 6
rt rt
00276*100400 MEG U » U
00314*177777 IMP0W
00315*006006- JSR ».G «
9 INP U T S P E C T R A L
00316*177777 R AY6 1 DENSITY
00317*000242* RAY 3
00320*020044- LDA G* AMES8
A A
003 2 1 * 100400 MEG V.' # %.«
229
0007 SIGAN
00322 1
1A AAA COM 0# 0# SZR s PLOTTING ?
UUJC J
a a q a o
000451 JMP N5
00324 A
w rtCw £ A A«
kj ¿1 JSR ©•PLOT
00325 uuoJlD RAY 6
00326 u
AA ÌO
l ' D u o c
_
JSR § . N EX TP
00327 ».»wOuO 1
A ft £ ft O 1 _
JS R ©.TXTT
A A AJuU
U U
QA 17 7 7 7 7 0UTP0W
00331 A
oA ur OAuAuZ
u_ * JSR #• G i OUTPUT SPECTRAL
A A ^JJC
nn 177777 RAY 7 • DENSITY
#
u «3*3O
A A Q A O AAAOA O • RAY4
00334 020044- LD A 0# A M E S 8
00335 1
1A A /iA A NE G A A
V» KJ o# o
00336 1
1U AU AAA
U U N>. COM 0*0# SZ R •
PLOTTING ?
00337 000435 JMP N5
O
v.* c* 3
O /i O
v.* ycOwcl" JSR 3» P L O T
00341 000332« RAY 7
00342 uoOuOd
ft ft £ ft n _
JSR S.NEXTP
f t f t £ftft/i_
00343 o u Ouu*! JSR f.RR
00344 U U U C fJ f t *7 ft «
RAY 1
00345 0002 7 2 * RAY 2
00346 A A A A A A
0 0 0 3 3 3
f
RAY4
00347 A A A A O 1
o o O o 31
.
JSR © • TXTT
A A O C A
u u o De# 177777 XP0W
00351 A A A A A AT
O O O O V.» f JSR ©.GXY # CROSS-SPECTRAL
00352 A A A A A A
o o o Oo
t
RAY5 9 DENSITY
A A A C O
u u O O J 0 0 0 346* R AY4
0 0 3 5 4 A A A 3A A C « O O O 3 RAY 6
00355 000341 * RAY7
0 0 3 5 6 AO O A A A C A
O 3 3 3
•
RAY 4
0 0 3 5 7 177777 RAY 8
A A A
020044-
A /:
v.' v.' 3 O 3 LDA 0# AM ESS
00361 1 oA oA * /4i Ao oA «
NEG ft
KJ 0
/>
lw<
0 0 3 6 2 i A A Ao Ao *4 ]
O O
/.
COM 0# 0# S Z R •
9 PLOTTING ?
0 0 3 6 3 000411 JMP N5
00364 A A A A A _ r
V. ' U v w / O c JSR ©.NEXTP
0 0 3 6 5 Ao Ao ^O Ao Ao 1 I .
JSR © •PLOT 9 PLOT MAGNITUDE
0 0 3 6 6 AO AO AV.* 3A C Co t 3 RAY4
0 0 3 6 7 oA Ao CO oA A3 oO — .
JSR ©.NEXTP
A A
^ • ^ .»0
A A A
9 ^ *
A A C A A |
o o Oo e I
_.
JSR @.PL0T •
9 PLOT PHASE
0037 1 A A A A C A «
o o o 3 3 f RAY8
0037 2 A0 O^ A0 A0 /0 ' 3A ^ LDA O.M 9 RESET M COUNTER
00373 040062- STA 0* M 1
00374 020043-N5: LDA 0» A M E S 7 9 I M P U L S E F U N C T I O N CALL
0 0 3 7 5 i1 A0 0A 0A 0A 3C COM 0.0# S N R •
9 R E Q U I R E S RXX A N D RXY
00376 000413 JMP N6
0 0 3 7 7 Ao oA o Ao A C3 . C
JSR e . H
A A /i A A A A A A |A
o o o 3 1 9
«
RAY 3
A A
V./ V_#
/. A
N»' 1
1 rs r\ r\ *\
o o o 3 3 o
C O #
RAY 5
A A /, A O
V»>^HUC 0 0 0 354* RAY 6
A A / .A A
\J KJ *-* \J O 020043- LD A 0# A M E S 7
y. r\ /, 1 A A
1o
/. A A
o o NEG «
0
A
%.<
00405 | A A A A /,
1 o o o o *4 COM 0 . 0 . SZR •
9 PLOTTING ?
00406 A A A /, A A
o o o *4 o O JMP N6
A A £. A O 1 «
o o O o ¿0 1 JSR ©.PLOT
00410 f A A A /.A A
o o o M V.» c RAY 6
004 1 1 0 2 0 0 3 7 - N 6 : LDA 0.AMES3 m
9 CONVOLUTION AND
0 0 4 1 2 |I Ao Ao Ao Ao C3 COM 0 . 0 . SNR 9 I N V E R S E CON VOL. CALL
004 1 3 0 0 0 4 3 3 JMP N7
004 14 A As>.i \DC A A O—
\j O d JSR ©.NEXTP
230
SIGAN
00415*006031 - JS R «•TXTT
00416*177777 CON V
00417*006022- JSR ««CON ; CONVOLUTION ROUTINE
00420*000410* RAYÓ
00421*000344* RAY1
00422*000355* RAY7
00423*020037- LDA 0# A M E S 3
00424*100400 NEG 0# 0
00425*100004 C0M A
kJ KJ0
0
A
J PLOTTING ?
rs a a o /. JMP
/ | O A * o **
N7
00427*006021- JSR «•PLOT
t»UU*lCd RAY7
/, « •
A A A r t Í A A
JSR «.NEXTP
uvy*4oc u y O u J I
n A y i O A ' A A / - n < » |
JSR « . TXTT
00433*177777 INC0N
00434*006023- JSR «•ICON *•
I N V E R S E CON. ROUTINE
00435*000420* RAY6
00436*000345* RAY2
00437*000371* RAY8
00440*020037- LDA 0# AM ES 3
00441*100400 NEG A
k.' kj
0
A
00462*100004 C0M A
O #
A C 7 D
O L.Í1
00463*000407 JMP N8
00464*006032- JSR «•NEXTP
00465*006021- JSR «•PLOT ♦
0 PLOT MAGNITUDE
00466*000456* RAY7
00467*006032- JSR «•NEXTP
00470*006021- JSR «•PLOT £ PLOT PHASE
00471*000457* RAY8
0 0 4 7 2* 0 2 0 0 0 3 - N 8 S LDA 0 # TTIN•
•
RESET FOR TELETYPE
00473*040001- STA 0.GETC •
0 IlN P U T
00474*002401 JMP «• ♦ 1
A A /, T C
U V- » H f J
• A A A A A A
U V . * U U V . » V.»
t
STA R T
A A A A A A •
. EN D S T A R T
231
-produce parity
-perform pagination
-initialise teletype
TTYR:
-initialise high speed paper tape reader
HSPTR:
-initialise cartrifile mag. tape uniy
CARTF:
-initialise analogue to digital converter
A2DI:
0001 SERVI
.TI TL S E R V I C E
. EMT T T I N . T T W R # T X T T . C R L F » N E X T P # R E S L T » Q U E S
• EMT TTYR# CARTF # A 2 DI » GTCC» G T C P » G T C T # H S P T R
. E X T D . T X T T » •N E X T P » • R E S L T # N L » N L P 1 »MS K
. EX TD G ETC# PUTC# C 5» C 5 P 1# C 12# C 15# SCM TR
. E X T D A M E S O # I C H M # S R A T E # .G T C T # .G T C P # . G T C C
•M R E L
0 0 0 0 0 * 0 5 4 4 4 7 TTIN* STA 3# RI 3 9 INPUT A CHARACTER
t • /iA /
r\ r > r \ r> a q a
LDA 3# M S K
A A n A O * A / O /
c uu JO |lu
A
SKPDM TTI 9 I D L E T. T Y P E R E A D E R
A A A A A • A A A 9 7 7
V.» V.» KJ O L* / f f JMP . - 1 9 G E T CHARAC. A M D
/. 9 £
u v.»u u M
r * r\ r\
v D v O lv a a t a
DI AC 0» TTI
0 0 0 0 5 * 163401 AMD 3» 0 » SKP 9 CLEAR READER
0 0 0 0 6 * 054441 TTWR: STA 3# RI 3 99 OUTPUT A CHARACTER
00007 * 040444 STA 0» SAVO 99 S A V E AC0
00010*176400 SUB 3» 3 •
9 Z E R O AC3
00011*101222 M0VZR A
KJ 9
A
9 OC 6
%J
7 P
.tr 99 M O V E 1ST B I T INTO CAR R Y
00012*175400 IMC 3# 3 99 SET AC3 TO 1
00013*101024 M0VZ 0# 0# S Z R 99 ALL BITS CHECKED ?
u v -»v» 1 h rD 99
A A A | /. « A A A ^ ' j r
K J K jv .fi JMP •— 3 M0# GO B ACK
U U U l J ucv.»4*j O
A A A 1 e t A A A /,Q/-
LDA 0# SAVO 99 YES# L O A D THE M0.
A A A | r • 1 A fi A A A
V V V 1 U X X v* \J sJ KJ M0VS A A
C» # C* 99 S W A P IT TO UPP E R HALF
A A A 1 A f | A | 1 A A
KJKJK.>1 i 1V 1 1VU MO VL 0 »0 99 R E M O V E P A R I T Y BIT
VOUCU X
A A A A A f «*7fC.JCUU
O A A
M0VR 3# 3 99
M O V E P A R I T Y INTO CAR R Y
A A A A 1 « 1A | A A A 9
L'V» C lUl OUU1 M0VR *_»#
r\ a
KJ 9 A D D P A R I T Y TO WORD
A A A A A f J A |A A A 9
u u e e lui Ouv MO VS A
KJPA KJ 9 S WAP TO L O W E R HALF
00023*06351 1 SKPBZ TT0 99 I DLE T • T Y P E P R I N T E R
A A A O / i tA A A 7 7 7
fwf k j k j k j JMP .-1 ; G E T CHARAC. A M D
00025*06111i DO AS 0# TT0 •
9 S T A R T T Y P I N G IT
00026*034013$ LDA 3# C 12
00027*116415 SUB# 0# 3# S M R
v.»V»U O V» U Xuuun3>
A A A Q A f A l A A A / ,< £
ISZ ML 99 K E E P T R A C K OF
00031*002416 JM P @R13 99 N U M B E R OF L I M E S
A A A Q A f
U U U O u
A A / . A A C tt
U U ^ t U U J j ) LDA 3.MLPI 9 PRINTED
00033*054004S STA 3» ML
00034*034413 LD A 3# RI 3 9 SAVE REENTRANT
9
A A A C A f A A A A A A D 1 /■ f A
V.» U V^' J U V.» V» V» v.» U v.* A X H é V.»
A A A C | • A A A A A A D 1 C • A
V.» O I VJ %-* KJ v_» n X u i U
A A A C
v.*uu J o
O ® A A A A A A
u u u v.» u v.» n
D 1 7
I r
•
•
A
K.f
234
A A O O
U U U G SERVI
A A A C O A A A r t A A C A J J A , A
L*W Ut.'UV/ J H W u Z KJ
O A A C /
0 0 C C C 0 SAV2Z
i A
Ci
A ^^c:c
uu u DO 0 5 4 7 7 4 TXTT* STA 3,R15 9 TYPE A TEXT STRING
A A A C^
u u u D O 010773 ISZ Rl 5
A O A C 7
v » o u D / 031400 LDA 2,0,3
A A A £ A A A I A A A
u u u u u y c l C U U LDA 0,0,2 «
9 L O A D FIRST WORD
A A A il 1
o u u u I 034006S LDA 3,MSK «
» LOAD MASK
a a a / q
u u ^ ' O g 163400 AND 3,0 »
9 MASK LEFT BYTE
A A A £ A
u u u O o 004723 JSR TTWR «
9 P R I N T IT
A A A / /. A A 1A A A
v » u u 0*4 W C l o w LDA u, u* c 9 L O A D WORD
A A A ^C 1 A
1 11 Juw
A A A
M0VS 0,0 •
9 S W A P BYTEiS
A A A / /
u u u O O 034006$ LDA 3,MSK
A A A / 7
v.' U \.»0 f 1 63405 AND 3 , 0 , SNR i MASK RIGHT BYTE
A A A 7 A
v ^ v ^ v . * / V.* 00 2 7 6 1 JMP @R15
A A A 7
U U U f
1
1 004715 JSR TTWR «&
9 P R I N T IT
A A A 7 A
U U U /G 151400 INC 2,2 2 INCREMENT POINTER
A A A 7 O A A A 7 r r
v .»u v #O 4 O j JMP .-13
A A A 7 /|
u u v .» / *4 0 5 4 4 0 6 CRLF* STA 3, R3 2 CARRIAGE RETURN
A A A 7 C
u u u /D 020014$ LDA 0, C 15
A A A 7 / A A Z A A 1
u u v .» f O U U U O l u 3 JSR GPUTC
A A A 7 7
SJXJK J t i 020013$ LDA 0, Cl 2 m
0 LINE FEED
A A 1 A A A A £ A 1 A «
».»U 1 \.*O UV/ D U l u l JSR GPUTC
A A 1A «
UC# i V.» 1 00 2 4 0 1 JM P §R3
A A 1 A A A A A A A A A A . A
V U 1u c u u u u u u r \ j ! U
A A |A A •
\JU l u o 0 5 4 4 1 0 NEXTPJ STA 3,R2 , GO TO N E X T P A G E
A A | A i* A A A A 14
U J) LDA 0, C12
A A I A C
v. «c* 1 u D006010$ JSR 0PUTC
A A | A /“
u u 1 v/O 024004$ LDA 1,NL J> FILL CURRENT
A A 1A 7
u u l u f 1 24014 COM# 1 , 1 , S Z R ,• P A G E WITH
u u 1 l u 000774 JMP • -4 ; BLANK LINES
00111 0 0 4 7 6 3 JSR CRLF
0 0 1 1 2 0 02401 JMP 0R2
0 0 1 1 3 Au Au uA uA uA uA IDt Oc ** A
U
0004 SERVI
00241*041400 STA A
9
A
KJ 9
'ì
ó
•
» S T O R E IN B U F F E R
00242*010412 ISZ DABP «
9 INC* D A T A B Y T E PT.
00243*010410 ISZ NOBYT •
0 INC. NO. OF B Y T E S
00244*000761 JMP FI
0 0 2 4 5 * 0 3 4 4 0 5 OUTs LDA 3# B B U F A 9 R E S E T D A T A B Y T E PT.
00246*175120 M0VZL 3*3
00247*054405 STA 3 . DABP •
9 S T O R E IN DATA
U V J G J U v j v ^ v O LDA 2 *R6 9 BYTE POINTER
00251*002404 JMP ©R5
0 0 2 5 2 * 1 7 7 7 7 7 BBUFA5 RAY3 i P O I N T E R IO S T O R A G E
0 0 2 5 3 * 0 0 0 0 0 0 N OBYT: V
0 0 2 5 4 * 0 0 0 0 0 0 DABP: A
U
A A A C C I n C * A
U U C J J V.» O
v A A A A a a n/r * A
UV.*C J U n o « V»
0 0 2 5 7 * 0 2 4 0 2 1 $TTYR: LDA l* • G T C T * T E L E T Y P E I N I T I A L I S A T I O N
A A A / A 9 / > O A / | A O
JMP • +2
0 0 2 6 1 *024022SHSPTR: LDA 1* • G T C P *• H I G H S P E E D P A P E R
00262*044007S STA 1 * GETC ; T A P E R E A D E R INITIAL.
O l A A « /lA A
u o g O j U U 1 JMP 0*3
0 0 2 6 4 * 0 5 4 7 0 2 CARTF: STA 3* R16 * C A R T R I F I L E INITIAL.
0 0 2 6 5 *0 2 4 0 2 3 S LDA 1* • G T C C
00266* 044007$ STA 1 # GETC
00267*006406 JSR ©GUNIT * S E T T A P E NO.
A A 0 7
V V G f
A
U
V A A i l /i A £
LA.» JS R ©GFILE «
S E T F I L E NO.
9
0005 SERVI
00 334•1 P O O 00 FDFC A
V.» SAMPLES
00335*040410 FSTA 0* TEM
00336*074407 FFIX TEM
O O #
\JVJ
| o A A A A
f
l u u u u u FEXT
A A O / . A fA / | A £
LDA 0* T E M + 1
A A
00 341 * 10 1 4 0 0 INC L»# U
00342*040015$ STA 0# S C N T R S T O R E IN S A M P L E
00343*002401 JMP 0R7 COUNTER
A A Q /. /.
U U U U
t A A A A A A
U U U O R7i A
U
A A
U
A /|C fAU
U ^J
A A A A A
U U U U TEMJ A
KJ
A A A /| / » A A A A A A A
u u O h u \ .• c » V.» u u U
A A A /|7
^ ^ 0 ^ 4 /
fv.»
A A
u
A A A A
c C20: 20 16 C H A N N E L S
• END
Enter
( Subroutine^
I
Load plot pointer
I
Load first No.
Increment pointer
Decrement counter
Set range
to 6
Multiply by
10
Cal.1Inti]«| Ml / .F.|
CaliINT2| = |MAX/S.F.+.95|
©
F igure A5.5 PLOT- A u t o r a n g ing /a u t o s c a l i n g plotter
239
COCI PLOT
. TITL PLOT
• ENT P L O T
. EXTN Y A X I S , S P C E 3 , F E N T
•N R E L
r\
y
a
y
a
K.' y
a
V.» 0 5 4 5 3 3 PLOT: STA 3# R7 ; PLOTTING ROUTINE
A A A A «
V.» 1 0340125 LDA 3 , WSA ; CHANGE FORMAT
A A A A A A A A t ff
V.» O y <z y y y 1 3> LD A o , w ; SPECI F I F A T I O N
A A A A O
V.» v.' v.» v.» O 0 41521 STA 0, 121,3
A A A A /»
M
A A A A
y <Zy y
A A
y <Z
C
T LD A 0, D
A A A A C
U U L» U J 041522 STA 0,122,3
A A A A /*
u u u u O
A A A A
y cZy y
O'"* C
c ,y J LDA o ,m i ; SE T P L O T T I N G C O U N T E R
A A A A H
V.» V,' / 040525 STA 0, C2í ;
A A A | A
V.» U V.» 1 U 032523 LDA 2 , @R7 5 P O I N T E R TO P L O T
A A A
y y y
1 1
11 177777 FENT SOURCE 9
A A A 1 A A Q | A A A
vy 1 Gi u G 1y y y FLDA r\
y 9
c\
y 9
rs
¿
A A A 1 A 1 A C A A A
y O y 1O i y y y y y FM0V 0, 1
i
r> r> r\
y y.» y 1
/>
H 1 11000 FM0 V 0,2
i c
r\ r^ r\
y uv.» 1 J 0 2 1 0 0 0 PLl: FLDA A
y 9
A
y 9
A
€Z
A A A
y~* y y
\
1O
r
122416 FSUB# i ,o , f s l e ; F I N D M A X I M U M
A A A 1 7
y \J y X /
1 A
1y
C A A A
y y y y FM0V 0, 1
A A A
y.' y.» y
O A
(Z y 112416 FSUB# 0 , 2 ,f s l e ; F I N D M I N I M U M
A A A A 1
L» V.» U G 1
l i t A A A
1 X X y y y» FM0V 0,2
A A A A A
y y y C .< Z
I A /, A A A
X y H y y y FIC2
f\ r\ r\ c\
u u v.» c O 014511 FDSZ C2
o a a a /i
000771 F JMP PLl
A A A A C
U V.' C ^ 044456 FSTA 1 , MAX ; S A V E MAX.
A A A
u u
A £
g D 050457 FSTA 2,MIN ; S A V E MIN.
A A A A T
y* y j y j CZ. i 146403 FSUB 2, i , f s n r ; MAX. > MIN. ?
A A A A A
\w' O 145000 FM0V 2,1 í NO
A A A
y y y
A
O
1
I 125003 F M0V i, i , f s n r ; IS R A N G E = 0 ?
A A A O
y y y O
A
CZ 024477 FLDA i, s i x ; YES, SET IT = TO 6
A A A A A
W OO a
y
r y r s r s c\ /» a
<z y y y *4 3) FLDA 0, O N E ; CALCULATE SCALE FACTOR
r* r\ r\ o
v_» y O
/*
Hr 030475 FLDA 2 , SIX
^ ^ A O C
v»> u v.» O J 1 3 2412 FSUB# i , 2,f s l t ; RANGE < 6 ?
v.»uOO
A A A A £ A A A /i 1 A
F JMP PL3 J YES
O / 034005SPL2:
A A A A A
w ^ FLDA 3 , TWO ; NO
A A A /, A
y» y 004412
y *-* y F JSR PL4 ; S. F. *2
000/11 0 3 4 4 6 6 FLDA 3, T W O 5 ; S. F . * 2* 5
00042 004410 F JSR PL4
A
y
A A y, A
y 034005$
y *H O FLDA 3 , TWO ; S* F • *2
00044 004406 F JS R PL4
000/15 A A yA T T/ G A
y y i F JMP PL 2 ; RE T U R N
0 0 0 4 6 0 3 4 0 1 7 $PL3: FLDA 3, TEN J L O A D TEN
O
A A A /» *7
u v v y M / 1 O y ¿Z y y F DI V 0 y
A
9 9• S.F./10
A A A C a
O o 164100 FMPY 3,1 i R A N G E * 10
/A A\
y y»y
C
y X
I
000764 F JMP PL2-2 ; RETURN
A A A
V.» U L» J G
c O
160100 PL4: FMPY 3,0 J S . F . * AC3
\w»
A A A C Q
y y o O 164200 F DI V 3, i ; RANGE/AC3
a a
v.» v.'D
a c y,
*-* 132415 FSUB# i,2,f s g e ; RANGE = < 6 ?
A A A C C
v_* y y
A A
y y
1 /| A A
X *4 y y F JMP y90A
A
9• NO, RE T U R N
000 56 103400 F RND A
y
A
9 y 9 »
YES, R O U N D S.F.
240
0002 PLOT
•
00057*040430 FSTA 0# SF 9 S T O R E IN S. F.
•
00060*034425 FLDA 3,MIN # L O A D MIN.
00061*114200 FDIV 0# 3 •
# D I V I D E IT BY S. F.
•
00062*054427 FSTA 3 # INT1 # S T O R E IN I N T E G E R 1
•
00063*074426 FFIX INT1 # T R U N C A T E INT. 1
00064*060425 FFL0 INT1
00065*034416 FLDA 3# M A X L O A D MAX.
00066*114200 FDI V 0,3 # D I V I D E IT BY S. F.
•
00067*030436 FLDA 2» P95 L O A D .95
A A r t 7 A | 1 C 7 A A A •
U U U 1U 1 3 / UOV.» FADD 2# 3 # A D D .95 -T O IT
00071*054422 FSTA 3# INT2 î S T O R E IN I N T E G E R 2
00072*074421 FFIX INT2 J T R U N C A T E INT. 2
00073*060420 F FLO INT2
00074*030415 FLDA 2 » INTl 5 L O A D MIN.
•
00075*034416 FLDA 3 # INT2 * L O A D MAX.
00076*156400 F SUB 2,3 # MAX.-MIN.
00077*030432 FLDA 2# SIX * L O A D SIX
00100*156411 FSUB# 2,3,FSGT; M A X . - M I N . =< 6 ?
I A 1 f /. /. A •
A A
V.» 1U 1
A A A
OV.» FJMP PL 5 # YES, GO ON
I A A 1 A A « A A •
A A
u u l u e U U 1 H v.»U FJMP 0,3 # NO, R ETURN
« A •
A A
U
A
1 U O
A A A A A
V.» V.'
A
U MAX: A
A A | A /, f A A A A A A A
V.» 1^ M U U V.» U U L» SJ
A A
U U
1 A
lu J
C f A A A A A A
V.» L* S J \ J \ J ^ MIN: A
u
A A |A ¿ |A A A A A A A
L» 1 V.' O V* L»
8A A A A A A
A A
u u
1 A 7
1u / u u u u u u SF: A
KJ
A A 1 A
u u 1 1 u u u u u u KJ
A A 1 « | VA
^ » ^ 1 1 1
A A A A A
V.» V.» U V.» u INTI: A
\y
A A 1 |A « A A A A A A A
U U 1 i d. \ JK JK A \ J K J y
A A I | n | A A A A A A
<.•\ J 1 1 O \J \ J vy \ j INT2: A
V.»
A A 1 1 /. I A A A A A A A
UU 1 1 H \y V.» V.» U KJ V/
A A « i A*
U U 1 l J \ J KJ U U U ISTRT: U
A A f | / f A A A A A A A
U U l 1 U V.» V.» U V.» U V.1
7 V A A A A A A
A A
uu 1 1 #
1 i
uu u uu u D DA T A : A
U
A A 1 A A 9 A A A A A A A
U U 1 C U U U U V.' U V V»'
00121*037424 r \j y O •
A
KJ *
A A C
KJ KJ Zj
00122*075341
1 0 7 « A /.A
A A
«y ty 1 y u *4 u
A A A
P5: 0- 5
A A | A 7, | A A A A A A
V.' 1 O M *y *y y \y vy vy
0 0 1 2 5 * 0 4 0 3 6 3 P95: 0 .95
00126*031463
0 0 1 2 7 * 0 4 0 4 5 0 TW05 : 2. 5
/A A | A A f A A A A A A
V.' V.' 1 O \y W U U U U
A A « 7 /, I A A A A A A A
V»' 1 O H U U U V.' U v.' C2: U
A A 1 7 c « A A A A /. A
V.» u 1 J J u u u u Hu C40 : 40
A A t 7 ¿ f A A A A CO
U V.' 1 O Ü V.» V.' U V.» J o C52: 52
A A 1 77 • A A A A C 7
UU 1 O í V.» U U u J O C53: 53
| y, • r\ r\ r\ r\ r c
a a
V.' \y 1 M V-'
a
\y V.' \y *y y y C55: 55
00141 *0 2 0 7 5 0 PL5: FLD A 0, INTl
00142*0300175 FLDA 2, TEN
0 0 1 4 3 ’ 140100 FMPY 2,0
A A
00144*103400 FRND xj 9 vy
*
00145*040750 F S T A 0 , ISTRT S E T L E F T M O S T Y
•
00146*074747 F F IX I STRT J A X I S VALU E
00147*060746 F FLO I STRT
r > A I c r> 1 <r > n n n / i
v.»v.i 1 J \ ' 1 F EXT
00151*006006? JS R 0 . CRLF
241
CC O 3 PLOT
00152 A A O A v
r
u m o ü j JSR e.CRLF
0 0 1 5 3 126120 ADCZL 1» 1 9 G E N E R A T E -2
0 0 1 5 4 127120 ADDZL 1* 1
«
9 -2*4
0 0 1 5 5 A A A A I1 1| * KJKJKJKJ FENT
00156 C24733 FL DA WINT1
0 0 1 5 7 AU GO U 7# 7O AO a
FLDA 0# SF
/>/> 1
i Ou 104100 FMPY 0# 1
00161 0 3 0 7 4 0 P L 6 : FLDA 2# P C C 5 R O U N D OFF THE
•
P
0 0 1 6 6 146400 FSUB 2, 1 •
SCALE P
001 6 7 |1 A A A A A FEXT
U U i / v/ Ouu / 3)
A A £ A A »7
\J kj «JSR 3.TXTT
00171 177777 SPCE3
0 0 1 7 2 125405 INC 1»1 SNR 9
0 0 1 7 3 Au A A Oy , A £ JMP •♦ 6
0 0 1 7 4 AU AU AU 1| JC DC f FENT
00175 020712 FLDA 0» SF
0 0 1 7 6 1V 1uuu
f A 7 A A A
FADD 0» 1
001 7 7 r\r\ /Osr\ K »K J K J g F JMP PL6
A A A A A 1 A A A A A
U U Ü U U 1U U V . » V . » U FEXT
A A ¿A A £
uu cu 1 u u O u u ü J>
a a 1
JS R 3.CRLF
A A £ A A O <T « PRINT Y AXIS
O U OU U / «T JSR § • T XTT
A A O A O
UL» G L » C 9
UV.»CU O 177777
A A A A Q
Y AXIS
A A A A /
V V.' H. 0 3 2 7 2 7
KJ LDA 2» 9R7 •
P O I N T E R TO P L O T
rsr\ryr\c AU GA UA UA GA UA <T3>
K,* KJ U LDA 0* M 1 • SOURCE
P
A A o o o
o \^' 050675
a, cz C, FSTA 2 b DDATA
A A O O O
\J ^ 074674
CZ. CZ,sJ FFIX DDATA 9 T R U N C A T E IT
00224 060673 F FL0 DDATA
A £ 7A
0 0 2 2 5 Av.»A \_» O fu FLDA 0 # ISTRT •
* C A L C U L A T E SPA C E S TO
0 0 2 2 6 030671 FLDA 2, DDATA P PLOTTING POINT
0 0 2 2 7 1 12400 FSUB 0* 2
A A A Q A
V_» V_» O V.»
Cj 050667 FSTA 2 # DDATA
00231 0 7 4 6 6 6 FFIX DDA T A
00 232 1 v.*uuuu | A A A A A
FEXT
0 0 2 3 3 145000 MOV 2b 1
0 0 2 3 4 A AV.' O£ AV.» AU 7í ^*T JSR @. TXTT
9
P LEAVE 3 SPACES
0 0 2 3 5 000171 * SPCE3
00236 I O lv/uu
| O 1 A A A
MO V Ib 2
00237 A A A 7 A A
V.» í V* v_» LDA 0 # C 53
00 240 0 1 0 0 1 1 $ I SZ C9P1
0024 1 020677 LDA 0, C55 TYPE EITHER A
9
00 24 3 0 2 4 6 5 5 LDA 1# D D A T A + 1
0 0 2 4 4 125112 M0V L # 1 # Íb SZC
242
0004 PLOT
00245*126460 SUBC 1» 1
00246*124000 COM 1* 1
00247*020666 LDA 0# C40
00250*125405 INC 1 * 1*S N R
1 t A A A / . O Q
V.» U C J I V.* %J \J *4 \J O JM P • +3
00264*101400 INC A
\J »
A
KJ
00274*0200155 LD A 0# DD
00275*041522 STA C m 122# 3
7
k.
0 7
» kjcz* 1
A
O • A A
o u Or uA wA / Ch JSR 9 m TXTT
•
9 P R I N T Y A XIS
A A 0 7 7 • A A A O A
V ' V.« O / / U U U C U
7
J
•
YAXIS
A A O A A • A O /, / A A
v.'u u v j O ^ u O O LDA 3# R7
00301*001401 JMP 1» 3 •
9 RE T U R N
. END
243
-* Get
Store data
Noy
• TI TL A T O P C O N V E R T E R
• E X T D I CHN# SR A TE# • TX T T# N # M # N U M # N U M P 1
• E X T D SCNTR# TEN
•NREL
0 0 0 0 0 * 0 5 4 5 0 2 A 2 DC: STA 3# SAV3
A A A A ( f J A A Ç A A
SUBZL
y y X
y k
y y G y G y VJ# VJ
A A t
* SET WAIT C O U N T E R
A A A A A f <A / . A A A A
y y y y G. y y y y y STA 0# V C N T ; TO 1 (SYNC.)
A A P A A I 1 1
l 1
A | /.A
M0 VOL
y y y yl y y *4 y
A
KJ 9
A
KJ # SET C L O C K TO
•
A A A A / , « A
v.' v_> v j ¿0*4
A A A O .
v.» H u v.» vj O STA 0# BUFP+ 1
A A A A C
V.» V.» V.» Ci J
9 A /. A A A / i .
STA 0# BUFP+ 2
A A A A / ;
U v.» u u O
f f A A /, A A
1 u c M u v . ' SUB P
y »
P
y
C
A A A A H
VJ VJ V.» f
f A /. A A A A -
KJ KJ KJ KJ d STA 0# BUFP
J
A A A A A 9 A A A
V.» V.» V.» V.' JJ C T
LDA 0# TM SK
VJL*L*O 1
A A A A 1 9 A
VJ
£O fi
O A A A
V»
A
MSK0 «
V*
.>
' • M A S K OUT AL L INTER.
#
A A A A A
V» K.» K.1 sJ d
fV.'V.</;
A A A A A
O U U O Î JSR @ • TXTT ; E X C E P T A TO D
A A A A A 9 A A A A C A f
KJ KJ K.» \ J \D K.> KJ KJ d Z J K J M ESS » A N D CLOCK
00034*060210 NIOC tti ; C L E A R T* T Y P E R E A D E R
00035*063610 SKPDN TTI J IDLE T . T Y P E R E A D E R
v»v»v»o O
A A A /■ uuvj f f f
A 9 A A A A A A
JMP .-1
V»v'VJO f «j O/u c l u
A A A A 7 A A A I A
NIOC tti ; C L E A R IT
00040*060177 IN TEN •
9 E N A B LE IN T E R R U P T S
00041 *0601 14 N I 0 S RTC ; S T A R T CLO C K
0 0 0 4 2 * 0 2 4 0 0 2 - GET: LD A 1 # BUI* P J L O A D NO. OF WORDS
245
0002 ATO DC
A /i A
U 1 25005 MOV 1 # 1 * S N R »• LAST ONE ?
00044 000442 JM P EMPTY •
# YES# B U F F E R E M P T Y
00045 014002- DSZ BUFP » NO# DEC* NO. 0F
0 0 0 4 6 A A A /• A t1 KJKJ KJ * A v J JMP • -**1 * W O R D S IN B U F F E R
A A A /l *7 A 7 A A A 7 __
' J U U 4* f v y O vy v y vy o LDA 2*BUFP+1
A A A C A A A | A A A
U U Ji.» u c 1 vy vj vj LDA 0# 0# 2 #
•
L OAD WORD
vy vy v y D 1 155400 INC 2 »3 # INC. A D D R E S S
A A A C A
vy J c 024006- LDA 1# B U F P + 4
a a
u u u
a c
J O
o
136015 A DC# 1*3# S N R #« AT MAX. ADDRESS ?
A A A
u u u
c /i
J*4 034005- LDA 3* BUFP«*3
A A A C C
O U M D J 054003- STA 3# B U F P * 1 # NEXT WORD ADDRESS
A A A
v » u u
C £
J O 034425 LD A 3 # INAD
7
a
u
a a c
J 9 151222 M0VZR 2 * 2 # S Z C #• EVEN A D D R E S S ?
A A A r A
034424 LD A 3* 0 U T A D #* NO# O U T P U T
A A A / 1
u u u O 1 101113 M0VL# 0*0# SNC j YES# INPUT. N E G A T I V E ?
A A A /! A
v» u u D c 126461 SUBC 1*1* SKP ; NO* SET UP M O S T
A A A £ A
o j j D O 1c O u U KJ A DC 1# 1 #
•
S I G N I F I A N T 16 B I T S
00064 045400 STA 1*0*3 ; STORE DOUBLE
A A A £ C
v .* u u O 041401 j STA 0*1*3 # PRECISION NUMBER
A A A £ £
VJ VJ VJ O O 175400 INC 3*3 •
# NEXT ADDRESS
A A A £ 7
VJ VJ VJ O f 175400 INC 3*3
A A A 7 A | A f A A O
V j V j VJ 9 VJ 1 v y 1 k j \j c MOV 0 * 0 * S Z C #• EVEN ADDRESS ?
A A A 7 1 A A A /| A Q •
U U U f I v y u O H v y o JMP •♦ 3 # NO* O U T P U T
A A A 7 O
UJV.» f c 054 4 1 1 STA 3# I N A D •
# YES# I N P U T
A A A 7 O
vy vy vy i O
r t A A / ,/ > Q
JMP • * 2
A A A 7 /|
v j v j vy i *4 054410 STA 3# O U T A D
v.»v.»u f 3 000745 JMP GET
A A A 7 £
U ^ » U / O 177777 PRAYlS RAY1 # P O I N T E R S TO I N P U T A R R A Y S
A A A 7 7
vyvyvy r f 177777 PRAY2: RAY2
A A 1A A A A A A A A T lT M t A
o u 1 vy u vyvyvyvyvyvy I IT J J vy
A A « A 1 A A A A A A A
v y v y X vy X KJ KJ KJ S J KJ KJ vy
A A 1 A A A
JV.» i u c 0 0 0 0 0 0 SAV3: vy
1A
A A
vy vy 1u O
O
0 0 0 0 0 0 INAD: A
vy
A
A A | A / .
I vy »4 0 0 0 0 0 0 OUTADl vy
• A
A A
vy vy 1 v»
C
J 1 7 7 5 7 3 TMSK* 177573
•
A A 1 A £
vy vy i vy D 0 2 0 0 0 1-EMPTY: LD A 0# SCN T # SAMPLE COUNTER
«
A A t A 7
vy vy 1 vy 9 101004 MOV 0 * 0 # SZR # A L L S A M P L E S TAKEN
A A 1 I A •
vyvy X 1 v» A A A 7 7 0
vy v/ vy 9 J c JMP GE T # NO
•
00111 vy *4 Avy vy A
A A £ <T
vy O o> STA 0# N U M # YES
•
00112 0 3 0 7 6 4 LDA 2# P R A Y 1 # P O I N T E R TO X
A A A A 1 A C
0 0 1 1 3 v y c v y v y I v i 3) LDA 0# S C N T R
0 0 1 1 4 11 Avy Av y/*4. Avy Av y NEG A A
v y # vy
0 0 167 r\r\r\
a /: a
K j O K j < Z s K j KJ Nf'QC A
KJ
A A
\ J kj
1 7
X 9
A
kj
A A A /. O
JMP DISM+1
00171 0 3 0 4 0 7 DISMS LD A 2* SAV2 •
0 DISMISS ROUTINE
00172 024405 LDA 1* SAVI •
0 RESTORE A C S •
00177 A A A A A A C A TI 1 •
O H v l •
A
KJ
A A A A A A C A IIO * A
V.» V» C V.» U C» U V.» U V.» OH V C # V.*
A A A A
e; e» o e' X
|
0 6 0 2 0 0 CLERI NIOC A
KJ
A A
v.»u cu c
A A A
0 62721 ADCTl DI CP 0* ADC V
A A
u u cu O
A A Q
060121 NI OS A DCV
A A
e» V_'
A A /.
e» *4 0 2 4 0 0 2 - ADI s LDA 1 * BUFP 0 NO. OF WORDS IN BUF .
A A A A
e' V.' C U J
C A Q A A /> O _
e» 3 e' e» k j 9 LDA 2* B U F P + 5 ; MAX. NO OF WORDS
/-xa r^r\ s
v .»U Cu U 132415 SUB# 1*2* SNR 0« S I Z E OF E U F F E R ?
•
A A
V-'
A A 7
O V.» 9 063077 HALT 0 B U F F E R FULL* F A T A L
A A
UUG 1 U
A 1 A A 1 A A A O—
e» X <.> k j k j & I SZ BUFP
0 0 2 1 1
A O A A A / i ^
e» O e» e' e# *4 LDA 2* BU F P + 2 ; N E X T AVAIL. A D D R E S S
0 0 2 1 2
A /i 1 A A A
uH X u u o STA 0*0*2 0 STORE WORD
0 0 2 1 3 141400 INC 2*0 0 INC. A D D R E S S
00214 024006- LDA 1* B U F P + 4
00215 122015 ADC# 1*0* SNR 0 A T M A X . A D D R E S S ?
00216 A A A A A
e* o e» e' e* 3
C .
LD A 0* BUFP+ 3 0 Y E S * G ET M I N . A D D R E S S
00217 STA 0* B U F P + 2 0 S T O R E N E X T A V A I L . A D D .
A A
v.' e'
QO
CPG, e»
A 1 A A
e» 1 e»
1 A , .
X e' I SZ CHNF 0 I N P U T OR O U T P U T ?
•
00221 000406 JMP N EXT 0 INPUT
00 2 2 2 0 2 4 0 0 1$ LDA 1 * I Ch'N 0 O U T P U T
00 2 2 3 065021 D0A 1* A DCV 0 S E T A TO D M U L T I P L E X E R
00 224 024407 LDA 1 * CM 2 0 TO I N P U T C H A N N E L
00 2 2 5 044010- STA 1 * CHNF 0 RESE.T FLAG
00226 000743 JMP DI SM
00227 0 6 3 621 NFiXTl SKPDN ADC V 0 GET N EX 1 A TO L
0 0 2 30 000777 JMP . - 1 0 CON VERSI ON
247
0004 A TO r e
00 231 062621 * DI CC 0» ADCV
00232*000752 JM P ADI
00233* 177776 CM2: -2
0 0 2 3 4 * 0 6 0 1 1 4 TIME: NI0S RTC ♦
9TIMING ROUTINE
A A
V.1u c
R
O J
^
U
1 /i^
1 t v J u u
^
DSZ WCNT i DEC. WAIT C O U N T E R
A A 0 7 4
v.i v.i ¿o O D
t a
v.»
a a h
v.» f
o /:
O O JMP DISM+3 •
9N O T R E A D Y TO START
00237*060121 NI0S ADCV •
9S T A R T A TO D C O N V E R T E R
A A O/l
U V.»
A • A A A A A Aff
Ov.'v.'v.'cl) LDA 0 » SRATE •
9RESET WAIT COUNTER
A A A /| 1 f A /. A A A A _
v > C . *4 1 H V.» V.» V.» STA 0# W CNT
00242*014001- DSZ SCtNT i SAMPLING FINISHED ?
00243*000730 JMP DI S M + 2 •
9NO
00 2 4 4 * 0 60 221 NIOC ADCV 9YES» C L E A R A TO D
00245*060214 NIOC RTC 9CLEAR CLOCK
00246*020730 LDA 0# SAV0
A A O /, 7
f
• A A
V.»
A A A
CL K.9
A
JMP *»0 0 R ETURN
9
00251*052012
00252*020117
00253*052123
00254*051101
00255*026124
00256*051440
00257*051124
00260*045511
00261*020105
0 0 2 6 2 * 0 4 7 1 0 1 ANY KEY< 1 5>< l2> < 12>*'
00263*020131
00264*042513
00265*006531
00266*005012
/A O *7 t A A A A A A
V.' V.' o O / u U S.» V.» V.»
END
248
F igure A 5 ,7 DATAIN- Ge n e r a l d a t a e n t r y r o u t i n e
249
0001 DATAI
• TI TL DATAIN
• ENT DATAI N
• E X T D . C RL F * • T X T T , N # M , N U M , M U M P 1, TEMPO
•EXTD TEMPI,C15,WSA,AMES1,TEN
• EXTN F EM T, M E A N
-MREL
A A A
U U V »
1 7
1 / 000417 JMP D4 p YES, L A S T D ATA
A A A A A
V.1V.» V»* 035001 D32 LDA 3, 1,2 •
NO, IS IT A
p
A A A
U U U
O
G
1
1 175004 MOV 3, 3, SZR p• B R E A K CHARAC- ?
A A A O O
K.» KJ K J C Z CZ
A A A /. A A
O U U H u O JMP • +3 • YES, S T O R E NO-
p
a a
V.» L*
a
U
a
G O
Q A A A A
o o o o
1 O
X g
•
F EMT
r% r\
v.' v_»
o /i
O *4
A A A A ZI *7
KJ K J KJ f O f F JMP D2 p NO, GET N E X T CHARAC-
A A A Oc;
O U G J
A A A A Q A
KJ K J K.» \ J C Z O
•
F EMT
A A A O A
v. î u u g D 064010$ FLD3 TEMPI
A A A
K.’ K.* KJ <Z
0 7
f 045400 FSTA 1,0,3 p S T O R E IN ARR A Y
A A A
u u v/ O
A A
X
«
X 1 A A A A
O O O O F IC 3 •
p INC. P O I N T E R
A A A A
U V
.*G» O
1
1
f
X
A 7 A A
O f KJ KJ
A
O
ff
^ FST3 TEMPI
A A A A A
KJ KJ KJ O ( Z X
1 A A A A A
K.I KJ KJ KJ KJ F EXT
A A A O T
KJ KJ KJ O 024003$ LD A 1 M , p I N C R E M E N T DATA
A A A A /,
O v.' v j O H 125400 IMC 1, 1 •
p COUNTER
A A A
«J U U
A C
U J v u v » / O I JMP DI
A A A A /:
o Kj o O 0 0 0 0 2 5 ’D4: FEMT
A A A A 7
V.# V.» V»* O /
A Z A A A C Ç
o O o o o D FFL0 NUM p C A L C U L A T E NO- OF
A A A /.A
KJ KJ K j KJ 030005$ FLDA 2, NUM •
p CORRELATIONS REQUIRED
v.»*41 0 2 4 0 1 4 $
A A A /l
c/
1
FL DA 1 TEM
,
A A A
V.» V.» *4 X /iA
tZ
1 A A A A A
J u g u v .» F DI V 1,2 •
p N/1C
A A A /i A
v.» U H J 050410 FSTA 2, TEM
00044 074407 FF IX TEM
000 4 5 X 1 A A A A A
V,' G» V.' G* FEXT
00046 024406 LDA 1 , TEM+ 1
STA 1, M • S T O R E IN M
A A A
KJ K^**-*# 0 4 4 0 0 4 $
/i 7 p
V.'UV.»Jv.» 0 3 4 4 0 2
A A A C A
LDA 3, R5
JMP 1,3 • RETURN
V» Jc»J 1 001401
A A A
V
C 1
p
0 0 0 5 P v.» u u nj•
•
A A A A A A
V.
» v.*
D C
V.»
A
000 5 3 0 0 0 0 0 0 TEM: A
K.»
A A A A A A A
000 54 V.KH»U G' K.9
*
A A A C C
V
.V.» -J 000 1 15 Cl 152 1 15 p ASCII M
F igure A5.8 LIST(1)- No r m a l i s e da t a
251
0001 LIST
• T ITL L I S T
• ENT LI S T 1 # L I S T 2
. EXTD .CRLF#.TXTT#N#NUM*TEMPO#TEMPI#AMES1
• E X T D M# M 1* • R E S L T # •PLOT# C5* C 5 P 1
• EXTN F E N T # M E A N
•N R E L
r\ r\ r\ r\
y y y y 0 5 4 5 0 2 LISTI: STA 3# R5 ; CALCULATE MEAN
A A O A
y y y y
1
1
A A A A
y CZy y y O x
A A (T
LDA 0# N ; AND NORMALISE
A A A A A /| STA
y y y y CZ
A
U ^
A A A c ff
U U U j i 0 * TEMPO
%.»
A A A A
y y O 177777
Q
FENT
/»
A A A A
U V.» D» M 102400 FSUB 0#0
A A A
y y y j y Zj
A C
066475 FLD3 3R5 ; P O I N T E R TO A R R A Y
a r
v*'v.»
a a a
O 0 2 5 4 0 0 Dl: F LDA 1 * 0*3 ; LOAD NUMBER AND
A A A A 7 1
y uuu / 1 cOv»uu FADD 1 *0 ; A D D TO SUM
L'V.' 1 y 1I I y y y y
A A A 1 A | A A A A
FIC3 ♦ INC. P O I N T E R
f
A A A « 014005 S
1
y.»yty* 1 1 FDSZ TEMPO ; DEC. C O U N T E R
1
A A A
W % J I CZ 000774
A
FJMP Di
1
y y y 1 O
A A A
024004?
O
FLDA 1*NUM
i /. 1
X C Z y C Z y y^
y y y 1 *4 FDI V 1*0 ; CALC. M E A N
A A A I r
L»U U ID l UUUUV.»
1 A A A A A
F EXT
r \ r \ s rs r \ i ^
y y y 1 O y y y y y 1 J JS R @.CRLF
vJ y
A A A
y 17/
1
y y \ D y y cZ x
c JSR @ •TXTT ; PRINT HEADING
A A A A A
CZy \ m n MEAN
A A A A
y y cz 1
1 A A A A A
y y y y y sj
A § F ENT
A A A
y y y <Z CZ
A A
140001 FFDCF ^
A • PRINT MEAN
»
v.#y j y
A A A A A
r. O «
1 A A A A A
u u o v / u FEXT
A A A
y yyCZ**
A/i A A £
y y y y y 1X
A A | ff
JSR @.CRLF
A A A
y y y c z D
A C
uCuuu
A
LDA 0# N
L* V.»O
A A A il
CZ O v.#
A
M
A
^• A A C
O
C
X STA 0* TEMPO
r\ 7
A A A
v y y cZ t 032453 LDA 2#@R5 ; P O I N T E R TO ARR A Y
r\ r \ r \ a a
y y y O A A A |f
y y y y c z 1
A A
FENT
A A A
y y y
A|
o I A A C A A A
y CZ y y y y
H O«
U CZ * FLDA 1 * 0# 2 * LOAD NUMBER
A A A A A
y y y % j CZ 106400 FSUB o# i ; MEAN SUBTRACTED
V.»U V'O O
A A A A A
M Oy y y A C O A A
F STA i * o #2 ; STORE NORMALISED
v-'v' O /. 1
A A A A A / . A A A
FIC2 • DATA IN A R R A Y
p
r \ />
L'V .» O O 0 1 4 0 0 5 S
o c
FDSZ TEMPO ; DEC. C O U N T E R
u v.' O O / y y y f i y F JMP
A A A A
V.»
A A A 7 7 7
D2
v.' O 7/ Iy y y y y
A A A A 1 A A A A A
FEXT
L» /. 0 3 4 4 4 2
A A A A
LDA 3* R5
0004 1 001401 JMP i#3 ; RET U R N
000 4P 0 5 4 4 4 0 LIST2J STA 3* R5 ; LISTING OR PLOTTING
00043 031400 LDA 2*0*3 * P O I N T E R TO ARRAY
00044 A A A A
c z y y y y
A A
X LDA 0# N
0 0 045 /i / STA
A A A A
*-* y y y
ff
O X 0* TEMPI
0 0 046 0 2 4 0 14S LD A i* C5 ; RESET RESULT FORMAT
V'y y M/l*/
A A A 7 044015$ STA 1* C5P1
^ v.»a c
A
O ^ 024007S LDA 1* A M E S 1
aaac 1
y y y O 1 125004 MOV l* l# s z r ; DATA LISTING ?
^
y yn yr \ ZD
c cz
o 000414 JMP D4 ; NO* SKIP OVER L I S T I N G
t~> O <
"
cQ
rs r\
y y ZD y
r\
v^i v .' y y J v ,'
>»Lnot
' JO • FENT
000 54 r D < »' v.< FLDA 1*0*2 J LOAD NUMBER
00055 144001 FFDCF l ; L I S T IT
•
1( r.>iHUL'V.1
y,n n n F I C2
00056
253
0002 LIST
/“% /*\ /> C *7 • 1n n n o r t
L.' L/ o D / * U U U U v . ' FEXT
A AA
v .» v . » u O u •050423 ST A 2 # SAV2 •
9 S A V E AC 2
A r i
o v-f O 1 •006012$ JSR © • RESLT
A ^ A / n
u u u D c •030421 LDA 2* S A V 2 #
•
R E S T O R E AC2
u u o U O •014006$ DSZ TEMPI
A A A / /, • <•><■> r> *
7¿*7
U U U U H u u u / O f JMP D3
A A A / C
U U O J *034415 LDA 3# R5
A A A / /
O O • 1 2 4 5 1 3 D4J NEGL# l#1#SNC • DATA PLOTTING ?
9
A A A / A
v / u ^ O # •001401 JM P 1,3 *
•
NO# E X I T
A A A H A • r t O A> *A t « % C
O U U f V.» u c u u It* « LDA 0#M ; YES# SET C O U N T E R
a a a *7t
V.*L»U / I •040011$ STA 0#M 1 ,
*
T O P L O T M OR N V A L U E S
A A A 7 A
yj\-> s** i •026410
o LDA 1 • @R5
A A A 7 A
U S J U # O •044402 STA I s +2
*
A A A 7 /,
U V J U / •006013$
*4
JSR ©•PLOT ,• P L O T THE A R R A Y
0 C 0 7 5 1A A A A A A A
U
7 /
A A A
v.»u u /O • 0 3 4 4 0 4 LDA 3# R5
A A A 77 *V/ IU J>
A A A A 1 A ff
LDA
K J W i i 0, M
A A « A A
u u i v u •040011$ STA 0.M1 s RESET M COUNTER
uu
A A 1 A
l u
1
•001401
1 JMP 1,3 ,• R E T U R N
A A 1 A A 1A A A A A A DC# A
u u l u c \J Kt Kf KJKJ %J> f \ D # U
A A
V-iU IU J w u w u w u jCAXIO#
| A Q ■ A A A A A A
H V c «
A
KJ
END
Figure A5.10 CORREL- General correlating routine
001 CORRL
• T ITL C O R R L
• ENT RR
• EXTN F ENT
• NREL
A A A A
*054461 RR: STA 3 * R9 9 CORRELATION ROUTINE
A A A f
V
_» Kj V» I *0240055 LDA 1*C5
•
9 R E S E T R E S U L T FORMAT
A A A A A
CL * 0 4 4 0 0 6 $ STA
K J KJ XJ
1* C 5 P 1
A A A A A
O * 0 3 5 4 0 2 •
v.* xj LDA 3* 2* 3 9 POINTER F OR STORAGE
A A A A /.
UUUO H *054010$ STA 3 TEMPI
m
A A A A C
s»' v.' v.' O * 1 5 2 4 0 1 S UB 2# 2 * S K P M J=0
/\ s s r \ r
\~J x.* \ j c' D * 1 5 1 4 0 0 RR1: m
INC 2 m 2 9 SET COUNTERS
A A A A *7 t A C A A A /i
W#I STA 2 m J
AAA 1A A A
v
.#V#u 1V.I * 102400 S UB X.9 9 XJ
AAA 1 |
UL»U 1 1 *040446 STA C#NJF
AAA I A
L»V'U 1 G * 0 3 4 0 0 2 $ LDA 3#N
•
\
»>o 1O • 1 5 6 4 0 0 SUP 2 3m 9 N-J
AAA 1 /• •
1M * 0 5 4 4 4 2 STA 3* N J
x j x j \j
9 STORE IN FLOAT. LIM I
AAA 1 r
c' x j i O * 0 5 4 4 4 3 STA 3*NJF+1
AAA 1£ •
tw
*v.* 1O *141120 M0VZL 2 0m 9 2 * J
AA A 17
S
-*\
^*X J i / * 026442 LDA 1# SR9
A A A AA t 1A 7 A A A •
V.'V/U CL L' i u / u u u ADD 0# 1 9 SET P O IN T E R
A O 1
^ U k .iL 1 *044007S STA 1 TEMPO
m 9 TO X C I + J )
AA A A A •A A A A A A e
V.' k.' U G C U LU U U O T LDA 0*M
AAA A Q
V.' V
.»V.' G O *034436 LDA 3#R9
A A A A/1
v.» ^ g: *4 •112415 SUB# C 2m SNR m 9 J=M ?
A A A A r
U U*U G J *001403 J MP 3 m 3 9 YES* E X IT
A 0/1
k JU U G U *031401 LDA 2 m I m 3 9 P O I N T E R TO D E L A Y E D
AAAAH
UUUG / * 1 7 7 77 7 FENT 9 S I G N A L LOADED HERE
A A A OA A /A
u u u j u *102400 F SUB XJ 9 XJ
A A A A 1
V.1V
-*kJ O i *064007$ FLD3 T E MP O 9 P O I N T E R TO R R C J )
A A A A A « A A 1A A A ODO•
v.*\
>»\
^*O u o i u u u n r \ cl • FLDA 2*0# 2 9 LOAD X C I )
A A A A A
V' vj OO *035400 FLDA 3*0*3 9 LOAD X C I + J )
a a a O/ i f 1A A kA A •
V.I \.i o M 1 /V
.« 1 G# F MP Y 3* 2 9 X C I)*X C I+J)
A A A A C
V U U*O J * 143000 FADD 2*0 9 A DD TO SUM
f 1A / .A A A
000 3 6 1 k.'*lVUkU FIC2 9 INC. POINTERS
9 1 1 AAAA
00037 l 1 UV.' k
^»V.' FIC3
A A A /, A
V V H 'H U*014416 Frsz NJ 9 DEC. COUNTER
0004 1 * 0007 7 1 F J MP RR2
00042 *060415 F FLO N JF
0 0 0 4 3 * 0 2 4 4 14 F L DA 1* N J F
00044 V1 1A A A A A
f.k.»f .\.»V FDI V 1* 0 9 CAL. RRCJ)
0004 5 *0640105 r L 13 tempi 9 LOAD STORAGE P O I N T E R
A A A •
0004 6 •041400 F STA ^ • 9 XJ 9 O 9 STORE I N R R C J )
r\
0 0 0 4 7 * 140001 F F DC F 9 P R IN T RRCJ)
A A A C A
v *U U Jv .' * 1 1 0 0 0 0 F.I C3 9 INC. P O I N T E R
9 «AA A A k A A
00051 v' K ' o FST3 TEMPI 9 STORE I T
• | AAAAA
00052 1 U U U V
.*k
.» F EX T
A A A r A
v*^' k.* 3 3 * 0 0 6 0 0 1 ? JSR £ • RESLT 9 F 0 RM A T TH E R E S UL T
00054 *030004? LDA 2* J 9 L O A D COUNTER
00055 •000731 JMP RR 1
f AA AA A A \T 1t /A
0 0 0 5 6 V■O \J»V* V* V-* V X J • V»
**\
0 0 0 57 • 0 0 0 0 0 0 N J F :
256
V >>.■O 2 CORRL
• TI TL IMPUL S E
• ENT H
• EX T P C 5# C 5 P 1 # • N EX T P # . I X T 1 # RXXO * J » K # M- » S Uh
. EX T L SUr* P I * T EN. PC » T E M P I , . R E S L T * D E L û M
. N R EL
A A A A A
V'V.' V' V» 0 5 4 5 0 4 H: S T û 3» R I O P IMPULSE F UM C i l 0 M
A A A A 1
V.»V* V.» V.» 1 0 P400 15 LPA 1 »C5 P ROUTIME
A A A A A
V.»V» V.* V» C * 044002$ STA 1* C 5 P 1
A A A A Q A A /•A A A
V
.* \»* O L* vJ C* vj v.» O J ' JSR ® • N EXTP
/> r\ r\ ^ A A / A A /, (
V» U U U U H Î JSR e - TXTT •
P PRIMI HEADIMG
c
V.* V
. »*>
_» O 177777 IMPULSE
A A A A ^
177777 F EM T
A A A A T
V' V» V.» L' / 066475 F LT 3 «RIO P POINTER TO RXX
A A A I A
V. ' V' V.' 1 V' 021400 FLl A A A O
9 vj P vJ
A A A 1 1
V* V»V» 1 1 040005$ F STA 0» R X X 0
CCO 1 P 070014? F ST 3 TEMPI
0 0 0 1 3 J A A A A «A F EX T
0 0014 126400 SU E 1» 1
AA A 1 C
V.' V .» 1 D 044006? SIA 1» J P SET J=0
A A A 1 /■
\w ' V*V> 1 P 034466 L L 'A 3 » R 10
A A A 1»7
V.1V.' L' 1 / 025401 L P A 1» 1»3 P POINTER TO RXY
A A A A A
$»V» O r * v.» 044013? S TA 1» TEMPO
/A r\ r s o 1
o V* r . I 10 2 4 0 0 H 1 ï SUE P
A
VJ
A
VJ
000 2 P 0 3446 P L L A 3 » R 10
/•
0 0 0 P 3 V.' M A A A A 7
V.»v. »«j f J
<r
S T A 0» H P SET K =0
0 0024 0 4001 1$ S 1A 0 » S UM- P SET SUM TO 0
0 0 0 P 5 0 4 0 0 1 P? SIA 0» S U M P 1
000 P 6 0 3 1 4 0 2 M 'A 2» P» 3 P P O I N T E R TO Hirt)
0 0 0 P7 0 3 4 0 1 4 $ L PA 3» T E M P I P S E T P O I N T E R 10
A A A OA
V_* V..1\.* O V
.• 024006? L P A 1» J 1 R X X ( J - K ) W I T H n=
0 0 0 31 1 P 7 0 0 0 ADP 1 » 1
000 3 3 0 5 4 0 1 4 ? STA 3 » 1 EM P 1
0 0 0 34 0 P 4 0 0 6 ?H2: L P A 1 » J
0 0 0 3 6 .A /A »A n A
V.* $»v.» M O V» JM P H3 P YES
A A A A A /; «
0 0 0 37 V*VJ V' \J V# V? r EM T •
P K< = J - 1
A A A /. A
V»V.» *-* * 064014? FLP3 T EMPI
Q •
000 4 1 0 P 1400 FLPA A
VJ 9
A
VJ » O P LOAD RXX(J-K)
0 0 0 4 P /A A C A A A
U 6 J U VJ O rLDA 1 » 0» 2 P LOAD H( K )
1
0 0 0 4 3 1O V.‘ V« VJ
,-N /l A A A
F IC2 P INC. K POINTER
000 4 4 014014? F US Z TEMPI P DEC. J PO IMI E R
00045 014014? F PS Z I EM P 1
00046 104100 FM P Y 0 » 1
0 0 0 47 0 2 0 0 1 1 ? FLDA 0» S U M
V'
«A «A A C A
V* VJ J u 1P 3 0 0 0 F Apr 1»0 P ADD TO S U M
000 5 1 0 4 0 0 11? F STA 0» S U M
O
6> <A *A C O
v» v.» sj 1
« A A A A A
VJ Vj V.* VJ VJ f f:x t
0 0 0 5 3 »i «**
/A 1 A A A O <T
V V VJ V I SZ K •
P I MC * K
00054 CP0007 ? LDA 0 »K
00 0 5 5 000757 JM F HP
0 0 0 5 6 0 0 0 0 37 *M 3 : F FM Í
259
0002 IMPUL
00057 *0 2 0 0 0 5 * FLDA A
1 ilAAv/
V} V V A
G • RESLT
75*034010* LDA 3*M ; LOAD M c o u n t e r
76*0100065 I SZ J J INC- J
77'024006* LDA 1*U
00100*166414 SUB# 3*1* S Z R 5 J=M ? F I N I S H E D ?
O l o J W K . ’ I G yj JM P H1 i NO* C O N T I N U E
00102*034402 LDA 3* RIO J YES# EXIT
^»v.<K'O 1HuO
1 A Q I A A 1 / . A ' )
JM-P 3* 3
A A 1 A /, • A A A A A A n 1 A # A
u v 1 \.*h n l ^ t V.»
END
260
Enter
( Subroutine^
I
Save return address
Set result format
" '1 --------
Set pointers to
H. Y, X________
I
CalH(0)*X(0)*AX
2 1
¿X[H(l)*X(0)+H(0)»X(n
Add It to sum
ISet
Calc. Y (J)
Load X(K)
H(J-K)
1 -i
Store and print
Calc, summation part Y( J ) ___________ _
of Y(J)
1.... .... Inc. J
Increment K
0001 CONVO
. T ITL C O N U O L U T I O N
.ENT CON
. EX TN h EN T
. E X TI C 5 , C 5 P 1 , J , K , M , S U K
. EX TP T EK P O ,T EM P 1, DEL AK , • RES L T
. N H EL
r»
3* R 1 2
a a
V.» V./ *054561 CON : STA ; CON V/OLUTI ON ROUTINE
r\ r\ r\ i
V'V.»V.» V *1 * 0 2 4 0 0 1 ? LPA 1* C 5
A A A A A
V »V.» V » CL • 0 4 4 0 0 2 $ STA 1 ,C5P1
A A
V.»
A A O
V.» V-l O * 0 2 1 4 0 0 LPA w» 9 AV' 9A O
A •
9 P O I N T E R T0 H
A A A A
K.» *025401 LD A u 1,3 ; P O I N T E R T0 X
#a r\ r\ r\ c LDA
v» v.' v.' Z) * 0 3 1 4 0 2 2 ,2,3 » P O I N T E R T0 Y
A A A A f
V 1 1.1 ' O * 0 4 0 5 6 0
. V V V. STA 0 , TP0 5 (TEMP. POINTERS T0 )
r\ rs r\ r\ »7 STA
'044560 1, T P 1
A A A 1 AN
1 V' * 0 5 0 5 6 0 STA 2, TP2
A A A I 1
XX'1 1 1 0 0 0 K0 V 0, 2
P\»A.“\ 1 Q F EN T
v.n.a' X c- * 1 7 7 7 7 7
A A A 1O
v'Vw'v.' X O * 0 6 4 5 5 4 PLI 3 TPI ; P O I N T E R 10 X
fA A I AAA P L LA r\ r\ r\ •
LOAD H C0 )
0 0 0 1 4 v.' Xu v . ' U ü \J S \J> » d. »
A A A Í 7
v.» v.» v. ' X / * 0 4 4 5 4 3 ESTA 1,X0
A A A Q <*S
V» ¿ '0v^> * 0 3 0 0 1 1 $ P LDA 2 , DELAM
000 PI *1 04 1 0 0 F MPY 0,1 ; H <0 ) * X ( 0 )
00022
• 1 O A
l O u
1 A A
i u u FMPY 1,2
00023
« 1 C A
X 03 v.» O
00
A A
F HLD 2,2
•
A A A
V
A /| 9 1 A
X
/ .A A A
*-* v_* v.» v_* F I C2 9 INC. POINTERS
f 1 1 A A A A
F I C3
00025 1 X U L ' V.» V.»
000 26 * 0 2 1 0 0 0 PLL A r\
9 KJ 9
0
d
•
9 LOAD HCl)
0 0 0 27 * 0 2 4 5 3 3 F L DA 1,xo
A A A Q A
v^» Ov. » * 1 2 0 1 0 0 P MPY 1,0 ; H ( 1) * X (0 )
A A A A 1
V.* V.' O X * 0 2 5 4 0 0 P LD A 1 , 0,3 ; LOAD X ( 1)
A A A A A
X* V.» O *034532 P L DA 3 , H0
A A A O Q
V- ' UV. ' O J * 1 3 4 1 0 0 P MPY 1,3 ; H( 0 ) *X ( 1)
0 0 0 34 * 1 6 3 0 0 0 F A LD 3,0
9 1 A A A A A a r\
^ a a t c
v,/ v>* v.» O X O u FHLV K.> 9
A A A
V 1 Q /
v_> C *0 2401 1$ FLDA 1 , DEL AM
/*>A A Q
v.' v ' v.» O
7*
f 1 2 0 100 P MPY 1
1 A
* O
•
A A A ¿J • «
X
1 A A A A
X V.» V* V.» V..» P I C3 9 I N C . AN L STORE
000 4 1 *0 70 5 2 6 PST3 t p i ; POINTER
0004 2 *064526 F LD3 TP2 ; LOAD P T . T0 Y
000 4 3 * 051400 P STA 2 , 0,3 î STORE IN Y C0 )
9 1 1A A A A
0 0 0 4 4 1 1 uv. »v. »v. » F I C3
r\ Q •
00045 *041400 P STA K.' 9 k.* 9 sJ 9 STORE IN Y ( 1)
0 0 0 4 6 * 150001 F P DCP 2
« 1 1A A A A
00047 X X V.»V
.»V
»* F I C3
•AP» A C «A
v » v_» v* 03 v_» * 0 7 0 5 2 0 F SI 3 TP2
« 1 A A A A A
0 0 0 51 X V.' V» V.»V.»V* P EXT
0 00 52 * 0 5 0 5 1 4 STA 2, TP0
v.1
A A A C O
.v *v.1 03 O • 0 0 6 0 1 2 $ JSH (i. H ESLT
« AAAA 1 O •
000 5 4 v* v » X F EN T
A
0 00 5 5 ' 1 4 0001 P P PCF U'
• i A A »AA A
0 0 0 56 i V
v' V V' V
.»VJ r EXT
262
0002 CON VO
00057 006012? JSR @ •RESLT
^ rs r\ £- r\ rs rs
V-'V.*V.»O V»* 102520 SUEZL O 9 v.»
»*>AO/ 11 r\ rs
101400 INC KJ 9 \J
v.»v.*v.'A
O OO 0 2 4 0 0 5 S C 1s LDA 1 *M
A
v.»u o OA oQ
A A A /.A A A 7
u M vjv.»vjO J STA 0# J
00064 034475 LDA 3* R 1 2 ; LO A D EXIT ADDRESS
000 6 5 1 0 6 4 15 SUP* 0 * 1 * SNR ; J =M ?
00066 001403 JMP 3*3 * YES* EXIT
A
v•v A.'A A 7/
v.»O 030477 LDA 2 * TP0 ; N O * S E T PT. TO H ( J )
r\ r\ r\ n r\
V.'V'V f v.» 000054• FENT
A'A
L L*»>7# 1 1 104000 F ICP
0007 P 064475 FLD3 tp l ; SET POINTER TO X(J)
A
L'V A.•V
A*7/ VAJ FLDA A 9 A\J 9A d • L O A D HCJ)
0 2 1 0 0 0 P
A.'A
v v.»A 7# */*i4 025400 FL DA i*o*3 ; L O A D XCJ)
rs r\ r\ ^ j c
l»v.» / O 1 1 0 0 0 0 F IC 3
v*l'v**f7O/- 070471 FST3 TP 1
A A A 79 7 9
V.'V'O 030463 F L DA 2 *X 0
AA V.'11VA.»O
A r\ r\ • H ( J ) * X ( 0)
140100 FMP Y d P kj p
AA u U» 11
1 A FLDA
030463 2* H O
A
V.' 1A 1A A
v.fA 144100 FMPY 2 * i ; HC0)*X<J)
A A 1
V.' 1 v.»O A Q 1 o a nn
FADE:
1 U / v-» 0 » 1
JV-N»^U 11 ,
v-
.\
«^/l 124300 FHLV \> 1
. a /-\i c:
V.*V» I V.»J 044006$ FSTA 1 9 SUM
A
V.»A V1 11A
L'C / jAAAAA F EXT
A A 1
V'V.» k» 9 A 7
050457 STA 2» T P 0
A
O AV.»1 111A
050007$ STA 2 9 TEMPO ; HCJ)
0 0 1 1 1 030450 L DA 2# HI 2
0 0 1 IP A
v'TOt ZJ v
A."A I 1
LDA 3*1*2 ; POINTER TO X
00113 054010? STA 3* T E M P 1
0 0 1 14 SUPZL A.»9 A
V O
10 P 5 2 0
00 1 15 0 2 4 0 0 3 ? C 2 : LDA 1 *U
00 13 3 0 P 5 4 0 0 FLDA i*c*3 ; L O A D X (K )
0 0 1 34 104100 FMPY o*i ; X CK )* H ( J - K )
0 0 13 5 0 2 0 0 0 6 ? F L E'A 0* S U M
00 13 6 1 A 7 A A A
1U / V» F A DD 0 * 1
00 13 7 0 4 4 0 0 6 ? F STA 1* SUM
1#-># *\A A rs
0 0 1 40 .• V
I v .‘ K,> V.*Vw' F EX I
A A •
0 0 14 1 101400 TMC s^» P > INC. K
0 0 1 4 P #A «A'‘J ^ JMP C2
0 0 14 3 0 0 0 1 3 0 • 0 3 : •F F.N I
0 0 1 4 /i 0 6 4 4 P 4 FL1 3 TP2 ; POINTER 10 YCJ)
00 145 0 2 0 0 0 6 ? FLDA 0* S U M
oo i n e 024011? FL1 a 1* DELAM
0 0 1 47 1 0 4 1 0 0 FMPY 0 * 1
0001 IMVEH
.TI TL IN V E R S E C O N
• ENT ICON
. EX TW F ENT
. EX T P C 5 , C 5 P 1 , T E M P O , T E M P I , DEL AM, J , K , M , S U M
. FXTD SUM P I , . R E S L T
. MR F L
o o o o o 10 54533 icnws ST4 3 , R 12 J INVERSE COMVOLUTION
»-Nr\AA 1 r r\ /, »Arv t <T LPA
v.H» v.» v.» 1 v.» r 1 1 , C5 9 ROUTIW E
00002* 0^4002? STA 1 ,C5P1
00003*177777 F KW T
00004 *066477 •
FLD3 0 RI 2 9 P O I N T E R TO IMPULSE
r\ r\ 0
00005 *021400 FLPA \J 9 \j 9 O 9 FUNCTION H
r\ tj /*%
00006 *040473 F STA s.» * n ^
r\ r\ r\ r\ *7 » a n a a a y, ff
v•v»v_' v* / v.* / v» H r F ST3 TEMPI
A A A i A t 1A A A A A
V»V*V» 1 I F EXT
0 0 0 1 1 ’ 126400 SUF 1, 1 9 SET J=0
00012 *0 44006 $ STA 1, J
00013*034470 LPA 3 , R12
00014*025401 I.PA 1 ,L 3 9 P O I N T E R TO OUTPUT
000 15 ’ 0 4 4 0 0 3 ? STA 1, TEMPO 9 SIGNAL 1
A A
0 0 0 1 6 * 1 0 2 4 0 0 CO 1: SUF X
^>* V» 9 SET K=0
r\ r\ r\ \ <7 f r\ /x rs r\ n <r ST4
\_»v'v.' 1 1 Vw' f - 0>K
•
k»V. X.'Y0V_' 1 1- STA 0 , SUM 9 SET SUM TO 0
00021*040012? STA 0 , S UM P 1
000 22 * 0 344 6 1 L1A 3 , R 12
00023*031402 LPA 2,2, 3 9 P O I N T E R TO I N P Ü I
00024 *034004? LPA 3 * TEMP 1 9 SIGNAL X
00025* 024006? LPA li J 9 S E T P O I N T E R TO
•A Q « 1 Q7 A A A ALL 1* 1 H ( J - K ) wI TH K = 0
• .» u o C
jv 1 £6 / v.» x.» 9
A A A QH f 1 A O A A A ADF 1,3
VPV_*V*AA I L yJ f X.»
0 0 0 30 * 0 5 4 0 0 4 ? STA 3* T EMPI
0 0 0 3 1 * 0 2 4 0 0 6 ? CO2: LPA 1, J
00032*106415 SUF* 0*1 SWR , 9 K=J ?
A.A ,A Q Q f A A A /. Q A
#
•
„ V.'V.» J vJ X.' V' i-
-W* JKP CO 3 9 Y ES
*A/AA A A /. 9 A A A A «A A «
JH
v'L'v.' v* x^/v.» O F KW T J H< J - 1
00035 *064004? FLL 3 TEMPI
»A#AfA A r 9 <
AQ 1 /. #A#A r\ 0
ov.» ü C w CO i x^x
.* FL r A v.' 9 \,f 9 O 9 LOAD HCJ-K)
/A/A# AO *7 9 A A C »AO A
V»V. » ^ ' 0 I C
Z- D O V.' V
_* F L DA 1,0,2 9 LOAD X( K)
AAA y, .A 9 1 «A /, /A.A »A
V»X
^* X^' 1X ^
'*-•V .»V
.»V
** F I C2 9 INC. X POINTER
000 4 1 * 0 1 4 0 0 4 ? F LSZ T EMP 1 9 DEC. H POINTER
00042 *014004? FPSZ TEMPI
0 0 0 4 3 * 104 100 F MP Y 0* 1
0 0 0 4 4 * 0 2 0 0 1 1? F L PA 0 , S UM
00045 * 123000 F App L 0
00046*040011? F STA 0 , SUM
A A A /, T • X A A A A A
V •1»V » *-« f l v.»x.'v»\.'v.» F F.X T
0 0 0 5 0 ' 0 10 0 0 7 ? I SZ K 9 I NC • H
*
“x \f
200 5 1 * 0 2 0 0 0 7 ? LPA v.* 9 c\
0008 IN VER
^ a c r 1 1 A A A A
\_t V.« Zj O 1 F I C 3
/n C *7 A 7 A A A A ff
V.I V.IV.I J / \J • v.* \ j y,t %J J> F S T 3 T E M P O
0 0 0 7 4 0 2 4 0 0 6 $ L D A l i J
A A A 7 C •
kj y*' \ j 1 Zj 1 6 6 4 1 4 S U B # 3 # 1 » S Z R » J = M ? F I N I S H E D
a A A 7 / A A A 7 A A •
S*'S*»KJ i O V.» I c u J M P C 0 1 9 N O * C O N T I N U E
A A A 7 7 •
U\.K* f / 0 3 4 4 0 4 L D A 3 * R 1 2 9 Y E S * E X I T
A A 1 A A
K » V 0 0 1 4 0 3 J M P 3 * 3
A A 1 A 1 A A A A A A M A • A
1 V.' 1 o u u u u u nvy • VJ
A A 1 A A A A A A A A A
l u c u u u u u u \j
A A 1 A A A A A A A A n 1 A # A
v_' \^» K ' O V.»U U U V . » C » n 1 G « \j
END
Enter A
C
z
Subrout ine y
T z l
Save return address
Set result format
T
„ , t ttR \. . Load 2 DT
Calc, cosi— )+ 1 Multiply by sum
*m •
I
„ , /ttRK \
Calc, cos-- Load storage pointer
* m /
I 1>
Load RR (R) Store and print
P.S.D.
A A A «
SPECT
• TITL SPECTRAL
• ENT G
• EXTN P E N T # S P C E 3
•N R E L
A A AAA
V.' KJ V
.» • 0 5 4 5 3 6 G: STA 3» RI 1 •
* POWER SPECTRAL DENSITY
AAAA 1
V-M.' 1 •024001$ LDA 1» C5 * ROUTINE
A A A A A
V» V' V
.» O •044002$ STA 1» C 5 P 1
AA AA Q
u L» V
.»v.* O • 0 2 1400 LDA A A
U #U #
7
Ü i POINTER FOR STORAGE
A A A A J , 1 A /«A A A A ff
V
.»U O U H u H \ Jv .» ü O J> STA 0 » TEMPI
A A A A C
V
.»V.»V
.» J •006014$ JSR CRLF
r\ r\ f\ r\ r
V.»c» y j O •006015$ JSR © T XTT
.
A A A A 7
V
_*# • 177777 SPCE3
A A A 1A
V
.»k j V
.» X V
.» • 102400 SUB A
U J
A
L» 5 SET C O U N T E R S
0001 1 ' 0 4 0 5 3 3 STA 0»MF # M FLOATING
AAA 1A
V
.»k j V
-*1r. •024004$ LDA l» M
AAA
1
1 A
O*044532 STA 1 »M F + 1
^ ^ ^ i y»
v
_*V» <.• 1 * A ’0 31401 LDA 2» 1» 3 •
P O I N T E R TO A U T O C O R
AAA 1 C
v.*ov.' D 1 • 177777 FENT •
9 RELATION ARRAY
r
UV.'L»1 D
A A A I
•060526 FFL0 MF
AAA
V
.»
1
i_*v,*1 r 7 •A AAA A
U C U U uO/! <
T
T F LDA 0» DELAM
A A A AA
L' V
.»V.' CO • 0 2 4 0 1 1 $
V» FLDA 1» TWO
AA A A
K.» CZ 1
• X1 O'"' 1rsr\
1
c . kj FMPY 1»0 1 2*DT «
AAA AA
•040515
K .' K.» K.r C m C . F STA 0» DT2
A A A A O
V
.» kj r> O •024521 FLDA 1»MF
A A A A
(Z>
>1
• 1 04100 FMPY 0» 1 Î 2*DT*M
r \ r > r \ r^s tz
K.» U U G J 1L»^
•A A A A | A
UGUW* FLDA 0» O N E
.1 ^
AA A A^
v
t •
c:O loo oo oc\<\ c\r\ a
EDI V 1
l
A
KJ
9 Î FREQ. S T E P HZ.
A A A Q7
V
.»V.' V
.» f • 140001 FFDCF A
KJ P R I N T F I R S T STEP
9
A A A 7 A
V
-' U v) V
_*
fo o l u u u
A A | A A A
F LDA 0 » 0» 2 LOAD FIRST AUTO
•
A A A A | V 1 A / .A A A
V
.»u u O X 1 FI C2 •
CORRELATION VALUE
#
A A A 7 A t A /.A A a 7
L» 1 #O
<
O H V.» U U f ^
r
FSTA 0» RXXO
A A A 7 7 « 1A A A A A
v.' v
_*O O 1 u u u u u FEXT
•006014$
A A A 7 /j
V
./ JSR ©• C R L F
a a a o c
•006014$
v 'U U O J JSR ©.CRLF
AAA 7 /
uuO O •006014$ JSR ©•CRLF
AAA 7 7
v.' V • 102400
.»v.» O f SUB A
9
A
SJ
V
.»V_
»V.» • 101400 Gl:
A A A /. /A
INC A A
U# U 9 S E T TO 1ST. FREQ.
00041 f A /, A A A C tt
U H V .IU U Jλ ST A 0» K
A A A y, A
V' V •040505
.»C' ^ C STA 0» K F + 1
00043 •024004$ LDA 1» M
00044 •034472 LDA 3» RI 1 9 LOAD EXIT ADDRESS
0 0 0 4 5 • 122015 ADC# 1 » 0» SNR 9 K >M ?
«
00046 •001402 JMP 2» 3 * YES» EXIT
f\
0 0 0 4 7 • 10 2400 SUB A
V 9 C*
A A AC A
»
V V.» V
.» •040476 STA 0» KF 9 STORE IN FLOAT. K
00051 • 04047 1 STA 0» RF
r\
00052 • 10 1400 INC a
yj » yj
Ar\r\ a
v.» c í SPECT
u v.O O 040013$ FSTA C* SUM ; SET SUM = R X X C O
r\ r\ r\ c
'JU U J
*7
/ 060467 F FL0 KF
r\ r\ r\ g r/ \
KJKJKJ O
1 A n n o «
1 UU U L»U FEXT
A A A i! 1
V
.»V.»O 1 1 2 2 4 1 5 G2 J SUB# 1*0» S N R J R > M ?
/>/ ■> r\ ¿l o
K J KJ %*J U ¿ 000434 JMP G3 *«
YES* CAL. G (K)
r* r> o
v.»v.»v.' O O 040460 STA 0* R F + 1 ♦# L O W P A R T OF R
A A A ¿ /i
X.' v.» O *4 000054* FENT
*
A A A i! C
v* v.* O O 060455 FFL0 RF FLOAT R
AAA / /
O O 020454 FLDA 0# RF •
LOAD R
r \ r \ /r * 7
v->O f 024455 FLDA 1*MF 9
•
LOAD M
A AA 7
O
A
030016$ FLDA 2# PI 9 L O A D PI
AA A 7 1
V' V.' # 1 140100 FMPY 2*0 #• P I * R
A A A 7 A
f 0
1
1
r \ r \ c \ r \ r\
C , \ j C* k j y j FDI V 1*0 9
•
PI*R/M
A A A 7
V
.»V.' V
O
.' / O 110060 FC0S 0*2 # C0SC P I * R / M )
•
A A A 7 /i
v 'c» o / *4 0 3 4 0 IOS FLDA 3» O N E
A A A 7 C
v-»uu r j
1
1
7 0 AAA
/ O l »v^v.* FADD 3*2 •
CC0SCPI*R/M>+1)
r\ r> r\* 7
vj v.» i
¿L
O 024450 FLDA 1» K F
A A A7 7 1 AA 1A A
V-'V'V' / / i Cv.» l u u F MPY 1*0 •
9 PI*R*K/M
A A 1A A
V.»V* l u u I
I A A A /! A
U U O V
.» O FC0S A
U #
A
L» 9 COSC P I * R * K / M )
AA
IV.» 1
1A 1 1 1 ^ 1 AA
1 l u lu u F MPY 0*2 I (C0SXC0S+1 )
•
AA 1A A
V.» l u c
A A C A A A
U ü JU U U FLDA 1*0*2 9
•
LOAD NEXT AUTO
AA
V.' 1O
1A A
O
1 A /. A A A
1 UHUUU FIC2 *• C O R R E L A T I O N V A L U E
A A I A /,
1 44100 FMPY 2» 1
AA
1 I AC
U J I
A A A A 1 A ff
OO V _* O FLDA 2» SUM
\w
M w* 1^0
A A 1A £ 1O A A A A
FADD 1*2 •
9 A D D TO SUM
uu lu 7#
A A 1A
050013$ FSTA 2» SUM
11
AA 1 1A
V-'O VJ 1 V.»\^
f A A A A A
FEXT
00111 102400 SUB A . A
Vj >» O
00112 040430 ST A 0» RF
00113 010426 ISZ R
00114 020425 LDA 0* R
001 1 5 000744 JMP G2
001 1 6 0 0 0 0 6 4 *G3í FEMT
00117 020420 FLDA 0» DT2 #
DIGIT INTERVAL
V»»V' I C V CAL. GCK)
A A 4 AA
-» 110100 F MPY 0*2 •
00 14 1 UUV-'V., V-1'-1 nD •♦
A A A A A A A
*
END
270
OGGI X SPEC
• T I TL X S P E C T R A L
• EN T GXY# P 180
• EXTN FENT# S P C E 3
•NREL
AA A A A
y u L» c» V-» •054514 G X Y l STA 3# R20 #• C R O S S - S P E C T R A L
A A A A |
y %J U X•024001$ LDA i#C5 • DENSITY FUNCTION
#
A A A A A
K.I y y V.» C* • 0 4 4 0 0 2 $ STA 1# C 5 P 1
v.» u u u O • 0 2 5 4 0 4 LDA 1# 4# 3 # S T O R A G E PT. F O R
r \ c\ r\
A A A A /»
v.»o VJ v.' *4 • 0 3 1 4 0 5 LDA 2# 5# 3 #• S T O R A G E PT. F O R
A A A A r
v.» v.» u O •044003$ STA 1# T E M P O í MAGNITUDE
A A A A / f A C A A A /j <
o O P
STA 2 # TEMPI a
» PHASE
A A
A A A A 7
í• 102400 SUB y 9 y #a S E T H I G H M=C
* y Iy • 0 4 0 5 0 5 STA
A A A 1 A
u l 0#MF
r \ r s r\ \ «
u v.» u 1 i • 0 2 4 0 0 5 $ LDA 1#M
A A A 1O
v_#y \J X cL • 0 4 4 5 0 4 STA 1#MF+1 #a S E T L O W M
y j y j y j X O ’0 0 6 0 1 1 $
A A A 1 Q
JSR @.CRLF #a F O R M A T
r \ r \ c \ i/,
y y y X ^ •006011$ JSR § • CRLF
c» u u i J • 0 0 6 0 1 2 $
A A A i c
JSR 0.TXTT
A A A 1 r •17 77 7 7
y y y XO SPCE3
C» U U 1 í• 177777 FENT
A A A « *7
A A A AA
V-' V.»V.» <-* y
• A A A A A 7 <
V.»
r
1 41 FLDA 0# DEL AM p CAL. 1ST FREQ.
A A A A
y y
1
cí X • 0 2 4 0 1 4 $ FLDA 1 # TWO • A N D FREQ. STEP
p
2 * DT
A A A AO • 1 Q A 1A A
v.» y.* cz. r . X C*y l y y FMPY 1#0
•
p
A A A A A
V.» V_r O O •040511 FSTA 0# DT2
O a A O /i
O Nw» y d * * + •060471 FFL0 MF
n n / > nc:
y y y C* Zj •024470 FLDA 1#MF
a a a a A
O O U u D • 104100 FMPY 0# 1 p 2*DT*M
A A A 0 7
v.» U U G f
« A A A A
y cz y y
1A C
l u x F LDA 0# ON E
A A A O A • | A A A A A
V.' W J V.» l c u o y y FDI V 1#0 9 FREQ. S T E P HZ.
A A A O 1
y y y y X • 140001 FFDCF A
y
•
9 P R I N T F I R S T STEP
f
A A A 7 A
u u O c
$A A A A A
Xy y y y y F EXT
A A A 7 7
u v» u O O *006011$ JSR 0.CRLF
A A A o y,
•006011$ JSR @.CRLF
o v.» v.» O v) ’0 0 6 0 1 1 $
A A A 7 C
JSR @.CRLF
v.» y y y y • 0 3 4456
A A A Q /:
LDA 3# R20
A A A 7 7
^ V_r O 9 •021401 LDA 0 # 1# 3 •
9 P O I N T E R TO RYX
»V.' *-* V.' • 0 2 5 4 0 2
A A A /.
v.» Vw LDA 1» 2# 3 9 P O I N T E R TO ACR)
00041 •031403 L DA 2 # 3# 3 9 P O I N T E R TO BCR)
00042 • 0 3 5 4 0 0 LDA 3# 0# 3 0
9 P O I N T E R TO RXY
0 0 043 • 0 5 4 4 7 3 STA 3# TP0 9 RXY
0004 4 • 0 4 0 4 7 3 STA 0# TP 1 9 RYX
00045 • 04447 3 STA 1# TP2 0
9 A( R)
00046 •050473 STA 2# TP3 9 BCR)
00047 • 0 30467 L DA 2# TP0 9 LOAD POINTER TO
A A r
a
v.» V.» v .* O y • 1 0 2 4 0 0 SUB V.» # *.<
í*m t l
K » y y y 1 ’ 0 2 4 0 0 5 $ LDA 1#M
.A
v * O v.*
c O •106015 GX1 : A DC# Ó#1#SNR ; K=M?
00 0 5 3 •000^26 J MP GX 2
00054 •000017* ir E U T
00055 •021000 FLDA 2 ; RXYCK)
vuv
r>
' oG.
rv r\
XSPEC
v-' U V .O O 064461 FLD3 TP 1
c;7
v.»\.»U J f 025400 FLDA 1 0» 3
»
•
0 RYXCK)
r> r\ r\
v.»k.»v» O c*
1 1 AA AA
1 Xu u u u FI C3 •
# INC. A N D S T 0 R E
A A A /* <
k
,#V.» O 1 070456 FST3 TP1 •
POINTER
A A A / A
k. ' k j ^ i O d 1 uM uu u FIC2
A A A £ Q
v.»v' v.» O O
1 7 1A A A
l O 1 u u u FMOV 1,2
00064 1A 7 A A A
iu i u u V
.» FADD 0 1 , 0 2 * AC K )
142400
A A A/* c
uv'v^O J FSUB 2, 0 9
0 2*BC K )
v» v
~ 064452
r> A A / /:
»D O FLD3 TP2
00067 0 4 5 4 0 0 F STA 1,0,3 0 S T O R E IN ACK)
A A A 7 A | IA A A A
V
.»V
.»V.’ # 1 1 \ .» u u u FI C3
0007 1 0 7 0 4 4 7 FST3 TP2
7
A A A
V
.» V 064447
A
.' / G. FLD3 TP 3
A A A 7 A
041400
U V 'U f J FSTA A
U l
A
ux
Q
O 0 S T 0 R E IN BCK)
A A A 7 /. | 1A A A A
T
w»V
.»V
.# / X XV U U U FIC3
o / ~
l.' ^
>* 7 c:
v.*
/ O 070444 FST3 TP3
A A A 7 /■ 1A A A A A
O vJ / O 1 U U U U V» FEXT
A A A 7 7
i / Xu XH u u INC A
V#
A
U
AA | A A
V
.» 1V.»V
_*
A A A 7 CO
u u u f Og JMP GX l
A A « A 1
v.» K» J 0 3 4 4 1 3 GX2: LD A 3, R20
A A 1A A
k' lu e 031402 LD A 2,2,3 •
0 P O I N T E R T0 ACR)
A A 1A A A A A A C / ,»
v' 1 O UUUU FENT
| r> r \ r y
^ V» 1 U H
a
u e
q
X L» vue* FLDA A A A
U X U X G 0 L O A D 2 * A C0 )
A A 1A C 1A A 7 A A
V.*u I u j X u u O u u FHLU A A
UX U
m
0 CALC. ACC)
A A 1 A /"
V
.» 1 O 124300 FHLV 1,1 9
0 CALC. ACM)
s\ «
a i
u u lu
*7
/ 040415 FSTA 0, A0 •
0 S T O R E IN AC0)
A A 1 1A
^ u» 1 X 044416 FSTA 1, AM 0 S T O R E IN ACM)
| A A A A A
0 0 1 1 1 X v.»^ vuo FEXT
0 0 1 1 2 102400 SUB A
0
A
KJ
00115 MF:
A A A A A A
V.» V.' V
.»U
A
U
00116 A A A A A A
^ U»V.»V.' V_'
A
V
_*
001 17 KF : A A A A A A
V
-*Vj L' U U U
A
<.*
»A< A 4 or\ A A A A A A A
v.» 1 e v-
ì V
-» V
»*V
.»V.' O
00 1 2 1 A A A A A A
^ VJ V.I U»V
.» R ?
A
001 2 2 A A A A A A
V.' U U U V» RF: A
V
_»
00 123 A A A A A A
U U V.' V
.»V
-»
A
U
00124 A A A A A A
V.' V
.» V
.» A0: A
U
00 125 A A A A A A
VUV
.»V» V.»V
.»
A
V.'
00126 r\ r\ /
VJ V
a r\ rv
.*< J \ J VUv
_» AM: A
V.»
00127 A A A A A A
V
.»VUV.»V
-»V
.»^
A
VU
A A 1AA
y~* y..> 1 o v»
A A A A A A
V.»V
U SUM : 1
A
U
A A A A A A A
0 0 131 Vj W*V
.» VUV.»
00132 A A A A A A
V.' V' v
_» v„' V
.» SUM2: A
V.'
00 133 A A A A A A
V.' V
.»U» V
.»V
.»
A
V.*
00 34 1
A A A A A A
V.»VJ V.»V
.* DT2: A
V
-*
00 135 A A A A A A
V.» V.' K.J
A
U
00 36 1 OW .» U TPOs
A A A A A A
*V
.»O
A
V.* ; RXY ( TEMP. POINTER T0 )
00 1 37 OV.»V* TP :
A A A A A A
<.» V 1
A
V.* ; RYX
00140 -'K
-*C V'TP2:
A A A A A A
V U'
A
V
.» ; AC R)
00141 TP3:
A A A A A A
v_( vu y.* ^ v»
A
v
~ ' ; B(R)
0 0 142 101400 GX3: INC A
V.»0
A
KJ
•
0 SET T0 I S T . F R E Q .
00 143 040006$ STA 0, K 0 SET L O O P C 0 U N T E R S
00 144 040754 STA 0 K F+ 1
,
00 14 5 126400 SUB 1, 1
00 146 044754 STA 1, RF
00 14 7 044750 STA 1,KF
00 1 50 024005? LDA 1,M
273
0003 XSPEC
0 0 1 5 1 • 122015 ADC# 1# 0# SNR 9 K = M ? (FINI SHED)
00152*000541 JMP GX 6 9 Y E S , G 0 T0 P R I N T 0 U T
00153*000103* FENT i NO, C O N T I N U E
00154*060743 FFL0 KF
00155*126400 F SUB 1, 1
00156* 0 4 4 7 5 4 FSTA l* S U M 2 9
9 SET QCK) T0 0
00157* 0 2 0 7 4 5 FLDA 0# A0 9 L O A D ACO)
00160 * 034746 FL DA 3* AM •
9 L O A D ACM)
FEXT ♦
o o
9
XSPEC
0 0 2 4 4 0 0 0 204* GX5X FENT
00245 024667 F LDA 1 DT2
*
«
9 2*DT
00246 034662 FLDA 3 , SUM1
0 0 2 4 7 134100 F MPY 1*3 •
9 CALC. COO
/n
a a q
054660
c
FSTA 3# SUM 1
00251 03 0 6 6 1 FLDA 2 , SUM2
9
A A O^O
U U G JC 1o u Iv .iv FMPY 1,2 9 CAL. Q(K>
00253 050657 FSTA 2 , SUM2
0 0 2 5 4 150100 FMPY 2,2 •
9 G(K>*Q<K.)
0 0 2 5 5 174100 F MPY 3,3 •
9 C ( K)*C(K)
002 5 6 1C7 AAA
FADD 2,3
00257 174240 FSQR 3,3
174020 FALG 3,3 9 LNCMAG.)
00261 0 2 0 0 1 5 $ FLDA 0 , TWNTY
0 0 2 6 2 114100 FMPY 0,3
•
9 20+LNCMAG.)
00263 A A A A «
FLDA 0 , L N l 0G
0 0 2 6 4 114200 9
FDIV 0,3 9 20*L0G(MAG. )
00265 064003$ FLD3 9
TEMPO 9 P O I N T E R TO GXY
00266 055400 FSTA 3 ,0,3 •
9 STORE THERE
00267 FIC3
rs a o 7
UV.»C # V
A
.» f A 7 A A A A C
U U V.' u u J) F ST3 TEMPO
00271 0 2 0 6 3 7 FLDA 0 , SUM 1
00272 024640 FLDA 1» SUM2
0 0 2 7 3 104200 FDIV 0 , 1
e
9 Q(K)/CCK)
0 0 2 7 4 124040 FATN 1# 1
0027 5 124400 FN EG Ì, 1 9 LAG NETWORK
0027 6 0 3 4 0 1 7 $ FLDA 3 , DEG
00277 164100 FMPY 3, 1 9 P H A S E IN D E G R E E S
A A QA A
V
.»u o u u
A
e*O
O h/!u u v .*_ A A A
FLDA 3# P 180 0
9 180 DEG. S H I F T
a a on 1 1 1
l»v.»Ou 1 1 1 A AA A
UUG FM0V o ,o ,f s l t ; QUAD. TERM NEG. ?
AA AA A
U U OU G 176400 F SUB 3# 3 •9 NO
A A AA A
u u Ju O lOr FADD 3* 1 •
9 YES
A A Q /» 0 6 4 0 0 4 $
A
FLD3 TEMPI
UU JU H 9 L O A D P O I N T E R TO QXY
A A AA C
V' V
.» J V 045400
,' J F STA 1,0, 3 9 STORE THERE
A A QA/
o u O t 'O 1I 1 I AA AA
O U U L» FIC3
V.' O V.»f
A A AA 7 A 7 A A A /, ff
U f J FST3 TEMPI
AAA I A | A A A A A
U V.' J 1 I FEXT
0031 1 A AAA A / ■(T
u g u u v .»0 3
> LDA 0, K
00312 AAA
v j u v .*D O u JMP GX3
00313 034601 GX6: LDA 3, R2C 9
9 ANSWER PRINTING
00314 031404 LDA 2,4 , 3 9 P O I N T E R TO GXY
00315 025405 LDA 1 5, 3
, 9 POINTER TO QXY
00316 A C A A A Q ff
O J O V . 'U u i1 STA 2, TEMPO
00317 044004$ STA 1,TEMPI
00 320 102400 SUB A A
CGC5 XSPEC
00337*000762 JMP GX7
0 0 3 4 0 * 0 0 6 0 1 1 SGX8: JSR 3.CRLF
00341*006011$ JSR §•CRLF
00342*006011 $ JSR e.CRLF
00343*020001$ LDA 0*C5 ; RESET RESULT
00344*040002$ STA 0 » C5P1
»
» FORMAT
00345*102400 SUB 0# 0
0 0 3 4 6 * 0 4 0 0 0 6 $GX9: STA C , K
00347*024005$ LDA 1#M
00350*036416 LDA 3,@PR20 •
> LOAD EXIT ADDRESS
00351 * 122415 SUB# 1# 0* SNR •
» FINISHED ?
V.H.*O J c yu 1 JMP 6 » 3
r'A Tt\'î» rtAAAA /.<r LDA 2 » TEMPI
00354*000326* F ENT
00355*021000 FLDA A A
SJ* KJ* d
A
L O A D QXY
00356*140001 FFDCF A
KJ
•
9 P R I N T IT
AA ATA f I A / ,A A A
F IC2
yuOOv.i 1o y u u y F EXT
00361*050004$ STA 2 »TEMPI
00362*006013$ JSR 0.RESLT
00363*020006$ LDA C*K
00364*101400 INC 0# 0
00365*000761 JMP GX9
0 0 3 6 6 * 0 0 0 1 1 4 *P R 20: R20
•Z R E L
0 0 0 0 0 - 1 4 1 2 6 4 P 180s -180.0
UVJUV.» 1
A A A A 1. A A A A A A
END
276
• TI TL FRE S P
.ENT FRESP
. EXTN FENT
. EX T P C 5 * C 5 P 1* T E M P O * T E M P I 9 D E L A M * . R E S L T * M
. EXTD P E G * L N L 0 G * T W N T Y * M 1ß P 1 8 0
. MREL
/> a #■> #
-\
V»V.»V* V
»' *0 54 5 6 4 FRESP: STA 3 * R 18 ß FREQUENCY R E S P O N S E
A A A A 1
V'V»»V.« V • 1 * 0 2 4 0 0 1 $ LPA 1*C5 ß S A V E R E S U L T F O R MA T
r\ r\ r \ g
*044564 STA 1* CSAV
A A A A A
c' o v.»v»O *024562 LPA 1*C3
A A A A
*044001$ STA 1 * C5
A A A A C
\
w 'V
.» V
.* J *044002$ STA 1* C 5 P 1
A A A A £ 0
*024561 LPA 1* C50 ß SET PLOTTING COUNTER
A A A A ?
*044013$ STA 1* M 1
A A A ! A
V.»V _
*V » 1V » * 0 2 5 4 0 1 LPA L L 3 ß OUTPUT POINTER G
A A A « 1
’ 031402 LDA 2* 2 9 3 ß OUTPUT POINTER PHI
AAA 1A
S^» V.'V-' 1 C ' 0 4 4 0 0 3 $ STA 1* TEMPO
A A A « A t A C A A A /. < r
v^'v'v*' I O V» KJ K.* *A I 1 STA 2 , TEMPI
I y, A A
v
.»v.' * 102400
1 *4 S UB SJ ß
AAA 1C
V
.* A0*040554
3 STA 0 *X •
ß SET FREQ. COUNTER
AAA « £
V-'V'V.f 1 C* ' 0 4 0 5 5 4 STA 0 *X + 1
AAA 1A
1 / *M i l l i TRI : F LNJT ß CAL. FREQ. RADS./SEC
A A A pA 0
*060551 FF L O X ß IN LOG. DECADES
A A A A « • A A A r r ,- N 0
U U V1' 1 V CZ. v.» Z3O FLPA Ci X ß IE. 0 W = 0 . 1E X P ( . 2 3 0 2 5 9 X )
AAA A A
\.'V' U C C '024555 F LDA 1*P23
A A A p ß
'104100 FKPY 0* 1
A A A p^j
* 124220 F EXP Ì 9 1
«/
"N
A A r^
V
.*V.' O 03 V.» CL D
3v.» FLDA 0 >P 1
A A
V»\»*V
A
O ^
» C. P '104100 FMPY 0* 1
.-Nr \ r \ r \ *7
v» v_' v' < 01/ * 0 4 4 5 5 2 FSIA LOW
AAA AA
v»*V.» O *0 2 0 5 4 3 F L FA 0 *P I 2
A AA Q 1
^'V-'V.'O I * 1 0 4 2 0 0 F FI V 0» 1 ß F R E G« IN HZ.
000 3 2 ’ 144001 F F PC F 1
A A A A A
v.» v» O O
t1 A A A A A
I V
.» V
.» V
_# F EX T
0 0 0 3 4 V^' t /* V»O
•A A /A0^A A
JSR P.RESLT
0 0 0 3 5 *0 2 2 5 2 7 LIA 0* PR18 ß POINTER 10 H
AAA A / 1
v_> ^*v 'O C *0 4 0 5 5 3 STA 0 , TP0
AAA O
v'v.a*vi / * 10 2 4 0 0 SUE
.
'J
KH V.'
ß
r\r\r\^
j. * 0 4 0 5 4 3
V'V.’ V.i
a
V
-‘ STA 0* J ß SET J=0
0 0 0 4 1 *0 4 0 5 4 3 STA 0 J+l
9
0 0 0 4 2 *0 4 0 5 4 3 STA 0 * S UM 1
0 0 0 4 3 *0 4 0 5 4 3 STA 0 * S UM 1 + 1
00044 * 040543 STA 0 * S UM2
00045 ’ 040543 STA 0 * S UM2 + 1
« A A A A 1 -7 9 ' T O O *
00046 V.* V.» V
.» 1 # i l \ CZ, • F EN T
A A A /, O
V
_*V.»V.* M f ’060534 F F LO U
A A A CA
l.' V .»J V * 0 2 4 5 3 1
*V F L DA 1* oy
000 5 1 *0 20 5 3 2 F L FA 0, J
9 A AA A A C < T
000 5 ? O v' 3 — F L DA 2* DEL AM'
A A A C O
V* V.» \^r ^3* ' * 1 0 4 100 FMPY 0* 1 ß J * W
00 1 50 07 4421 • FFIX X
9
00 15 1
4A A A A A
1 V
.» V
.» V* F EXT
279
003 FHE S P
00150* 0 0 6 0 0 6 $ JSR @ •RESLT
00153*010417 I sz X+l •
» IMG. 1c REG* COUNTER
00154*020416 LDA 0>X+1
00155*024412 L DA 1 , C5C •
» 4 DECADES ?
00156*106503 SUPL 0 * SNC
00157*000640 JMP TR1
00160*024406 L DA l , CSAV 9 RESET RESULT FORMAT
00161*0440015 STA 1,C5
00 1 6 2 * 0 3 4 4 0 2 LDA 3* R 18
00163*001403 JMP 3, 3 9 RETURN
r\
1
V'V.'
| s /, t
U H R 18:
r\ r\ r\ rv
V.'U
»A
A A 7/1
l
A
C
| A A A A A A
V.» V,» L ' V/
A
KJ
0 0 1 7 3 * 0 4 0 5 4 4 PI 2: 6.283185 ; two pi
00174*103754
0 0 1 7 5 * 0 4 0 0 3 1 PI J 0. 1
00 176* 1 14631
0 0 1 7 7 * 0 4 0 0 7 2 P23 5 0 . 2 3 0 259
00200* 17 1 100
AAAA
L* 1 V.» 0W:
1 | A A A A A A
V.' ^ \^f V»
r\
V.»
AA AA Q | A A A A A A A
U v» ¿Sv.» V U U V.» V.' V.'
A A Q/S O
V « V» c v ' O
f A A A A A A
V,»V.»V.» V . ' V ' V J: .1
A
V
A A A A /, • A A A A A A r\
V.'V_' ^ U V » V ' U V.7U Vw*
A A A A CT • A A A A A A
t>-' V.» O V.* V.» v/ «W
* SUM 1J A
V.'
AAAA/; • A A A A A A
v .»u V.« V,» V_* L* V.'
A A AA 7
/
« A A AAA A
v . ' * w * v.» »^* v.* SUMP: A
V'
AAA 1 A • A A A A A A ,A
W.' V» ¿T- 1 V.' V-» Vj V.'
A A A I
V' V.' c . 1 1
) f A A A A A A
V.» V-» TP0: A
EM t
280
The principal algorithm used was that of Fletcher and Powell and
follows
to G|<(Xi , X 2 , •• X|\|) = 0
where k = 1, 2, 3 ... M
N M N
\ 9Gj 9Gk _ 3Gk 3F
E E E
A J3xj 3x; " 3xj 3x j
¡-1 j=1
for k = 1, 2, 3 .. . M
achieved and the programme package was tested on some problems whose
SLMrO.
Z-YCOR**:
DC 20 I =1»N
SUM=SUH«XII)*Z/CI«1)
20 Z -Z* YCOR
T R U C z CXCOR-S UM 1 • t N V M l / Z
X I NV Ï -TRL' C
W R I T E C l W R I T E * 1 2 0 ) t 1 1 X C l ) » 1 = 1 # NVJ
I SWI CH - 0
C
CALL SCALER C N » H « X 1 1 . D - l 2 1 1 . D6 • 5 . D - 6 1 1 . D 2 » G0 » Gn I
C
c
W R I T E ( I W R I T E * 1 2 0 ) C I • X C l ) » 1 =1 « NV)
C A L L L I N K C 1 »X »F tC RADF .G t G R A D G )
C
CA L L CONK I N CL I NK t N * M» X . F t GR ADF t G t GRAD G 1 1 .D 2 * 1 . D 1 tH tF L A Mt DE L X* E P S •
1 DE F X * D E L G , 2 0 » 6 1 » 0 » I E R » S » G 0 » G N « W » A )
C
W R I T E C I W R I T E »801
I S W I CH -1
WRITE Cl W RI TE »110)
C A L L L I N K C 1 »X tF »G RADF tG » G R A D 6 I
VAL-YRCm
C
W R I T E C I W R I T E «1201 CI»XXCI)»I-1*NV)
WRITE C I W R I T E t l O O ) V I N I T t VAL t T F
I S Wl C H -0
G C TC 10
C
30 S T OP
C
40 FORMAT C * 1 T Y P E I N THE F O L L O W I N G DATA S E P A R A T E D BY A C O M M A * » / * N V ’ .
1 *» M t M C . H S » V I N I T » M U . T A U t XCOR t Y C O R V / )
50 F CRMAT C )
60 FORMAT C* T Y P E I N THE E S T I M A T E S OF THE N C O E F F I C I E N T S * t / * X U ) * . * #
I X 1 2 !...ETC*//)
70 F 0 RMA T t * * t / / / R X » * I Ni T I A L V A L U E S * )
80 F ORMAT C * 1 * » q X t *F I N AL V A L U E S * )
30 FORMAT C / / » * NV - * t I 2 t 5 X f * M - * » I 2 » 5 X » * M C z * . I 2 t 5 X f * H B z * t F 8 - 5 . * VI
l N l T z * t F 8 . 5 » 5 X . *MUz * » F 8 . 5 t 5 X t •T A U= * t F 8 . 5 t / * C l z * »F 1 1 . 8 . E X • * C2
2 - * .F 1 1 . 3» 5X» *XCOR- * F 8 . 5 . 5X t * Y C 0R= * t F 8. 5 . 5 X • *T HE TA F - * t F 1 0 « 3 / / )
100 FORMAT <//* I N I T I A L V E L O C I T Y z * , F l C U 5 / / * M I N . D E S C E N T TI ME z « , l p
1 E 1 5 . 3 / / * THETA F z * » 0 P F 1 0 . 2 )
110 F CRMAT ( / / t 3 X t *Y *t 15X t *W * t 15Xt * X * # 1 5 X *T * t 13X» *T HE TA * » 1 3 X t * V E L * t l 2 X
1 »*PAR.X*t8Xt*D IF.IN X *)
120 FORMAT ( / / * C O E F F I C I E N T S A RE * t q f 2 X • *X t • 1 1 1 » • ) z * , 1 PE 1 5 . 8 ) )
E ND
* 284
S L B R PUT I fi E C CMKI N i F N S # N * M * X * F t G R A D F » G * G R A D G ? B * R * H * F L A M f D E L X * E P S » D
1EF X » D E L G t l T » I P R . l N D * I E R t S t G O . G N * W » A 1
M A I N O P T I M I S A T I O N S U B R O U T I N E : P E R F O R MS OR C O O R D I N A T E S ALL CALCULAT
I C N S TO D E T E R M I N E THE S C A L E D C O E F F I C I E N T S OF THE
INDEPENDENT V A R I A B L E S .
I M P L I C I T DOU B L E P R E C I S I O N ( A - H t O - Z I
D I M E N S I O N X C 2 5 ) * G R A D F ( 2 5 Ï t G ( 2 5 1 t G R A D G C 2 5 ) t H ( 5 0 ) » F L A M ( 2 5 ) t St 2
151 t GO f 25 ) » GN ( 2 5 1 1 W t 5 0 ) » A I 2 5 1
CCMMON l R E A D » I W R I T E t I S W I C H t N N t M M * X S C C l O O ) » F S C » G S C t l O O Î t D E L D I F t 1 0 0 )
IF fIN D -1 ) 10t20t20
10 C » L L M A I ’J P l H « 1 « N t S * GO» G Nt A LF A » D E L X » DE LG • W )
20 IP-IPR/10
I R=I P P - 1 C *1 P
I F I I P ) 4 0 t <i 0 « 3 0
30 W R I T E t I W R I TE 15 0 0 ) N , M »1 T • I P R» I ND »B t R t DE LX t E P S » OE F X »DE LG
W R I T E t l W R I T E *57 01 t I . X 11 ) » 1 =1 • N )
I F 1 1 . G T . M ) GO TO GO
DO 50 U l . M
50 F L A M t l ) =0 . 0
GO C O N T I N U E
70 C A L L LAM3 t N * M • G R ADF * SR ADG * G» S t F LA Mt RtD E LG » W t A »1 Cv )
I E L - I E L ♦1
L =1
80 P H I = F
I F l l . G T . M ) GO TO 100
DO 90 I =1 tM
90 P 3 I = P H I * G t I I * I B * G f I l - F L A M ! I ) )
100 I C K C G T = L
I - Ü C K C 3 T . L E . 1 ) GO TO 110
I F I T C K C G T . G E . 2 ) GO TO 310
03 TO 1110 » 3 1 0 ) • I CKCGT
1 1 0 . G M = 0.
285
I F I l . G T . N ) GO TO 150
DO m o 1 =1 , N
T=GRADF I I J
K=I
I F | l . G T . M ) GO TO 130
DO 120 J =1 . M
T =T* GRAD G ( K I *1 2 . * B « C ( J l - F L A M I J l l
120 K =K t N
130 G N I I ) =T
1 C| 0 G N A =GNA « T * T
150 GNA=SGRT!GNA)
GS - 0«
IF » l . G T . M ) GO TO 170
DO 160 1=1 •M
160 G S=GS«B « G I I I * G I I I
170 I C K CG T =L
IF I I C K C G T . L E . I I GO TO 180
IF I I C K C G T . G E . 2 I GO TO 350
GO TO 1 1 8 0 . 3 5 0 1 . I C K C G T
C
C PRINT OUT DATA FOR LAMBDA DETERMINATION
C .
180 I F I I P » 200 » 2 0 0 * 1 9 0
1 9 0 W R I T E 11 W R I T E »5201 I EL » I E T * I C V
W R I T E ( I WRI TE »5701 ( 1 » X( I I »1 =1 «N I
W R I T E ( I W R I T E « 5 1 0 ) < F L A Mi l l » 1 = 1 « Ml
W R I T E I I W R I T E * 530 I ( 6 1 1 1 . 1 = 1 . Ml
W R I T E f I W R I T E t 54 0 ) P H I « F « G S » G N A
W R I T E f I W R I T E *5701 11 « XSC11 I * 1 =1 «N I
C .
C D E T E R M I N E WHETHER CONV E R GE NC E HAS B E E N OB T A I N E D
C
200 I F I I C V - 1 ) 2 1 0 . 2 2 0 »230
210 I F I T N L ) 2 3 0 * 2 2 0 . 2 3 0
220 I E R= 0
GO TO 450
C
C LINEAR ITERATION
C
2 3 0 L =2
I NL = 1
I T N=0
I F I l . G T . N ) GO TO 250
00 240 1 = 1 «N
240 G O U I =GNI I I
250 CONTINUE
260 I F I I E T —I T J 2 8 0 . 2 7 0 * 2 7 0
2 7 0 I E R =1
GO TO 4 5 0
230 I F l l . G T . N ) GO TO 300
DC 2 90 1 = 1 . N
290 Wl I ) = x m
300 I E T = I E T * 1
ITN=TTN«1
I N L F = I NL
T NL =0
A LF A =1. 0
286
310 C fi LL LI NUN ( I N L . N . K F •W» S * AL F A * X » PH I . D E L X * 0 E F X . E P S )
l- 1 T N L - 1 I 34 0 . 330* 320
370 IER r?
GO TO 450
330 C A L L F N S t 1 • X t F . G R ADF * G » G R A O G )
G 3 TT 80
c
C 03TATN GN.GNA* AND GS
C
340 CA L L F NS ( 2 . X * F • G R A D F . G . GRADG)
G D TC 11C
C
C FIND INDEX INL FOR LINEAR ITERATIONS
C
350 INL=INL«2
IF (ALFA) 370.370.360
360 I F IGKA-DELG! 370.370.330
370 I N L -0
C
C PRINT OUT DATA FOR LINEAR ITERATION
C
330 IF (IP*IR ) 41 0 * i l l 0* 390
390 IF I M O D ( I E T t l R M 410.400•410
400 W R I T E ( I W R I T E . 5 5 0 ) I E T . I N L . A LF A , NF
W R I T E ( I W R I T E . 5701 ( I . X ( I I . I =1. N )
W R I T E ( I W R I T E . 5401 P H I . F . G S . G N A
C
C U P DAT E H AND 03TAIN NEW S
C
410 CALL MATUP t H . 3 . N •S t G Ot G N. A LF A » D EL X . DE L G . W J
C
C C O NV E R GE NC E QBTAINEO IF INLP=0 A ND I NL ” 0 ( I T N = 2 >
C
1= I I N L F * I N L ) 420.420.430
420 IER=0
GO TO 450
C
C 03TAIN NEW SET OF LAMOA VALUES OR CARRY OUT A FURTHER ITERATION
C
430 IF (IT N -1 ) 260.2S0.440
440 IF I I N L • ( K-M«1-ITN1f 70t70»260
C
C CCNCLUSICN
C
450 IF ( I P ) 490.490.460
460 W R I T E ( I W R I T E . 5 G 0 ) 1ER
K=1
1= t l . G T . N ) GO TO 480
D C 4 70 I - 1 . N
Wt I î - H ( K 1
470 K = K * N « 1 - I
4 3 0 C ONT I NUE
WRITE ( I W R I T E . 5 8 0 ) tW tI!«I= l*N >
490 P E T U R N
C
500 FORMAT ( / / ' I N I T I A L DATA FOR C O N S T R A I N E D MI N i M l Z A T 1 0 Nf » //• N = * . I 3 .
1 3 X . , P =* t I 3 . 3 X . » I T r * . I 4 . 3 X » t l P R = t . I 2 . 3 X . , l N D = » . I l / » D z ' . l p E l O ^ . E x *
28^
Y A RR A Y ( ! ♦1 ! = -Y R ( 1 )
50 C ONTINUE
IFn sw iC H .C G .il call °LT c XARPAYf YARRA Y J
FüNFzYRfdl
RE TURN
C
60 F CRVAT ( 1X.1P8E16.8 1
C
E ND
SL'BP.CUTIhE FUNC (XtFtG)
I M P L I C I T DOUB L E P R E C I S I O N t A - H t O - Z l
DIMENSION X (2 5 ) t G(25 )
COMMON / C H U T E / NV » N » NSTE P »MC t XC OR t YCOR t T A L* t V I NI T * S T E P . MUt C I » C2 t YR t
1 9) t X X ! 251• TF
M I N I I F t 0 . DO 1
FrFUNF(X}-1.D20*U
CONSTRAINT EQUATIONS
THET-FTHETACYRC1 I I
THET-ATAMTHET )
G IDzTHET-TFA
RETURN
291
F fNCTÏQN F THE TA CY I
C
C EVALUATES THE F U N C T I C N FOR T HE TA AS A FUNCTION OF Y
C
I H P L I C I T DOUBLE P R E C I S I O N C A - H t Q - 2 )
CQMMON / : H U T E / N V . N f NST E P t HCt XCO Rt Y C O R» T A U t V I Nl T • S TE P • MU • C l * C 2# Y RI
I9Î.XX(25)tTF
S-O.
A X-Y
D3 10 I - l t N V
S =S * XX( I l * A X
10 A X - A X * Y
F THE TA - S
RETURN
C
END
292
SUBROUTINE DS l MQ ( A • B1 • 32 . N. KS I
C
C STIVES A SYSTEM OF SIMULTANEOUS LINEAR EQUATIONS
C
I M P L I C I T DO U B L E P R E C I S I O N I A - H . O - Z )
D I ME NS I C ^ A ( 2 5 ) . B 1 C 5 0 I * B 2 C 5 0 »
C
C F ORWARD S O L U T I O N
C
T C L - l . D -20
KS z 0
J*Jz-N
I F I I . G T - N 1 GO TO 120
DC 110 J z 1» N
JY z J* 1
Jv'z J 1
RIGAzO
I TzJJ-J
IF I J . G T . N I GO TO 30
DC 20 I z J , N
C
C SEARCH F OR MAXI MUM C O E F I C I E N T IN COLUMN
C
I J z l T ♦I
I c IA BSt3 IG A 1 -A BSlA (U m 10.20.20
10 B I G A z A ( I 1
IMAXzI
20 C O N T I N U E
30 C O N T I N U E
C
C TEST : 0 R PIVOT LESS THAN T O L E R A N C E (SINGULAR MATRIX!
C
I F ( A3 S t 91 GA l - T O L 1 <40.«40.50
«10 K S z l
RETURN
C
C INTERCHANGE ROWS IF NECESSARY
C
50 I I z J * N* t J - 2)
I T z I MAX-2
IF I J . G T . N J GO TO 70
DC GO K z J *N
I I z T 1 ♦N
I2 z Il.IT
SAVEzA i l l I
A I I I J z ( 121
AII2JZSAVE
C
C DIVIDE EQUATION 9Y LEADING COEFFICIENT
C
GO A! 1 1 1 z A C I I l / B I G A
70 SAVE z B I C I MAX 1
B1 I I M A X J z B l I J )
B 3IJ>=SA VE/ BIGA
S A V E z B 2 ! I MAX 1
0 211 rA X ) - B 2 ( J )
B 2 I J 1 - S A V E /BIGA
C
C ELIMINATE NE X T VARIABLE
C
I F ! J - Nl 8 0 * 1 3 0 * 8 0
80 I G S = N * < J - l J
00 13 0 I X = J Y * N
1 XJ - I QS« I X
IT= J-IX
I F t J Y . C T . N I G O T O 100
DO 90 J X r J Y . N
I X J X - N <I ¿ X - 1 M I X
JJX - IX JX*IT
SO A I I X J X ! = A ( I X J X ) - f A ( I X j > « A ( J J X ) >
1 00 B l I I X I = B 1 ( I X I - l Q l ( J ) « A l I X J l )
110 B 2 I I X > = B 2 ( Z X I - C B 2 ( J H A ( I X J ) I
120 C O N T I N U E
C
c BACK 5 OL UT I ON
c
130 IF IN-1) 1 7 0 . 1 7 0 » m o
140 NY -N -1
IT iN*N
IF 1 l . G T . K Y I GO TO 160
00 150 J =1 •NY
I A =1 T - J
1 3 -N- J
I C -N
00 150 K z 1 * J
BllIB)= £ l(IB!-A fIA > «B l(T C )
B2 I 19 I —3 2 C I B I - A f I A J * B 2 ( I C I
I A -I A-N
150 I C - I C - 1
160 C C N T I N U F
170 R E T U R N
C
END
294
SUBROUTINE MDW I XP • DY t Y 1
C
C E V A L U A T E S DW/OY t T A N ( T H E T A I t I C O S « THETA 1•V E L ) • •-1
C R E Q U I R E D FOR F I N D I N G Y C O R r W t X C G R t T
C I-E. D Y t l l t DY12 ) t D Y l 3)
C
R E A L MU
I M P L I C I T DOUBL E P R E C I S I O N l A - H t O - Z l
COMMON I R E A D t l W P l T E . I S W l C H t N N t M M . X S C I l O O ) f F s C t G S C l I O O ) f D E L D I F ( l O O )
COMMON / C H U T E / NV t N » NSTE Pt MC t XCOR r YCOR t TAUf V I Nl T* S TE P t MU» C 1 * C2 * YR <
1 9) t X X ( 2 5 ) * TF
DI MT N S I OF X P l M C 1 t D Y t N V 1
THET-FTHETAtY )
THET-ATAMTHET )
C S - C O S t THE T 1
SS-SINtTHET»
V E L 2 - V I N I T • •2 •2 • * ( 9 « 8 1 * Y « X P C 1 I I
C
C V - l V l N I T * * 2 * 2 t G Y « V I I I • • 0 •5
C
VEL-SQRTt ABSt VEL21I
s=o.
Z= 1 •
DC 10 I r l . N V
S - S * X X I I 1* 2
10 Z-Z< Y
TE R M - 9 . 8 1 « S S - V E L 2 * t C S * * 3 ) * S
TE RM-AB S t T E R M ) « T A U * t 1 . « C 1 / V E L « C 2 * V E L )
C
C C A L C U L A T E DW/DY
C T E RM = T A U ( l 4 C l / V « C 2 * V ) t ! G * S l N ( T H E T A ) - V * * 2 * C 0 S * * 3 ( T H E T A ) ! )
C DW/DY z - ( MU•V « T E R M » / C O S t T H E T A )
C
D Y t 1 ) - - t T E R M« MU * V E L ) / C S
DY t 2 ) -TAMt THF T )
D YI 3 >—1• / ( C S * V F L )
RETURN
C
END
295
SGIrW(3 1
IF C1 . GT• M ) GO TO 100
D C 90 J r 1
SSIrSGI-FLAMt J M A ( K )
90 K zK « 1
100 S G I r A B S f S G I l
AGIrSQRTfAC L I )
Q-Q*A C L )
LrL#M«l
T r A G I t DEL Q
TrSGI/T/DELG
1 10 R D z A M A X l t R D . T )
C
C COMPUTE S L I I I AND NEW V A L U E S FOR F L A M t I I
C - D r M A G N I T U D E OF CHANGE I N FLAM
C
120 Kzl
Q- Q *Z
I F f l . G T . M I GO TO l d O
DO 130 I r l i H
A CK ? r A ( K )« Q
SH Ilr-G V tll
130 K=K« K#1
1 «10 C O N T I N U E
CALL D S I F G ( AtW»SL t M t I E R )
IF C I E R ) 160*160*250
1 5 0 Z r 10. * Z * 1 . 0 - 1 2
GO TO 10
160 D r O . 0
IF f 1 . G T . M I GO TO 180
0 C 170 I = 1 * M
Zr Ul I J - FLAM! I I
FLA M fl)rw (I)
W CD rSLCIl
170 D rD ♦Z * Z
180 DrS QRT( D1
IF f l . G T . N ) GO TO 200
DO 190 I r 1 * N
KrI
s m iro .o
DC 190 J r 1 * M
SL I D r S L l I K U ( J I * GMt K )
190 K: K«N
C
C MODI FY SL C I ) I F NE C E S S A R Y
C - W r N O R M A L l S S D G R A D I E N T OF SUM OF G V ( I ) * » 2 1 I r l TO M)
C -SWrW»SL
C - S P r M A G N l T U D E OF COMPONENT OS SL P E R P E N D I C U L A R TO W
C - R r M A X I F U M ALLOWED VALUE CF S P / A 3 S C S W )
C
200 TrO. O
I F I 1 . GT. N ) GO TO 2<l0
DC 230 I : l , N
WI I )= 0.0
KrI
IF C1 . GT. M ) GO TO 220
DC 210 J -1 * M
298
W IIlzW d M G V lJM G M IK )
210 KzK*U
220 Z= W «I)♦ SlG N tl.D -20.U (Iil
WI I > - I
230 T =T * 2 • Z
240 T-SQRTCT)
SW = 0 . 0
I F t 1 •G T • K I 60 TG 260
DO 250 I - 1 t N
W I I l =WI I 1f l
250 SW=SW» W H ) » S L I I J
260 SW=SW*SÏGN( 1 . D - 2 0 . S W )
S P -0.0
I F I 1 . 6 T . K l GO TO 280
DO 270 1 = 1 tN
Z =S L 11 ) - S W * W C I I
SLITI=Z
270 SP=SP«Z*Z
280 SP=S3RTISP 1*1.0-20
T = A M I N 1 ( P * A B S ( S W I / S P » 1. »
IF 1 1 . G T . N I GO TO 300
D C 2 90 I = 1 . N
290 SL II1 = T *SL( I)«SW *W <I1
300 CONTI MUE
C
C CONVERGE »CE p A RAME TER
C
IF I I C V I 310*330*330
31 0 I I 00=0
D C 320 I =1*4
320 H I S I I I = 1 . D20
330 I C V =2
T = H I S < 11
DC 340 1 = 1 . 3
T =T * H Ï S C H 1 )
340 H I S I I I = H I S I 1* 1 I
HI S I 4 1=D
T=T/4.
IF IR D -1 0 0 .1 360*360*350
350 I 100=0
RE TURN
360 I 1 0 0 = I 1 0 0 « 1
I F I R D - 1 . 1 3 7 0 . 3 7 0 * 389
370 I C V = 0
R : TURN
380 I F ( 1 1 0 0 - 5 ) 4 1 0 . 3 9 0 . 3 9 0
390 I r 1 0 - T ) 4 1 0 * 4 0 0 * 4 0 0
400 I C V=1
410 R E T U R N
C
END
299
DC ¿4 10 I =1» NO
' n o E D I I ) =3 E G I N ( I ) « A LF A • S ( I )
'420 C O N T I N U E
C
C TEST WHETHER IMPROVEMENT IS SIGNIFICANT
C
I F f F K E E P - F N - E P S I *130 « *1*10 » *1*10
*130 I N D r - 2
RETURN
*4*40 1= I A B S l AL F A 1 - DE L I *130 • *450» *150
*150 I NDr -1
RETURN
C
END
StBPOUTIFE LINK C I F •X * F •G R A O F * 6 t G R A D G J
C
C APPLIES SCALING FACTORS TO F U N C T I O N AND INDEPENDENT VARIABLES
C
I M P L I C I T DOUBLE P R E C I S I O N C A - H . O - Z 1
D I M E N S I O N X I 2 51 . GR ADF I 25 11 -G 125 I * G R A D G I 2 5 1 * X R A W l l O O l
COMMON I R E A D t l W R U E t I S W I C H t N N » M M * X S C ( 1 0 0 J »F SC .G SC ( 1 00 J t DE L O I F ( 1 0 0 )
N=NN
K=HH
C
I F f l . G T . N ) GO TO 20
D3 10 1 = 1 . N
10 X R AWCI ) = m M X S C C I I
20 CONTINUE
IF IMJ EC.G0.30
30 C A L L GRADF G f I F « X R AW. c * G R A O F • G . G R A D G 1
I F I l . G T . M l GO TO 50
00 *10 J=1 v M
40 G l Jl=G( Jl/GSCC J )
50 CONTINUE
GO T P 70
SO CA L L GRADF G « I F t X R A W . F v G R A D F « G t G R A D G 1
70 F =F /F SC
I c IIF-11 140.140.80
80 IF f l . G T . K 1 GO TO 130
DO 120 1 = 1 . N
G R A D F ( I 1 =GRADF C I 1 / F S C * X S C f I 1
IF I M 1 120.1 2 0 .9 0
90 K=I
I" 11 . GT. M ) GC TO 110
DC 100 J =1 . M
GRADG ( K 1 = G R A D G I K I / G S C I J M X S C I I I
100 K =K * N
110 CONTINUE
120 CONTINUE
1 30 CONTINUE
140 RETURN
C
E NO
SUBROUTINE MA TUP I H* T ND * N*S » G Q • G N* A L F A * D î L X • DE LG * W I
C
C COMP UT E S SEARCH DIRECTIONS
C
I M P L I C I T DOUBLE P R E C I S I O N I A - H . O - Z J
D I HE N S I OK HC 50 11 S I 2 5 I * G Q C 2 5 I * G N Î 2 5 J * WC50)
GE N E R A T E I D E N T I T Y MATRIX
o n
icn k =k ^ : - j
Gn TC l a o
170 K =K« 1
130 C 3 NTT NUE
190 CONTINUE
200 W II)=T
C
C -COMP UTE POSITION CHANGE Z AND S T O R E IN G ( N* 1 ) TO Gl2*N)
C
210 K=N
IF l l . G T . K ) GO TO 230
00 220 1 =1 »N
K =K» 1
22 0 Wl K 1=ALF A * S t I 1
C
C -COMPUTE Q-YtH.Yt P =Y . Z AND R=SQRTI7*Z)
C
230 R-0.0
Q-0. 0
P =0 . 0
K -N
IF l l . G T . N ) GO TO 250
DO 240 I =1 . N
K- K ♦1
R=R« WCK ) * W( K )
Q= Q «G N (I ) * W( I )
240 P = P « G N i I 1 *WC K )
250 R - S Q R T I R )
C
C -LEAVE H UNCHANGED IF P OR R ARE VERY SMALL
C
IF 1 R- DE L X ! 3 0 0 . 3 0 0 . 2 E 0
260 IF I P - R ♦DE LG J 300* 3 0 0 . 2 7 0
C
C - C C MP U T F NEW H
C
270 K:l
I e « 1 . GT. N ) GO TO 290
DO 280 1 =1. N
L I =N* I
DC 280 J =1 . N
L J =N» J
H t K J =HC K 1 * W ( L I ) * Wt L J l / P - K C I > « W t J * / Q
230 K= K ♦1
290 CONTINUE
C
C COMPUTE NEW D I R E C T I C N
C
300 I F ( I ND -4 ) 3 1 0 . 3 8 0 . 3 1 0
310 I F I 1 « GT• N ) GO TO 370
DC 3 E 0 1 =1 . N
T =0 • 0
K=I
I F f l . G T . N ) GO TO 350
DC 340 J =1 * N
T=T*G0(J1«H(KI
IF I J - I » 320.330.330
306
3 20 Kr K♦ N- J
GC TP 3 0
330 K-KM
340 CONTINUE
350 CONTINUE
360 SII*r-T
370 CONTINUE
380 RETURN
E ND
SUBROUTINE PLT (XARRA YtY AR RA YJ
I M P L I C I T OOUBLE P R E C I S I O N ( A - H . C - Z )
REAL XXRAY.YYRAY
C
C SUBROUTINE TO P L O T RESULTS
C
COMMON I R E AD» I U R I T E » I S W I C H . N N t M M f X S C l l O O ) f F s C f G S C t l O O ) » D E L D I r i l O O )
CCKMON / C H U T E / NV » N r N S T E P » MC »XCOR t Y C O R » T A U » V I N l T » S T E P » M U * C 1» C 2 • Y R i
1 9 ) • X X 1 251 * TF
0 I HE K S I O K X A R R A Y I 5 0 ) . YARR AY f 50 I * X X p A Y ( 5 0 ) t YY RA Y i 50 I
1= I I N T C Q . G T • 0) GO TO 10
CALL START
NSTFPlzNSTEP«!
XXRA Y l N S T E P 1 * I ) z 0
XX RA Y f NST E P l * 2 l = «8<l6 6
Y Y R A Y f N S T E P 1 « 1 1=0
YYRAY<NSTEP1«2I=.8*|66
10 I K 7 E Q=I NTE Q* 1
DO 20 I -1» N S T E P 1
XXRA Y ( I ) =SNGLC X A R R A Y ( I I )
20 Y Y R A Y t l ) - S N G L l Y A R R A Y ( I ) « 5 « )
CALL L I N E ( X X R A Y » Y Y R A Y * N S T E P 1 » 1 » I N T E Q 1
I F < l M T E 3 . E Q . i l CALL AX E N 1 5 . 9 05 5» . 59 OS • 1 . 1 8 1 2 *2 t 0. . 1 . 1 1 5»
1»X A X I S ( HE TRE S ) * » 5 . 9 0 55 * . 590 6 » 1 . 1 8 1 2 » 2 » - 5 . » I •»15 t
2 *Y A X I S ( M E T R E S ) t t - . 0 5 » - . l l
RETURN
C
E NO
308
I M P L I C I T DOUBLE P R E C I S I O N C A - H . C - Z J
D I M E N S I O N DY t 3 1 t Y R I M C I . X T f l O O ) . Y M i l O O Î . Y K t l O O Î
D I M E N S I O N X P C I O O Ï » Z K l l O O Î t WK CI O O ) * X K t l O O î
NQ - M C - 1
TWrYRi1 J
S T E P H - O . 5•STE P
T ? r T W ♦S TEPH
C
C EVALUATE RUNGE KUTTA COEFFICIENTS Z Kt WK * XK » YK
C
DC 10 I r l , K Q
10 XTIT)= YRII«1]
DO 20 I = 1 »NQ
20 XPlT)= XTfI 1
CALL F 4 DW t X P . D Y . T W I
DO 30 I Mr 1 »NQ
30 ZKIIM IrO Y(IM )
DO 4 0 I -1 » N9
40 X P tl I r X T C I M STEPH*ZK( T )
C A L L F4DW t X P » D Y * T Z I
DC SO I r l »NQ
50 WKI I J r D Y t I I
DC S 0 I r l »NQ
GO X P i I l r X T I I MSTEPH*WK!I1
C A L L F4DW I X P » D Y » T Z 1
DD 70 I r l » NQ
70 X K 11 î r D Yf I î
DO 80 I r l » NQ
80 XFtI > rXTfIM ST EP*XK fI 1
C A L L F 4DW CXP »D Y »TW« S T E P J
OC 90 I r 1 *NQ
90 YK! T 1 r DY ( I I
C
C SET NEW V A L U E S F OR W t X C O R » T T H E N Y CH R
C
DO 100 1 = 1 »NQ
T E R M r Z K ( I Î « 2 . * W K ( I )« 2 . « X K ( I ) *YKCI ï
100 YM IIîrXTCIî«STEP*TERM /6.
Y R 11 J r T W« STE P
DO 110 I r 1 »NQ
110 Y R 11 «11 r YMt I >
RETURN
C
END
309
c
I M P L I C I T DO U B L E P R E C I S I O N f A - H t Q- 2 1
D I F E N S I O A X I 25 1 • G R A D G 1 C 2 5 ) * G R A 0 G 2 I 2 5 ) * G I l O O ) » G p U O O i * G M l O O l ,
1 3RAD F1I100)* GRADF2I100)* ABGDllOO)
C OMMON I READ » I W R I T E t I S V J I C H » NNt MHt XSC ( 1 00 ) *F SC *G SC Cl 00 ) *DE L D I F ( 100 I
C
NN=N
MM=M
C
c FIND FIRST AND SECONO DERIVATIVES
C
DELIAS8RTC A B S ID E L II
DE L3 =DE L • 1. / 3 . )
FSC rl.0
IF I l . G T . N ) GO TO 20
DO 10 1 =1 t N
10 X S C I I ) =1 . 0
20 CONTINUE
IF I K ) 6 C *6 0 * 3 0
30 IF l l . G T . M ) GO TO 50
DC 1«0 J =1 *M
qo G scu izi.o
50 CONTINUE
60 C A L L L I N K 11 » X •F •GRAD=1 * G*GRADG1 1
IF l l . S T . N ) GO TO 110
DO 100 I = 1 tN
X I I ) =Xt I l « D E L 3
CALL L I N K t l * X . F P * G R A D F l * G P * G R A D G l )
X I I >=XCI ) - 2 . * D E L 3
C A L L L I N K I 1 # X t F M t G R A D F 1 t GMt GRADG11
X I I >=X( I ) « D E L 3
GRADF1 I I ) = I F P - F M ) / 2 . / D E L 3
G R A D F 21 I 1=1F P « F M - 2 . *F ) / D E l 3 / D E L 3
IF I M 1 100*100*70
70 K:I
I F t l . G T . M ) GO TO 90
DO 80 J : l » M
G R A O S H K l =«GPt J I - G M I J1 1-/2./DEL3
G R A D C 2 I K ) -C GP I J M G M C J ) - 2 « * G ( J ) ) / D E L 3 / D E L 3
80 K=K«N N '
90 CONTINUE
100 CONTINUE
110 CONTINUE
C
C FIND F 5C AND F I R S T A P P R C X I M A T I ON TO XSC (SCALE FACTORS)
r
F SC =ADS l F J
IF IM) 260*260*120
120 K =1
I F l l . G T . H ) GO TO 160
DC I S O J =1 * M
s =0 . o
I F ( l . G T . K l GO TO 140
03 130 1 = 1 *N
S=S*GRADG1< K ) * * 2
1 3 0 X= K ♦1
140 C O N T I N U E
150 A 3 G D I J l z S Q R T I S l
I GO A N=N
R N = 1 0 •• S3 R T t ANI
IF t l . G T . K ) GO TO 210
Ou 200 1 = 1 tN
F S1=ABSCGRADFIIII/F »
FS2=SQRTt A B S C G R A D F 2 C I l / F I )
S =0. 0
K=I
IF t l . G T . M GO TO 190
DO 1 SO J =1 t N
G F A T = A B S ( G R A D G 1 I K I I / ABGD( J I • TOP*RN/BOT
IF IGRAT-1.1 1 8 0 t l 8 0 t l 7 0
170 T = A 9 S ( G F A D G 2 C K l / G R A D G 1 C K I I
S=ANAX1IS*T)
180 K =K♦N
1 90 S =A N A X H F S 1 # F S 2 * S * 1 . / T 0 P )
S =A NI N I t S t l . / B O T I
2 0 0 x s c n 1 = 1 . /s
FIND G SC <G S C A L E )
210 K=1
I F I l . G T . M ) GO TO 250
DO 24 0 J =1 * M
S =0. 0
IF ( 1 . G T . N I GO TO 230
DO 220 1 =1 * N
T =G RADG1 ( K I * X S C <I I
S=S« T*T
220 K= K ♦1
230 CONTINUE
240 GSCIJ1=S3RTfS I
250 CONTINUF
260 I F 1 1 . G T . N I GO TO 320
DO 310 1 = 1 *N
S =0 . 0
ic tN) 300*300*270
270 K=l
I e t l . G T . M ) GO TO 290
DO 280 J = 1 * M
TrGRADG1( K ) / G S C ( J )
S =S ♦T * T
280 K= K ♦N
293 S =2. *B * S
300 S=S*CRA0F2CIl/FSC
S=ANAX1( S t l . / T O P / T O P )
S=ANAXltSQRTtl./S)*BOT)
XSC II)= S
31 1
DE L DI F ( I l = X S C 11 1« D E L I
3 10 X I I ) : X ( I l/XSCCI )
3?0 CONTINUE
RF TURN
E ND
312
there arises the need to compute the step response of the system.
included here.
vdv______
= f-
(E - £ - Fv ) (Ev - c - Fv2 )
v v
vdv
-F
(v2 - I v + £ )
-i r(2v - ! + & dv
2F , 2 E v .c
(V -
F + F}
313
9 9
2 2E v ç E_. ,c__E_
(completing the square)
v - 2F F 2F + F 4F2
0 E
2V - F +
f-
2F v 2 - —E v ,fc— 2 Ev c
F F - T F
1 I / 2 E v , cN E r dv
Ô F ln (v - -p + — ) - f (4.5)
2F F F 2F2 2 E v .c
V F + F
dv
2F‘ 2 2E v c
V - 2F + F
Let (jü = v -
2F
E _ c_
a2 -
4F2 F
dv
/-.
2F2 J E_ 2 _ ,_Ef_
(v ‘ 2F 4F2
Upon substitution
E f_ dv
2F2 w2 - a2
1
E < '
2F2 2a a) - a u) + a
[ In (a) — a ) - In (a) + a )]
4F2a
314
E v2
’1n (oo - a ) = AT
2F ln (V? " f V + f 4F2a
vi
Ev £ gv2
1 E 1n (oo - a )
In (v2- -f v + ~) -
2F F F 4F2a b (oo + a)
Vl
2 Ev . c. E ^1 n (oo + a )
2 p In (v - -p + — ) +
4F2a (oo - a)
Vl
1
2 p"
I
In
, 2 Ev
(v — ~
c.
+ —) + —
E In (a) + a )- v2
F 4F2a (oo - a)
"v 1
v2
,n (V 2 - | v + § t F + In
F F (co - a)
Vl
1v 2
co + a 4F2a 1
In
oo - a (v2 - Ev + £)J f AT
F F
vi
E Ev c\
2F
ico2 + a\ 4F2a v i 2 - F + F
In = AT
CM
oo2 - aj ooi + aJ - Ev + c
>
CM
F F/
(4.6)
315
E
0)1 - Vi -
2F
E
0)2 = V2 -
2F
2 _ E2 £ _ E2 4F
4F2 " F 4F2 4F2
a = /E2 - 4Fc
Cons ider
0)i0)2-VV - EP-^V1+v2^+^E2p"
1 2 2 2
and
E
0)2 ” 0)1 — V 2 “ - Vi - 2 p = v2 - Vi
2F
VV - p-^V1+v2^+"ÿp
1 2 2
T
Y+a ^v2“ Vi) --L-
T1 + £
4F2 F
e
/ /
And final Iy
E
4F2a
2FVlv2 - E (vi + V 2 ) - 2Fa (v 2 - Vj) + 2c
2F v i v 2 - E (vi + V 2 ) - 2Fa (V2 - vi) + 2c
_L
2F
Evi2 - Evi + c (t2 - t i )
e (4.7)
^Fv!2 - Ev 2 + c
316
(.,2 , c Ev Ev c. ,
(v + p - — + — - F ) dv
( V 2 + £. -
Ev
F
= 1 + dv
— (2v - i) - £ + I L \
2F v f’ F 2F2
1 + dv
(v2 + c _ Ex
’ /
E E c F‘2
1 + — {2v - pr) - p" + 2 F 2"
dv = -Fds (4.10)
Ev
(v2 + j
integrat ing
^ E i /v2 + c V2 /
EV^ + dv
V 2F ln F - T 2F2 /
\ J _ Vl Ev
_ (v2 +
F F '
=-F(s2-si) (4.11)
c ly_____ dv
2F2 F J, o .c Ev which is similar to R.H.S. of
(vz + ^ -
Equation 4.5.
Iett ing
a =
4F2
317
and
a) = v -
2F
afa)
2F2 F i (a)2 - a 2 )
E2 c dfa)
2F2 F 1 1 1
2a o)-a o)+a
u
1. /.E2 c \
2a 2p2 ^r"l (fa)-a ) — In ( (aH-a ) ]
1 / E2 c\ (fa)-a )
2a [ 2F2 " F In
(fa)+a )
V + ^ ln(v2 + f -
V2 c i . /fa)-a v2
2a 2F2 FÌ
F / 1niI
fa)+a,
Vl Vl
= -F( s 2 - Si)
i.e.,
c I . /faj-a V2
+ 1 il
v + w ,n r + f F 2a 2F2 F lrW
VI
=-F(s2-si) (4.12)
_ e_F(s2 - si)
e(V2“ Vi)
i.e
L J_ E2 _ c
2p / , x\ 2 a 2F2 F
( V 2 “ V i ) | FV 2 2 + C - EV2 o)2-a /o)i+a
Fv i 2 + c - 0)2+a (a)i-a
Now
OJi—OJ2 - V i “ V2
/ /
E E2 E2
viv2 - 2F (vi+V2) + /i?2- a (Vl-v2) - w
F^ c
V1V2 - 2 f ivi+v2) + “ a(vi-v2 ) - ,JpT + f
_ 2Fv i v 2 - E ( v i + v 2 ) - 2F a ( v i - v 2 ) + 2c
2FV1V2 - E ( v 1+ v 2) + 2Fa( vi ~V2) + 2c
written.
MONTAGNER, G.J. and TROTT, G.W. "A General Programme for System
Chute System” .
Symp. on Solids and Slurry Flow and Handling in the Chemical Process
System” .
IFAC Symp. on Auto. Control for Agricul. Uni. Saskat. June 1974
pp.D6 1-15.
System” .