Apparatus and method for synchronizing nonlinear systems using filtered signals5655022Abstract A filtered cascaded synchronized nonlinear system includes a nonlinear transmitter having stable first and second subsystems. The first subsystem produces a first transmitter signal for driving the second subsystem, and the second subsystem produces a second transmitter signal for driving the first subsystem. A first filter filters the second transmitter signal to produce a filter output signal. A subtractor subtracts the filter output signal from the second transmitter signal to produce a transmitter output signal which is transmitted to a nonlinear cascaded receiver. The receiver includes an adder for summing the received transmitter output signal with a receiver filter output signal to restore frequencies that were subtracted from the second transmitter signal in order to produce a first receiver drive signal. The receiver includes cascaded third and fourth subsystems that are respective duplicates of the first and second subsystems. The third subsystem is driven by the first receiver drive signal to produce a first receiver signal in synchronization with the first transmitter signal. The fourth subsystem is driven by the first receiver signal to produce a second receiver signal in synchronization with the second transmitter signal. A second filter filters the second receiver signal to produce the receiver filter output signal. Claims What is claimed and desired to be secured by Letter Patent of the United States is: Description CROSS REFERENCE TO RELATED PATENTS AND APPLICATIONS
______________________________________
Resistor 42 = 10K.OMEGA.
Resistor 62 = 220K.OMEGA.
Resistor 44 = 10K.OMEGA.
Resistor 64 = 150K.OMEGA.
Resistor 46 = 10K.OMEGA.
Resistor 66 = 150K.OMEGA.
Resistor 48 = 20K.OMEGA.
Resistor 68 = 330K.OMEGA.
Resistor 50 = 100K.OMEGA.
Resistor 70 = 100K.OMEGA.
Resistor 52 = 50K.OMEGA.
Resistor 72 = 100K.OMEGA.
Resistor 54 = 3K.OMEGA.
Potentiometer 74 = 10K.OMEGA.
Resistor 56 = 20K.OMEGA.
Capacitor 76 = 0.01 .mu.F
Resistor 58 = 100K.OMEGA.
Capacitor 78 = 0.01 .mu.F
Resistor 60 = 100K.OMEGA.
Capacitor 80 = 0.001 .mu.F
______________________________________
Drive circuit 22 can be subdivided into two subparts 14 and 12. Although the illustrative subparts 14 and 12 shown in FIG. 3 correspond to the two circuits forming drive circuit 22, this is not necessary, and the division of a given drive circuit into subparts in order to determine the proper configuration for a synchronized response circuit is made in accordance with the analysis described herein. Subpart 14 corresponds to the w subsystem (subsystem 14 in FIG. 1), subpart 12 corresponds to the v subsystem described above. Those parts of subpart 14 which affect the signal at X4 and those parts of subpart 12 responsive thereto, respectively, constitute response part 15 (FIG. 2) and drive part 17 (FIG. 9.), to provide feedback. Response circuit 16 is substantially a duplicate of subpart 14 of drive circuit 22 (the specifications for primed components, such as resistor 50', are the same as the specification for unprimed components, such as resistor 50) and corresponds to subsystem w' (subsystem 16) described hereinabove. Signals X.sub.1, X.sub.2, X.sub.3, and X.sub.4 are characteristic voltages of drive circuit 22. The signal X.sub.4 is connected as drive signal S.sub.d through a buffer amplifier 25, which ideally is an operational amplifier having linear characteristics such as an AD381 manufactured by Analog Devices, to response circuit 16 at the junction in circuit 16 corresponding to the junction in circuit 22 at which the signal X.sub.4 is generated. Signal X.sub.4 replaces the circuitry (subpart 12) of drive circuit 22 which is missing in response circuit 16. The subsystem of buffer amplifier 25 is the secondary means 19. Drive circuit 22 is an autonomous system and behaves chaotically. It can be modeled by the following equations of motion for the three voltages X.sub.1, X.sub.2 and X.sub.3 shown in FIG. 3. X.sub.1 =X.sub.2 +.gamma.X.sub.1 +cX.sub.3 X.sub.2 =-.omega..sub.2 X.sub.1 -.delta..sub.2 X.sub.2 .epsilon.X.sub.3 =(1-X.sub.3).sup.2 (sX.sub.1 -r+X.sub.3)-.delta..sub.3 X.sub.3, (9) where .gamma.=0.12, C=2.2, .omega..sub.2 =10.0, .delta..sub.2 =.delta..sub.3 =0.001, .epsilon.=0.001, s=1/6, and r=0.0. An analysis of the sub-Lyapunov exponents for the response circuit 16 requires a transformation of the equations of motion from the (X.sub.1, X.sub.2, X.sub.3) system to the (X.sub.1, X.sub.2, X.sub.4) system. This is done by analyzing the circuit, and finding that X.sub.3 =.alpha.X.sub.4 -.beta.X.sub.1 where .alpha.=6.6 and .beta.=7.9. This gives the following equations of motion: X.sub.1 =X.sub.2 +.gamma.X.sub.1 +c(.alpha.X.sub.4 -.beta.X.sub.1) X.sub.2 =-.omega..sub.2 X.sub.1 -.delta..sub.2 X.sub.2 .epsilon.X.sub.4 =(1/.alpha.){(1-(.alpha.X.sub.4 -.beta.X.sub.1).sup.2)(sX.sub.1 -r+.alpha.X.sub.4 -.beta.X.sub.1)-.delta..sub.3 .alpha.X.sub.4 -.beta.X.sub.1 -.beta.X.sub.2 -.beta..gamma.X.sub.1 -.beta.c(.alpha.X.sub.4 -.beta.X.sub.1)}(10) The equations of motion for the response are just the X.sub.1 and X.sub.2 equations. The sub-Lyapunov exponents are calculated directly from the Jacobian of the X.sub.1 and X.sub.2 equations, which is a constant in this case. It will be appreciated that conventional methods for calculating Lyapunov exponents, as analytical, measurement, numerical and otherwise can be used, such as, for example, those described by Eckmann et al., Rev. Mod. Phys., Vol. 57, p.617 et seq. (1985); Lichtenberg et al., Regular and Stochastic Motion, Springer-Verlag, New York (1983); Rashband, Chaotic Dynamics of Nonlinear Systems, John Wiley and Sons, New York (1990); and Wolf et al., Physica, Vol. 16D, p. 285 et seq. (1985). The sub-Lyapunov exponents in this case are -16.587 and -0.603, implying that synchronization of the two electrical circuits 22 and 16 will occur. X.sub.4 is the drive signal S.sub.d for the response subsystems and (X.sub.1, X.sub.2) and (X.sub.1 ', X.sub.2 ') are the synchronized signals S.sub.o and S.sub.o '. Circuit 22 itself runs in the realm of a few hundred Hz. Response circuit 16 synchronizes with drive circuit 22 within about two milliseconds. It has been observed experimentally that small changes (.about.10%) of the circuit parameters do not affect synchronization greatly, in that the response voltages still remain close to their counterparts in drive circuit 22; but larger changes (.about.50%) do. Even though the sub-Lyapunov exponents for the larger changes both remain negative, the response voltages no longer remain close to their drive counterparts. The circuit of FIG. 3 has been used to transmit a pure frequency signal hidden in a chaotic signal as follows. With circuits 22 and 16 operating in a synchronized mode, a sine wave of a few hundred Hz was added to the X.sub.2 signal from the drive circuit and sent to the response circuit. The X'.sub.2 signal produced by response circuit 16 was then subtracted from the sum of the X.sub.2 signal and the sine wave, thereby extracting the sine wave from the chaotic signal. Spectral analysis of the (X.sub.2 +sine wave) combination signal showed that the sine wave could not be detected in the chaos of the X.sub.2 signal. The smallest sine wave that could be extracted this way was approximately 40 millivolts peak to peak compared to a two volt peak to peak X.sub.2 signal, or a 50:1 ratio of chaotic signal to sine wave. Many other possible choices for the drive circuit are possible and may require transformation of the circuit equations to model them. This can be determined as described hereinabove for nonlinear circuits by analyzing the circuit dynamics in terms of the sub-Lyapunov exponents to determine which signal(s) to choose as a drive signal or signals, and which subcircuit is to be used as a model for the response circuit. It will also be appreciated that the prior art as described previously is applicable to any system which requires synchronization of remote signals and/or their low correlation with each other. For example, the prior art is particularly suited for use in control devices relying on wide-frequency-band synchronized signals. Similar principles as discussed previously can be applied to cascaded subsystems which allow the multiple signals to be synchronized. In the following discussion the previously discussed design of synchronized subsystems is built on by cascading two or more subsystem responses. The objective here is to get a synchronization of the response with its counterpart in the drive system but to build a response setup which produces signals in synchronization with one or more of the original input drive signals. The new synchronization signal may be used to process the original input drive, to detect parameter changes between the drive system and the responses, and to detect other information transmitted along with the output of the drive system. FIG. 4 illustrates a cascaded system having a drive system 1400 and a response system 1500. The drive system 1400 includes two subsystems A and B which are interdependent and may or may not overlap. Subsystem A drives subsystem B with signal S.sub.A and subsystem B drives subsystem A with signal S.sub.B. The response system 1500 produces a signal S.sub.B" which is to be synchronized with a signal S.sub.B produced in the drive system 1400. The subsystem B of drive system 1400 transmits a drive signal S.sub.B to the response system 1500. The response system 1500 includes two subsystems A' and B" that are cascade connected. As in the single stage subsystems discussed earlier (see FIG. 1), subsystems A' and B" are duplicates of subsystems A and B, respectively, which have all-negative sub-Lyapunov exponents. The subsystem A' receives the drive signal S.sub.B and provides a response signal S.sub.A, to the subsystem B". The subsystem B" in turn produces signal S.sub.B" in synchronization with signal S.sub.B. Unlike the single stage synchronization systems discussed earlier (see FIG. 1), in the cascade system shown in FIG. 4, the same signal S.sub.B which the response system 1500 synchronizes with respect to is also used to drive the response system 1500. In the single stage synchronization system 10 of FIG. 1, the synchronized signal S.sub.o may be different than the drive signal S.sub.d. The response system 1500 with the cascaded subsystems A' and B" is not only capable of producing the signal S.sub.B" synchronized with the signal S.sub.B but is also capable of producing the signal S.sub.A' which is in synchronization with the signal S.sub.A. Because the signal S.sub.B" can be compared to the signal S.sub.B, the fact of synchronization can be clearly determined allowing those on the response system side to rely on the synchronization of the S.sub.B and S.sub.B" signals in concluding that signal S.sub.A' is in synchronization with signal S.sub.A. Because of the nature of nonlinear dynamical systems driven in the chaotic regime, properties of one chaotic system do not necessarily carry over to another chaotic system. Nevertheless, the prior art applies to any chaotic system in general, so long as the chaotic system includes at least two stable subsystems. The two response signals or outputs S.sub.A' and S.sub.B" are produced as follows. The first subsystem A' accepts the input signal S.sub.B and produces its response signal S.sub.A' in synchronization with its counterpart (S.sub.A) in the drive system 1400. The second subsystem B" is driven by signal S.sub.A' from the first subsystem A'. The second subsystem response S.sub.B" produces signal S.sub.B" in synchronization with its counterpart S.sub.B in the drive system 1400, which in this case is the original drive signal S.sub.B coming from the element B. The subsystems A' and B" are selected so that all of the essential elements of the drive system 1400 that are not present in the first subsystem A' are present in the second subsystem B" and vice-versa. In other words, the logical union of subsystems A' and B" includes all of the essential elements of the drive system 1400. It is to be noted that each subsystem A' and B" in the response system 1500 is driven by a signal which supplies information in the drive system 1400 which is lacking in the drive subsystem. Thus, subsystem A' in the response system 1500 is driven by the same signal S.sub.B that drives subsystem A in the drive system 1400. Subsystem B" in the response system 1500 is driven by signal S.sub.A' produced by subsystem A', just as subsystem B in the drive system 1400 is driven by signal S.sub.A produced by subsystem A. As discussed earlier, subsystems A, A', B and B" must have all-negative sub-Lyapunov exponents. In other words, subsystems A, A', B and B" are stable subsystems. The same principles-discussed above concerning cascaded systems with 2 subsystems apply equally well to cascaded systems with more than 2 subsystems. In particular, each of the cascaded subsystems in the response system 1500 is a duplicate of a stable subsystem in the drive system 1400. Each subsystem in the response system 1500 is driven by a signal which supplies information from the complete system that is lacking in the drive subsystem, in particular, by a signal corresponding to the signal which drives the corresponding subsystem in the drive system 1400. To understand the theory behind the cascaded system of FIG. 4 it is necessary to build on the previous discussion of equations 1-3. Once the first subsystem of the response system 1500 is created a second system is created, say modeled by the set of differential equations r=a(r,s) and s=b(r,s), where r and s are subsets of variables of u in the same way that v and w are subsets of variables of u. The r variables are the drives for the second subsystem just as the v variables were for the first subsystem. The functions a and b are the corresponding vector field components. If this second subsystem is a stable subsystem (See Pecora et al., Synchronization in Chaotic Systems, Physical Review Letters, Vol. 64, No. 8, February 1990 and Pecora et al., Driving Subsystems With Chaotic Signals, Physical Review A, Vol. 44, No. 4, August 1991, both incorporated by reference herein, for a discussion of how to determine whether stability exists), the s variables synchronize with their corresponding variables in the first system and with the drive signal. This then provides a signal in synchronization with the input drive (one or more of the variables). For any two dynamical systems to become synchronized, they must start in the same basin of attraction. That is, their starting points (initial conditions) must be in the same set of points which will converge to the same attractor. Since many dynamical systems can have more than one attractor, it is possible for two such systems to start in different basins. If the response system 1500 has somewhat different parameters than the drive system 1400, the synchronized signals will not be exactly equal and in general will have a difference which at small parameter changes will be proportional to the derivative of the vector fields with respect to the parameters. As discussed below, this effect along with others in the dynamical system allows communication using signals from nonlinear systems, including chaotic ones. The details of cascaded synchronized systems and the circuit design, construction, and operation thereof will now be discussed. FIG. 5 functionally illustrates an example of a cascaded system. It includes a drive system 800 which includes elements X.sub.1, X.sub.2, X.sub.3, and X.sub.4 characterized by state variables x.sub.1, x.sub.2, x.sub.3 and x.sub.4, respectively, (FIG. 6). Element X.sub.4 constitutes subsystem A, and elements X.sub.1, X.sub.2, X.sub.3 constitutes subsystem B. Both subsystems A and B are stable, that is they have all negative sub-Lyapunov exponents. Subsystem A drives subsystem B with signal S.sub.A and subsystem B drives subsystem A with signal S.sub.B. The response system 900 is a cascade of two subsystems A' and B" where the first subsystem system A' includes a single element X.sub.4' and the second subsystem B" includes three elements X.sub.1", X.sub.2" and X.sub.3". The first subsystem A' is a duplicate of the subsystem A in the drive system 800 and the second subsystem B" is a duplicate of the subsystem B in the drive system 800. The drive system 800 drives the first subsystem A' with signal S.sub.B and the first subsystem A' drives the second subsystem B" with signal S.sub.A'. The second subsystem B" produces an output signal S.sub.B" in synchronization with drive signal S.sub.B. The operation of the elements in this example is modeled by the following equations: dx.sub.1 /dt=-.alpha..sub.1 [.beta..sub.1 A.sub.1 x.sub.1 -.gamma..sub.1 x.sub.2 +x.sub.3 -x.sub.4 +g.sub.1 (x.sub.4)+.delta.x.sub.1 ],(11) dx.sub.2 /dt=-.alpha..sub.2 (x.sub.1 +.delta.x.sub.2), (12) dx.sub.3 /dt=-.alpha..sub.3 (x.sub.2 +.delta.x.sub.3), (13) dx.sub.4 /dt=-.alpha..sub.4 ((-.beta..sub.4 /R.sub.v)x.sub.1 +.gamma..sub.4 A.sub.4 x.sub.4 +g.sub.2 (x.sub.4)), (14) dx".sub.1 /dt=-.alpha..sub.1 [.beta..sub.1 A".sub.1 x".sub.1 -.gamma..sub.1 x".sub.2 +x".sub.3 -x'.sub.4 +g.sub.1 (x'.sub.4)+.delta.x".sub.1 ],(15) dx".sub.2 /dt=-.alpha..sub.2 (x".sub.1 +.delta.x".sub.2), (16) dx".sub.3 /dt=-.alpha..sub.3 (x".sub.2 +.delta.x".sub.3), (17) dx'.sub.4 /dt=-.alpha..sub.4 ((-.beta..sub.4 /R.sub.v)x.sub.1 +.gamma..sub.4 A'.sub.4 x'.sub.4 +g.sub.2 (x'.sub.4)), (18) where the g.sub.1 and g.sub.2 functions are defined as: g.sub.1 (x)=.beta..sub.5 (.vertline.x-2.5.vertline.-.vertline.x+2.5.vertline.), (19) g.sub.2 (x)=.beta..sub.6 x+.gamma..sub.6 (.vertline.x-1.3.vertline.-.vertline.x+1.3.vertline.)+.epsilon.(.vertline. x-2.6.vertline.-.vertline.x+2.6.vertline.) (20) and the constants are .alpha..sub.1 =1098, .alpha..sub.2 =10980, .alpha..sub.3 =4972, .alpha..sub.4 =10980, .beta..sub.1 =1.466, .gamma..sub.1 =2.466, .gamma..sub.4 =10.sup.5, .gamma..sub.4 =0.5, .beta..sub.5 =0.5, .beta..sub.6 =0.5, .gamma..sub.6 =0.164, and .epsilon..sub.6 =0.361. The constant .delta., set at 0.2, is a phenomenological damping constant used to account for leakage current in the capacitors. Its value was set to make the stability of eqns. (11)-(20) match the stability of the actual circuit. A.sub.1 and A.sub.4 are variable parameters normally set at 1.0. As R.sub.v is decreased from 50,000 ohms to 46,000 ohms, the circuit goes from a limit cycle through a period doubling to a one-well chaotic attractor to a two-well chaotic attractor. With R.sub.v held constant the drive system 800 and response system 900 can produce a number of synchronized signals with the output S.sub.B" of the element B" being used to confirm synchronicity as previously discussed. If R.sub.v is varied information can be transferred. FIGS. 6-9 illustrate the circuit details of an example of a system of FIG. 5 where multiple synchronized signals can be produced and synchronization verified. FIG. 6 depicts the details of the drive system 800. This circuit 800 includes the following particular circuit elements:
______________________________________
Resistor R1 = 100k.OMEGA.
Resistor R11 = 221k.OMEGA.
Resistor R2 = 100k.OMEGA.
Resistor R12 = R.sub.v
Resistor R3 = 100k.OMEGA.
Resistor R13 = 100k.OMEGA.
Resistor R4 = 100k.OMEGA.
Resistor R14 = 200k.OMEGA.
Resistor R5 = 68.2k.OMEGA.
Resistor R15 = 100k.OMEGA.
Resistor R6 = 100k.OMEGA.
Resistor R16 = 100k.OMEGA.
Resistor R7 = 100k.OMEGA.
Resistor R17 = 100k.OMEGA.
Resistor R8 = 68.2k.OMEGA.
Resistor R18 = 100.OMEGA.
Resistor R9 = 1m.OMEGA.
Resistor R10 = 100k.OMEGA.
Capacitor C1 = 910pf
Capacitor C3 = 910pf
Capacitor C2 = 910pf
Capacitor C4 = 910pf.
______________________________________
R.sub.v is selected from among 47.8 k.OMEGA. and 46.9 k.OMEGA. with 47.8 k.OMEGA. preferable. Resistor tolerances are preferably 1% and all capacitors are preferably 5% mica capacitors. The system also includes operational amplifiers 01, 02, 03, 04, 05, 06 and 07 all of which are 741 type amplifiers and diode DO which is an IN485B type. The circuit details of the functions g1(x) (eqn. 19) and g2(x) (eqn. 20) are depicted in the circuit diagrams of FIGS. 7 and 8, respectively. Returning now to the example shown in FIG. 6, if one cuts the circuit at points a and b, the resulting systems A and B are stable. Subsystem B consisting of X.sub.1, X.sub.2, and X.sub.3 can be driven with the S.sub.A signal from the full system. Subsystem A consisting of X.sub.4 may be driven with the S.sub.B signal from the full circuit. When driving the B subsystem including elements x.sub.1, x.sub.2, and x.sub.3, it does not actually matter whether the S.sub.A driving signal is coming from the full circuit or from an A (or A') subsystem synchronized to the full circuit. Conversely, when driving the A subsystem, it does not actually matter whether the S.sub.B driving signal is coming from the full circuit or from a B or B" subsystem synchronized to the full circuit. This arrangement, in which the stable subsystems are driven by signals from subsystems and not necessarily the full circuit, is called "cascaded synchronization". FIG. 7 depicts a circuit with response g1(x) (eqn. 19). In this circuit the resistors R=10 k.OMEGA., the operational amplifiers 08, 09, 010 and 011 are 741 types and the diodes D1, D2, D3 and D4 are preferably type IN485B. FIG. 8 depicts a circuit with response g2(x) (eqn. 20) where operational amplifiers 012 and 013 are 741 type amplifiers and
______________________________________
Resistors R21 = 27.4k.OMEGA.
Resistors R29 = 50.1.OMEGA.
Resistors R22 = 27.4k.OMEGA.
Resistors R30 = 50.1.OMEGA.
Resistors R23 = 49.9k.OMEGA.
Resistors R31 = 50.1.OMEGA.
Resistors R24 = 49.9k.OMEGA.
Resistors R32 = 50.1.OMEGA.
Resistors R25 = 200k.OMEGA.
Resistors R33 = 20k.OMEGA.
Resistors R26 = 200k.OMEGA.
Resistors R34 = 178k.OMEGA.
Resistors R27 = 825k.OMEGA.
Resistors R35 = 156.2k.OMEGA.
Resistors R28 = 825k.OMEGA.
Resistors R36 = 100k.OMEGA.
______________________________________
Diodes D5, D6, D7 and D8 are type IN485B. FIG. 9 depicts the circuit details of the response system 900 of FIG. 5. In this circuit 900 the resistor, capacitor, amplifier and function components are the same as previously discussed regarding FIGS. 6-8. In FIG. 6 any of the nodes can be used as the source of the signals to be synchronized. However, the drive signal must come from a particular cut point as discussed above. The above discussion of separated system synchronization is performed with electronic hardware components or other equivalent devices. Using the same concepts, systems can also be synchronized using software. FIG. 10 shows a filtered cascaded synchronized nonlinear system having a transmitter 1100 and a receiver 1200. The transmitter 1100 includes subsystems A and B which are independent and may or may not overlap. Subsystems A and B each contain 1 or more variables, that is, each subsystem is at least 1-dimensional but may contain more than 1 dimension. Neither subsystem A nor subsystem B is contained within the other subsystem. At least part of subsystem A is external to subsystem B and at least part of subsystem B is external to subsystem A. Subsystem A drives subsystem B with signal S.sub.A and subsystem B drives subsystem A with signal S.sub.B. The signal S.sub.B is the input to filter 1110, and the output of filter 1110 is the signal S.sub.f. The subtractor 1120 subtracts the filter output signal S.sub.f from S.sub.B to produce the broadcast signal S.sub.t, which is transmitted to the receiver 1200. The receiver 1200 is responsive to the broadcast signal S.sub.t and produces an output signal S.sub.B" in synchronization with signal S.sub.B. The receiver 1200 consists of subsystem A' which is a duplicate of A and B" which is a duplicate of B. The broadcast signal S.sub.t is used as one input to the adder 1220. The output of the adder 1220 is the signal S.sub.d. The signal S.sub.d is used to drive the subsystem A', and the output signal S.sub.' from subsystem A' is used to drive the subsystem B". Subsystem B" does not directly drive subsystem A', and the sub-Lyapunov exponents (as defined in U.S. Pat. Nos. 5,379,346 and 5,245,660) for subsystems A' and B" are all negative. The signal S.sub.B" from subsystem B" is used as an input for the filter 1210, which is identical to the filter 1110 in the transmitter. The filter 1210 produces an output signal S.sub.o which is used as an input for the adder 1220. When the receiver 1200 is synchronized to the transmitter 1100, then the signals in subsystem A' reproduce the signals in subsystem A and the signals in subsystem B" reproduce the signals in subsystem A. If the subsystems A' and B" and the filter 1210 are not exact replicas of the subsystems A,B and the filter 1100 (which will be the case in an electronic circuit implementation of the present invention), then the signals in A' and B" can be made arbitrarily close to the signals in A and B by making the differences between A and A', B and B", and filters 1110 and 1210 arbitrarily small. In order to determine if the receiver 1200 will synchronize to the transmitter 1100, it is necessary to determine the stability of the receiver 1200 in the synchronized state. Techniques for determining the stability of such a system are well known; see, for example, J. M. T Thompson and H. B. Stewart, "Nonlinear Dynamics and Chaos", (Wiley, New York, 1986) or F. C. Moon, "Chaotic Vibrations", (Wiley, New York, 1987). The synchronization of the receiver 1200 to the transmitter 1100 may be confirmed by comparing the receiver output signal S.sub.B" to the receiver driving signal S.sub.d. When the receiver is synchronized to the transmitter, then signal S.sub.B" will match signal S.sub.d. FIG. 11 shows a filtered cascaded synchronized nonlinear system having a transmitter 2100 and a receiver 2200 when the nonlinear systems are nonautonomous, that is, they have a periodic forcing part F (2130). The description of the transmitter 2100 is the same as the description of the transmitter 1100 in FIG. 10 except that the periodic forcing source F provides periodic forcing signals F.sub.A and F.sub.B to subsystems A and B. Either F.sub.A or F.sub.B may be zero, but they may not both be zero. The receiver 2200 in FIG. 11 contains a periodic forcing source F' (2240) which provides the periodic forcing signal F'.sub.A to subsystem A' and periodic forcing signal F'.sub.B to subsystem B". If F.sub.A in transmitter 2100 is zero, then F'.sub.A in receiver 2200 is zero, and if F.sub.B in the transmitter 2100 is zero, then F'.sub.B in receiver 2200 is zero. The receiver 2200 in FIG. 11 operates in the same manner as the receiver 1200 in FIG. 10, except that it is necessary to match the phase of the periodic forcing source F' (2240) in the receiver 2200 to the phase of the periodic forcing source F (2130) in the transmitter 2100. The receiver 2200 contains a phase control system 2230 responsive to signals S.sub.B" and S.sub.d. The phase control system 2230 generates an error signal .DELTA. proportional to the phase difference between F and F'. The periodic forcing source F' uses the error signal .DELTA. to match the phase of F' to the phase of F. The procedures for producing the error signal .DELTA. are described in U.S. patent application Ser. No. 08/267,696 (Navy Case No. 75,496), entitled: "SYNCHRONIZATION OF NONAUTONOMOUS CHAOTIC SYSTEMS", filed Jun. 29, 1994, Inventors: Thomas L. Carroll et al. The systems in FIGS. 10 and 11 may be any nonlinear dynamical system or combination of systems, provided that they may be subdivided into stable subsystems. The systems may be electronic circuits, they may be sets of differential equations or recursion relations (maps) to be solved on a computer, they may be implemented in digital signal processing systems or other physical or electronic systems, or they may be any other physical system that can be broken into stable subsystems. The filters, adders and subtracters may also be electronic devices or they may be implemented as computer algorithms. It is also possible for part of the system to be of one type (such as a computer algorithm) and the other part of the system to be of some other type (such as an electronic circuit). FIG. 12-18 show how the present invention may be built as an electronic circuit. The circuit details of an electronic example of a chaotic circuit 3000 are shown in FIGS. 12, 13 and 14. This circuit 3000 includes the following particular circuit elements:
______________________________________
Resistor R1 = 10k.OMEGA.
Resistor R2 = 39.2k.OMEGA.
Resistor R3 = 10k.OMEGA.
Resistor R4 = 10k.OMEGA.
Resistor R5 = 10k.OMEGA.
Resistor R6 = 10k.OMEGA.
Resistor R7 = 100k.OMEGA.
Resistor R8 = 1M.OMEGA.
Resistor R9 = 1M.OMEGA.
Resistor R10 = 100k.OMEGA.
Resistor R11 = 1M.OMEGA.
Resistor R12 = 100k.OMEGA.
Resistor R13 = 100k.OMEGA.
Resistor R14 = 100k.OMEGA.
Resistor R15 = 5.2k.OMEGA.
Resistor R16 = 100k.OMEGA.
Resistor R17 = 100k.OMEGA.
Resistor R18 = 1M.OMEGA.
Capacitor C1 = 1nF Capacitor C2 = 1nF
Capacitor C3 = 1nF
______________________________________
Resistor tolerances are preferably 1% or better and all capacitors are preferably 5% mica capacitors. The system also includes operational amplifiers Op1, Op2, Op3, Op4, Op5, Op6, Op7, Op8 and Op9, all of which are 741 type amplifiers. The circuit details of the circuit g of FIG. 12 is depicted in the circuit diagram of FIG. 13, having the following particular elements:
______________________________________
Resistor R19 = 100k.OMEGA.
Resistor R20 = 100k.OMEGA. Resistor
R21 = 100k.OMEGA. Resistor R22 = 100k.OMEGA. Resistor
R23 = 680k.OMEGA. Resistor R24 = 2M.OMEGA. Resistor
R25 = 680k.OMEGA. Resistor R26 = 2M.OMEGA. Resistor
R27 = 100k.OMEGA.
Potentiometer P1 = 20k.OMEGA.
Potentiometer P2 = 50k.OMEGA.
Potentiometer P3 = 20k.OMEGA.
Potentiometer P4 = 50k.OMEGA.
Potentiometer P5 = 20k.OMEGA. in
Potentiometer P6 = 20k.OMEGA. in
parallel with a 100.OMEGA.
parallel with a 100.OMEGA.
resistor (not shown)
resistor (not shown)
Potentiometer P7 = 20k.OMEGA. in
Potentiometer P8 = 20k.OMEGA. in
parallel with a 100.OMEGA.
parallel with a 100.OMEGA.
resistor (not shown)
resistor (not shown)
______________________________________
Resistor tolerances are preferably 1% or better. The system also includes operational amplifier Op10 which is a type 741, and diodes D1, D2, D3, and D4 which are of type 1N485B. As explained further below, the potentiometers P.sub.1 -P.sub.8 are used to match different circuits g to each other. The circuit details of the circuit f of FIG. 12 are depicted in the circuit diagram of FIG. 14, having the following particular elements:
______________________________________
Resistor R28 = 10k.OMEGA.
Resistor R29 = 10k.OMEGA.
Resistor R30 = 490k.OMEGA.
Resistor R31 = 490k.OMEGA.
Resistor R32 = 50k.OMEGA.
Resistor R33 = 50k.OMEGA.
Resistor R34 = 20k.OMEGA.
Resistor R35 = 100k.OMEGA.
Resistor R36 = 100k.OMEGA.
Resistor R37 = 100k.OMEGA..
______________________________________
Resistor tolerances are preferably 1% or better. The system also includes operational amplifiers Op11 and Op12 which are of type 741, and diodes D5, and D6 which are of type 1N485B. The circuit shown in FIGS. 12-14 is modeled by the following equations: dx/dt=.beta.[y-z] (21) dy/dt=.beta.[-.GAMMA..sub.y .multidot.y-g(x)+.alpha..multidot. cos (.omega..sub.t .multidot.t)] (22) dz/dt=.beta.[f(x)-.GAMMA..sub.z .multidot.z] (23) g(x)=-3.8+0.5*(.vertline.x+2.6.vertline.+.vertline.x-2.6.vertline.+.vertlin e.x+1.2.vertline.+.vertline.x-1.2.vertline.) (24) f(x)=0.5*x+.vertline.x-1.vertline.+.vertline.x+1.vertline.,(25) where .alpha.=2.0, .GAMMA..sub.y =0.2, .GAMMA..sub.z =0.1, the time factor .beta.=10.sup.4 /sec, and the angular frequency .omega..sub.t =2.pi.f.sub.t, where the transmitter forcing frequency f.sub.t =769 Hz. The cosine term in Equation (22) is provided by a signal S.sub.1 supplied by an HP 3300A function generator 2130. The functions g(x) (equation (24)) and f(x) (equation (25)) are piecewise linear functions produced by the circuits shown in FIGS. 13 and 14, respectively. Equations (21-25) model the B subsystem (FIG. 11) of the transmitter 2100, and equation (23) models the A subsystem (FIG. 11) of the transmitter 2100. This circuit is designed so that it is possible to create a synchronizing subsystem. Equations (21-22) (with z treated as a parameter) constitute the well known 1-well Duffing equations. For the parameter settings used here, the behavior of such a subsystem is periodic, indicating that the largest Lyapunov exponent for this subsystem is zero. Equation (23) was added to the Duffing system of Equations 21-22 to provide an instability for certain values of x, thereby leading to chaos. If the feedback loop between equations (21) and (23) were not completed, i.e. if the subsystem of equations (21-22) were not dependent on the z-variable produced by the subsystem of equation (23), or if the subsystem of equation (23) were not dependent on the x-variable produced by the subsystem of equations (21-22), then the system of variables x, y and z would be periodic. In other words variables x, y and z would each be periodic or a fixed point. The largest conditional Lyapunov exponent with respect to the signal S.sub.1 would be less than or equal to zero. The feedback loop between equations (21) and (23) can be disconnected by cutting the system at node T.sub.1 and grounding the input (x) to the circuit f, or by cutting the system at node T.sub.2. Such disconnection would remove the dependence of equation (23) on the variable x, or the dependence of equation (21) on the variable z, respectively. The conditional Lyapunov exponents for the transmitter system of FIGS. 12-14 calculated from equations (21-25) with the above parameters are 284s.sup.-1, -1433s.sup.-1 and -1854s.sup.-1. The sinusoidal forcing term cos(.omega..sub.t .multidot.t) of equation (22) is treated as a parameter in this calculation, so the zero exponent attributable to signal S.sub.1 does not show up here. Since one of the conditional Lyapunov exponents is positive, therefore the system modeled by equations (21-25) and shown in FIGS. 12-14 operates in the chaotic regime. In FIG. 15, R38=10,000 ohms, R39=10,000 ohms, R40=10,000 ohms, R41=10,000 ohms, R42=10,000 ohms, R43=5000 ohms, R44=10,000 ohms, R45=10,000 ohms and R46=10,000 ohms. In FIG. 16, resistor R47=31,380 ohms and C4=10.sup.-8 F. All operational amplifiers are type 741. The transmitter 6000 in FIG. 15 may be described by the following differential equations when the filter 4000 of FIG. 16 is used (equations 26-30): dx/dt=.beta.(y-z) (26) dy/dt=.beta.-.GAMMA..sub.y y-g(x)+.alpha. cos (.omega.t)+A (27) dz/dt=.beta.f(x)-.GAMMA..sub.z z (28) du/dt=dx/dt-u/RC (29) s.sub.t =x-u (30) dv/dt=dx"/dt-v/RC (31) s.sub.d =s.sub.t +v (32) dz'/dt=.beta.f(s.sub.d)-.GAMMA..sub.z z' (33) dx"/dt=.beta.y"-z' (34) dy"/dt=.beta.-.GAMMA..sub.y y"-g(x")+.alpha. cos .omega..sub.r t+.PHI.r+A(35) .alpha.=1.9, .GAMMA..sub.y =0.2, .GAMMA..sub.z =0.1, A=0, .beta.=10.sup.4 s.sup.-1, .omega.=2.pi..times.780 Hz. The receiver 7000 (FIG. 15) is shown in FIGS. 17 and 18. The receiver 7000 consists of a cascaded response circuit 5100 (FIG. 17) and a phase control circuit 5200 (FIG. 18), along with the filter 4000 (FIG. 16) and the adder formed by operational amplifier op14 (FIG. 15). The cascaded response circuit 5100 is described by equations 33-35. In FIG. 17, R51=10,000 ohms, R52=39,200 ohms, R53=10,000 ohms, R54=10,000 ohms, R55=10,000 ohms, R56=10,000 ohms, R57=100,000 ohms, R58=1,000,000 ohms, R59=1,000,000 ohms, R60=100,000 ohms, R61=1,000,000 ohms, R62=100,000 ohms, R63=100,000 ohms, R64=100,000 ohms, R65=5,200 ohms, R66=100,000 ohms, R67=100,000 ohms, R68=1,000,000 ohms, C6=1 nF, C7=1 nF, C8=1 nF. The phase control circuit 5200 is shown in FIG. 18. Referring now to FIG. 18, the details of a phase-detector/controller 5200 is shown. This phase-detector/controller 5200 is responsive to the receiver drive signal S.sub.d and to the receiver output signal S.sub.B" for producing a correction signal .DELTA. responsive to the phase difference between the transmitter forcing signal F1 and the receiver forcing signal F.sup.1. The strobe input signal S.sub.B" generated by the response system 5100 is applied to an amplifier 5220 with a high gain such as a 741 type amplifier with a gain of -100. The output of the amplifier 5220 is applied to a conventional comparator 5221, such as an AD 790. The comparator 5221 produces an output when the input signal S.sub.B" is less than zero. The positive-going signal from the comparator 5221 triggers a conventional Schmitt trigger circuit 5222, such as an SN 74121 monostable multi-vibrator. As a result, the Schmitt trigger circuit 5222 produces a pulse of about 1 microsecond (.mu.s) duration when the strobe input signal S.sub.B" crosses 0 in the negative direction. A difference device 5225, such as a 741 operational amplifier, generates the difference signal S.sub.d -S.sub.B" between the receiver drive signal S.sub.d and the strobe signal S.sub.B". The difference signal S.sub.d -S.sub.B" produced by the difference device 5225 is applied to the signal input of a conventional sample and hold circuit 5226, such as an LM 398, and the output of the Schmitt trigger circuit 5222 is applied to the logic input of the sample and hold circuit 5226. In other words, the difference S.sub.d -S.sub.B" between the receiver drive signal S.sub.d and the strobe signal S.sub.B" is applied to the sample and hold circuit 5226, which holds the difference seen when the strobe signal S.sub.B" passes through 0 going negative. The sampled signal produced by the sample and hold circuit 5226 is applied to the negative terminal of a 741 type amplifier 5228, and the correction signal .DELTA. is applied to the positive terminal of the amplifier 5228 thereby providing negative feedback. The amplifier 5228 thus accumulates the sampled difference signal and the correction signal .DELTA.. The correction signal .DELTA. is produced by a conventional integrator 5229, having a long time constant preferably of about 10 seconds (s), such as type 741 amplifier with a mica capacitor used for feedback, that averages the output of the amplifier 5228. In other words, the output of the sample and hold circuit 5226 is applied to an integrator to produce a correction signal .DELTA. proportional to the average phase difference transmitter forcing signal F and receiver forcing signal F.sup.1. Referring back to FIG. 11, a signal generator 2130 responsive to the correction signal .DELTA. produced by the phase-detector/controller 2230 of FIG. 18 which is itself responsive to a receive drive signal S.sub.d produced by the circuit shown in FIG. 15 preferably utilizes an HP 8116A function generator (not shown). Such a signal generator 2130 multiplies the correction signal .DELTA. produced by the phase-detector/controller 5200 of FIG. 18 by a factor of 1/100 and uses the resulting signal to modulate the frequency of the HP8116A function generator. The transmitter 2100 of eqs. 26-30 and the response system 2200 of eqs. 31-35 are not identical. The transmitter 2100 and receiver 2200 are effectively identical when they are synchronized. It is necessary for the synchronized state to be stable. A Lyapunov exponent calculation from the equations shows that the largest Lyapunov exponent in the response system is -319 s-1, indicating that the response system is stable. An alternate filter 4000 is shown in FIG. 19. The resistor values were given by RA.sub.ij =:
______________________________________
j = 1 j = 2 j = 3
______________________________________
i = 1 204,000 .OMEGA.
408,000 .OMEGA.
1026 .OMEGA.
i = 2 102,000 .OMEGA.
204,000 .OMEGA.
513 .OMEGA.
i = 3 68,000 .OMEGA.
136,000 .OMEGA.
342 .OMEGA.
i = 4 51,000 .OMEGA.
102,000 .OMEGA.
256 .OMEGA.
i = 5 40,800 .OMEGA.
82,000 .OMEGA.
205 .OMEGA.
______________________________________
and the capacitor CA was 10.sup.-8 F. This filter was described by the equations: w=dx/dt (36) u.sub.i /dt=-(2/R.sub.i2 C)u.sub.i -(1/R.sub.i2 C)(1/(R.sub.i3 C+1/R.sub.i1 C)v.sub.i (1/R.sub.i1 C) (37) dv.sub.1 /dt=u.sub.i (38) s.sub.t =x+.SIGMA.v.sub.i (39) s.sub.d =s.sub.t -.SIGMA.r.sub.i (40) dq.sub.i /dt=-(2/R.sub.i2 C)q.sub.i -(1/R.sub.i2 C)(1/(R.sub.i3 C)+1/(R.sub.i1 C))r.sub.i -(1/R.sub.i1 C)dx"/dt (41) dr.sub.i /dt=q.sub.i (42) where the resistor values are given by the above table. The largest Lyapunov exponent for the response system when alternate filter 4000 was used was found to be -10 s-1, indicating that the response system was stable. The transmitter filter output signal v was added to the drive output signal x in equation 39 because the filter of equations 37 and 38 inverted the input signal. For the same reason, the receiver filter output signal was subtracted from the transmitted signal in equation 40. FIG. 20(a) shows the power spectrum of the drive system output signal x from equation 26, while FIG. 20(b) shows the power spectrum of the transmitted signal s.sub.t described in equation 30, demonstrating the change in the power spectrum caused by the filtering. It may be shown that this technique also allows phase synchronization of the periodic forcing parts of nonautonomous synchronized nonlinear systems. The controller 5200 of FIG. 18 was used to control the phase of the response circuit periodic forcing to match that of the drive circuit. The controller 5200 generated a series of voltages that corresponded to the value of the response system output signal x" when the input signal s.sub.d crossed zero. If the drive and response circuits were synchronized, these voltages would all be zero. An integrator 5229 with a time constant of 1 s averaged the series of voltages to produce an error signal .DELTA., which was used to vary the frequency of the response periodic forcing 2240 to bring the phase into sync with the drive periodic forcing 2130. FIG. 21 shows the periodic forcing F" for the response vs. the periodic forcing F for the drive. There is some fluctuation of the response phase and a constant phase offset which is an artifact of the control circuit, but the basic principle works. This demonstrates that the nonperiodic part of the chaotic signal carries information about the phase of the periodic part. Most of the phase fluctuation is believed to be caused by component mismatch between the two circuits. There is also a phase flip caused by a sign change in the filters. Several authors have demonstrated communication between cascaded chaotic circuits via parameter switching in the sending circuit [U. Parlitz, L. O. Chua, L. Kocarev, K. S. Halle, K. Shang, Transmission of Digital Signals by Chaotic Synchronization, International Journal of Bifurcations and Chaos, vol. 2, p. 973 (1992)]. Parameter switching may also be used with the filtered nonautonomous chaotic circuits. The forcing offset A in eq. 27 was switched between .+-.1.0 V, and the parameter switching was detected by monitoring the error signal .DELTA. generated by the controller 5200. FIG. 22 shows the offset signal A as a time series and the resulting error signal .DELTA. coming from the response system controller. As may be seen by the sharp edges on the error signal transitions, the offset signal could be switched up to a factor of about 4 faster. The switching speed is limited by the time constant of the integrator that produces .DELTA., about 1 s for this system. Autonomous nonlinear systems may also be synchronized when the driving signal is filtered. The Piecewise Linear Rossler (PLR) system is a nonlinear system that may be synchronized in a cascaded fashion (T. L. Carroll, "A simple circuit for demonstrating regular and synchronized chaos", American Journal of Physics, vol 63, #4, pp. 377-379, April 1995). A bandpass filter was used to isolate a large periodic component in the output of the PLR circuit and reduce its presence in the transmitted signal. Reducing the size of the periodic component reduced the power contained in the transmitted signal by a large amount, so that the transmitted signal could be sent with less power. The PLR system and the filter were described by the equations: dx/dt=-500(x+10y+20z) dy/dt=-10.sup.4 (-x-0.13y+0.02y) dz/dt=-10.sup.4 (z-g(x)) (43) du/dt=-800u-5.times.10.sup.7 v-400(dy/dx) dv/dt=u (44) y.sub.t =y+1.5v (45) dw/dt=-800w-5.times.10.sup.7 r-400(dy'/dt) dr/dt=w (46) y.sub.d =y.sub.t -1.5r (47) dx'/dt=-500(x'10y'+20z') dy'/dt=-10.sup.4 (-x'-0.13y.sub.d +0.02y') dz/dt=-10.sup.4 (z'-g(x')) (48) g(x)=:0 if x<3, 15(x-3) otherwise (49) Equations 43 are the drive system and equations 44 and 45 are the drive system filter. Equations 48 are the response system and equations 46 and 47 are the response system filter. The transmitter filter output signal v was added to the drive output signal y in equation 45 because the filter of equations 44 inverted the input signal. For the same reason, the receiver filter output signal was subtracted from the transmitted signal in equation 47. Equation 49 is the nonlinear function g(x). Equations 43-49 form an embodiment of the present invention as a complete algorithm. FIG. 23 shows the output signal (the y signal) from the drive system of equations 43 (the solid line in FIG. 23). This signal has a large periodic component, as may be seen in the power spectrum of the y signal in FIG. 24. The bandpass filter of equations 44 is tuned to this periodic component. The filter output signal v is then subtracted from the drive system output signal y to give the transmitted signal y.sub.t. The transmitted signal y.sub.t is shown as a dotted line in FIG. 23. The power spectrum of the transmitted signal y.sub.t from equation 45 is shown in FIG. 25. The numerical integration routine generated 20,000 point output time series of y and y.sub.t were squared and integrated to give an estimate of the power in each signal. The power in the y signal was 155,371 (arbitrary units), while the power in the y.sub.t signal was 15,941 (arbitrary units). Filtering of the y signal to produce the y.sub.t signal reduced the power contained in the signal by a factor of approximately 10, reducing the power that must be transmitted. The receiver is composed of equations 46, 47, 48 and 49. The filter output signal r is subtracted from the transmitted signal y.sub.t to produce the receiver driving signal y.sub.d (equation 47). The driving signal y.sub.d is then used to drive the cascaded response system of equations 48 to produce the response system output signal y'. The derivative of the response system output signal y' is used to drive the receiver filter of equations 46 to produce the filter output signal r. FIG. 26 shows the response system output signal y' vs. the drive system output signal y to demonstrate that the drive and response systems are indeed synchronized. The filtered synchronized communications system may also be used to correct for the effects of filtering by the communications channel. FIG. 27 shows a transmitter 6100 which sends a signal to a receiver 6300 through a communications channel 6200. If the effect of the communications channel is to filter the signal S.sub.B with a filter of the form (1-F) to produce a signal S.sub.t, then the effect of the communications channel filtering may be removed by using a filter F for filter 6320 in receiver 6300. Therefore, what has been described in a preferred embodiment is a filtered cascaded synchronized nonlinear system which includes a nonlinear transmitter having stable first and second subsystems. The first subsystem produces a first transmitter signal for driving the second subsystem, and the second subsystem produces a second transmitter signal for driving the first subsystem. A first filter filters the second transmitter signal to produce a filter output signal. A subtractor subtracts the filter output signal from the second transmitter signal to produce a transmitter output signal which is transmitted to a nonlinear cascaded receiver. The receiver includes an adder for summing the received transmitter output signal with a receiver filter output signal to restore frequencies that were subtracted from the second transmitter signal in order to produce a first receiver drive signal. The receiver includes cascaded third and fourth subsystems that are respective duplicates of the first and second subsystems. The third subsystem is driven by the first receiver drive signal to produce a first receiver signal in synchronization with the first transmitter signal. The fourth subsystem is driven by the first receiver signal to produce a second receiver signal in synchronization with the second transmitter signal. A second filter filters the second receiver signal to produce the receiver filter output signal. It should therefore readily be understood that many modifications and variations of the present invention are possible within the purview of the claimed invention. It is therefore to be understood that, within the scope of the appended claims, the invention may be practiced otherwise than as specifically described.
|
Same subclass Same class Consider this |
||||||||||
