Macro and micro view on steady states in state space

03/06/2010 ∙ by Branislav Sobota, et al. ∙ Technical University of Kosice 0

This paper describes visualization of chaotic attractor and elements of the singularities in 3D space. 3D view of these effects enables to create a demonstrative projection about relations of chaos generated by physical circuit, the Chua's circuit. Via macro views on chaotic attractor is obtained not only visual space illustration of representative point motion in state space, but also its relation to planes of singularity elements. Our created program enables view on chaotic attractor both in 2D and 3D space together with plane objects visualization -- elements of singularities.



There are no comments yet.


page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

The visualization is good idea to show imagines, ideas, design, construction, realization or effects. It is also one way of verification before realization our goals. Computer based visualization brings utilization of physical, or simulated electric parameters course graphical interpretation in non-linear circuit theory together with other fields. From beginning it was used 2D visualization with possibility of color utilization to be more illustrative or to explain actions proceeding in non-linear circuits [1, 9, 13, 15]. High-performance or parallel computers enable to take advantage of 3D state space axonometry [14]. Actual available solutions provide high performance visualization suitable for 3D interactive presentation of processes and effects with support for over million saturated colors in hi-resolution mode and for use in all graphics-intensive applications. This paper describes actual possibilities of PC for visualization of steady states of chaos generating circuit.

2 Methods used for trajectory visualization

In last 24 years there was intensive interest of scientific community to analyse and applied Chua’s circuit generating chaos. Presentation of trajectories needs to solve system (2) describing physical Chua’s circuit.




Next we consider control pulse , the resistance of the inductance .

For parameters (2) in [5] there were found chaotic attractors showed in Fig. 1 in Monge projection.


Computer program was designed in C language by author of [2] and used for clarifying of place in state space where chaos originates [3]. It is only short segment of intersection two surfaces related to circuit singularities and , or and . In despite of explaining with help of tables and 2D presentation was definite, 3D visualization provides faster and lighter illustration of actions which proceed in specific non-linear circuit. Therefore 3D visualization is valued as from scientific as from edifying point of view [4].

Figure 1: Monge projection of chaotic attractor system (2), for parameters (2) to plane: a) , b)

3 Visualization in 3-dimensional space

Chaos visualizing system was designed for visualization Chua’s attractor in 3D space in real time and it is based on visualizing kernel developed on DCI FEEI TU Košice [6]. An application is implemented in C++ language using OpenGL graphics library. The application can work with I-V characteristics, Chua’s attractor trajectory, or limit cycle and it can visualize elements of the singularity planes. Additionally this visualization depicts representative point movement and it creates chaotic attractor using two basic modes (continuous and sequential). This application can be used not only for concrete Chua’s circuit. It is usable also for Chua’s circuit like structures analyzed in [7, 8].

Chaos visualizing system provides three basic visualizing modes: continuous mode, sequential mode and I-V characteristics visualization mode. System allows using four projection types (3D projection (see Fig. 5), 2D projection, 2D projection and 2D projection) for better-examined circuit understanding. Settable basic visualizing parameters for chaotic attractor visualization are: drawing speed, points omission, chaotic attractor point size, comet length and attractor colour. The combination of these parameters defines final visualization form of chaotic attractor. In case of 3D graphic accelerator supported acceleration of graphic interface OpenGL is drawing of 3D primitives accelerated by graphic card. It enables increasing performance and using some graphical improvements cannot be used without acceleration in real time. Second way is use of parallel computational system [10, 11] for faster or better visualization.

3.1 Visualization program utilization

3D visualization of chaotic attractor for parameters (2) is shown in Fig. 2. To next manipulation with chaotic attractor as 3D object is necessary to fulfill the following steps:

  • Load input file – chaotic attractor. Input file size is above 500 MB, therefore program enables to choose number of trajectory points, which is loaded from file and consequently displayed.

  • Load I-V characteristics and from files.

  • Set background color for scene displaying and

  • Set colors and line-width of I-V characteristics and chaotic attractor.

Figure 2: Chaotic attractor depicted in continuous mode

Views on chaotic attractor from various sides can be obtained by rotation of camera position horizontally and also vertically around visualized object. Fig. 2a and b show horizontal camera swing out, while Fig. 2c, d show vertical camera swing out to chaotic attractor. In this way top view on observed object can be obtained.

Figure 3: Chaotic attractor depicted in sequential mode with comet effect

On Fig. 3 are displayed different looks to the same chaotic attractor. It is visualizing in sequential mode with comet effect [12]. The comet length is adjustable from 512 to 16384 bits. Fig. 2 and also Fig. 3 show that left and right discs of chaotic attractor are situated in plane. It is possible to display this plane in the program. Mathematical description of this plane presented in [3] is outlined by the following equation:


It is available to define in application menu input parameters for appropriate planes as e.g.: singularity coordinates (

), eigenvectors (

), width, length and planes colors. The elements of the singularity planes are displayed by application using of selected colors. Fig. 4a shows these elements of the singularities and representation as parallel planes getting across singularities and . Fig. 4b shows also the third singularity element , called . Fig. 4b shows, that plane is not parallel with and . Via macro views on chaotic attractor showed in Fig. 2– 4 we obtain not only visual space illustration of representative point motion in state space, but also its relation to planes of singularity elements.

Figure 4: Visualization of singularity elements in state space: a) and b) and

Vertical camera swing round with regard to Fig. 4 enables to see such a part of state space, where chaos arises in Chua’s circuit. It is intersection of two planes and or and . This situation is displayed in macro view on Fig. 5a. Marking of conjunction region of mentioned elements and enables to obtain micro view to just site of state space, where chaos originates.

Figure 5: 3D view on intersection of singularity elements and (view on place of chaos origination in Chua’s circuit)

Website contains dynamic versions of some pictures mentioned in this article. Video appears jerky. This is caused by camera position rotation around object (usually 5°). Program enables to set angle step of camera position rotation. It was used when Fig. 5 was generated.

4 Conclusion

The 3D visualization brings new dimensions to the visualization of physical or electric effects. Visualization of chaotic attractor and elements of the singularities provides better understanding of representative point movement in state space, what was still possible only with help of representation in projection planes. From computer graphics point of view are produced big data-sets. Possible parallel processing (e.g. on multi-core or multi-computers platform) shortens computational time. Big-screen display solutions increase quality and ability of immersion into 3D space. It is possible to use finer integration step for better quality and more detailed visualization output (continuous trajectory displaying) in this case. Generally, there are two important application areas in physical or electric effects visualization where big-display environments are used: displaying images at very high resolution in real time exceeding those of available screens (monitors) and/or graphics cards and providing a larger field-of-view and better immersion into the explored attractor‘s space. Whole system enables visualization for user defined Chua‘s attractor where user has in standard visualizing mode (3D projection) 6 degrees of freedom for motion in explored attractor‘s space, which mean translation in 3 axis and rotating around them. It means that system is capable to visualize the Chua’s attractor from any point and from any angle. Actual available solution provide so unique big wide-screen high performance visualization solution suitable for 3D interactive presentation of attractor‘s space. Input devices are standard keyboard and mouse generally, but space mouse or other specialized input device can be also used.


This work is supported by VEGA grant project No. 1/0646/09: “Tasks solution for large graphical data processing in the environment of parallel, distributed and network computer systems”, by KEGA no. 3/6386/08: “E-learning and web oriented education technology of subjects from the field of electrical measurement for presentation and distance form of study” and by Agency of the Ministry of Education of the Slovak Republic for the Structural Funds of the EU under the project Center of Information and Communication Technologies for Knowledge Systems (project number: 26220120020).


  • [1] P. Galajda, V. Špány, M. Guzan, The state space mystery with negative load in multiple-valued logic, Radioengineering, 8, 2 (1999) 2–7.
  • [2] M. Guzan, Multifunctionality of Chua‘s circuit at R changes, International Conference of Electrotechnic Teachers, 16–18 September 2008, Košice-Herľany, Slovakia, TU Košice, 2008, pp. 61–66 (in Slovak).
  • [3] M. Guzan, P. Galajda, L. Pivka, V. Špány, Element of singularity is a key to laws of chaos, 15th International Czech-Slovak Scientific Conference RADIOELEKTRONIKA 2005, 2005, pp. 33–36.
  • [4] R. Kreheľ, Sensoren in der Prozessautomation und Prozessinformatik, in: CO-MAT-TECH 2004: 12th International Scientific Conference, Trnava, Slovakia, 14–15 October 2004, STU Bratislava, 2004, pp. 665–672.
  • [5] T. Matsumoto, L. Chua, M. Komuro, The double scroll, IEEE Trans. Circuits Syst., Vol. CAS-32, No. 8, 1985.
  • [6] R. Mocnár, Chaotic attractor visualization in state space, Diploma Thesis, DCI FEEI TU Košice, 2006 (in Slovak).
  • [7] J. Petržela, Modeling of the strange behavior of the selected nonlinear dynamical systems, Part I.: Oscillators, PhD Thesis Edition, Vol. 502.
  • [8] J. Petržela, Z. Kolka, S. Hanus, Simple chaotic oscillator: from mathematical model to practical experiment, Radioengineering, 15, 1 (2006) 6–11.
  • [9] L. Pivka, V. Špány, Boundary surfaces and basin bifurcations in Chua‘s bircuit, J. Circuits Syst. Comput., 3, 2 (1993) 441–470.
  • [10] B. Sobota, Parallel hierarchical model of visualization computing, J. Inf. Control Manag. Syst., 5 (2007) 345–350.
  • [11] B. Sobota, J. Perháč, Cs. Szabó, An application of parallel, distributed and network computer systems to solve computational processes in an area of large graphical data volumes processing, Comput. Sci. Technol. Res. Survey, KPI FEI TU Košice, 2008, 3, pp. 37–42.
  • [12] V. Špány, Personal communication, Department of Electronics and Multimedia Communications, Faculty of Electrical Engineering and Informatics, TU Košice, January 2006.
  • [13] V. Špány, P. Galajda, M. Guzan, Boundary surfaces of one-port memories, in: Tesla III Milenium: Proc. 5th Int. Conf., Beograd, 1996, pp. IV. 130–IV. 137.
  • [14] V. Špány, P. Galajda, M. Guzan, The state space mystery in multiple-valued logic circuit with load plane – part I., Acta Electrotechnica Inform., 1, 1 (2001) 17–22.
  • [15] V. Špány, L. Pivka, Boundary surfaces in sequential circuits, Internat. J. Circuit Theory Appl., 18, 4 (1990) 349–360.