US4819200A - Real-time N-dimensional analog signal rotator - Google Patents
Real-time N-dimensional analog signal rotator Download PDFInfo
- Publication number
- US4819200A US4819200A US06/845,855 US84585586A US4819200A US 4819200 A US4819200 A US 4819200A US 84585586 A US84585586 A US 84585586A US 4819200 A US4819200 A US 4819200A
- Authority
- US
- United States
- Prior art keywords
- matrix
- signals
- signal
- output
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R13/00—Arrangements for displaying electric variables or waveforms
- G01R13/20—Cathode-ray oscilloscopes
- G01R13/206—Arrangements for obtaining a 3- dimensional representation
Definitions
- the apparatus can be an electronic instrument which, when used with a standard oscilloscope, enables display of a two-dimensional projection from any desired perspective of an n-dimensional surface S evolving in time, where n is three or greater.
- a three-dimensional surface S for example, represents the dynamic loci of three electrical signals x(t), y(t) and z(t) in the X-Y-Z space over some observation time interval.
- the rotation in this invention is applicable for N distinct electrical signals, and the signal can be rotated in real time via electronic circuits, rather than being merely rotated by digital computations.
- the invention applies to any electronic circuit which implements an "orthogonal matrix transformation" in any dimension.
- the invention is a practical, and presently the only, method for observing the cross-section of any strange attractor (or any n-dimensional Lissajous figure generated by N electrical signals) from any perspective, and in forward or reverse time.
- Such a surface can be described as the collection of points (x,y,z) in X-Y-Z space generated by three electrical signals x(t), y(t) and z(t). These signals may arise from electrical systems, or they may be obtained from sensors monitoring such complicated dynamical phenomena as turbulence encountered in airplane flights and stress patterns in structures.
- This invention allows such a surface to be analyzed in real time (instead of by computer software which is often too slow and is too bulky to be used as a portable instrument) by projecting it onto an oscilloscope screen. With this invention, the surface can be rotated to any desired orientation, and the rotation continued during observation. Moreover, any cross-section of the surface can be dissected and displayed instantaneously on an oscilloscope.
- this rotator can be implemented by converting the analog signals first into a digital signal via an A/D converter, implementing the orthogonal transformation digitally, and finally reconverting the output digital signal into an analog signal via a D/A converter.
- the orthogonal transformation in this case may be implemented by a dedicated microprocessor, or a special purpose digital signal processor.
- (x(t o ), y(t o ), z(t o )) can be thought of as the position of a particle, or a planet, in three-dimensional space. As time evolves, this particle, or planet, traces out a continuous trajectory. If the motion represented by (x(t), y(t), z(t)) is periodic, as is the case of a planet, this trajectory will eventually repeat itself. However, if the motion is not periodic, as is the case in many practical applications, this trajectory never repeats itself and over a long period of time, fills some volume in space whose overall envelope is a three-dimensional surface S.
- This surface is usually extremely complicated, with cross-sections which are fractals, and to analyze its geometric structure would require a detailed analysis of different perspectives and cross-sections of S on the oscilloscope screen. It is necessary in such analyses to take perspectives and cross-sections because the oscilloscope can display only two dimensional objects.
- Another important operation needed in studying the structure of S is to pass a plane surface D through S at any desired position in space and look at the intersection points between S and S o .
- This set of intersection points is called a cross-section of S with respect to D.
- such cross-sections are fractals.
- One possible approach for implementing the above task is to sample the signals with a high-speed analog-to-digital (A/D) converter and then process the data mathematically by a digital computer before outputting the transformed data into a digital-to-analog (D/A) converter for final display on the oscilloscope.
- A/D analog-to-digital
- D/A digital-to-analog
- one major problem with this approach is that accurate and high-speed data acquisition and processing, currently require very expensive hardwares in addition to a dedicated digital computer. Even then, the displayed signal is no longer in real time.
- Another limitation of the digital approach is that some of the detailed nature of the signals may be lost in the course of data acquisition, unless the signals are properly scaled and the resolution of the analog-to-digital converter is sufficiently high.
- An even more serious limitation of the digital approach becomes obvious in the case of complex (non-periodic) dynamics where the sequence of data to be taken is extremely large and may easily exceed the available computer memory.
- the present invention involves a novel electronic circuit for implementing the above-described mathematical operations instantaneously.
- the apparatus may be a plug-in unit added to an existing oscilloscope.
- the apparatus of the invention is completely automatic, in the sense that the user can either request a specific perspective or cross-section, or he can request that the perspectives be obtained over all possible orientations automatically.
- three or more electrical signals that have some relation to each other are fed to a three-dimensional or n-dimensional rotator.
- a three-dimensional rotator may comprise a series of three two-dimensional rotators linked together.
- the outputs from the three-dimensional rotator are then ready for application to a cathode ray oscilloscope.
- the invention can provide either a stationary image, which can be manually changed to achieve stepped rotation, or there can be automatic rotation.
- a reference plane generator to provide a reference plane that can be used as valuable study aid.
- plane cutting means are provided, which can be used to cut the rotating figure or stationary figure in two at any point along the plane to enable further study of such a figure. If the reference plane or cutting plane forms part of the display, it is a good idea to provide a multiplexor for alternating the signals in such a way that the eye reacts to them as if they were being simultaneously generated.
- preamplifiers may be used prior to sending the signals to the three-dimensional rotator. Also, amplification may be used in any degree where desired.
- three or more signals are linked together in such a way that any two of them can be fed to a cathode ray oscilloscope, and rotation of the configuration is provided to give images on the oscilloscope that correspond to either stationary or automatically rotating images.
- FIG. 1 is a block diagram of a generalized circuit embodying the generalized principles of the invention.
- FIG. 2 is a block diagram of the matrix multiplier of FIG. 1, as made up of three sub-matrix multipliers, for application to a three-dimensional system.
- FIG. 3 is a similar block diagram for a four-dimensional system, necessarily employing six sub-matrix multipliers.
- FIG. 4 is a block diagram corresponding to FIG. 1 with the elements of the controller shown.
- FIG. 5 is a detailed block diagram of a circuit embodying the principles of the invention.
- FIG. 6 is a block diagram of the bus structure for the rotators of FIG. 5, showing that a data bus A is shared by the three rotation angle counters, and that data from the bus A is converted into corresponding multiplying constants by a ROM.
- a data bus B supplies the multiplying constants to the three rotators, while a third bus provides address and controlling signals.
- FIG. 7 is a bus timing chart: each counter outputs its data into the bus A according to its address as determined by signals m2 and m3. The data on the bus A is then converted into four multiplying constants by a ROM, where address signals, m0 and m1 are used to select the corresponding function table in the ROM. These data are placed on bus B so that the multipliers can read their multiplying constants.
- FIG. 8 is a diagram of a controlling circuit and reference generator for a device embodying the invention.
- FIG. 9 is a circuit diagram of one of the three identical rotation angle counters of FIG. 5.
- FIG. 10 is a circuit diagram of a function ROM suitable for use in the circuit of FIG. 5.
- FIG. 11 is a detailed two-dimensional multiplier circuit diagram for use in the circuit of FIG. 5.
- FIG. 12 is a diagram of a sample-and-hold and comparator circuit for use in the circuit of FIG. 5.
- FIG. 13 is a set of four diagrams for the geometrical interpretation of a three-dimensional rotation operation, showing at:
- FIG. 14 is a geometrical interpretation of the perpendicular projection and cross-sectioning operation.
- a plane D is defined by a normal direction vector U, and the distance of the plane D from the origin is denoted by r, FIG. 14 shows:
- FIG. 15 is a simplified block diagram of a 3-D rotation instrument embodying the principles of this invention. This view may be considered a simplified version of FIG. 5, for it lacks continuous rotation.
- FIG. 16 is a diagram of a scalar multiplier circuit for use in this invention, the gain of this circuit is determined by a set of digitally-controlled switches S-1 to S-8 (the switches S-4 to S-7 and related components being omitted).
- FIG. 17 is a diagram of a simplified two-dimensional matrix multiplier circuit, the modules m1, m2, m3 and m4 denoting scalar multipliers like that of FIG. 16.
- FIG. 18 is a diagram of sample signals illustrating the display of the projection of the upper and lower surfaces of S.
- X 0 (t), Y 0 (t), and Z 0 (t) are sample waveforms
- C(t) a comparator output signal obtained by comparing Z 0 (t) with a dc signal of magnitude r
- U and L denote the time intervals corresponding to the upper and lower surfaces, respectively.
- FIG. 19 comprises a pair of views showing sample signals and illustrating the display of the cross-section of S in both backward and forward directions, the solid curves denoting the output waveforms of the sample-and-hold circuits.
- FIG. 19 shows at
- the broken lines indicate the original input signals, which were clamped to a constant level during its entire hold "H" time interval.
- FIG. 20 comprises a pair of views for a driven series R-L-diode circuit, showing at
- FIG. 22 comprises a set of four related views, three of which are projections of a rotated period-three limit cycle, with prescaling the same as for FIG. 21.
- the views are
- FIG. 23 comprises a set of four related views illustrating the cross-section of the period-three limit cycle, the prescaling being the same as for FIG. 21 and the rotation the same as for FIG. 22.
- the views are:
- FIG. 24 comprises a set of four related views showing the projection of the attractor located at the first chaotic band in FIG. 26, the prescaling being the same as in FIG. 21, the rotation the same as in FIG. 22.
- the input signal V s is a 8.6 KHz triangular waveform. The views are
- FIG. 25 comprises a set of six related views showing the cross-section of the attractor at the first chaotic band in FIG. 26. These pictures denote cross-sections of the rotated attractor in FIG. 24(a). The views show:
- FIG. 26 comprises three related views of a bifurcation tree and magnified backward cross-sections. The views are
- FIG. 27 comprises a pair of views relating to Chua's circuit. They show at
- the horizontal scale is 1 msec per divisions, and the vertical scale is: 5 V per division.
- FIG. 29 is a set of three related projections of the unrotated double scroll attractor, the prescaling being the same as in FIG. 28, the horizontal scale being 2 V per division and the vertical scale 2 V per division.
- the views are:
- the views are:
- FIG. 33 comprises a set of three related views showing the forward and the backward cross-sections of the double scroll attractor. Prescaling: same as in FIG. 31. Horizontal scale: 2 V per division. Vertical scale: 2 V per division.
- the views (a), (b), and (c) are the forward and the backward cross-sections at S1, S2 and S3 in FIG. 31(b), respectively.
- the views (a) and (b) are projections of the forward cross-sections S1 and S2 onto the X o -Y o plane, respectively.
- FIG. 36 comprises a set of three CRT displays, being double exposures of the double scroll attractor and a cross-section: Prescaling: same as FIG. 31.
- the views are:
- the invention is an electronic circuit for implementing orthogonal matrix transformations in any dimension. It may be termed an analog matrix multiplier.
- FIG. 1 there are at least two analog input signals X 1 (t), X 2 (t) and so on to X N-1j (t) and X N (t), fed into a matrix multiplier 50.
- the total analog input X(t) over the total time T is
- At least two analog output signals are obtained from the matrix multiplier 50. They may be termed Y 1 (t), Y 2 (t), etc., and the total output signal over the time T is
- I is the identity matrix
- Matrix multiplier constant controlling means control each element value of each matrix multiplication operator q 1 , q 2 , etc.
- the output signals result, as said, from the matrix multiplier's operation on the input signals.
- Each multiplier constant of the matrix multiplier is controlled so that the matrix operator is always orthogonal.
- the matrix multiplier 50 of FIG. 1 is controlled by a controller 51 that is supplied with a rotation angle input 52, and the matrix multiplier 50 is adapted to the number of dimensions, i.e., number of input signals. In each instance there are several sub-matrix or two-dimensional multipliers.
- N 3, as shown in FIG. 2, there are three sub-matrix multipliers 53, 54, and 55, because there are three possible combinations: [1,2], [1,3], and [2,3].
- six sub-matrix multipliers 53, 54, 55, 56, 57, and 58 must be connected, because six distinct combinations [1,2], [1,3], [1,4], [2,3], [2,4] and [3,4] are possible.
- the matrix multiplier controller 51 controls each sub-matrix multiplier 53, 54, 55, etc., so that each sub-matrix R k is an orthogonal two-dimensional matrix.
- the multiplication constant controller 51 of FIG. 1 comprises a micro-controller 60 that executes a sequence of instructions which are stored in a memory 61 and generates all multiplication constants of the matrix multiplier 50, such that the matrix Q becomes orthogonal with the given rotation angles which are entered by a data entry means 62.
- the instructions of the micro-controller 60 are stored in the memory 61 which performs the following transformation digitally:
- control signals from the micro-controller 60 are fed to the matrix multiplier array 50 by suitable interfacing means 63.
- the matrix Q is broken down into a sequence of matrix multiplications. What is looked for are some convenient parameters which uniquely determine the matrix Q. Usually these parameters are rotation angles, because these parameters are well understood if the dimensions are fewer than four.
- the transform is a combination of rotations around each axis. In the three-dimensional case, this observation is true, because after choosing two coordinates, only one remains, namely, the rotation axis. However, it is not true in the general case. There, the transform is a combination of rotations in two-dimensional sub-spaces. While a two-dimensional sub-space can be defined, an axis of rotation, when the dimension exceeds three, cannot be defined. Since all combinations to generate an orthogonal matrix must be performed n!/2(N-2)! combinations (or sub-spaces) must be included.
- any orthogonal matrix Q can be generated.
- the orthogonal matrix can be fully controlled by six parameters, i.e., six rotation angles for each two-dimensional sub-space.
- the apparatus of this invention is an electronic instrument for displaying any perspective of a three dimensional surface S generated by three time-varying (not necessarily periodic) signals.
- the surface S is a three dimensional Lissajous figure which need not be a closed curve; it is typical of all strange attractors.
- This analog (not digital) instrument is designed to rotate S along any axis (not just the X, Y, Z-axis) through any prescribed solid angles (0°-360°) in the three dimensional coordinate system in real time.
- the particular instrument herein described works as a preprocessor for a standard oscilloscope and is built with components capable of displaying time-varying signals with a frequency spectrum from 0 to 20 K Hz. Different frequency spectra may be used by modifying the instrument.
- three electrical signals X i , Y i , and Z i are fed into terminals 65, 66, and 67. These terminals preferably feed these signals to respective preamplifiers 68, 69, and 70. It is not always necessary to use preamplifiers, but often their use is helpful to ensure uniformity of the input signal levels. From the preamplifiers these signals respectively go to ganged switches 71, 72, and 73. When these switches are in their upper position, as they are shown in FIG. 5, the Z signal from the switch 73 goes to a two-dimensional rotator 74, which also receives the Y signal from the switch 72.
- the rotator 74 then becomes a YZ two-dimensional rotator which provides signals corresponding to a projection onto the YZ plane.
- the X signal from the switch 71 goes to another two dimensional rotator 75, which also receives a Z signal from the rotator 72, switch 59.
- the X signal from the rotator 75 and the Y signal from the rotator 74 are fed to a third two-dimensional rotator 76.
- timing generator 77 which includes also a reference plane generator, that feeds a Y reference staircase signal to the switch 72 and an X reference staircase signal to the switch 71, so that when the ganged switches 71, 72, and 73 are thrown to their other position, the timing generator signals X ref and Y ref pass to the rotators 75 and 74, respectively, instead of the respective X and Y signals from the pre-amplifiers 69 and 70.
- the timing generator 77 (as will be explained below in reference to FIG. 8) also sends a signal to three counters 78, 79, and 80 in sequence. These counters will be further explained below with reference to FIG. 9.
- the counter 78 provides the angles ⁇ x ; the counter 79 provides the angles ⁇ y ; and the counter ⁇ z provides the angles ⁇ z .
- the outputs ⁇ x , ⁇ y , and ⁇ z from these three counters 78, 79, and 80 are fed to a bus A which feeds them to a trigonometric function ROM 81.
- the ROM 81 (as will be explained below with reference to FIG. 10) itself sends multiplication signals to a bus B which takes them to the three two-dimensional rotators 74, 75, and 76, as is explained below with reference to FIGS. 6 and 7. These are the three controlling signals for the two-dimensional rotators.
- the X output from the rotator 76 goes to a sample and hold unit 82, while the Y signal from the rotator 76 goes to a sample and hold unit 83, and the Z signal from the rotator 75 goes to a sample and hold unit 84.
- sample and hold units 82, 83, and 84, are controlled by a sample-and-hold control circuit 85.
- the control unit 85 receives two signals. A current is fed to a terminal 86 and goes via a potentiometer 87, which determines intensity, to the control unit 85.
- the control unit 85 also receives a signal from a comparator 88.
- a Z reference signal may be taken just ahead of the switch 73 via a line 89 and sent to a normally open terminal of a switch S D .
- a Z signal may be taken between the rotator 75 and the sample and hold unit 84 via a line 90 and sent to a normally closed terminal of the switch S D .
- the switch 73 has a terminal with a lead to a switch 91 which normally goes to ground; however, the switch 91 is ganged with the switch S D , so that either the Z signal in the lead 90 or the Z reference signal from the lead 89 goes to the plus side of the comparator 88.
- the switch 91 When the switch S D is thrown to open the connection between the Z signal 67 and the rotator 74, the switch 91 then connects the switch 73 to a point r, which connects a potentiometer to the negative side of the comparator 88.
- the alternation of the ganged switches S D and 80 acts to display the main X, Y, Z signals in alternation with a reference plane. This can be done so quickly that, to the eye, they appear to be displayed simultaneously on the oscilloscope.
- data bus A and data bus B are used to control the rotators 74, 75, and 76.
- the data bus A is shared by the three rotation-angle counters 78, 79, and 80. Data from the data bus A are converted into corresponding multiplying constants by the ROM 81. Data from a bus B supplies these multiplying constants to the rotators 74, 75, and 76 respectively.
- address and controlling signals 92 are provided to all of the counters 78, 79, and 80, the rotators 74, 75, and 76 and the ROM 81. Because of the use of the bus structure, i.e., buses A, B, and 92, extensions to higher dimensional rotators can be done easily without changing the interface; what is given here is simply an example.
- FIG. 7 is a bus timing chart.
- Each counter 78, 79, and 80 outputs its data into the bus A according to the addresses determined by the rotator element address signals m2 (LSB) and m3 (MSB).
- the data on the bus A is converted into four multiplying constants by the ROM 81, where rotator address signals m0 (LSB) and m1 (MSB) are used to select the corresponding function table in the ROM.
- the bus A sets the rotation angle at
- the signal wr means "write-enable when wr is false". These data are placed on the bus B, so that the multiplier can read their multiplying constants.
- FIG. 7 shows when the address signals m0, m1, m2, and m3 are on and off.
- the bus A shows the angle of ⁇ y , which is on all the time in this example, as compared with the angles for ⁇ x and ⁇ z and another for ⁇ y while that says in effect, "don't care”.
- Bus B changes from the cosine of ⁇ y to +sine of ⁇ y to -sine of ⁇ y to the cosine of ⁇ y again. These changes result in rotation of the main signals.
- FIG. 8 is a diagram of the controlling circuit and reference generator of this device. It includes four LS 191 chips and several logic elements.
- the first LS 191 chip 93 generates the four address signals m0, m1, m2, and m4 through four inverters 94, 95, 96, and 97 and a write-enable signal wr.
- the final LS 191 chip 98 provides the clock signal RCK for the angle counters 78, 79, and 80.
- the third LS 191 chip 99 generates a multiplex signal MPX and cooperates with the second LS 191 chip 100 to generate two staircase signals X ref and Y ref for displaying a reference plane.
- FIG. 9 is a circuit diagram of one of the rotation angle counters and gives an example of the rotation angle entry circuit. It comprises two LS 191 chips 102 and 103 (which are binary counters, a synchronous up-down counter, a standard chip which is widely available in different power ratings), a preset data array 104 of switches and resistors, and an LS 244 chip 105 with eight-bits of output. (This is a standard component, widely available; sometimes called an octal buffer or a line driver.) This is used when continuous rotation is desired. For continuous rotation, to change the angle continuously in time, the LS 191 counters count the pulses from the clock RCK and increment or decrement the angle data ⁇ . If continuous rotation is not required, the circuit of FIG. 9 may be replaced by an eight-bit switch.
- FIG. 10 is a circuit diagram of an eight-bit function ROM 81 that can be used in the invention.
- the address of the ROM 81 is classified by an A0 to A7 terminal for input for the data bus A and D0 to D7 output to the data bus B. There is also input for the signals m0 and m1.
- the data of bus B is simply output from the ROM memory.
- the actual data of the ROM is shown in Table 1.
- Table II explains the relationship between the address signals m0, m1 and the output data.
- FIG. 11 is a detailed two-dimensional multiplier circuit diagram, i.e., a diagram of one of the three identical two-dimensional rotators 74, 75, 76. These have an inherent accuracy of one part in 256, better than 0.4%.
- the circuit is controlled by a through bus ("DATA") similar to a common microcomputer interface. This is applied to a set of four multiplying DAC, such as AD7524 made by Analog Devices.
- the address decoder may vary from application to application. All resistors shown in the circuit diagram should have at least 0.5% accuracy so that multiplication is sufficiently accurate and commensurate with that of the rotator.
- FIG. 12 is a diagram of a typical and usable sample-and-hold comparator circuit, including units 82, 83, 84, 85, and 88.
- Sample-and-hold units are, of course, obtained off the shelf, e.g., LF 398 for units 82, 83, 84, by National Semiconductor Corporation, as are comparators, that may be a LM111, also by National Semiconductor Corporation.
- the circuit is very simple the corresponding integrated circuit reference should be checked to avoid unexpected effects.
- polystyrene capacitors are preferably used for the 0.001 ⁇ f connection from each of the units 82, 83, 84 to ground. Additional caution should be employed during the wiring process to minimize the stray capacitors from the units 82, 83, and 84 to the 0.001 ⁇ f capacitors.
- a small portable electronic instrument When used with a standard oscilloscope, it enables display of a two-dimensional projection from any desired direction of a three-dimensional surface S evolving in time.
- the three-dimensional surface S represents the dynamic loci of three electrical signal x(t), y(t), and z(t) in the X-Y-Z space over some observation time interval, and can be interpreted as the three-dimensional generalization of the well-known Lissajous figures.
- the circuit herein described is at present the least expensive approach capable of yielding good accuracy for low-frequency applications, i.e., below 10 KHz.
- the instrument described here is an analog instrument that implements both rotation and cross-section operations in real time.
- Rotation of a vector S(t) is simply an orthogonal transform of S(t), namely,
- Q is an orthogonal three-dimensional matrix and S o (t) is the rotated vector y[x o (t),y o (t),z o (t)] T .
- the orthogonal matrix Q implemented in the design of this invention is defined by its decomposition or breakdown with respect to the X, Y, and Z axis; namely,
- FIG. 13 shows the geometrical transformation of this operation.
- the output coordinates rotate around the X o -axis by ⁇ x in a counter-clockwise direction. This is equivalent to rotating the brick B by ⁇ x in a clockwise direction in the output coordinates, as shown in FIG. 13(b).
- R y ( ⁇ y ) Applying the operation R y ( ⁇ y ) to the ⁇ x -rotated brick in FIG.
- FIG. 13(b) is equivalent to further rotating the coordinate along the Y o -axis by ⁇ y degree (clockwise), and hence the brick actually rotates in a counterclockwise direction as shown in FIG. 13(c).
- R z ( ⁇ z ) is equivalent to further rotating the coordinate along the Y o -axis by ⁇ y degree (clockwise), and hence the brick actually rotates in a counterclockwise direction as shown in FIG. 13(c).
- the vector U indicates the normal direction of the plane D, and r is the distance of D from the origin.
- the cross-section of a surface S on a plane D is the set of intersection points such that
- the forward cross-section of S(t) on a plane D is the set of intersection points such that S(t) crosses D from the lower surface to the upper surface for all t>0.
- a backward cross-section is similarly defined (See FIG. 14(b)).
- FIG. 15 shows a simplified block diagram of the instrument.
- Three identical two-dimensional rotators 124, 125, and 126 are used to perform the operations R x ( ⁇ x ),R y ( ⁇ y ), and R z ( ⁇ z ). Observe that the connections of the three two-dimensional rotators 124, 125, and 126 in FIG. 15 realizes the orthogonal transformation Q in Equation 2.
- U the unit vector [1,0,0] t , [0,1,0] T , or [0,0,1] T and then applying the perpendicular projection (2.4) to S o (t). Since S o (t) can be rotated by any desired angle, the projection of S(t) from any direction can be displayed on the oscilloscope screen.
- the three sample-and-hold (S&H) circuits 128, 129, and 130 and a comparator 131 on the right side of FIG. 15 are all that is needed for displaying any desired cross-section on D.
- the plane D is specified by choosing U in Equation 3 to be the unit vector [1,0,0] T , [0,1,0] t , or [0,0,1] T .
- the plane D can be defined, with the help of a switch S D , in terms of either the input coordinates (X,Y,Z,) or the output (rotated) coordinates (X o , Y o , Z o ). (The switch S D in FIG.
- U [0,0,1] T for convenience.
- a six-position switch (three for input coordinates, and three for output coordinates) is recommended.) Since U has only one nonzero element, the times when S(t) intersects D can be easily and accurately detected by comparing one of the signals, say Z(t), and a constant r. Since the signals are time-varying, their values at the intersection point should be held for a short period of time. Three sampled-and-hold circuits are used for this purpose.
- FIG. 15 Most of the functional blocks in FIG. 15 are implemented by using off-the-shelf integrated circuit modules. The detailed circuit for the design depends on the choice of integrated circuit modules.
- any rotation circuit for implementing the transformation of Equation 2 is a multiplier.
- the rotation angle is either a constant, or is changing continuously but slowly, if the surfaces is to be rotated continuously through all angles.
- This simple design is inaccurate and cumbersome in practice because the adjustments are tedious and error prone.
- a multiplying digital-to-analog converter is employed as the basic multiplier.
- the MDAC is a resistor array with electronically controllable switches (for example, MOS transistors). Hence, it can be used either as a multiplier, or as a variable-gain amplifier when incorporated with a summing operational amplifier. Eight- to twelve-bit MDACs are now widely available.
- FIG. 16 shows a detailed circuit diagram of this scalar multiplier 132 using an MDAC.
- an eight-bit R-2R ladder network generates a set of binary-weighted currents.
- Digitally controlled switches, S-1 to S-8, are used to switch these currents between ground and the operational amplifier input which is at virtual ground. Since these currents are proportional to the input voltage V i , the output voltage is given by
- the gain between the input and the output varies from 0/256 to 255/256.
- the multiplier constant can be adjusted to any value between 0.0 to 1.0, and is accurate to within the discretization error.
- FIG. 17 is a circuit diagram for a two dimensional matrix multiplier consisting of four scalar multipliers M1, M2, M3, and M4, each like the multiplier 120.
- Two operational amplifiers 133 and 134 are connected as adders to implement the summing operation required in the matrix multiplication.
- the relation between the inputs and the outputs of this circuit is ##EQU3## where r ij can assume one of 512 discrete values between -1.0 to 1.0.
- Conversion from a rotation angle ⁇ to the multiplication constants r ij can be implemented by either using a microcomputer or other digital circuit techniques.
- a Read-Only Memory (ROM) is used to map ⁇ into r ij . This approach is both simple and efficient because it is equivalent to looking up a trigonometric function table electronically.
- the mode of the sample-and-hold circuits 127, 128 and 129 is controlled by a comparator 131, which detects the time t k such that one of signals, say Z o (t k ), is equal to the reference voltage r.
- a comparator 131 By holding the value of the signals t and the time t k , the intersection point on the oscilloscope screen can be observed. If the signals are sufficiently fast, the intersections, although occurring at different times, will appear to the eye as a set of points on the cross-section of S on the reference plane D.
- the comparator output can also be used to control the beam intensity of the oscilloscope so that some part of the trajectory is blanked out from the screen.
- FIG. 18 is an example of such a selective display.
- X(t), Y(t), and Z(t) are sample waveforms. If the comparator is used to compare Z(t) and the constant r, its output C(t) will be as shown in FIG. 18. By applying C(t) or inverted C(t) to the intensity input of an oscilloscope, only the upper or the lower surface of S relative to D in Equation 3 will be displayed.
- FIG. 19 shows an example of a cross-section display. If the sample-and-hold circuits are controlled by C(t), their output waveforms, X k (t), Y k (t), and Z k (t) corresponding to the reference value r In FIG. 18 will appear as shown by the solid curves on FIG. 19. The dotted curves denote the original waveform. In this case, the values of X o (t), Y o (t), and Z o (t) are held whenever Z o (t) crosses the constant value r from the upper side. Hence, when the signals during the holding periods are displayed on the screen, a backward cross-section on D is obtained. In FIG. 19(b), the signals are held whenever Z o (t) crosses the constant value r from the lower side. The corresponding signals will display the forward cross-section with respect to D.
- the other circuit functions that are not shown in the block diagram in FIG. 15 include a controlling circuit for rotation angles, a reference plane generator, and a signal-reference plane multiplexer block, all shown in FIG. 5. These blocks are included for the user's convenience.
- Strange attractors of a second-order nonautonomous circuit can be displayed using our 3-D rotation instrument. They represent snap shots when the strange attractor S--which is being viewed when S is rotated continuously and periodically through all angles--is frozen momentarily at different selected angles of rotation.
- a conventional method to characterize the periodic orbit or strange attractor in such a system is to sample the waveforms at the input frequency for different phase angles of the input signal. Such a cross-section is called a Poincare map.
- the instrument can be used not only to display these conventional Poincar/e maps, but also to look at Poincare maps at different perspectives.
- the excitation signal is chosen as the third variable, although it is not a state variable. Since the excitation signal is a periodic function of time, this third axis can be interpreted as a periodically-folded map of time.
- FIG. 20 (b) shows a coordinate system for the nonautonomous case. This simple set-up provides an extremely flexible and powerful tool for analyzing the fine geometrical structure of strange attractors that was not possible before.
- a cutting plane can be defined in any desired position in the three-dimensional space.
- One advantage of this approach is that by observing the three-dimensional surface S from only a few appropriate directions, the local structure of S can usually be determined without taking discrete Poincare maps at many phases of the input signal.
- the cutting plane is chosen to be orthogonal to the excitation axis, then the cross-section becomes the conventional Poincare maps. In other words, the instrument is much more versatile than just taking Poincare maps.
- Reference plane--- As a frame of reference and reminder of the relative position of the projection plane and the input (unrotated) X-Y plane,--henceforth called the reference plane--, which need not coincide with the cutting plane D, this reference plane can be displayed simultaneously with the projected surface S or the cross-section. The first quadrant of this reference plane will always appear (when requested) as a uniform grid of points on the oscilloscope screen.
- the periodic waveform I D (t), V D (t), and V S (t) (at a fixed input frequency) are displayed in FIG. 21(a), along with the signal C(t) from the comparator. Note that the period of V D (t) is three times that of V s (t).
- the projection of the surface S generated by I D , V D , and V S onto the V D -I D plane is shown in FIG. 21(b). Note the "barely visible" 3-loop Lissajous figure indicates a period-3 relation between V D and I D . Here, the output signal is at the input of the sample-and-hold circuits and hence C(t) is irrelevant.
- FIG. 22 shows the position of the rotated V D -I D plane.
- the cutting plane D in this case is a plane parallel to the reference plane but passing through the V S -axis at V S -r.
- FIG. 22(b) shows the projection of the rotated surface S onto the X o -Y o plane (the oscilloscope screen), along with the position of the rotated reference plane. Note that the period-3 nature of V D (t) is much clearer than that of FIG. 21(b). This demonstrates the power of rotation provided by this instrument: with an appropriate choice of rotation angles ( ⁇ x , ⁇ y , ⁇ z ), the projection reveals many more details of the attractor that were otherwise hidden (as in FIG. 21(b) prior to rotation.
- FIG. 22(c) and 22(d) show the projection of the upper and lower surface achieved by blanking out the other portion with the held of the blanking signal C(t). Note that if FIG. 22(d) is superimposed on FIG. 22(c), a complete projection per FIG. 22(b), would be recovered, as expected.
- a second application of the comparator signal C(t), is to display the backward and forward cross-sections of the above period-3 attractor relative to the cutting plane D through V s -r, as depicted in FIG. 23(a).
- the signal I D (t), V D (t), and V S (t) at the output of the sample-and-hold circuits are shown in FIG. 23(b) along with C(t).
- the backward and forward cross-section shown in FIG. 23(c) and 23(d) comprises three isolated points indicating V D (t) is period-3 relative to V S (t).
- FIGS. 24(a), 24(b), and 24(c) show the projection of the rotated chaotic attractor (in the first chaotic band in FIG. 26(a)) on the three orthogonal planes (X o -Y o , Y o -Z o , or Z o -X o planes), where the input signal V S is a triangular waveform.
- V S is a triangular waveform.
- FIG. 26(c) Another cross-section measured with the input frequency chosen to lie within the second chaotic band in FIG. 26(a) is shown in FIG. 26(c). This three-leg attractor corresponds to a similar attractor elsewhere reported.
- the waveforms of the three state variables, -I L , V C1 and V C2 are shown in FIG. 28 where the horizontal scale is 1 msec per division. Note that all variables oscillate at an approximate rate of 10 Khz and that the waveforms appear to be chaotic.
- FIGS. 29(a), (b) and (c) show the "unrotated" projection of the double scroll onto the (-I L )-V C1 , (-I L )-V C2 , and the V c1 -V c2 planes, respectively.
- the relative position of the reference plane is depicted in FIG. 30(d).
- FIG. 32 Several cross-sections of the double scroll on D for different values of r are shown in FIG. 32.
- FIG. 31(b) shows the projection of the double scroll and its four forward cross-sections onto the (-I L )-V c1 plane.
- FIGS. 32(a), (b), (c) and (d) give the projection of the forward cross-sections S1, S2, S3, and S4 onto the (-I L )-V c2 , plane.
- a print image of the projection of the double scroll onto the (-I L )-V c2 plane is also displayed to show the relative location of these cross-sections and that of the double scroll. Note that the near-circular contour corresponds to the projection of the lower "hole" of the double scroll attractor.
- FIGS. 35(a) and (b) The projection of the two forward cross-sections S1 and S2 of the rotated double scroll on the cutting planes depicted in FIG. 34(a) are shown in FIGS. 35(a) and (b), respectively. Note that they are merely rotations by the same angle ⁇ z of the corresponding cross-sections in FIGS. 32(a) and (b), respectively, as expected.
- FIG. 36 shows the superposition of a projection and a cross-section of the double scroll after it has been rotated by some angle ⁇ x , ⁇ y , and ⁇ z .
- the cutting plane in this figure coincides with the projection plane (X o -Y o plane).
- the invention thus provides a real-time preprocessor for the oscilloscope for rotating three-dimensional Lissajous figures through any angle, either fixed at some discrete value, or varying continuously in a periodic manner.
- Applications of this three-dimensional rotation instrument for observing strange attractors are illustrated for both autonomous and non-autonomous circuits. Many of the projections and cross-sections of strange attractors have never been seen before and reveal many new insights into the local structure of strange attractors. Many of these cross-sections confirmed observations which were obtained previously through hundreds of hours of brute-force computer simulation.
- the simplified instrument described above can be built using only low-frequency inexpensive integrated circuit modules. More careful circuit design and choice of components can result in a significantly increased operating frequency range and accuracy. Also, the controlling circuit may be replaced by a micro-computer or interfaced with the host controller for increased flexibility.
- comparators may be helpful in dividing the surface S into several regions simultaneously. This feature increases the flexibility of our measurement process.
Abstract
N!/2 (N-2)!
Description
X(t)=[X.sub.1 (t), X.sub.2 (t) . . . X.sub.N (t)].sup.T.
Y(t)=[Y.sub.1 (t), Y.sub.2 (t) . . . Y.sub.N (t)].sup.T.
Y1=r.sub.11 X1+r.sub.12 X2 and
Y2=r.sub.21 X1+r.sub.22 X2
(θ=Dxπ/128)
(r=D/128-1.0).
TABLE 1 __________________________________________________________________________ 000H 0FFH COS(θ)(SAME FOR 300H 3FFH) O: FF FF FF FF FE FE FE FD FD FC FB FA FA F9 F8 F6 10 F5 F4 F3 F1 F0 EE ED EB EA E8 E6 E4 E2 E0 DE DC 20: DA D7 D5 D3 D0 CE DB C9 C6 C4 C1 BE BC B9 B6 B3 30: B0 AD AA A7 A5 A2 9E 9B 98 95 92 8F 8C 89 86 83 40: 80 7C 79 76 73 70 6D 6A 67 64 61 5D 5A 58 55 52 50: 4F 4C 49 46 43 41 3E 3B 39 36 34 31 2F 2C 2A 28 60: 25 23 21 1F 1D 1B 19 17 15 14 12 11 F E C B 70: A 9 7 6 5 5 4 3 2 2 1 1 1 0 0 0 80: 0 0 0 0 1 1 1 2 2 3 4 5 5 6 7 9 90: A B C E F 11 12 14 15 17 19 1B 1D 1F 21 23 A0: 25 28 2A 2C 2F 31 34 36 39 3B 3E 41 43 46 49 4C B0: 4F 52 55 58 5A 5D 61 64 67 6A 6D 70 73 76 79 7C C0: 80 83 86 89 8C 8F 92 95 98 9B 9E A2 A5 A7 AA AD D0: B0 B3 B6 B9 BC BE C1 C4 C6 C9 CB CE D0 D3 D5 D7 E0: DA DC DE E0 E2 E4 E6 E8 EA FB ED EE F0 F1 F3 F4 F0: F5 F6 F8 F9 FA FA FB FC FD FD FE FE FE FF FF FF 100H 1FFH SIN(θ) 0: 80 83 86 89 8C 8F 92 95 98 9B 9E A2 A5 A7 AA AD 10: B0 B3 B6 B9 BC BE C1 C4 C6 C9 CB CE D0 D3 D5 D7 20: DA DC DE E0 E2 E4 E6 E8 EA EB ED EE F0 F1 F3 F4 30: F5 F6 F8 F9 FA FA FB FC FD FD FE FE FE FF FF FF 40: FF FF FF FF FE FE FE FD FD FC FB FA FA F9 F8 F6 5O: F5 F4 F3 F1 F0 EE ED EB EA E8 E6 E4 E2 E0 DE DC 60: DA D7 D5 D3 D0 CE CB C9 C6 C4 C1 BE BC B9 B6 B3 70: B0 AD AA A7 A5 A2 9E 9B 98 95 92 8F 8C 89 86 83 80: 80 7C 79 76 73 70 6D 6A 67 64 61 5D 5A 58 55 52 90: 4F 4C 49 46 43 41 3E 3B 39 36 34 31 2F 2C 2A 28 A0: 25 23 21 1F 1D 1B 19 17 15 14 12 11 F E C B B0: A 9 7 6 5 5 4 3 2 2 1 1 1 0 0 0 C0: 0 0 0 0 1 1 1 2 2 3 4 5 5 6 7 9 D0: A B C E F 11 12 14 15 17 19 1B 1D 1F 21 23 E0: 25 28 2A 2C 2F 31 34 36 39 3B 3E 41 43 46 49 4C F0: 4F 52 55 58 5A 5D 61 64 67 6A 6D 70 73 76 79 7C 200H 2FFH -SIN(θ) 0: 80 7C 79 76 73 70 6D 6A 67 64 61 5D 5A 58 55 52 10: 4F 4C 49 46 43 41 3E 3B 39 36 34 31 2F 2C 2A 28 20: 25 23 21 1F 1D 1B 19 17 15 14 12 11 F E C B 30: A 9 7 6 5 5 4 3 2 2 1 1 1 0 0 0 40: 0 0 0 0 1 1 1 2 2 3 4 5 5 6 7 9 50: A B C E F 11 12 14 15 17 19 1B 1D 1F 21 23 60: 25 28 2A 2C 2F 31 34 36 39 3B 3E 41 43 46 49 4C 70: 4F 52 55 58 5A 5D 61 64 67 6A 6D 70 73 76 79 7C 80: 80 83 86 89 8C 8F 92 95 98 9B 9E A2 A5 A7 AA AD 90: B0 B3 B6 B9 BC BE C1 C4 C6 C9 CB CE D0 D3 D5 D7 A0: DA DC DE E0 E2 E4 E6 E8 EA EB ED EE F0 F1 F3 F4 B0: F5 F6 F8 F9 FA FA FB FC FD FD FE FE FE FF FF FF C0: FF FF FF FF FE FE FE FD FD FC FB FA FA F9 F8 F6 D0: F5 F4 F3 F1 F0 EE ED EB EA E8 E6 E4 E2 E0 DE DC E0: DA D7 D5 D3 D0 CE CB C9 C6 C4 C1 BE BC B9 B6 B3 F0: B0 AD AA A7 A5 A2 9E 9B 98 95 92 8F 8C 89 86 83 __________________________________________________________________________
TABLE II ______________________________________ m1 m0 Meaning Output data ______________________________________ 0 0r11 cosθ 0 1r12 sinθ 1 0 r21 -sinθ 1 1 r22 cosθ. ______________________________________
S.sub.o (t)=QS(t),Equation 1
Q=R.sub.z (θ.sub.z)R.sub.y (θ.sub.y)R.sub.x (θ.sub.x),Equation 2
D={D(x,y,z)|DU.sup.T =r},Equation 3
S.sub.P =S+(r-SU.sup.T)U.
S.sup.+ ≡{S|S(t)U.sup.T -r>0}.Equation 5
S.sup.- ≡{S|S(t)U.sup.T -r<0}.Equation 6
{S|S(t)U.sup.T -r=0}.Equation 7
{S|S(t)U.sup.T -r=0 and S(t)U.sup.T >0}.Equation 8
V.sub.o =-(D.sub.7 /2.sup.1 +D.sub.6 /2.sup.2 +. . . +D.sub.0 /2.sup.8)V Equation 10
Claims (20)
N!/2(N-2)!
Y1=r.sub.11 X1+r.sub.12 X2,
Y2=r.sub.21 X1+r.sub.22 X2,
Y1=r.sub.11 X1+r.sub.12 X2,
Y2=r.sub.21 X1+r.sub.22 X2,
Y1=r.sub.11 X1+r.sub.12 X2,
Y2=r.sub.21 X1+r.sub.22 X2,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US06/845,855 US4819200A (en) | 1986-03-28 | 1986-03-28 | Real-time N-dimensional analog signal rotator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US06/845,855 US4819200A (en) | 1986-03-28 | 1986-03-28 | Real-time N-dimensional analog signal rotator |
Publications (1)
Publication Number | Publication Date |
---|---|
US4819200A true US4819200A (en) | 1989-04-04 |
Family
ID=25296246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US06/845,855 Expired - Fee Related US4819200A (en) | 1986-03-28 | 1986-03-28 | Real-time N-dimensional analog signal rotator |
Country Status (1)
Country | Link |
---|---|
US (1) | US4819200A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101826839A (en) * | 2010-04-19 | 2010-09-08 | 浙江大学 | Inverter-based chaotic oscillating circuit |
US10830545B2 (en) | 2016-07-12 | 2020-11-10 | Fractal Heatsink Technologies, LLC | System and method for maintaining efficiency of a heat sink |
US11598593B2 (en) | 2010-05-04 | 2023-03-07 | Fractal Heatsink Technologies LLC | Fractal heat transfer device |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4084136A (en) * | 1976-10-21 | 1978-04-11 | Battelle Memorial Institute | Eddy current nondestructive testing device for measuring variable characteristics of a sample utilizing Walsh functions |
-
1986
- 1986-03-28 US US06/845,855 patent/US4819200A/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4084136A (en) * | 1976-10-21 | 1978-04-11 | Battelle Memorial Institute | Eddy current nondestructive testing device for measuring variable characteristics of a sample utilizing Walsh functions |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101826839A (en) * | 2010-04-19 | 2010-09-08 | 浙江大学 | Inverter-based chaotic oscillating circuit |
US11598593B2 (en) | 2010-05-04 | 2023-03-07 | Fractal Heatsink Technologies LLC | Fractal heat transfer device |
US10830545B2 (en) | 2016-07-12 | 2020-11-10 | Fractal Heatsink Technologies, LLC | System and method for maintaining efficiency of a heat sink |
US11346620B2 (en) | 2016-07-12 | 2022-05-31 | Fractal Heatsink Technologies, LLC | System and method for maintaining efficiency of a heat sink |
US11609053B2 (en) | 2016-07-12 | 2023-03-21 | Fractal Heatsink Technologies LLC | System and method for maintaining efficiency of a heat sink |
US11913737B2 (en) | 2016-07-12 | 2024-02-27 | Fractal Heatsink Technologies LLC | System and method for maintaining efficiency of a heat sink |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4283765A (en) | Graphics matrix multiplier | |
US4467275A (en) | DC characteristics measuring system | |
US3809868A (en) | System for generating orthogonal control signals to produce curvilinear motion | |
US4818931A (en) | Vector analyzer with display markers and linear transform capability | |
US4819200A (en) | Real-time N-dimensional analog signal rotator | |
US2899550A (en) | meissinger etal | |
DE60033794T2 (en) | CALIBRATION OF AN OPTICAL TRANSMITTER FOR POSITION MEASUREMENT SYSTEMS | |
Chua et al. | Three-dimensional rotation instrument for displaying strange attractors | |
US4090244A (en) | Method and circuit arrangement for determining the vector components of an oscillation | |
WO2000032018A1 (en) | Digital phase measuring system and method | |
US5034745A (en) | Data acquisition with vernier control | |
IL33610A (en) | Coordinate converter | |
RU2129953C1 (en) | Apparatus for program control of manipulator | |
US3027485A (en) | Cathode ray tube display | |
Lukes | Oscillographic Examination of the Operation of Function Generators | |
CN114646440B (en) | Three-axis six-degree-of-freedom vibration control input/output coordinate transformation matrix construction method | |
Real | Graphic data compressor for plotting instruments | |
RU2120850C1 (en) | Apparatus for program control of manipulator | |
RU2054349C1 (en) | Device for manipulator programmed control | |
DE3315809A1 (en) | 3-D Sensor cube | |
GB2096367A (en) | Digital signal processing apparatus | |
Ananthanarayanan et al. | Identification of kinematic parameters of multiple closed chain robotic manipulators working in coordination | |
JPS6055837B2 (en) | graphic display device | |
SU769579A1 (en) | Device for displaying coordinate system on crt screen | |
US3280317A (en) | Root locus plotter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: REGENTS OF THE UNIVERSITY OF CALIFORNIA THE, 2199 Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNORS:CHUA, LEON O.;SUGAWARA, TSUTOMU;REEL/FRAME:004558/0517 Effective date: 19860310 Owner name: REGENTS OF THE UNIVERSITY OF CALIFORNIA THE,CALIFO Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:CHUA, LEON O.;SUGAWARA, TSUTOMU;REEL/FRAME:004558/0517 Effective date: 19860310 |
|
CC | Certificate of correction | ||
FEPP | Fee payment procedure |
Free format text: PAYER NUMBER DE-ASSIGNED (ORIGINAL EVENT CODE: RMPN); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 19970409 |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |