CN1011624B - Sonic system for monitoring the trajectory of a projectile - Google Patents

Abstract

The present invention relates to an acoustic method for testing the trajectory parameters of a pill, which belongs to the technology for testing target trajectory parameters. A plurality of sound measurement systems which are arranged at different distance parts on a straight line are utilized to automatically measure the data of coordinate values of the pill running through the sound measurement systems and the flight time of the pill between the sound measurement systems, etc. The projective curves of pill tracks on a vertical plane and a horizontal plane, a vacuum extrapolated curve, a speed-distance curve, a diffusion-distance curve, etc., a regression equation and more than ten movement parameters are automatically drawn after the data is processed by a computer with special software. The present invention realizes the trajectory test requirement for one pill with multiple parameters, and overcomes various defects in single measurement.

Description

Sonic system for monitoring the trajectory of a projectile

The inventive method is to belong to shooting lane trajectory technology.

The testing apparatus and the method for the small-bore bullet flat path of shooting lane parameter are a lot, as piezoelectric quartz sensor, velocity measurement target, initial velocity radar, video camera, acoustic-electric precision target etc.But majority is used for the individual event parameter testing, as measuring pressure, temperature, speed etc.Acoustic-electric precision target development in recent years is very fast, and it is simple in structure, precision is high, easy operating.Britain's " defence " magazine in May, 1973 number reported that the projectile position indicator of Britain MS company is exactly this kind equipment, it is to replace original paper target with the area that a vertical bars and horizon bar are limited, and piezoelectric sensor is equipped with at the two ends of every rod.Rod caused vibration when the bullet impact ripple passed to excellent going up, and induction vibration is at the uniform velocity propagated along rod.Thereafter piezoelectric sensor sends electric pulse because of vibration.Concerning every rod, and the time interval between the output of the corresponding sensor of single shock wave of bullet generation, measuring just become along the rod axial bullet position of axle.Two rods are measured simultaneously, have just determined the position of bullet on measurement plane.Four signals are input to main amplifier, be input to the coordinate computing machine and display device gets on after detection, amplification, elimination echo, the logical circuit of computing machine becomes the conversion of signals at two intervals that four pulses limit along the displacement of two axes of coordinates apart from the face of measurement center initial point.

Method for above-mentioned these test bullets pass the target plane position all can not solve the problem of bullet test multiparameter and robotization.

The objective of the invention is to propose a kind of method problem that solves a trajectory robotization and a bullet multiparameter with acoustic method.

The bullet of supersonic flight is measured by sound target (being acoustic-electric precision target or projectile position indicator) at the coordinate in a certain cross section.

The sound target is made up of target plane, channel amplifier, transmission line, time difference measurement instrument, microcomputer, plotting apparatus etc.Target plane is by orthogonal in same plane and vertical respectively and be parallel to two planes that the sounding rod is formed of surface level, and this plane is perpendicular to surface level, and its normal is parallel with the projectile flight direction.During test, make bullet pass the effective coverage of target plane.

Sounding rod (1), (2) are that aluminum round rod (9) two ends after turning are loaded onto piezoelectric sensor and made.The length of sounding rod is generally 1.5 meters to 5 meters scopes.Its diameter is generally at 25 to 30 millimeters.After shock wave affacts the sounding rod, have an effect with rod, produce acoustic vibration, vibration wave is axially at the uniform velocity propagated along rod, the piezoelectric sensor at rod two ends is transformed into electric signal with vibration signal, these two electric signal have time difference △ t, and the size of △ t and the shock wave active position on the sounding rod is one to one, as long as measure the position that time difference △ t just can obtain application point.The inventive method coverage is 400 meters.

In two sounding rods of target plane, one is horizontal positioned, claims x rod (1), and one is that vertical the placement claims Y rod (2).With x rod (1) is example, and toward target plane, sounding rod (1) application point is projected in excellent central point left side x distance, then △ t on rod from muzzle _xBear, on the contrary if the excellent central point of application point being projected on rod the right x distance, then △ t _xBe positive.Y rod (2) situation is with x rod (1), and the coordinate that bullet passes target plane is:

x＝ 1/2 c·△tx……（1′）

y＝ 1/2 c·△ty……（2′）

Wherein c is sound wave axial velocity of propagation in the sounding rod.

Time difference △ t measures by time difference tester, if the sprocket pulse cycle is T, time difference polarity is K, and then (1 '), (2 ') two formulas can be transformed into:

x＝ (c)/2 △tx＝ 1/2 c·K·T·Nx……（3′）

x＝ (c)/2 △ty＝ 1/2 C·K·T·Nx……（4′）

N wherein _xAnd N _yIt is the count pulse number.As long as the input of the parameter in (3 '), (4 ') Computer Processing just can obtain x and the y value (is the coordinate values on a certain section on the projectile flight track.）

In 200 meters scopes of distance muzzle, near muzzle, settle a piezoelectric sensor, simultaneously with the same rectilinear direction of projectile flight on, be provided with and be a plurality of vertical target planes with this direction.Its quantity can be determined according to concrete content measurement, generally is 4 to 16.Spacing distance between every adjacent two sound target planes is also decided by the requirement of content measurement, generally is to be no more than 30 meters.First target plane also is 25 meters from muzzle.As shown in Figure 1, each target plane is to regulate target holder (3) by three adjusting screw(rod)s (5), make this target plane perpendicular to surface level and make the sounding rod (1) of horizontal positioned be parallel to surface level, the vertical sounding rod of placing (2) is perpendicular to surface level.4 measured signals of four piezoelectric sensors that will be installed in sounding rod (1), (2) two ends are received boxcar (4) with lead, transfer to time difference tester (30) through channel amplifier (6), the long line of transmission (7) again.

Give microcomputer with the time difference umber of pulse that time difference tester measures.Carry out data processing with the program of giving establishment earlier, just can reach purpose of the present invention.

The hyperchannel time difference of the present invention gathers and digital processing system is made up of following each several part: magnetic tape station or disk (17) display (18), microcomputer (19), plotting apparatus (20), computer interface (21), channel address generator (22), time mark generator (23), clock signal bus (24), status bus (25), reset bus (26), internal address bus (27), internal data bus (28), muzzle signal (29), and time difference tester compositions such as (30).X among Fig. 3 ₁～x _nBe the 1st to n sound target x time difference tester; y ₁～y _nBe the 1st to n sound target Y time difference tester; f ₁～f _nBe projectile time of flight tester x ^-, x ⁺, y ^-, y ⁺Represent that respectively x is to sounding rod (1) and Y negative terminal and the positive end signal to sounding rod (2).Wherein muzzle signal (29) be loaded on the rifle or muzzle near the signal that sonic transducer received.Through with the same processing of the signal of sounding rod (1), (2) after signal, be used to survey the enabling signal of flight time.Here the span of n is 4～16.Promptly sound target plane quantity.

Sounding rod (1), (2) two ends are except that piezoelectric sensor (10), also has an impedance transformer (11) (see figure 5), its power supply is supplied with by channel amplifier (6), four signals of each target plane are guided to interconnecting device (4), carry out circuit by a split conductor with channel amplifier (6) again and be connected, this root connecting line has four signal wires from sounding rod (1) (2), the power lead that channel amplifier (6) is arranged and the common ground of housing shroud.

Channel amplifier (6) behind wave filter (12), amplifier (13), voltage comparator (14) and timing rejector (15), drives the output (see figure 2) to signal by line driver (16).

The present invention requires the channel amplifier gain greater than 70db, band bandwidth 10～50KHz, and the amplifier input noise is less than 0.4mv, and input impedance is greater than 10 ⁶Ω, difference output impedance 50 Ω, reference voltage is adjusted at 2.5v, suppresses circuit and adopts monostable circuit, and the restriction of frequency is shot in the adjustment of metastable state time, promptly must be less than the largest interval time of per two bullets under the running fire state.The performance requirement unanimity of No. four amplifiers in each channel amplifier (6), gain should be able to be regulated with switch, to be applicable to single-shot, running fire, indoor and outdoors various situations such as (being shooting lane and target range).

In the inventive method there being the building block of used hyperchannel time difference tester: opens earlier and afterwards close circuit (31), channel address code translator (32), status register (33), time difference counter (34), alarm (35), polarity detector (36), output buffer (37) and line receiver (39) etc.And the flight time tester is similar to the hyperchannel time difference tester of above-mentioned introduction, does not just open earlier afterwards to close circuit (31), input indication (38) and line receiver (39).Its enabling signal is from the signal of the sounding rod (2) of last target plane, if first flight time tester is then used muzzle signal (29), shutdown signal is from the signal of the sounding rod (2) of back one target plane.

Line receiver (38) receives the differential pulse signal from long line, and it is transformed into standard TTL signal.

Open earlier that afterwards to close circuit (31) be the gate-control signal of time difference impulse meter (34), open counter gate from first comer in the signal of same sounding rod, begin counting, then the counting door is closed to the person, stop counting, make channel status trigger (33) set (effective status) simultaneously.

The output of output buffer (37) is subjected to the control of channel status register (33), channel address code translator (32) and alarm three signals such as (35).Could output signal when having only these three signals all effective.

Channel address code translator (32) is a kind of coincidence circuit, and when channel address was identical with the address of giving setting earlier, it was output as " effectively ", otherwise is invalid.Alarm (35) is the equipment that a kind of state can be provided when the counter that causes for various reasons is out of control, if counter (34) is out of control, then going out of alarm (35) is engineering noise, just often is output as " effectively ".

It is that of transmitting earlier in two signals that detect from same sounding rod is negative terminal signal or positive end signal that polarity detects (36), if what transmit earlier is the negative terminal signal, and high level of polarity detector (36) output then; On the contrary, be positive end signal, then output low level as if what transmit earlier.The output that polarity detects (36) is the part of data.

Input indication (38) is the drive signal that provides for panel display LED, when input, input pointer (38) just export a broadening pulse signal, the driving LED display.

Channel address generator and controller (22), it produces channel address, the start and stop that channel address scans are subjected to the control of outer signals, when the state of channel address passage pointed when being effective, scan and stop, at this moment the state of this passage is served status bus, put forward data, simultaneously the output buffer gate of this passage is opened by computer interface circuit (21) prompting computing machine; Channel address is taken away by computing machine as the part of data.The structure of channel data is 17 word lengths, and wherein 11 is time difference data, and 1 is polarity, and 5 is channel address.The length of channel address is decided on used port number, each target plane needs two placed channels, need a flight time passage between per two target planes, eight target planes are totally ten six placed channels and eight flight time passages totally two ten four passages, and then channel address needs 5.After computing machine is taken data away, provide " answer " signal, channel status is resetted, recover channel address simultaneously and scan.

It produces a high-frequency impulse by crystal oscillator time mark generator (23), after obtain the pulse of multiple frequency behind a series of frequency division, target is selected when various to adapt to.

Computer interface circuit (21) is the bridge of time difference tester and intercomputer, by it, computing machine gets the status information of arrival time difference tester, takes time difference data away by its computing machine, that channel reset that computing machine utilizes controller to make data are provided by it.

Some parameters of hyperchannel time difference tester of the present invention are:

Internal data bus: (D ' _0～11) (28) 12;

Internal address bus (A ' _0～4) (27): 5;

Clock signal bus (24): 1;

Status bus (25): 1;

Reset bus (26): 1.

Time difference measurement port number: 24

Computing machine and external unit:

Microcomputer (19) is general general eight micro computers.That the inventive method is used is TRS-80, also can use the micro computer of PS-80, MDR-Z80, DPS-80A, DJS-043, LWY-80A type.

Magnetic tape station or disk (17) are the external units of computing machine, and it is to be used for the equipment of stored routine.

Plotting apparatus (20) can automatically be drawn out various curves and print various parameters according to the instruction and the result of calculation of microprocessor (19).

Fig. 1 is sound target plane and the channel amplifier arrangement plan in the inventive method.Sounding rod (1), (2) are arranged among the figure, target holder (3), interconnecting device (4), adjusting screw(rod) (5), channel amplifier (6) transmits long line (7), prime amplifier (8) and time difference tester (30).

Fig. 2 is the channel amplifier block diagram in the inventive method.Block diagram is by wave filter (12), amplifier (13), voltage comparator (14), regularly rejector (15), line driver (16) and the long line of transmission (7) composition.

Fig. 3 gathers and digital processing system the hyperchannel time difference in the inventive method.Include among the figure: magnetic tape station or disk (17), display (18), microcomputer (19), plotting apparatus (20), computer interface (21), channel address generator (22), time mark generator (23), clock signal bus (24), status bus (25), reset bus (26), internal address bus (27), internal data bus (28), each parts of muzzle signal (29) and time difference tester (30).

Fig. 4 is the block diagram of the time difference tester in the inventive method.Include among the figure as lower member: open earlier and afterwards close circuit (31), channel address code translator (32), status register (33), time difference counter (34), alarm (35), polarity detector (36), output buffer (37) and line receiver (39), input pointer (38).

Fig. 5 is the sounding bar structure figure in the inventive method.Include metal garden rod (9) among the figure, piezoelectric sensor (10) and impedance transformer three parts such as (11).

Drop shadow curve and the target surface deflection figure of ballistic trajectory on vertical plane that Fig. 6 obtains for field measurement in the inventive method.K among the figure ₁To K ₈Be the sequence number of eight sound target planes, K ₈Numeral 1 to 4 on the sound target plane is meant the pairing curve of bullet preface of projectile, does not indicate that wherein a digital curve is meant the mean trajectory curve of 1 to 4 bullet.

Fig. 7 is actual measurement vertical projection center of impact track vaccum extrapolation curve.

Fig. 8 is that the field measurement ballistic trajectory is at surface level Projective Curve and sound target plane deflection figure.L is target center the distance between surface level on projection join line (claim target center line) of muzzle point at the projection on the surface level and each target plane among the figure.

Fig. 9 is a measured level projection average point of impact track vaccum extrapolation curve.

Figure 10 is actual measurement velocity of shot (V) and distance (S) relation curve and regression equation thereof.V among the figure ₀, V ₅₀, V ₁₀₀To V ₂₀₀Represent bullet respectively at muzzle, from 50 meters in muzzle, 100 meters speed to 200 meters.R is the degree of fitting of curve.

Figure 11 is relation curve and the regression equation thereof of distribution with distance.E among the figure _xBe the relation curve of deflection probable error and distance, E _yBe the relation curve of height probable error and distance, R ₅₀For half must in the relation curve of garden and distance, R ₁₀₀Relation curve for extreme spread garden and distance.

A specific embodiment of the inventive method is: on 200 meters straight line shooting lane, 25 meters at every interval is placed one can be adjusted to the sound target plane that falls in ground by adjusting screw(rod) (5), and makes the normal direction on its plane be parallel to the direction of projectile flight.Be total to 8 such sound target planes of device.

On each target plane, all include two orthogonal, and vertical respectively again and be parallel to the sounding rod (1) of ground level, (2).Each sounding rod be by one long 2 meters, diameter is the long rod of 30 millimeters aluminium matter, and at its two ends piezoelectric sensor (10) and impedance conversion (11) is housed respectively.

The signal that four piezoelectric sensors (10) on each target plane are surveyed arrives boxcar (4) through many stamens cable transmission, and channel amplifier (6) amplifies.Signal after the amplification is transferred to time difference tester by the many stamens long lines of transmission (7) again, thereby on the acquisition diverse location, piezoelectric sensor records respectively organizes time difference value.

Each group time difference signal is sent into digital processing system again, be little meter machine and peripheral each parts thereof, make it to carry out various calculating by giving the program of setting earlier, at last again by the plotting apparatus various parametric lines that draw automatically, and print the related parameter that has of each curve automatically, thereby realized the automatic operation of a whole set of proving installation.

Adopting the inventive method, is 7.62 millimeters small-bore gun when carrying out single-shot or burst-firing in actual field with the bullet diameter, and the shooting of each bullet can be drawn out the various parametric lines of Fig. 6 to Figure 11 automatically.And print the various parameters of table 1 automatically to table 3.

K wherein: target number;

S: the target spacing, first is the distance of muzzle to the first target plane;

Xo, Yo: the center of impact coordinate figure of one group of bullet;

Ex, Ey: direction and height probable error;

R ₅₀, R ₁₀₀: half must in garden and extreme spread garden;

T/V: average velocity between flight time or target between target

The root-mean-square error of RMS:T/V item

Table 2:

K/NR：8/8 DIS.B.T：25

M X（MM） Y（MM） T（MS） SCORE

1 -139.7 108.8 44.7 0

2 -193.8 32.2 44.475 7

3 -135.2 9.7 44.8 8

4 -171.2 27.7 44.925 7

TOTAL＝22

MEAN：-160 44.6 44.725

RMS： +/-.1904

Ex＝18.7 EY＝29.6

R ₅₀＝39.4 R ₁₀₀＝67.3

Wherein: K/NR:8/8: the data that are the 8th target.If table 1 is not drawn the form that then can all be depicted as the data of preceding 7 targets table 2,, then can only draw table 2 if selected table 1;

DIS.B.T:25 represents that the 7th target to the 8 ranges are from being 25 meters;

M: play preface

X, Y: coordinate values

T: from the flight time of the 7th target to the 8 targets

SCORE: the scoring ring number;

TOTAL=: total number of rings.

Table 3:

DIS：200.6

M T（MS）

1 322.45

2 321.5

3 323.175

4 323.95

MEAN：322.768

RMS：1.04

DIS wherein: be to cut a section length (being 200.6 meters here)

T: total flight time.

Show that from test result the absolute error that tests the speed (comparing with canopy of the heavens target) is that the bullet shown in ± 1.0 meter per seconds and the table 4 precision a little.

Table 4. precision (with mm is unit, every group of 20 bullets)

The accurate deviation of the distance location intensive boundary mark of probable error

Single-shot ± 1.0 ± 1.0

100 meters

Running fire ± 2.0 ± 1.5

Single-shot ± 2.0 ± 1.5

200 meters

Running fire ± 3.0 ± 2.0

From results measured as can be seen, adopt the inventive method, can be effectively applied to the modern military training, to the development and the quality indicator of small arms rifle, bullet, shortened the cycle of weapon development and the application that can supply the military science research aspect in the war production.Avoid the danger of artificial indication of shots fully, saved starting material and funds, improved efficient and measuring accuracy greatly.

00010 ORC 0D6EFH

00020 SELECT DEFS 2

00030 DELAY DEFS 2

00040 IT DEFS 1

00050 NUM DEFS 1

00060 NF DEFS 2

00070 NF3 DEFS 2

00080 NM DEFS 1

00090 TOTAL DEFS 2

00100 ACNF DEFS 2

00110 MEMORY DEFS 1440D

00120 DSPAC1 EQU 3FB8H

00130 DSPAC2 EQU 3FF8H

00140 AUTDIS EQU 3F78H

00150 TOTDIS EQU 3F7EH

00160 COUNT1 DEFS 1

00170 IP1 DEFS 1

00180 COUNT2 DEFS 1

00190 IP2 DEFS 1

00200 LD A，（OB829H）

00210 CP OCH

00220 JP NZ，BADBAS

00230 CALL ERACE

00240 LD A，L

00250 CP 2H

00260 JP Z，ENTRYO

00270 LD A，（SELECT）

00280 CP 1

00290 JR Z，NFMODE

00300 LD BC，（NF3）

00310 JR SEDISP

00320 NFMODE LD BC，（NF）

00330 SEDISP LD A，（SELECT）

00340 ADD A，2FH

00350 LD （3F7BH），A

00360 LD A，26H

00370 LD （3FFEH），A

00380 LD （3FFAH），A

00390 LD （3FBAH），A

00400 LD （3FBEH），A

00410 LD A，41H

00420 LD （3FFFH），A

00430 LD A，93H

00440 OUT （03H），A

00450 LD L，20H

00460 TOERAC LD A，7FH

00470 OUT （02H），A

00480 LD A，OFFH

00490 OUT （02H），A

00492 LD A，7FH

00494 OUT （02H），A

00500 DEC L

00510 JP NZ，TOERAC

00520 LD A，30H

00530 LD E，0

00540 LD （AUTDIS），A

00550 LD （AUTDIS+1），A

00560 PUSH BC

00570 LD HL，TOTDIS

00580 LD A，（TOTAL）

00590 CALL BINH2

00600 LD A，（SELECT）

00610 CP 1

00620 JR Z，START

00630 LD HL，DSPAC1

00640 LD A，OH

00650 LD （COUNT1），A

00660 CALL BINH2

00670 INC HL

00680 LD A，（IP1）

00690 CALL BINH2

00700 LD A，（SELECT）

00710 CP 2H

00720 JR Z，START

00730 LD HL，DSPAC2

00740 LD A，OH

00750 LD （COUNT2），A

00760 CALL BINH2

00770 INC HL

00780 LD A，（IP2）

00790 CALL BINH2

00800 START PUSH DE

00810 CALL OA7FH

00820 DEC L

00830 JP Z，ADJUST

00840 POP DE

00850 POP BC

00860 LD HL，MEMORY

00870 ST LD A，（DELAY）

00880 LD D，A

00890 ST1 LD A，（4301H）

00900 LD E，30H

00910 CALL ERACE

00920 IN A，（04H）

00930 SRL A

00940 JR C，ST1

00950 FIRST LD （HL），A

00960 LD A，（IP1）

00970 INC HL

00980 LD A，30H

00990 CALL ERACE

01000 FIRST1 INC HL

01010 IN A，（01H）

01020 LD （HL），A

01030 INC HL

01040 IN A，（00H）

01050 LD （HL），A

01060 INC HL

01070 CALL ERACE

01080 LD A，（IT）

01090 SRL A

01100 CALL ERACE

01110 LD A，31H

01120 LD （3FFFH），A

01130 ND2 IN A，（02H）

01140 SRL A

01150 JR NC，SECOND

01160 DEC D

01170 JP Z，DISPL1

01180 JR ND2

01190 SECOND LD （HL），A

01200 INC HL

01210 IN A，（01H）

01220 LD （HL），A

01230 INC HL

01232 IN A，（00H）

01234 LD （HL），A

01236 INC HL

01240 CALL ERACE

01250 LD A，（IP2）

01260 SRL A

01270 CALL ERACE

01280 LD A，32H

01290 LD （3FFFH），A

01300 TH LD A，（DELAY）

01310 LD D，A

01320 TH3 IN A，（02H）

01330 SRL A

01340 JR NC，THIRD

01360 DEC D

01370 JR Z，DISPL2

01380 JR TH3

01390 THIRD LD （HL），A

01400 INC HL

01410 IN A，（01H）

01420 LD （HL），A

01430 INC HL

01440 IN A，（00H）

01450 LD （HL），A

01460 INC HL

01470 CALL ERACE

01480 LD A，33H

01490 LD （3FFFH），A

01500 EADD1 INC E

01510 LD A，E

01520 PUSH BC

01530 PUSH HL

01540 LD HL，AUTDIS

01550 CALL BINH2

01560 POP HL

01570 POP BC

01580 DEC BC

01590 LD A，B

01600 OR C

01610 JP NZ，ST

01620 TRANS LD HL，MEMORY

01630 LD BC，（NF）

01640 LD （ACNF），BC

01650 TRANS1 LD A，（HL）

01660 AND 1H

01670 LD （HL），A

01680 INC HL

01690 LD A，（HL）

01700 AND 248D

01710 SRL A

01720 SRL A

01730 DEC HL

01740 ADD A，（HL）

01750 LD （HL），A

01760 INC HL

01770 LD A，（HL）

01780 AND 7H

01790 LD （HL），A

01800 DEC HL

01810 SRL （HL）

01820 INC HL

01830 JP NC，NOSET

01840 SET 3，（HL）

01850 NOSET INC HL

01860 INC HL

01870 DEC BC

01880 LD A，B

01890 OR C

01900 JR NZ，TRANS1

01910 LD A，42H

01920 LD （3FFFH），A

01930 LD A，（TOTAL）

01940 ADD A，E

01950 LD （TOTAL），A

01960 BASIC JP OA9AH

01970 ERACE LD A，7FH

01980 OUT （02H），A

01990 LD A，OFFH

02000 OUT （02H），A

02010 LD A，7FH

02020 OUT （02H），A

02030 RET

02040 DISPL2 DEC HL

02050 DEC HL

02060 LD A，（HL）

02070 SRL A

02080 SRL A

02090 SRA A

02100 PUSH AF

02110 DEC HL

02120 SRL （HL）

02130 JR C，$+9

02140 LD A，2BH

02150 LD （3FFDH），A

02160 JR $+7

02170 LD A，2DH

02180 LD （3FFDH），A

02190 POP AF

02200 LD （IP2），A

02210 LD A，（COUNT2）

02220 INC A

02230 LD （COUNT2），A

02240 PUSH BC

02250 PUSH HL

02260 LD HL，DSPAC2

02270 CALL PUTHEX

02280 INC HL

02290 LD A，（IP2）

02300 CALL BINH2

02310 INC HL

02320 INC HL

02330 LD （HL），2AH

02340 POP HL

02350 POP BC

02360 DISPL1 DEC HL

02370 DEC HL

02380 LD A，（HL）

02390 SRL A

02400 SRL A

02410 SRL A

02420 PUSH AF

02430 DEC HL

02440 SRL （HL）

02450 JR C，$+9

02460 LD A，2BH

02470 LD （3FBDH），A

02480 JR $+7

02490 LD A，2DH

02500 LD （3FBDH），A

02510 POP AF

02520 LD （IP1），A

02530 LD A，（COUNT1）

02540 INC A

02550 LD （COUNT1），A

02560 PUSH BC

02570 PUSH HL

02580 LD HL，DSPAC1

02590 CALL PUTHEX

02600 INC HL

02610 LD A，（IP1）

02620 CALL BINH2

02630 INC HL

02640 INC HL

02650 LD （HL），2AH

02660 POP HL

02670 POP BC

02680 JP ST

02690 PUTHEX LD B，A

02700 CALL BINH2

02710 RET

02720 BINH4 PUSH BC

02730 LD A，B

02740 CALL BINH2

02750 POP BC

02760 LD A，C

02770 BINH2 LD B，A

02780 RRA

02790 RRA

02800 RRA

02810 RRA

02820 CALL BINH1

02830 LD A，B

02840 BINH1 AND OFH

02850 ADD A，30H

02860 CP 3AH

02870 LD （HL），A

02880 INC HL

02890 RET C

02900 DEC HL

02910 ADD A，07H

02920 LD （HL），A

02930 INC HL

02940 RET

02950 INRUPT LD A，（SELECT）

02960 LD B，A

02970 LD D，O

02980 LD HL，O

02990 LOOP1 ADD HL，DE

03000 DJNZ LOOP1

03010 LD （ACNF），HL

03020 LD A，（NUM）

03030 CP 1H

03040 JP Z，SINGLE

03050 COUNTI LD HL，MEMORY

03060 LD BC，（ACNF）

03070 LD A，43H

03080 LD （3FFFH），A

03090 LD A，B

03100 OR C

03110 JP NZ，TRANS1

03120 JP 0A9AH

03130 SINGLE LD HL，O

03140 LD A，53H

03150 LD （3FFFH），A

03160 JP 0A9AH

03170 DEL160 EXX

03180 LD BC，0101H

03190 LOOP2 DEC C

03200 JR NZ，LOOP2

03210 DEC B

03220 JR NZ，LOOP2

03230 EXX

03240 RET

03250 ADJUST LD A，4AH

03260 LD （3FFFH），A

03270 POP DE

03280 POP BC

03290 LD A，3FH

03300 OUT （02H），A

03310 EXX

03330 LD BC，012AH

03340 LOOP3 DEC C

03350 JR NZ，LOOP3

03360 DEC B

03370 JR NZ，LOOP3

03380 EXX

03390 LD A，2FH

03400 OUT （02H），A

03402 LD A，6FH

03404 OUT （02H），A

03410 LD A，7FH

03420 OUT （02H），A

03430 AJACQU LD BC，（NF）

03440 LD HL，MEMORY

03450 LD D，OFFH

03460 AJASK IN A，（02H）

03470 SRL A

03480 JR NC，AJBYTE

03490 DEC D

03500 JP Z，AJBASI

03510 JR AJASK

03520 AJBYTE LD （HL），A

03530 INC HL

03540 IN A，（01H）

03550 LD （HL），A

03560 INC HL

03570 IN A，（00H）

03580 LD （HL），A

03590 INC HL

03600 CALL DEL160

03610 CALL ERACE

03620 INC E

03630 LD A，E

03640 PUSH BC

03650 PUSH HL

03660 LD HL，AUTDIS

03670 CALL BINH2

03680 POP HL

03690 POP BC

03700 DEC BC

03710 LD A，B

03720 OR C

03730 JP NZ，AJASK

03740 LD A，26H

03750 LD （3FFEH），A

03760 JP TRANS

03770 AJBASI LD B，O

03780 LD C，E

03790 LD （ACNF），BC

03800 LD HL，MEMORY

03810 LD A，42H

03820 LD （3FFFH），A

03830 LD A，B

03840 OR C

03850 JP NZ，TRANS1

03860 JP OA9AH

03870 BADBAS LD A，0AH

03880 LD （42E8H），A

03890 LD A，O

03900 LD （42E9H），A

03910 LD （42EAH），A

03920 RET

03940 ENTRYO LD A，O

03950 LD （43BCH），A

03960 RET

03970 END OD6EFH

10 CLEAR（200）∶CLS∶AJ$＝"Y"

15 POKE16526，178∶POKE16527，245∶POKE-2464，1∶POKE-2459，2

20 AJ＝2∶J＝USR（AJ）∶IFAJ$＝"Y"AJ$＝"N"∶GOTO15

25 PRINT"MOMENT 10 S＊"∶DEFINT J，N，P∶DIM FT$（9），S（8），X（9），Y（9），NR（8），RI（20），NI（8，25），XT（8，25），YT（8，25），T（25，8），N（30），XO（8），YO（8），C（3），D（3），L（3），IP（20，50），PF（25），NF（8），TS（8），R50（8），RH（8），DT$（8），EX（8），EY（8），NA（8），S1（8），S2（8），TT（25），A（3，4），TV（8），XX（9），YY（9）

30 DIMTM（8），SS（9），ZZ（9），ST（9）∶A1＝25.38∶B1＝30.38∶B2＝-20.

38∶YY＝3.75∶L（1）＝11.9∶L（2）＝17∶L（3）＝26.9∶X2＝11.38∶X3＝X2∶X4＝10.48∶C（1）＝-0.582∶C（2）＝10.2∶C（3）＝-7.071∶D（1）＝18.3∶D（2）＝-100.3∶D（3）＝98.3∶D5＝12.5∶R1＝17.94

35 SE＝-4097∶DY＝SE+2∶F1＝DY+2∶F2＝F1+1∶F3＝F2+1∶F4＝F3+2∶F5＝F4+2∶TL＝F5+1∶AF＝TL+2∶MY＝AF+2

40 PL$＝"N"∶LP$＝"N"∶INPUT"＊ DO YOU HAVE WX4675 PLOTTER（<EN>＝NO W＝YES）";PL$∶IF PL$<>"W"THEN GOTO 50

45 INPUT"TEST DATE∶";DA$∶NA$＝"CPLA 59185"∶PRINT" TESTER∶"NA$;∶INPUTNA$∶GT$＝"7.62"∶PRINT" TYPE OF GUN∶"GT$;∶INPUTGT$∶INPUT" GUN′S NO.";GN$

50 DP＝1∶INPUT" ＊ DIVIDING OF POSITION PORTS <EN>＝4MHZ";

DP∶F＝4/DP

55 DT＝255∶POKEDY，DT

60 U＝5150/（F＊2E4）

65 IT＝1∶INPUT"＊ HOW MANY TARGETS （1-8<EN>＝1）";IT∶POK EF1，IT

70 FORK＝1TO8∶Y（K）＝2＊K∶X（K）＝Y（K）-1∶NEXTK

75 FORK＝1TOIT∶PRINT" X（"K"）＝"X（K）;∶INPUTX（K）∶PRINT"

Y（"K"）＝"Y（K）;∶INPUTY（K）∶NR（K）＝INT（X（K）/2+.5）∶NF（K）＝NR（K）+20∶NT（K）＝X（K）+Y（K）+NF（K）∶PRINTTAB（15）"NR（"K"）＝";NR（K），"NF（"K"）＝"NF（K）∶PRINTTWO-TARGET′S DISTANCE S（"K"）＝";∶INPUT S（K）

80 IF（X（K）+1）/INT（X（K）/2+1）＝2ANDX（K）<16ANDY（K）<17ANDY（K）-X（K）＝1THENNEXTKELSEPRINT"AGAIN INPUT THE PORTS"∶K＝K-1∶NEXTK∶SS＝0∶FORK＝1TOIT∶SS＝SS+S（K）∶NEXTK∶TT（O）＝SS

85 NUM＝1∶INPUT"＊ CONTINUAL FIRE NUMBERS （1-20<EN>＝1）";

NUM∶POKEF2，NUM∶NM＝NUM∶POKEF5，NM∶AU＝NUM

90 TV$＝"N"∶INPUT"＊DO YOU FIND TIME & SPEED（<EN>＝NO Y＝YES）";TV$∶IF TV$＝"Y" THEN TU＝3 ELSE TU＝2∶FT＝1 GOTO105＊ DO YOU FIND TIME & SPEED（<EN>＝NO Y＝YES）";TV$∶IFTV$＝"Y" THENTU＝3ELSETU＝2∶FT＝1∶GOTO58

95 DF＝1∶INPUT"DIVIDING OF FLYING PORTS<EN>＝4MHZ";DF∶FT＝DF/4/1E3∶SS＝0∶FORK＝1TOIT∶FT$（K）＝"Y"∶PRINT" K＝";K;∶IN

PUT" DO YOU FIND THE TIME（<EN>＝TIME N＝SPEED）";FT$（K）∶IFFT$（K）<>"Y"THENPRINT"TWO-TARGET′S DISTANCE S（";K;"）＝（M）";∶INPUTS（K）

100 SS＝SS+S（K）∶NEXTK∶TT（0）＝SS

105 NF＝NUM＊TU＊IT∶IFNF<256J＝0∶I＝NFELSEJ＝1∶I＝NF-256

110 IFTV$＝"Y"THENST＝2∶POKESE，ST+1∶GOTO115ELSEST＝1∶POKESE，ST+1∶GOTO115

115 POKEF3，I∶POKEF3+1，J∶POKEF4，NF/TU ∶POKEF4+ 1，0∶IFNUM＝1 THENNM＝20∶INPUT"＊ TOTAL-FIRE NUMBERS IS THIS TEST（1-20<EN>＝20）";NM∶POKEF5，NM

120 CLS∶PRINT@5，"IT＝"IT" NUM＝"PEEK（F2）" NM＝"PEEK（F5）"TV$＝"TV$" DP＝"4/DP"MHZ";∶IFTV$＝"Y"THENPRINT" DF＝"4000/DF"KHZ";∶GOSUB130∶GOSUB2085ELSEGOSUB130∶GOSUB2085

125 NT＝NT+1∶PRINT@384，"＊＊ TEST NO."NT;∶INPUT"＊＊";NT∶PRINT@530，"＊＊＊ JUST A MOMENT，PLEASE ＊＊＊"∶FORI＝iTO28∶N（I）＝0∶NEXTI∶FORI＝1TO8∶FORJ＝1TO20∶XT（I，J）＝0∶YT（I，J）＝0∶T（J，I）＝0∶NI（I，J）＝0∶IP（I，J）＝0∶PF（J）＝0∶NEXTJ，I∶IW＝1∶ID＝MY+1∶LN＝0∶GOTO145

130 IFTV$＝"N"THENRETURN

135 PRINT∶A＝75∶B＝A∶N＝0∶FORK＝1TOIT∶N＝N+1∶PRINT@B，"FT$（"K"）＝"FT$（K）;∶B＝B+12∶PRINT@B，"S（"K"）＝"S（K）;∶B＝B+12∶IFN＝2THENN＝0∶A＝A+64∶B＝A∶NEXTKELSENEXTK

140 RETURN

145 IFAJ$＝"Y"THENGOSUB1635ELSEIFIT<>1ANDNUM＝1THENA＝80∶GOSUB2060ELSEA＝.04∶GOSUB2060∶GOTO150

150 NUM＝PEEK（F2）∶NM＝PEEK（F5）∶NF＝PEEK（F3）∶IFNUM<>1THEN23OELSEFORK＝1TOIT∶IFPL$＝"W"ANDK<>1THEN170ELSECLS∶GOSUB270∶GOSUB600

155 PF＝0∶XO＝XO（K）∶YO＝YO（K）∶NA（K）＝0∶S1（K）＝0∶S2（K）＝0∶FORM＝1TONM∶IFK＝1THENGOSUB850∶PF＝0∶IFJ＝0THEN165ELSEGOTO160ELSEGOSUB270∶G126＝0∶GOSUB305∶GOTO165

160 GOSUB270∶GOSUB210∶GOSUB185∶IFLN＝1THENIFG85＝1THENLN＝0∶G85＝OELSELN＝0∶G85＝0∶M＝M-1ELSEGOTO165

165 NEXTM

170 IFNM<>0THENGOSUB450∶GOSUB1205ELSENT＝NT-1∶GOTO120

175 NEXTK∶CT$＝"N"∶IFIT＝1THEN180ELSEGOSUB545∶GOSUB530

180 IFPL$＝"W"THENGOSUB1420∶GOTO120ELSEGOTO120

185 IFIT<>1ANDNUM<>1G126＝1∶FORK＝1TOIT∶GOSUB270∶GOSUB305∶NEXTK∶RETURN

190 IFAJ$＝"Y"G126＝0∶FORK＝1TOIT∶GOSUB270∶GOSUB305∶NEXTK∶K＝1∶RETURN

195 IFIT＝1THENIFNUM＝1THENP＝M∶PF＝1∶GOTO205ELSEP＝1∶PF＝PF+

1∶GOTO200ELSEFORK＝1TOIT∶IW＝2＊K-1∶P＝M∶PF＝PF+1∶GOSUB270∶GOSUB420∶NEXTK∶K＝1∶RETURN

200 IFM＝1THEN205ELSEIW＝IW+2∶ID＝ID+6

205 GOSUB420∶RETURN

210 IFTV$＝"N"ORAJ$＝"Y"THENRETURN

215 J＝MY∶B＝0∶N＝0∶K＝0∶FORI＝1TONF∶N＝N+1∶A＝PEEK（J）∶B＝B+A∶J＝J+3

220 IFN＝3THENN＝0∶IFIT＝1THENK＝1ELSEK＝K+1ELSENEXTI

225 IFB＝NT（K）THENPRINT@362，" ";∶B＝0∶NEXTI∶PRINT@362，" ";∶RETURNELSEPRINT@362，"＊＊";∶GOSUB435∶RETURN

230 CLS∶A＝.04∶GOSUB2060∶CLS∶K＝1∶M＝1∶GOSUB270∶GOSUB600∶GOSUB850∶GOSUB210∶PF＝0∶IFNUM<>0GOTO235ELSENT＝NT-1∶GOTO120

235 FORK＝1TOIT∶NA（K）＝0∶NEXTK∶K＝1∶FORM＝1TONUM∶GOSUB270∶GOSUB185∶NEXTM∶FORK＝1TOIT∶XO＝XO（K）∶YO＝YO（K）∶S1（K）＝0∶S2（K）＝0∶GOSUB450∶GOSUB1205

240 NEXTK∶CT$＝"N"∶IFIT＝1THEN245ELSEGOSUB545∶GOSUB530

245 IFPL$＝"W"GOSUB1420

250 IFCC＝0GOTO120ELSEGOSUB255∶CC＝0∶GOTO120

255 CLS∶PRINT"CC＝";CC∶FORI＝MYTOMY+NF＊3-1∶PRINTPEEK（I）;∶IFI＝-3883THENINPUT"PRESS THE <EN> TO CONTIUE <C>";G$∶NEXTIELSENEXTI

260 IFLP$＝"Y"LPRINT"CC＝";CC∶FORI＝MYTOMY+NF＊3-1∶LPRINTPEEK（I）;∶NEXTI

265 PRINT∶INPUT"PRESS TNE <EN> TO CONTIUE";I∶RETURN

270 PRINT@0，"T.NO.（K/NR）";K;"/";NR（K）;

275 PRINT@109，STRING$（18，32）;∶PRINT@173，STRING$（18，32）;

280 PRINT@128，STRING$（21，32）;∶PRINT@237，STRING$（18，32）;

285 PRINT@192，STRING$（21，32）;∶PRINT@301，STRING$（18，32）;

∶RETURN

290 IFTV$＝"Y"THENIF PF（PF）<>NF（K）THENPRINT@301，"IP＝"PF（PF）" ";∶GOSUB400ELSEGOTO295ELSEGOTO300

295 PF（PF）＝0

300 RETURN

305 PRINT@45，"FIRE NO.";M;" ";∶X（9）＝XT（NR（K），M）∶Y（9）＝YT（NR（K），M）

310 X（9）＝X（9）-XO/10∶Y（9）＝Y（9）-YO/10∶X（9）＝INT（X（9）＊100+.5）/100∶Y（9）＝INT（Y（9）＊100+.5）/100∶XT（NR（K），M）＝X（9）∶YT（NR（K），M）＝Y（9）

315 R＝SQR（X（9）＊X（9）+Y（9）＊Y（9））∶R＝INT（R＊100+.5）/100∶RI（M）＝R

320 PRINT@109，"X＝";X（9）＊10;"MM";

325 PRINT@173，"Y＝";Y（9）＊10;"MM";

330 IF R>0.4 THEN 350

335 N＝10∶NA（K）＝NA（K）+N∶PRINT@237，"N＝10";

340 PRINT@128，"JUST CENTRE ！！";∶PRINT@256，" ";

345 GOTO 360

350 GOSUB 1090

355 IFTV$＝"Y"THENGOTO360ELSEGOTO370

360 IFFT$（K）＝"Y"THENT（M，NR（K））＝INT（T（M，NR（K））＊1E4+.5）/1E4∶PRINT@301，"T＝"T（M，NR（K））"MS";ELSEIFT（M，NR（K））＝0THENPRINT@301，"V＝"S（K）"/0";ELSEIFG126＝1THENT（M，NR（K））＝INT（S（K）/T（M，NR（K））＊1E5+.5）/100∶PRINT@301，"V＝"T（M，NR（K））"M/S";ELSEPRINT@301，"V＝"T（M，NR（K））"M/S";

365 IFIT>＝2ANDPL$＝"W" Y（9）＝Y（9）/2

370 GOSUB 710

375 NI（NR（K），M）＝N

380 RETURN

385 PRINT@128，IP（P，IW）;"S.-PORT IS FAULT！";∶GOSUB390∶RETURN

390 IFIT＝1THENPRINT@172，IP（P，IW）;∶IFTV$＝"N"THENPRINTPEEK（ID）;PEEK（ID+1）;ELSEGOTO395ELSEPRINT@172，IP（P，IW）;∶PRINT@236，IP（P，IW+1）;∶GOTO400

395 PRINT`236，IP（P，IW+1）;∶IFTV$＝"N"PRINTPEEK（ID+3）;PEEK（ID+4）;

400 IFAU＝1THENM＝M-1∶G85＝1∶IFIT＝1THENGOTO405ELSEK＝ITELSECC＝CC+1∶PRINT"-";CC;∶INPUT"<EN>";G$

405 GOSUB435∶RETURN

410 PRINT@128，DE;"EX-PORT IS FAULT！";∶DE＝0∶GOSUB390∶RETURN

415 PRINT@128，IP（P，IW）;"EX-PORT IS FAULT！";∶PRINT@192，IP（P，IW+1）;"EX-PORT IS FAULT！";∶GOSUB390∶RETURN

420 REM

425 IF（IP（P，IW）＝X（K）ORIP（P，IW）＝Y（K））AND（IP（P，IW+1）＝X（K）ORIP（P，IW+1）＝Y（K））THENIFIP（P，IW）＝IP（P，IW+1）THENGOSUB385∶RETURNELSEG126＝1∶GOSUB305∶GOSUB290∶RETURNELSEIF（IP（P，IW）＝X（K）ORIP（P，IW）＝Y（K））OR（IP（P，IW+1）＝X（K）ORIP（P，W+1）＝Y（K））THENGOTO430ELSEGOSUB415∶RETURN

430 IFIP（P，IW）＝X（K）ORIP（P，IW）＝Y（K）THENDE＝IP（P，IW+1）∶GOSUB410∶RETURNELSEDE＝IP（P，IW）∶GOSUB410∶RETURN

435 IFTV$＝"Y"THENB＝3ELSEB＝2

440 J＝MY∶N＝0∶A＝361∶C＝A∶FORI＝1TONF∶N＝N+1∶C＝C+3∶PRINT@C，PEEK（J）;∶J＝J+3∶IFN＝BTHENC＝C+OELSEIFN＝6THENN＝0∶A＝A+64∶IFA

＝937THENA＝361∶C＝A∶INPUTG$∶PRINT@C，"CO.";ELSEC＝AELSEGOTO445

445 NEXTI∶C＝C+3∶PRINT@C，" ";∶RETURN

450 CLS∶N＝0

455 PRINT"TARGET′S NO.（K/NR）∶";K;"/";NR（K）

460 PRINT"NO";" X（MM）"，" Y（MM）"，" R（MM）"，∶IFTV$＝"Y"THENIFFT$（K）＝"Y"THENPRINT" T（MS）";ELSEPRINT"V（M/S）";

465 PRINTTAB（58）"SCORE"∶PRINTSTRING$（63，"-"）

470 I＝N＊10+1 ∶J＝N＊10+10

475 FOR M＝I TO J

480 PRINT M;" ";XT（NR（K），M）＊10，YT（NR（K），M）＊10，RI（M）＊10，∶IFTV$＝"Y"PRINTT（M，NR（K））;

485 PRINTTAB（59）NI（NR（K），M）

490 IF M＝NM THEN 520

495 NEXT M

500 IFPL$<>"W"INPUT"PRESS THE ′ENTER′ TO CONTINUE";G$

505 N＝N+1

510 CLS

515 GOTO 460

520 PRINTTAB（52）"TOTAL∶";NA（K）∶PRINTSTRING$（63，"-"）∶RETURN

525 RETURN

530 CLS∶PRINTTAB（0）"K";TAB（4）"EX";TAB（11）"EY";TAB（18）"R50";TAB（25）"R100";TAB（33）"XO";TAB（40）"YO";TAB（47）"MEAN（T/V）";TAB（57）"RMS"∶PRINTSTRING$（63，"-"）

535 FORK＝1TOIT∶PRINTTAB（0）K;" ";TAB（4）EX（K）;TAB（11）EY（K）TAB（18）R50（K）;TAB（25）RH（K）;TAB（33）S1（K）;TAB（40）S2（K）;TAB（47）TM（K）;TAB（57）INT（TS（K）＊100）/100∶NEXTK∶PRINTSTRING$（63，"-"）

540 IFPL$<>"W"THENINPUT"PRESS <EN> TO CONTIUE";G$∶PC$＝"V"∶XY＝0∶GOSUB2455∶PC$＝"T"∶XY＝1∶GOSUB3225∶PC$＝"Y"∶GOSUB2720∶GOSUB2445∶PC$＝"X"∶GOSUB2720∶GOSUB2450∶PRINT∶INPUT"PRESS <EN> TO CONTIUE";G$∶RETURNELSERETURN

545 IFTV$＝"N"RETURN

550 CLS∶PRINT"＊＊ TOTAL FLYING TIME ＊＊＊"∶FORM＝1TONM∶T＝0∶FORK＝1TOIT∶IFFT$（K）＝"Y"THENTT＝T（M，NR（K））ELSEIFT（M，NR（K））＝0THENGOTO560ELSETT＝INT（S（K）/T（M，NR（K））＊1E6）/1E3

555 T＝T+TT

560 NEXTK∶TT（M）＝T∶NEXTM

565 IFNM<＝1MT＝TT（NM）∶ST＝0∶GOTO575

570 S1＝0∶SX＝0∶FORM＝1TONM∶S1＝S1+TT（M）∶SX＝SX+TT（M）＊TT（M）∶

NEXTM∶SX＝NM＊SX-S1＊S1∶MT＝INT（S1/NM＊1E3）/1E3∶ST＝INT（SQR（SX/NM/（NM-1））＊1E3）/1E3

575 TT＝0∶FORK＝1TOIT∶TT＝TT+S（K）∶NEXTK∶TT（0）＝TT

580 PRINT"TOTAL DISTANCE∶";TT;"（M）"∶PRINT"M"，"T（MS）"∶PRINT STRING$834，"-"）

585 FORM＝1TONM∶PRINTM，TT（M）∶NEXTM∶PRINTSTRING$（34，"-"）∶PRINT"MEAN∶"，MT∶PRINT"RMS∶"，STMEAN∶"，MT∶PRINT"RMS∶"，ST

590 IFPL$＝"N"INPUT"PRESS <EN> TO CONTINUE"G$∶RETURN

595 RETURN

600 REM

605 N1＝63

610 SET（N1，26）∶SET（N1，27）∶SET（N1，28）

615 SET（N1-1，27）∶SET（N1+1，27）

620 FOR I＝N1-1 TO N1+1

625 SET（I，O）

630 NEXT I

635 FOR I＝1 TO 24

640 IF I<＝3 THEN X＝INT（I＊17/3）+1∶GOTO 660

645 IF I>3 AND I<12 THEN X＝INT（（I-3）/3）+18 ∶GOTO 660

650 IF I>＝12 AND I<＝16 THEN X＝INT（（16-I）/2）+19 ∶GOTO

660

655 X＝INT（29＊（I-16）/8）+16

660 GOSUB 830

665 NEXT

670 X＝45

675 FOR I＝25 TO 46

680 GOSUB 830

685 NEXT I

690 FOR I＝N1-45 TO N1+44

695 SET（I，46）

700 NEXT I

705 RETURN

710 X（9）＝FIX（X（9）＊1.8）+N1

715 Y（9）＝27-FIX（Y（9）＊0.9）

720 IF X（9）<0 OR X（9）>127 THEN 775

725 IF Y（9）<0 OR Y（9）>47 THEN 775

730 IF POINT（X（9），Y（9）） THEN 795

735 FOR I＝1 TO 3

740 SET（X（9），Y（9））

745 FOR J＝1 TO 30∶ NEXT J

750 RESET（X（9），Y（9））

755 FOR J＝1 TO 30∶NEXT J

760 NEXT I

765 SET（X（9），Y（9））

770 RETURN

775 FOR I＝1 TO 100

780 PRINT@192，"EXCEED CRT RANGE！！";

785 NEXT I

790 PRINT@256，" ";∶RETURN

795 FOR I＝1 TO 3

800 RESET（X（9），Y（9））

805 FOR J＝1 TO 30∶NEXT J

810 SET（X（9），Y（9））

815 FOR J＝1 TO 30∶NEXT J

820 NEXT I

825 RETURN

830 Y＝N1-X∶XX＝N1+X

835 SET（Y，I）

840 SET（XX，I）

845 RETURN

850 IFAJ$＝"Y"J＝1∶RETURN∶REM

855 AJ＝0∶J＝USR（AJ）∶IFNUM＝1THENIFJ＝OTHENNM＝M-1∶M＝PEEK（F5）∶RETURNELSEGOTO865ELSENF＝PEEK（AF）+PEEK（AF+1）＊256∶NUM＝INT（NF/TU/IT）∶NM＝NUM

860 IFNF＝ORETURN

865 J＝MY∶PN＝O∶PF＝O∶CC＝O∶FORI＝1TONF∶GOSUB1045

870 IF IP>＝21 GOTO 985

875 IFIP>16TTHENPRINT@128，"IP＝";IP;"BAD！！（1-16）"∶INPUT"<EN>";G$

880 IFIT<>1ANDNUM<>1THENGOTO885ELSEPN＝PN+1∶IP（M，PN）＝IP∶IFNUM＝1THENN（IP）＝M∶GOTO890ELSE885

885 N（IP）＝N（IP）+1

890 ON IP GOTO900，905，910，915，920，925，930，935，940，945，950，955，960，965，970，975

895 GOTO1040

900 XT（1，N（1））＝IX＊U∶GOTO1040

905 YT（1，N（2））＝IX＊U∶GOTO1040

910 XT（2，N（3））＝IX＊U∶GOTO1040

915 YT（2，N（4））＝IX＊U∶GOTO1040

920 XT（3，N（5））＝IX＊U∶GOTO1040

925 YT（3，N（6））＝IX＊U∶GOTO1040

930 XT（4，N（7））＝IX＊U∶GOTO1040

935 YT（4，N（8））＝IX＊U∶GOTO1040

940 XT（5，N（9））＝IX＊U∶GOTO1040

945 YT（5，N（10））＝IX＊U∶GOTO1040

950 XT（6，N（11））＝IX＊U∶GOTO1040

955 YT（6，N（12））＝IX＊U∶GOTO1040

960 XT（7，N（13））＝IX＊U∶GOTO1040

965 YT（7，N（14））＝IX＊U∶GOTO1040

970 XT（8，N（15））＝IX＊U∶GOTO1040

975 YT（8，N（16））＝IX＊U∶GOTO1040

980 GOTO1040

985 IFIP<29THENIFNUM＝1THENN（IP）＝MELSEN（IP）＝N（IP）+1ELSEP RINT@128，"IP＝";IP;"BAD！！（21-28）"∶INPUT"<EN>";G$

990 IFIT<>1ANDNUM<>1THEN995ELSEIFIT＝1ANDNUM＝1THENPF＝1∶PF（PF）＝IPELSEPF＝PF+1∶PF（PF）＝IP

995 ON IP-20 GOTO 1000，1005，1010，1015，1020，1025，1030，1035

1000 T（N（21），1）＝FT＊IX∶GOTO1040

1005 T（N（22），2）＝FT＊IX∶GOTO1040

1010 T（N（23），3）＝FT＊IX∶GOTO1040

1015 T（N（24），4）＝FT＊IX∶GOTO1040

1020 T（N（25），5）＝FT＊IX∶GOTO1040

1025 T（N（26），6）＝FT＊IX∶GOTO1040

1030 T（N（27），7）＝FT＊IX∶GOTO1040

1035 T（N（28），8）＝FT＊IX∶GOTO1040

1040 NEXT I∶RETURN

1045 REM

1050 IP＝PEEK（J）∶J＝J+1∶N＝PEEK（J）∶J＝J+1∶NL＝PEEK（J）∶J＝J+1

1055 IFIP<17THENIFN<8THENIX＝N＊256+NL∶GOTO1065ELSEIX＝-1＊（256＊（N-8）+NL）∶GOTO1070ELSEIX＝N＊256+NL∶GOTO1080

1060 IX＝N＊256+NL∶GOTO1075

1065 IFIX>＝1792THENGOSUB1085∶IX＝9999.9/U∶RETURNELSERETURN

1070 IFIX<＝-1792THENGOSUB1085∶IX＝-9999.9/U∶RETURNELSERETURN

1075 REM

1080 IFIX>＝3840THENGOSUB1085∶IX＝99999/FT∶RETURNELSERETURN

1085 LN＝1∶PRINT@192，IP;IX∶RETURN

1090 REM

1095 X＝ABS（X（9））

1100 IF X>A1 OR Y（9）>B1 OR Y（9）<B2 THEN 1180

1105 IF Y（9）<＝YY THEN 1165

1110 IF Y（9）<＝L（3） AND X<＝X4 THEN 1165

1115 FOR I＝1 TO 3.1

1120 IF Y（9）<＝L（I） THEN 1140

1125 NEXT I

1130 IF X>X4 THEN 1180

1135 IF SQR（X＊X+（Y（9）-D5）＊（Y（9）-D5））>R1 THEN 1180 ELSE

1165

1140 C1＝C（I）∶D1＝D（I）∶Y＝C1＊X+D1

1145 ON I GOTO 1150，1155，1160

1150 IF Y（9）<＝Y OR X<＝X2 THEN 1165 ELSE 1180

1155 IF Y（9）>＝Y OR X<＝X2 THEN 1165 ELSE 1180

1160 IF Y（9）<＝Y OR X<＝X4 THEN 1165 ELSE 1180

1165 N＝10-INT（（R-0.4）/5）∶NA（K）＝NA（K）+N

1170 PRINT@237，"N＝";N;∶PRINT@256，" ";

1175 RETURN

1180 N＝0∶PRINT@237，"N＝";N;

1185 FOR I＝1 TO 10

1190 PRINT@128，"OUT OF TARGET ！！";

1195 NEXT I

1200 PRINT@256，" ";∶RETURN

1205 REM

1210 IF NM＝1 THEN R＝0∶RI（NM）＝0∶S1＝XT（NR（K），1）∶S2＝YT（NR（K），1）∶EX（K）＝0∶EY（K）＝0∶GOTO1370

1215 S1＝0∶S2＝0∶SX＝0∶SY＝0∶CA＝0∶FORM＝1TONM

1220 IFNUM＝1THENGOTO1230ELSEIFLN<>1THENGOTO1230ELSEGOTO1225

1225 IFABS（XT（NR（K），M））<>9999.9ANDABS（YT（NR（K），M））<>9999.9THENGOTO1230ELSECA＝CA+1∶RI（M）＝0∶GOTO1235

1230 S1＝S1+XT（NR（K），M）∶SX＝SX+XT（NR（K），M）＊XT（NR（K），M）∶S2＝S2+YT（NR（K），M）∶SY＝SY+YT（NR（K），M）＊YT（NR（K），M）

1235 NEXTM∶NM＝NM-CA∶SX＝NM＊SX-S1＊S1∶IFSX<OTHENEX（K）＝0∶GOTO1245

1240 EX（K）＝.6745＊SQR（SX/NM/（NM-1））∶EX（K）＝INT（EX（K）＊100+.5）/10

1245 SY＝NM＊SY-S2＊S2∶IFSY<OTHENEY（K）＝0∶GOTO1255

1250 EY（K）＝.6745＊SQR（SY/NM/（NM-1））∶EY（K）＝INT（EY（K）＊100+.5）/10

1255 S1＝S1/NM ∶S2＝S2/NM∶NM＝NM+CA

1260 CA＝0∶FOR M＝1 TO NM

1265 IFNUM＝1THENGOTO1275ELSEIFLN<>1THENGOTO1275ELSEGOTO1270

1270 IFABS（XT（NR（K），M））<>9999.9ANDABS（YT（NR（K），M））<>9999.9THENGOTO1275ELSECA＝CA+1∶RI（M）＝0∶GOTO1280

1275 RI（M）＝SQR（（XT（NR（K），M）-S1）＊（XT（NR（K），M）-S1）+（YT（NR

（K），M）-S2）＊（YT（NR（K），M）-S2））

1280 NEXT M

1285 FOR L＝1 TO INT（NM/2）

1290 N＝NM-L+1

1295 FOR M＝L+1 TO N

1300 J＝M-1

1305 IF RI（L）<RI（M） THEN 1315

1310 X＝RI（L）∶RI（L）＝RI（M）∶RI（M）＝X

1315 IF RI（N）>RI（J） THEN 1325

1320 X＝RI（N）∶RI（N）＝RI（J）∶RI（J）＝X

1325 NEXT M

1330 NEXT L

1335 NM＝NM-CA∶N＝NM/2-INT（NM/2）

1340 IF N<>0 THEN 1360

1345 N＝NM/2 ∶J＝N+1

1350 R＝（RI（N）+RI（J））/2

1355 GOTO 1370

1360 N＝INT（NM/2）+1

1365 R＝RI（N）

1370 R＝INT（R＊100+.5）/10∶R50（K）＝R

1375 NM＝NM+CA∶RI（NM）＝INT（RI（NM）＊100+.5）/10∶RH（K）＝RI（NM）

1380 S1（K）＝INT（S1＊100+.5）/10∶S2（K）＝INT（S2＊100+.5）/10∶IFCT$＝"Y"XO（K）＝S1∶YO（K）＝S2

1385 IFTV$＝"Y"THENTM＝0∶SM＝0∶FORM＝1TONM∶TM＝TM+T（M，NR（K））∶NEXTM∶TM（K）＝TM/NM∶FORM＝1TONM∶SM＝SM+（T（M，NR（K））-TM（K））[2∶NEXTM∶IFNM>1THENIFSM<>0THENTS（K）＝SQR（SM/（NM-1））ELSETS（K）＝0ELSETS（K）＝0ELSEGOTO1390

1390 PRINT"MEAN∶";S1（K）;"MM"，S2（K）;"MM"，，∶IFTV$＝"Y"THENGOSUB1415∶IFFT$（K）＝"Y"THENPRINTTM（K）;"MS"ELSEPRINTTM（K）;"M/S"ELSEPRINT" "∶GOTO1400

1395 IFTV$＝"Y"PRINT"RMS∶";，，，"+/-";TS（K）;"MS"

1400 PRINTSTRING$（63，"-"）∶PRINT"X ∶"EX（K）"MM"，，"EY∶"EY（K）"MM"

1405 PRINT"R50∶";R50（K）;"MM"，，"R100∶";RH（K）;"MM"

1410 IFPL$＝"N"THENINPUT"PRESS <EN> TO CONTINUE";G$∶RETURNELSERETURN

1415 TM（K）＝INT（TM（K）＊1E4+.5）/1E4∶TS（K）＝INT（TS（K）＊1E4+.5）/1E4∶RETURN

1420 CLS∶PRINT" WX4675 PLOT.SUB."∶IFIT＝1TB$＝"N"∶GOTO1430∶REM

1425 TB$＝"Y"∶INPUT"＊＊ DO YOU WANT ONLY TOTAL TABLE（<EN>＝YES N＝NO）";TB$∶IFTB$＝"Y"OP$＝"N"∶LB$＝"N"∶LT$＝"Y"∶INPUT

" ＊DO YOU WANT LAST TARGET′S DATA TABLE（<EN>＝YES N＝NO）";LT$∶GOTO1440

1430 OP$＝"Y"∶INPUT"＊＊ DO YOU WANT ONLY THE POINT（<EN>＝YES N＝NO）";OP$

1435 LB$＝"Y"∶INPUT"＊＊ DO YOU WANT THE POINT MARK（<EN>＝YES N＝NO）";LB$

1440 WP$＝"Y"∶INPUT"＊＊ DO YOU WANT THE WAVE-PACKED（<EN>＝YES N＝NO）";WP$∶IF WP$＝"Y"AND IT<>1 TC$＝"Y"∶INPUT"＊＊ DO YOU PLOTER WITH TWO-CURVE（<EN>＝YES N＝NO）";TC$

1445 H＝0∶INPUT"I＊TARGET CENTRE-MUZZLE HIGHT DIFFERENCE（+/-） （MM）";H

1450 LPRINT"J"1∶GOSUB1535

1455 FORK＝1TOIT∶SS＝0∶FORM＝1TOK∶SS＝SS+S（M）∶NEXTM∶XX（K）＝BC＊SS∶GOSUB2285∶GOSUB2250∶GOSUB2300

1460 IFIT<＝2SN＝3∶ZZ＝0.9∶LPRINT"S2"

1465 IFTB$＝"Y"THENIFK＝ITTHENGOTO1525ELSEX＝0∶Y＝1300∶GOSUB2250∶GOTO1475ELSEGOSUB1660∶IFOP$＝"N"THENGOSUB1825ELSEGOTO1470

1470 GOSUB1935

1475 GOSUB1575∶GOSUB1735∶NEXTK

1480 IFTB$＝"N"ORLT$＝"Y"GOSUB2300∶GOSUB1885

1485 PC$＝"Y"∶XY＝1∶GOSUB1765∶IFWP$＝"Y"K＝K+1∶GOSUB2300∶PC$＝"Y"∶XY＝1∶GOSUB1790

1490 IFTB$＝"Y"PC$＝"Y"∶XY＝1∶GOSUB2720∶GOSUB2315

1495 GOSUB1590∶GOSUB1615

1500 GOSUB2395∶IFPP$＝"Y"THENGOSUB1590ELSEGOTO1505

1505 XY＝0∶PC$＝"V"∶GOSUB2455∶IFSD$<>"Y"GOTO1520

1510 XY＝1∶PC$＝"T"∶GOSUB3225

1515 GOSUB2990∶GOSUB1590

1520 GOSUB2280∶GOSUB2385∶RETURN

1525 IFLT$＝"Y"N＝1∶GOSUB2290∶GOSUB1660∶GOSUB1825∶GOSUB1935∶OP$＝"N"∶GOTO1475

1530 OP$＝"Y"∶GOTO1475

1535 AA＝INT（345/（IT+.5））＊10∶BC＝0∶FORK＝1TOIT∶BC＝BC+S（K）∶NEXTK∶BC＝AA＊IT/BC∶X＝0∶Y＝1300∶GOSUB 2250∶LPRINT"B15"∶FOR K＝1TOIT∶SS＝INT（BC＊S（K））∶IFK＝1THENSL＝0∶SC＝SS∶GOSUB1565∶GOTO1545

1540 IFSS>＝（AA/2-25）THENSL＝1∶SC＝AA/2-25∶GOSUB1565∶SL＝0∶SC＝SS-AA/2+25∶GOSUB1565ELSESL＝1∶SC＝SS∶GOSUB1565.

1545 GOSUB1570∶NEXTK∶LPRINT"L0"∶IFABS（S1（IT））>ABS（S2（IT））THENRH＝ABS（S1（IT））ELSERH＝ABS（S2（IT））

1550 TW＝130/2.2/（RH（IT）+RH）

1555 IFAA>＝1300THENKK＝1∶AA＝1300∶BB＝0ELSEKK＝（AA-50）/1300∶BB＝80

1560 RETURN

1565 LPRINT"L"SL∶Y＝0∶X＝SC∶GOSUB2255∶RETURN

1570 X＝0∶Y＝20∶GOSUB2255∶Y＝-40∶GOSUB2255∶Y＝20∶GOSUB2255∶RETURN

1575 X＝XX（K）-AA/2+25∶Y＝1950∶GOSUB2250

1580 Y＝650+BB∶GOSUB2245∶X＝XX（K）+AA/2-25∶Y＝650∶GOSUB2245∶Y＝1950-BB∶GOSUB2245∶X＝XX（K）-AA/2+25∶Y＝1950∶GOSUB2245

1585 RETURN

1590 LPRINT"S3"∶X＝0∶Y＝2555∶GOSUB2250∶LPRINT"J1"∶P$＝"TEST DATE∶"+DA$∶GOSUB2265

1595 Y＝Y-60∶GOSUB2250∶P$＝"TEST NO.∶"+STR$（NT）∶GOSUB2265

1600 Y＝Y-60∶GOSUB2250∶P$＝"GUN NAME∶"+GT$∶GOSUB2265

1605 Y＝Y-60∶GOSUB2250∶P$＝"GUN NO.∶"+GN$∶GOSUB2265

1610 RETURN

1615 LPRINT"S2"∶Y＝Y-60∶GOSUB2250∶P$＝"TW＝"+STR$（TW）+"KK＝"+STR$（KK）∶GOSUB2265

1620 Y＝Y-100∶GOSUB2250∶IFNUM＝1THENP$＝"SINGLE FIRE！！"ELSEP$＝"CONTINUE FIRE！！"

1625 GOSUB2265∶Y＝2050∶GOSUB2250∶P$＝"VERTICAL PLOT"∶GOSUB2265

1630 X＝10∶Y＝1250∶GOSUB2250∶P$＝"H＝"+STR$（H）+"MM"∶GOSUB2265∶RETURN

1635 PRINT""∶PRINT"＊＊＊ RANDOM GENERATE X Y T，JUST A MOMENT ＊＊＊"

1640 XX＝1792/6＊2∶YY＝3840

1645 FORM＝1TONM∶XT（NR（IT），M）＝（RND（XX）-XX/2）＊U

∶YT（NR（IT），M）＝（RND（XX）-XX/2）＊U

∶T（M，NR（IT））＝RND（YY）＊

FT ∶NEXTM

1650 IFIT＝1RETURN∶REM

1655 FORK＝1TO（IT-1）∶KP＝K/IT∶FORM＝1TONM∶XT（NR（K），M）＝INT（XT（NR（IT），M）＊KP＊100）/100∶YT（NR（K），M）＝INT（YT（NR（IT），M）＊KP＊100）/100∶T（M，NR（K））＝INT（T（M，NR（IT））＊KP＊1E4）/1E4∶NEXTM，K∶RETURN

1660 REM

1665 X＝XO∶Y＝YO∶GOSUB 2250

1670 P$＝"K/NR∶"+STR$（K）+"/"+STR$（NR（K））+" DIS.B.T∶"+STR$（S（K））

1675 GOSUB 2265

1680 X＝X∶Y＝Y-15＊ZZ∶GOSUB 2250

1685 LPRINT"L1"

1690 LPRINT"B20"

1695 X＝X+.7＊SN＊100∶∶GOSUB 2245

1700 LPRINT"L0"

1705 X＝XO∶Y＝YO-70＊ZZ∶GOSUB 2250∶IF TV$＝"N" THEN T$＝" "ELSEIF FT$（K）＝"Y" THEN T$＝"T（MS）"ELSET$＝"V（M/S）"

1710 P$＝"M X（MM） Y（MM） "+T$+"SCORE"

1715 GOSUB 2265

1720 X＝XO∶Y＝YO-85＊ZZ∶GOSUB 2250

1725 X＝XO+.7＊SN＊340∶GOSUB 2245

1730 RETURN

1735 REM

1740 LPRINT"S1"∶C$＝"N"∶FOR M＝1TONM

1745 X＝XT（NR（K），M）＊100＊TW＊KK+XX（K）∶Y＝YT（NR（K），M）＊100＊TW+1300

1750 GOSUB 2250∶GOSUB 2275∶IFLB$＝"Y"THENGOTO1755ELSEIFK＝ITTHENGOTO1755ELSEGOTO1760

1755 P$＝STR$（M）∶GOSUB 2265

1760 NEXTM∶LPRINT"S5"∶X＝S1（K）＊10＊TW＊KK+XX（K）∶Y＝S2（K）＊10＊TW+1300∶GOSUB2250∶GOSUB2275∶LPRINT"S"SN-1∶RETURN

1765 IFIT＝1RETURN∶REM

1770 PRINT"＊＊＊ TO PLOTER MEAN "PC$"-S CURVE ＊＊＊"∶FORK＝1 TOIT∶ZZ（K）＝S1（K）＊10＊TW＊KK+XX（K）∶YY（K）＝S2（K）＊10＊TW+1300

1775 NEXTK∶ZZ（IT+1）＝0∶YY（IT+1）＝-H＊10＊TW+1300

1780 GOSUB2595∶IFUS$<>"Y"GOSUB2945

1785 US$＝"N"∶RETURN

1790 REM

1795 IF IT＝1THENGOSUB1810∶RETURN

1800 IF TC$＝"Y"THENGOSUB 2905 ELSEGOSUB 1810

1805 RETURN

1810 Y＝1300-H＊TW＊10∶YO＝Y∶K＝1∶FORM＝1TONM∶PRINT"＊＊＊ TO PLOTER "PC$"-S CURVE "M" ＊＊＊"∶IFYO<0ORYO>2600THENX＝XT（NR（1），M）∶Y＝YT（NR（1），M）∶K＝1∶GOSUB1820∶GOSUB2250ELSEX＝0∶Y＝YO∶GOSUB2250

1815 FORK＝1TOIT∶X＝XT（NR（K），M）∶Y＝YT（NR（K），M）∶GOSUB1820∶GOSUB2245∶NEXTK，M∶RETURN

1820 X＝X＊100＊TW＊KK+XX（K）∶Y＝Y＊100＊TW+1300∶RETURN

1825 REM

1830 FORM＝1TONM∶X＝XO∶Y＝YO-130＊ZZ-35＊ZZ＊（M-1）∶GOSUB 2250

1835 P$＝STR$（M）∶GOSUB 2265

1840 X＝XO+0.7＊SN＊40∶GOSUB 2250

1845 P$＝STR$（XT（NR（K），M）＊10）∶GOSUB 2265

1850 X＝X+0.7＊SN＊80∶GOSUB 2250

1855 P$＝STR$（YT（NR（K），M）＊10）∶GOSUB 2265

1860 X＝X+0.7＊SN＊80∶GOSUB 2250∶IF TV$＝"Y"THENGOTO1865ELSEGOTO1870

1865 P$＝STR$（T（M，NR（K）））∶GOSUB 2265

1870 X＝X+0.7＊SN＊90∶GOSUB 2250

1875 P$＝STR$（NI（NR（K），M））∶GOSUB 2265

1880 NEXTM∶RETURN

1885 IFTV$＝"N"ORNM＝1ORIT＝1RETURN

1890 K＝1∶GOSUB2285∶XO＝XO+520∶X＝XO∶Y＝YO∶GOSUB2250

1895 P$＝"DIS.∶"+STR$（TT（0））∶GOSUB2265

1900 X＝XO∶Y＝Y-70＊ZZ∶GOSUB2250∶P$＝"M T（MS）"∶GOSUB2265

1905 X＝XO∶Y＝YO-85＊ZZ∶GOSUB2250∶X＝X+.7＊SN＊110∶GOSUB2245

1910 IFOP$＝"Y"GOTO1925

1915 FORM＝1TONM∶X＝XO∶Y＝YO-130＊ZZ-35＊ZZ＊（M-1）∶GOSUB2250∶P$＝STR$（M）∶GOSUB2265∶X＝XO+.7＊SN＊40∶GOSUB2250

1920 P$＝STR$（TT（M））∶GOSUB2265∶NEXTM

1925 X＝XO∶Y＝YO-130＊ZZ-35＊ZZ＊NM∶Y＝Y-20＊ZZ∶GOSUB2250∶X＝X+.7＊SN＊110∶GOSUB2245∶X＝XO∶Y＝Y-35＊ZZ∶GOSUB2250∶P$＝"MEAN∶"∶GOSUB2265∶X＝XO+.7＊SN＊50∶GOSUB2250

1930 P$＝STR$（MT）∶GOSUB2265∶Y＝Y-35＊ZZ∶X＝XO∶GOSUB2250∶P$＝"RMS∶"∶GOSUB2265∶X＝XO+.7＊SN＊50∶GOSUB2250∶P$＝STR$（ST）∶GOSUB2265∶RETURN

1935 REM

1940 X＝XO∶Y＝YO-130＊ZZ-35＊ZZ＊NM∶GOSUB2250∶P$＝STRING$（23，32）+"TOTAL＝"+STR$（NA（K））∶GOSUB2265

1945 Y＝Y-20＊ZZ∶GOSUB2250∶X＝X+.7＊SN＊340∶GOSUB2245

1950 X＝XO∶Y＝Y-35＊ZZ∶GOSUB2250∶P$＝"MEAN∶"∶GOSUB2265

1955 X＝XO+.7＊SN＊50∶GOSUB2250∶P$＝STR$（S1（K））∶GOSUB2265

1960 X＝X+.7＊SN＊70∶GOSUB2250∶P$＝STR$（S2（K））∶GOSUB2265

1965 IF TV$＝"Y"THEN GOTO1970 ELSE X＝XO∶Y＝Y-35＊ZZ∶GOTO1980

1970 X＝X+.7＊SN＊80∶GOSUB2250∶P$＝STR$（TM（K））∶GOSUB2265

1975 Y＝Y-35＊ZZ∶X＝XO∶GOSUB2250∶P$＝"RMS∶"+STRING$（12，32）+"+/-"+STR$（TS（K））∶GOSUB2265

1980 Y＝Y-35＊ZZ∶GOSUB2250∶P$＝"EX＝"+STR$（EX（K））∶GOSUB2265

1985 X＝XO+.7＊SN＊150∶GOSUB2250∶P$＝"EY＝"+STR$（EY（K））∶GOSUB2265

1990 X＝XO∶Y＝Y-35＊ZZ∶GOSUB2250∶P$＝"R50＝"+STR$（R50（K））∶GOSUB2265

1995 X＝XO+.7＊SN＊150∶GOSUB2250∶P$＝"R100＝"+STR$（RH（K））∶GOSUB2265

2000 LPRINT"S"SN-1∶RETURN

2005 CLS∶PRINT"NF＝";PEEK（F3）+PEEK（F3+1）＊256，"DELAY＝";PEEK（DY），"TOTAL＝";PEEK（TL），"U＝";U∶PRINT"NF3＝";PEEK（F4）+PEEK（F4+1）＊256，"IT＝";PEEK（F1），"NUM＝";PEEK（F2），"NM＝";PEEK（F5），"ACNF＝";PEEK（AF）+PEEK（AF+1）＊256，"O′CLOCK＝";PEEK（SE）-1，"TV$＝";TV$

2010 PRINT∶IFTV$＝"Y"THENB＝3ELSEB＝2

2015 J＝MY∶N＝0∶A＝256∶C＝A∶FORI＝1TONF∶N＝N+1∶PRINT@C，PEEK（J）;∶J＝J+1∶C＝C+3∶PRINT@C，PEEK（J）;∶J＝J+1∶C＝C+3∶PRINT@C，PEEK（J）;∶J＝J+1∶C＝C+4∶IFN＝BTHENC＝C+1ELSEIFN＝6THENN＝0∶A＝A+64∶IFA＝960THENA＝256∶C＝A∶INPUTG$∶PRINT@C，"CO.";ELSEC＝AELSEGOTO2020

2020 NEXTI∶PRINT@C，"<END>";

2025 PRINT∶FORIP＝0TO28∶IFN（IP）＝0THENNEXTIPELSEPRINT"N（"IP"）＝"N（IP）;∶NEXTIP

2030 PRINT∶INPUT"COUNTINE DISPLY XT，YT PRESS <ENTER>";G$

2035 PRINT∶FORK＝1TOIT∶PRINT"NM＝";NM;" K＝";K;" NR（K）＝";NR（K）∶PRINT"M"，"XT（NR（K），M）"，"YT（NR（K），M）（CM）"，"T（M，NR（K））（MS）"∶FOR M＝1TONM∶PRINTM，XT（NR（K），M），YT（NR（K），M），T（M，NR（K））∶NEXTM∶NEXTK∶END

2040 CLS∶PRINT"DISOLAY & ADJUST MS6M12′S DELAY TIME SUB."∶REM

2045 B＝PEEK（-2011）∶IF B<>217 THENPRINT"NOT MS6M12！！"∶RETURNELSEGOSUB2070

2050 INPUT"INPUT DELAY TIME（0.04-516）";A

2055 IFA>5160RA<.04 THENPRINT"AGAIN";∶GOTO2050

2060 IFA<2 THENB＝1∶C＝INT（A＊1000/7.84-22/7.84-1）ELSEC＝0∶B＝INT（（A＊1000-22）/2014.9）∶IFB＝256THENB＝OELSEGOTO2065

2065 POKE-2042，C∶POKE-2041，B∶B＝PEEK（-2011）∶IFB<>217THENPOKE17130，0∶POKE17129，0∶RETURNELSERETURN

2070 C＝PEEK（-2042）∶IFC＝OTHENC＝256

2075 B＝PEEK（-2041）∶IFB＝OTHENB＝256

2080 A＝（7.84＊（C+1）＊B+22）/1000∶PRINT" DELAY TIME（MS）＝"A∶RETURN

2085 IFAJ$＝"Y"RETURN∶REM

2090 XX＝0∶PRINT@384，"＊＊ SELF DIAGNOSTIC ＊＊";∶PRINT@530，"＊＊＊ IN PROGRESS ＊＊＊"∶NF＝TU＊IT∶POKEF3，NF∶POKE

F3+1，0

2095 AJ＝1∶J＝USR（AJ）∶J＝PEEK（AF）∶IFJ<>NFTHENGOSUB2190∶GOTO2180

2100 M＝MY∶IP＝0∶FORN＝1TONF∶IP＝IP+PEEK（M）∶M＝M+3∶NEXTN

2105 IE＝0∶FORN＝1TOIT∶IE＝IE+X（N）+Y（N）∶NEXTN

2110 IFTV$＝"Y"THENFORN＝1TOIT∶IE＝IE+NF（N）∶NEXTN

2115 IFIP<>IETHENGOSUB2205∶GOTO2180

2120 IFIT>1GOTO2165

2125 T＝408.75∶PC＝INT（4＊T/DP＊100）/100∶IFTV$＝"Y"FC＝INT（4＊T/DF＊100）/100

2130 FP＝642∶PP＝770∶J＝MY∶FORM＝1TONF∶GOSUB1045

2135 IFIP>16THENGOTO2145

2140 IFIX>PC+6ORIX<PC-6THENGOSUB2235∶GOTO2155ELSENEXTM∶GOTO2160

2145 IFIP＝NF（1）THENIFIX>FC+60RIX<FC-6THENGOSUB2225∶GOTO2155ELSEGOTO2155ELSEGOTO2150

2150 IFIX>1GOSUB2225∶GOTO2180

2155 NEXTM

2160 IF FP<>6420RPP<>770THENGOTO2180

2165 PRINT@530，"＊＊＊ OK！ ＊＊＊";

2170 NF＝AU＊TU＊IT∶IFNF<256THENJ＝0∶I＝NFELSEJ＝1∶I＝NF-256

2175 POKEF3，I∶POKEF3+1，J∶RETURN

2180 XX＝XX+1∶IFXX<＝0GOSUB2095

2185 PRINT@530，"＊＊＊ BAD！ ＊＊＊";∶PRINT@915，"＊ CHECK IN DETAIL，PLEASE！＊";∶PRINT@851，"";∶END

2190 PRINT@578，"TOTAL ERROR";∶PRINT@645，"CORRECT∶";NF;∶PRINT@709，"REAL ∶";J;∶IFJ＝ORETURN

2195 M＝MY∶PRINT@725，"";∶FORN＝1TOJ∶PRINTPEEK（M）;∶M＝M+3∶IFN＝10THENPRINT"";∶PRINT@775，" ";

2200 NEXTN∶PRINT"";∶RETURN

2205 PRINT@578，"CHANNEL ERROR";∶PRINT@645，"CORRECT∶";∶FORN＝1TOIT∶PRINTX（N）;Y（N）;∶NEXTN

2210 IFTV$＝"Y"THENFORN＝1TOIT∶PRINTNF（N）;∶NEXTN

2215 PRINT∶PRINT"REAL∶";∶M＝MY∶FORN＝1TONF∶PRINTPEEK（M）;∶IFN＝11THENPRINT

2220 M＝M+3∶NEXTN∶RETURN

2225 PRINT@578，"FLYING VALUE ERROR"，"DF＝"4000/DF"KHZ""FC＝"FC;

2230 PRINT@FP，IP;IX;"，";∶FP＝FP+9∶RETURN

2235 PRINT@706，"POSITION VALUE ERROR"，"DP＝"4/DP"MHZ""PC＝"PC;

2240 PRINT@PP，IP;IX;"，";∶PP＝PP+9∶RETURN

2245 LPRINT"D" INT（X）"，"INT（Y）∶RETURN

2250 LPRINT"M" INT（X）"，"INT（Y）∶RETURN

2255 LPRINT"I"INT（X）"，"INT（Y）∶RETURN

2260 LPRINT"R"INT（X）"，"INT（Y）∶RETURN

2265 LPRINT"P"P$∶RETURN

2270 LPRINT"X"INT（X）"，"INT（Y）"，"INT（Z）∶RETURN

2275 LPRINTC$ INT（Z）∶RETURN

2280 LPRINT"H"∶RETURN

2285 ZZ＝.6∶SN＝2∶LPRINT"S"SN-1∶N＝K

2290 IFN<＝4THENYO＝2575∶XO＝660+720＊（N-1）ELSEYO＝635∶XO＝660+720＊（N-5）

2295 RETURN

2300 C$＝"J"

2305 IF K<＝5 THEN Z＝K+1∶GOSUB 2275ELSE Z＝K-5∶GOSUB 2275∶Z＝K-4

2310 C$＝"N"∶RETURN

2315 K＝2∶GOSUB2285∶XO＝XO+50∶X＝XO∶Y＝YO∶GOSUB2250∶LPRINT"S2"

2320 P$＝"K S（M） XO（MM） YO EX EY

R50 R100 T/V RMS "∶GOSUB2265

2325 X＝XO∶Y＝YO-35∶GOSUB2250∶X＝XO+1900∶GOSUB2245

2330 FORK＝1TOIT∶X＝XO∶Y＝YO-80-35＊（K-1）∶GOSUB2250∶P$＝STR$（K）+" "+STR$（S（K））∶GOSUB2265

2335 X＝XO+21＊13∶P$＝STR$（S1（K））∶GOSUB2250∶GOSUB2265∶P$＝STR$（S2（K））∶GOSUB2380∶P$＝STR$（EX（K））∶GOSUB2380∶P$＝STR$（EY（K））∶GOSUB2380∶P$＝STR$（R50（K））∶GOSUB2380∶P$＝STR$（RH（K））∶GOSUB2380

2340 IFTV$＝"Y"P$＝STR$（TM（K））∶GOSUB2380

2345 IFTV$＝"Y"THENIFFT$（K）＝"Y"THENP$＝"MS"∶GOSUB2380ELSEP$＝"M/S"∶GOSUB2380ELSEGOTO2350

2350 IFTV$＝"Y"P$＝STR$（TS（K））∶GOSUB2380

2355 NEXTK

2360 X＝XO∶Y＝YO-80-35＊IT∶GOSUB2250∶X＝XO+1900∶GOSUB2245

2365 X＝XO∶Y＝Y-35∶GOSUB2250∶P$＝"TOTAL∶"+STR$（TT（0））∶GOSUB2265∶IFTV$＝"Y"THENX＝XO+21＊58∶P$＝"TIME（MS）∶"∶GOSUB2250∶GOSUB2265ELSEGOTO2375

2370 IFTV$＝"Y" P$＝STR$（MT）∶GOSUB2380∶X＝X+21＊9∶P$＝STR$（ST）∶GOSUB2380

2375 RETURN

2380 X＝X+21＊9∶GOSUB2250∶GOSUB2265∶RETURN

2385 PRINT∶PRINT∶T$＝"N"∶INPUT"Y＝COPY<EN>＝END";T$∶IFT$＝

"Y"THENPRINT"PUT A NEW PAPER，PLEASE"∶INPUT"PRESS<EN>TO START COPY";G$∶GOTO1420ELSERETURN

2390 RETURN

2395 IFIT<>1THENPRINT""∶PP$＝"Y"∶PRINT∶INPUT"＊＊ DO YOU WANT THE PARTIAL FLOW PLOT（<EN>＝YES N＝NO）";PP$ELSERETURN

2400 IFPP$<>"Y"RETURN

2405 L＝0∶INPUT"＊ TARGET CENTRE-MUZZLE LEVEL DIFFERENCE（+/-） （MM）";L∶H＝L

2410 INPUT"＊＊ PUT A NEW PAPER，TO PRESS <EN> START PLOTER ＊＊";G$

2415 GOSUB2445

2420 X＝10∶Y＝1250∶GOSUB2250∶P$＝"L＝"+STR$（L）+"MM"∶GOSUB2265∶LPRINT"J1"∶GOSUB1535

2425 FORK＝1TOIT∶GOSUB2285∶GOSUB2300∶GOSUB1575∶GOSUB1735∶NEXTK∶PC$＝"X"∶XY＝1∶GOSUB1765∶IFWP$＝"Y"THENK＝K+1∶GOSUB2300∶PC$＝"X"∶XY＝1∶GOSUB1790

2430 PC$＝"X"∶XY＝1∶GOSUB2720

2435 X＝50∶Y＝2320∶GOSUB2250∶P$＝"TW＝"+STR$（TW）+"KK＝"+STR$（KK）∶GOSUB2265∶Y＝2050∶GOSUB2250∶P$＝"HORIZONTAL PLOT"∶GOSUB2265

2440 GOSUB2450∶RETURN

2445 FORK＝1TOIT∶A＝S1（K）∶S1（K）＝S2（K）∶S2（K）＝-A∶FORM＝1TONM∶A＝XT（NR（K），M）∶XT（NR（K），M）＝YT（NR（K），M）∶YT（NR（K），M）＝-A∶NEXTM，K∶RETURN

2450 FORK＝1TOIT∶A＝S2（K）∶S2（K）＝S1（K）∶S1（K）＝-A∶FORM＝1TONM∶A＝YT（NR（K），M）∶YT（NR（K），M）＝XT（NR（K），M）∶XT（NR（K），M）＝-A∶NEXTM，K∶RETURN

2455 REM

2460 IFTV$<>"Y"ORIT<4RETURN

2465 PRINT" "∶IFPL$＝"W"SD$＝"Y"∶INPUT"＊＊ DO YOU WANT THE SPEED AND DISTRIBUTION PLOT（<EN>＝YES N＝NO）";SD$∶IFSD$<>"Y"RETURN

2470 IFPL$＝"W"INPUT"＊＊＊ PUT A NEW PAPER，TO PRESS<EN>START PLOTER ＊＊＊";G$

2475 PRINT"＊＊＊ TO COMPUTE MEAN "PC$"-S CURVE ＊＊＊"

2480 FORK＝1TOIT∶IFFT$（K）＝"Y"THENIFTM（K）＝OTHENTV（K）＝0∶PRINT"SPEED＝0，DO NOT PLOTER V-S CURVE"∶RETURNELSETV（K）＝INT（S（K）/TM（K）＊1E5+.5）/100ELSETV（K）＝TM（K）

2485 NEXTK∶SS（1）＝S（1）/2∶FORK＝2TOIT∶SS＝0∶FORI＝1TOK-1∶SS＝SS+S（I）∶NEXTI∶SS（K）＝SS+S（K）/2∶NEXTK

2490 IFSS（1）<25THENFORK＝2TOIT∶ZZ（K-1）＝-SS（K）∶YY（K-1）＝LO

G（TV（K））∶NEXTK∶XY＝-1ELSEFORK＝1TOIT∶ZZ（K）＝-SS（K）∶YY（K）＝LOG（TV（K））∶NEXTK∶XY＝0

2495 GOSUB2755∶RETURN

2500 IFPL$<>"W"RETURN∶REM

2505 A＝1∶B＝2∶XO＝250∶YO＝300∶X＝XO∶Y＝YO∶LPRINT"S2"∶LPRINT"J1"∶GOSUB2250

2510 LPRINT"X1，100，12"∶X＝XO-50∶Y＝YO-40∶GOSUB2250∶FORI＝1TO5∶X＝X+B＊100∶GOSUB2250∶P$＝STR$（100＊I）∶GOSUB2265∶NEXTI

2515 X＝X+B＊150+50∶GOSUB2250∶P$＝"S（M）"∶GOSUB2265

2520 X＝XO∶Y＝YO∶GOSUB2250∶LPRINT"XO，100，12"∶X＝XO-120∶Y＝YO-10∶GOSUB2250

2525 FORI＝1TO5∶Y＝Y+A＊200∶GOSUB2250∶P$＝STR$（200＊I）∶GOSUB2265∶NEXTI

2530 Y＝Y+A＊150+110∶GOSUB2250∶P$＝"V（M/S）"∶GOSUB2265∶LPRINT"S3"∶X＝600∶Y＝140∶P$＝"SPEED"∶GOSUB2250∶GOSUB2265∶Y＝100∶P$＝"TIME"∶GOSUB2250∶GOSUB2265∶X＝750∶Y＝120∶P$＝"-DISTANCE PLOT"∶GOSUB2250∶GOSUB2265∶LPRINT"S2"

2535 GOSUB3215∶GOSUB3205

2540 LPRINT"J2"∶MM$＝"N"∶FORS＝0TO500STEP25∶V＝BA＊EXP（-BB＊S）∶Y＝V＊A+YO∶X＝S＊B+XO∶IFX>＝3450ORY>＝2600ORY<YOTHENMM$＝"Y"∶GOTO2550

2545 IFS＝0ORMM$＝"Y"THENGOSUB2250ELSEGOSUB2245∶MM$＝"N"

2550 NEXTS∶RETURN

2555 REM

2560 IFPL$<>"W"PRINTPC$"＝"A（1，4）"+"A（2，4）"＊S+"A（3，4）"＊S＊S"∶PRINT"R＝"RR，"RMS＝"RS∶FORS＝0TOTT（0）STEP50∶PRINTPC$;S;"＝";INT（A（3，4）＊S[2）＊10+.5/10，∶NEXTS∶PRINT""∶RETURN

2565 GOSUB2830

2570 IFPC$＝"V"THENK＝7∶GOSUB2285∶X＝XO-450∶XO＝X∶Y＝YO-50ELSEK＝8∶GOSUB2285∶X＝XO-300∶XO＝X∶Y＝YO-50

2575 GOSUB2300∶LPRINT"S2"∶P$＝PC$+"＝"+STR$（A（1，4））+"+"+STR$（A（2，4））+"＊S+"+STR$（A（3，4））+"＊S＊S"∶GOSUB2250∶GOSUB2265

2580 Y＝Y-50∶X＝XO+200∶P$＝"R＝"+STR$（RR）+"RMS＝"+STR$（RS）∶GOSUB2250∶GOSUB2265∶Y＝Y-35∶FORS＝0TOTT（0）STEP50∶P$＝PC$+STR$（S）+"＝"+STR$（INT（A（3，4）＊S[2＊10+.5）/10）∶Y＝Y-35∶GOSUB2250∶GOSUB2265

2585 IFS＝500X＝X+340∶Y＝Y+35＊11∶GOSUB2250

2590 NEXTS∶RETURN

2595 REM

2600 US$＝"N"∶X5＝0∶X6＝0∶X7＝0∶X8＝0∶Y1＝0∶Y2＝0∶Y3＝0∶Y4＝0

2605 FORK＝1TOIT+XY∶X5＝X5+ZZ（K）∶X6＝X6+ZZ（K）[2∶X7＝X7+ZZ（K

）[3∶X8＝X8+ZZ（K）[4∶Y1＝Y1+YY（K）∶Y2＝Y2+YY（K）＊ZZ（K）∶Y3＝Y3+YY（K）＊ZZ（K）＊ZZ（K）∶Y4＝Y4+YY（K）[2∶NEXTK

2610 A（1，1）＝IT+XY∶A（1，2）＝X5∶A（1，3）＝X6∶A（1，4）＝Y1∶A（2，1）＝X5∶A（2，2）＝X6∶A（2，3）＝X7∶A（2，4）＝Y2∶A（3，1）＝X6∶A（3，2）＝X7∶A（3，3）＝X8∶A（3，4）＝Y3∶Q1＝A（2，4）-A（1，2）＊A（1，4）/A（1，1）∶Q2＝A（3，4）-A（1，3）＊A（1，4）/A（1，1）∶QY＝Y4-A（1，4）[2/A（1，1）

2615 FORJ＝1TO3

2620 IFA（J，J）<>0GOTO2655

2625 FORI＝JTO3

2630 IFA（I，J）<>0GOTO2650

2635 NEXTI

2640 PRINT"NO UNIQUE SOLUTION ！"∶US$＝"Y"∶FORI＝1TO3∶PRINT"A（"I"，"J"）＝"A（I，J），∶NEXTI

2645 RETURN

2650 GOSUB2705

2655 C＝1/A（J，J）

2660 GOSUB2710

2665 FORI＝1TO3

2670 IFI＝JGOTO2685

2675 C＝-A（I，J）

2680 GOSUB2715

2685 NEXTI，J

2690 UU＝A（2，4）＊Q1+A（3，4）＊Q2∶Q＝QY-UU∶IFQ<OTHENRS＝QELSERS＝SQR（Q/（IT-3+XY））

2695 IFQY＝0THENRR＝1ELSERR＝SQR（UU/QY）

2700 RETURN

2705 FORK＝JTO4∶B＝A（J，K）∶A（J，K）＝A（I，K）∶A（I，K）＝B∶NEXTK∶RETURN

2710 FORI＝JTO4∶A（J，I）＝A（J，I）＊C∶NEXTI∶RETURN

2715 FORS＝JTO4∶A（I，S）＝A（I，S）+C＊A（J，S）∶NEXTS∶RETURN

2720 IFIT<3RETURN∶REM

2725 PRINT"＊＊＊ TO COMPUTE MEAN"PC$"-S CURVE ＊＊＊"

2730 FORK＝1TOIT∶SS＝0∶FORI＝1TOK∶SS＝SS+S（I）∶NEXTI∶ZZ（K）＝SS

2735 IFPC$＝"Y"THENYY（K）＝S2（K）ELSEYY（K）＝-S2（K）

2740 NEXTK∶ZZ（IT+1）＝0∶IFPC$＝"Y"THENYY（IT+1）＝-HELSEYY（IT+1）＝H

2745 GOSUB2595∶IFUS$<>"Y"GOSUB2555

2750 US$＝"N"∶RETURN

2755 X5＝0∶X6＝0∶X7＝0∶X8＝0∶Y1＝0∶Y2＝0∶Y3＝0∶REM

2760 FORK＝1TOIT+XY∶X5＝X5+ZZ（K）∶X6＝X6+ZZ（K）[2∶X8＝X8+ZZ（K）＊YY（K）∶Y1＝Y1+YY（K）∶Y2＝Y2+YY（K）[2∶NEXTK∶X7＝X5/（IT+XY）∶Y

3＝Y1/（IT+XY）

2765 BB＝（X5＊Y1-（IT+XY）＊X8）/（X5[2-（IT+XY）＊X6）

2770 AA＝Y3-BB＊X7

2775 B＝Y2-（IT+XY）＊Y3[2

2780 C＝0∶FORK＝1TOIT+XY∶C＝C+（ZZ（K）-X7）＊（YY（K）-Y3）∶NEXTK∶C＝C＊BB

2785 IFB<CTHENCC＝0∶CV＝0ELSECC＝SQR（（B-C）/（IT+XY-2））∶CV＝CC

2790 DD＝SQR（C/B）∶BA＝EXP（AA）

2795 IFPL$<>"W"PRINTPC$;"＝";BA;"＊EXP（";-BB;"＊S）"∶PRINT"R＝"DD，"RMS＝"CV∶FORS＝0TOTT（0）STEP50∶PRINTPC$;S;"＝";INT（BA＊EXP（-BB＊S）＊10+.5）/10，∶NEXTS∶PRINT" "∶RETURN

2800 GOSUB2500

2805 XO＝700∶YO＝2100∶X＝XO∶Y＝YO∶LPRINT"S2"

2810 P$＝PC$+"＝"+STR$（BA）+"＊EXP（"+STR$（-BB）+"＊S）"∶GOSUB2250∶GOSUB2265

2815 Y＝Y-50∶X＝XO+200∶P$＝"R＝"+STR$（DD）+" RMS V＝"+STR$（CV）∶GOSUB2250∶GOSUB2265∶Y＝Y-35∶FORS＝0TOTT（0）+100STEP50∶P$＝PC$+STR$（S）+"＝"+STR$（INT（BA＊EXP（-BB＊S）＊10+.5）/10）∶Y＝Y-35∶GOSUB2250∶GOSUB2265

2820 IFS＝300X＝X+340∶Y＝Y+35＊7∶GOSUB2250

2825 NEXTS∶RETURN

2830 REM

2835 A＝1.25∶B＝1.4∶XO＝1050∶YO＝70∶X＝XO∶Y＝YO∶LPRINT"J1"∶LPRINT"S2"

2840 GOSUB2250∶LPRINT"X1，70，20"∶X＝XO-50∶Y＝YO-40∶GOSUB2250

2845 FORI＝1TO5∶X＝X+200＊B∶GOSUB2250∶P$＝STR$（200＊I）∶GOSUB2265∶NEXTI

2850 X＝X+B＊150∶GOSUB2250∶P$＝"S（M）"∶GOSUB2265∶X＝XO∶Y＝YO∶GOSUB2250∶LPRINT"XO，63，12"∶X＝XO-120∶Y＝YO-10∶GOSUB2250

2855 FORI＝1TO5∶Y＝Y+100＊A∶GOSUB2250∶IFPC$＝"Y"THENP$＝STR$（-300+100＊I）ELSEP$＝STR$（300-100＊I）

2860 GOSUB2265∶NEXTI∶Y＝YO+A＊300-10∶X＝XO-250∶GOSUB2250∶P$＝PC$+"（CM）"∶GOSUB2265

2865 LPRINT"J2"∶YO＝445∶FORK＝1TOIT+XY∶IFPC$＝"Y"THENY＝YY（K）/10＊A+YOELSEY＝YO-YY（K）/10＊A

2870 X＝ZZ（K）＊B+XO∶GOSUB2250∶LPRINT"N2"∶NEXTK∶LPRINT"J2"∶MM$＝"N"∶FORS＝0TO1000STEP50∶V＝A（1，4）+A（2，4）＊S+A（3，4）＊S＊S∶Y＝V/10＊A+YO∶X＝S＊B+XO∶IFPC$＝"Y"THENY＝V/10＊A+YOELSEY＝YO-V/10＊A

2875 X＝S＊B+XO∶IFX>3450ORY>2600ORY<70THENMM$＝"Y"∶GOTO288

5

2880 IFS＝0ORMM$＝"Y"THENGOSUB2250ELSEGOSUB2245∶MM$＝"N"

2885 NEXTS

2890 RETURN

2895 INPUT"A＝";A

2900 A＝A+1∶PRINTA"＝"PEEK（A）;∶INPUT"B＝";B∶POKEA，B∶PRINTA"＝"PEEK（A）∶GOTO2900

2905 REM

2910 FOR M＝1 TO NM∶PRINT"＊＊＊ TO PLOTER "PC$"-S CURVE"M"＊＊＊"

2915 FOR K＝1 TO IT

2920 ZZ（K）＝XT（NR（K），M）＊100＊TW＊KK+XX（K）

2925 YY（K）＝YT（NR（K），M）＊100＊TW+1300

2930 NEXTK∶ZZ（IT+1）＝0∶YY（IT+1）＝-H＊10＊TW+1300

2935 GOSUB2595∶IFUS$<>"Y"GOSUB2945

2940 US$＝"N"∶NEXTM∶RETURN

2945 A＝ZZ（IT）∶B＝A/30∶C＝XX（1）/2∶REM

2950 FORS＝0TOASTEPB∶GOSUB2985

2955 IFS<CTHEN2960ELSEGOTO2965

2960 IFY>2600ORY<0THENS＝S+B∶GOSUB2985∶MM$＝"Y"∶GOTO2960

2965 IFS＝0ORMM$＝"Y"THENGOSUB2250∶MM$＝"N"ELSEGOSUB2245

2970 NEXTS

2975 X＝ZZ（IT）∶Y＝YY（IT）∶GOSUB2245

2980 RETURN

2985 Y＝A（1，4）+A（2，4）＊S+A（3，4）＊S＊S∶X＝S∶RETURN

2990 IFIT<4RETURN∶REM

2995 A＝5∶B＝2∶XO＝1950∶YO＝300∶X＝XO∶Y＝YO∶GOSUB2250∶LPRINT"S2"∶LPRINT"J1"

3000 LPRINT"X1，100，12"∶X＝XO-50∶Y＝YO-40∶GOSUB2250∶FORI＝1TO5∶X＝X+B＊100∶GOSUB2250∶P$＝STR$（100＊I）∶GOSUB2265∶NEXTI

3005 X＝X+B＊150+50∶GOSUB2250∶P$＝"S（M）"∶GOSUB2265

3010 X＝XO∶Y＝YO∶GOSUB2250∶LPRINT"XO，100，12"∶X＝XO-90∶Y＝YO-10∶GOSUB2250

3015 FORI＝1TO5∶Y＝Y+A＊40∶GOSUB2250∶P$＝STR$（40＊I）∶GOSUB2265∶NEXTI∶Y＝Y+A＊40∶GOSUB2250∶P$＝"（MM）"∶GOSUB2265∶LPRINT"$3"∶X＝2200∶Y＝120∶GOSUB2250∶P$＝"DISTRIBUTION-DISTANCE PLOT"∶GOSUB2265∶LPRINT"S2"

3020 GOSUB3215

3025 XY＝0∶PC$＝"EX"∶LPRINT"J2"∶GOSUB3180∶GOSUB3050

3030 XY＝0∶PC$＝"EY"∶LPRINT"J3"∶GOSUB3185∶GOSUB3050

3035 XY＝0∶PC$＝"R50"∶LPRINT"J4"∶GOSUB3190∶GOSUB3050

3040 XY＝0∶PC$＝"R100"∶LPRINT"J5"∶GOSUB3195∶GOSUB3050

3045 RETURN

3050 IFIT<4RETURN

3055 PRINT"＊＊＊ TO COMPUTE "PC$"-S CURVE ＊＊＊"

3060 FORK＝1TOIT∶SS＝0∶FORI＝1TOK∶SS＝SS+S（I）∶NEXTI∶ZZ（K）＝SS

3065 IFPC$＝"EX"YY（K）＝EX（K）∶GOTO3085

3070 IFPC$＝"EY"YY（K）＝EY（K）∶GOTO3085

3075 IFPC$＝"R50"YY（K）＝R50（K）∶GOTO3085

3080 YY（K）＝RH（K）

3085 NEXTK

3090 GOSUB2595∶IFUS$<>"Y"GOSUB3100

3095 RETURN

3100 XO＝2400∶Y＝1600∶P$＝PC$

3105 IFPC$＝"EX"YO＝2500∶LPRINT"J2"∶X＝1800∶GOTO3125

3110 IFPC$＝"EY"YO＝2100∶LPRINT"J3"∶X＝1900∶GOTO3125

3115 IFPC$＝"R50"YO＝1700∶LPRINT"J4"∶X＝2000∶GOTO3125

3120 LPRINT"J5"∶X＝2100∶GOSUB2250∶GOSUB2265∶GOTO3150

3125 GOSUB2250∶GOSUB2265

3130 X＝XO∶Y＝YO∶P$＝PC$+"＝"+STR$（A（1，4））+"+"+STR$（A（2，4））+"＊S+"+STR$（A（3，4））+"＊S＊S"∶GOSUB2250∶GOSUB2265

3135 Y＝Y-50∶X＝XO+200∶P$＝"R＝"+STR$（RR）+"RMS＝"+STR$（RS）∶GOSUB2250∶GOSUB2265∶Y＝Y-35∶FORS＝50TOTT（0）+100STEP50∶P$＝PC$+STR$（S）+"＝"+STR$（INT（（A（1，4）+A（2，4）＊S+A（3，4）＊S＊S）＊10+.5）/10）∶Y＝Y-35∶GOSUB2250∶GOSUB2265

3140 IFS＝300X＝X+340∶Y＝Y+35＊7∶GOSUB2250

3145 NEXTS

3150 IFPC$＝"T"ORPC$＝"V"THENXO＝250∶YO＝300∶X5＝0ELSEXO＝1950∶YO＝300∶X5＝S（1）

3155 MM$＝"N"∶FORS＝X5TO520STEP25∶V＝A（1，4）+A（2，4）＊S+A（3，4）＊S＊S∶Y＝V＊A+YO∶X＝S＊B+XO

3160 IFX>＝3450ORY>＝2600ORY<YOTHENMM$＝"Y"∶GOTO3170

3165 IFS＝X5ORMM$＝"Y"THENGOSUB2250ELSEGOSUB2245∶MM$＝"N"

3170 NEXTS

3175 RETURN

3180 FORK＝1TOIT∶Y＝EX（K）＊A+YO∶GOSUB3200∶NEXTK∶RETURN

3185 FORK＝1TOIT∶Y＝EY（K）＊A+YO∶GOSUB3200∶NEXTK∶RETURN

3190 FORK＝1TOIT∶Y＝R50（K）＊A+YO∶GOSUB3200∶NEXTK∶RETURN

3195 FORK＝1TOIT∶Y＝RH（K）＊A+YO∶GOSUB3200∶NEXTK∶RETURN

3200 SS＝0∶FORI＝1TOK∶SS＝SS+S（I）∶NEXTI∶X＝SS＊B+XO∶GOSUB2250∶LPRINT"N1"∶RETURN

3205 FORK＝2TOIT∶Y＝TV（K）＊A+YO∶GOSUB3210∶NEXTK∶RETURN

3210 X＝SS（K）＊B+XO∶GOSUB2250∶LPRINT"N5"∶RETURN

3215 X＝TT（0）＊B+XO∶Y＝1100∶GOSUB2250∶Y＝Y-400∶GOSUB2245∶X＝X-25∶Y＝Y+25∶GOSUB2250∶X＝X+25∶Y＝Y-25∶GOSUB2245∶X＝X+25∶Y＝Y+25∶GOSUB2250∶X＝X-25∶Y＝Y-25∶GOSUB2245∶RETURN

3220 FORK＝1TOIT+1∶Y＝ST（K）＊A+YO∶GOSUB3210∶NEXTK∶RETURN

3225 IFTV$<>"Y"ORIT<5RETURN∶REM

3230 PRINT"＊＊＊ TO COMPUTE MEAN "PC$"-S CURVE ＊＊＊"∶ST（0）＝0∶SS（0）＝0∶FORK＝1TOIT∶IFFT$（K）＝"Y"THENST（K）＝TM（K）ELSEST（K）＝INT（S（K）/TM（K）＊1E6）/1E6

3235 ST（K）＝ST（K-1）+ST（K）∶SS（K）＝SS（K-1）+S（K）∶NEXTK∶ST（IT+1）＝0∶SS（IT+1）＝0

3240 IFSS（1）<50THENFORK＝2TOIT+XY∶ZZ（K-1）＝SS（K）∶YY（K-1）＝ST（K）∶NEXTK∶XY＝0ELSEFORK＝1TOIT+XY∶ZZ（K）＝SS（K）∶YY（K）＝ST（K）∶NEXTK∶XY＝1

3245 GOSUB2595∶IFUS$＝"Y"RETURN

3250 IFPL$<>"W"THENPRINTPC$"＝"A（1，4）"+"A（2，4）"＊S+"A（3，4）"＊S＊S"∶PRINT"R＝"RR，"RMS＝"RS∶FORS＝0TOTT（0）STEP50∶PRINTPC$;S;"＝";INT（（A（1，4）+A（2，4）＊S+A（3，4）＊S＊S）＊1E6+.5）/1E6，∶NEXTS∶RETURN

3255 A＝1∶B＝2∶XO＝250∶YO＝300∶X＝0∶Y＝1300∶GOSUB2250∶LPRINT"J3"

3260 GOSUB3220

3265 XO＝700∶YO＝1700∶X＝330∶Y＝1555∶P$＝"T（MS）"∶GOSUB3125∶RETURN

3270 END∶REM

Claims (6)

1, a kind of method that sound wave arrives definite each parameter of pill impacting of mistiming of different acceptance points of measuring of utilizing, be to adopt near the sound target plane of muzzle, placing an acoustoelectric sensor and forming by two orthogonal sounding rods, and an acoustoelectric sensor is installed respectively at the two ends of every rod, wherein a sounding rod is to place horizontal level, another root is to be in position, plane vertically, this target plane placed on the trajectory and with ballistic path direction in vertical state, when bullet when passing through target plane, four acoustoelectric sensors measure the time that sound wave arrives each sensor respectively, after these electric signal are transported to time difference tester, can obtain several mistimings of horizontal direction and vertical direction, again with these mistiming electric signal input computing machines, by analysis, just can record each parameter of pill impacting after the calculating, the invention is characterized in, flight bullet trajectory on the different distance of muzzle, can place simultaneously a plurality of same sound target planes, like this, when bullet passes through each target plane, several acoustoelectric sensors on each target plane can measure the time that acoustic signals arrives constantly in difference, these acoustical signals convert electric signal to after impedance transformer (11) through sensor, send into channel amplifier (6), each electric signal after processing and amplifying is by Chuan Shu Long line (7) feed-in time difference tester (30), can obtain each time difference value thus, again these time difference signals are sent into microprocessor (19), under predefined programmed control, carry out various analytical calculations, its result can be stored by magnetic tape station (17), or is shown and drawn various parametric lines by display (18) and plotting apparatus (20).

2,, it is characterized in that said sound target plane can be 4 to 16 scopes according to the said method of claim 1.

3, according to the said method of claim 2, it is characterized in that said each target plane, their front and back positions distance to each other can decide according to different content measurements, generally is no more than 30 meters.

4, according to the said method of claim 1, it is characterized in that said sounding rod, its length is generally 1.5 to 5 meters scopes.

5,, it is characterized in that said channel amplifier (6) is to have comprised wave filter (12), amplifier (13), voltage comparator (14), regularly rejector (15) and line driver (16) according to the said method of claim 1.

6, according to the said method of claim 1, it has comprised to open earlier and has afterwards closed circuit (31) to it is characterized in that said time difference tester (30), channel address code translator (32), status register (33), time difference counter (34), alarm (35), polarity detector (36), output buffer (37), line receiver (39) and input pointer (38).

Images