Lagrangian for Circuits with Higher-Order Elements

The necessary and sufficient conditions of the validity of Hamilton’s variational principle for circuits consisting of (α,β) elements from Chua’s periodical table are derived. It is shown that the principle holds if and only if all the circuit elements lie on the so-called Σ-diagonal with a constant sum of the indices α and β. In this case, the Lagrangian is the sum of the state functions of the elements of the L or +R types minus the sum of the state functions of the elements of the C or −R types. The equations of motion generated by this Lagrangian are always of even-order. If all the elements are linear, the equations of motion contain only even-order derivatives of the independent variable. Conclusions are illustrated on an example of the synthesis of the Pais–Uhlenbeck oscillator via the elements from Chua’s table.


Introduction
It is well known that the Lagrange formalism adopted from classical mechanics [1] can be used for describing phenomena in electric circuits. The equations of motion of a conservative system can be derived from the Lagrange function L = L x 1 , .., x n , .
x n , t according to the Euler-Lagrange equation Here, x i , for i = 1,..., n, are the generalized coordinates of the system. The dot above x i denotes the derivative with respect to time, or, in other words, the corresponding component of the generalized velocity, and n∈ Z+ is the number of degrees of freedom. The Lagrange function is usually derived as a difference between the kinetic and the potential energy of the system, or the energy of inductors and capacitors in electrical engineering. In the latter case, the electric charges or integrals of the electric voltages, the fluxes, usually supersede the role of the generalized coordinates. The coordinates are selected depending on the method of the circuit analysis. The generalized velocities are the derivatives of the coordinates with respect to time, thus electric currents and voltages. Systematic differentiation of the Lagrange function (1) with respect to coordinates, velocities, and time, yields the equations of motion (2) expressing the equilibrium of the generalized forces. For electrical engineering, these equations represent Kirchhoff's voltage or current laws (the KVL or KCL).
The so-called Hamilton variational principle holds for the conservative systems. Let the system go through the states x i (t 1 ) and x i (t 2 ), i = 1,..., n, at instants t 1 and t 2 , t 2 > t 1 . The principle states that, within the time interval (t 1 ,t 2 ), the system "selects" for its motion, between the initial and the final state, such a trajectory that is a stationary point of the action functional x 1 , .., .
x n , t dt The functional (3) is a definite integral, which generates a number for each trajectory. From all the conceivable trajectories that might pass through the terminal states at instants t 1 and t 2 , it is only for the existing trajectory that the action (3) takes the stationary value (the minimum, maximum, or saddle point). That is why Hamilton's principle is denoted as the principle of stationary action. Any virtual variation (δx 1 ,..., δx n ) of the existing trajectory must give rise to a zero variation of the first order δA of the action. The existing trajectory is, therefore, the solution of the optimization task It can be easily proven that the conditions (2) and (4) are equivalent [1]. If the system fulfils the Hamilton's variational principle (4), then its equations of motion are in the form of the Euler-Lagrange Equation (2) and vice versa.
External excitation can be included in the modified Lagrange function L' as where Q 1 ,..., Q n are components of the generalized external force. Then Hamilton's variational principle also holds for systems represented by this modified Lagrange function. The presence of arbitrary dissipative forces in the system, dependent on generalized velocities, means a violation of the necessary condition for the validity of Hamilton's variational principle. Additional terms, incompatible with the desired form (2), appear in the equations of motion. The dissipative forces are usually derived as gradients of the dissipative function D, which depends on generalized velocities. The equations of motion are then in the form It is evident from (6) that not all the generalized forces can be generated from a single scalar function L as Hamilton's principle states. The dissipative forces must be generated via another scalar function, the dissipative function D.
The utilization of classical Lagrange formalism leads to equations of motion (2), which are of the 2 nd order. It is a logical consequence of the fact that the equations of motion of classical mechanics represent an equilibrium of inertial, potential, and dissipative forces. However, some existing basic higher-order processes are not governed by this simple rule. The interest in studying systems whose dynamics explicitly depends on the accelerations and higher-order derivatives of generalized coordinates with respect to time gave rise to Ostrogradsky's work [2] (1850), which generalizes the Lagrange and Hamilton formalism in this sense. Other works such as [3] by Whittaker (1904) followed [2], introducing the Lagrange function in the form L = L x 1 , .., x n , x where m is the maximum order of the derivative of a generalized coordinate with respect to time. The equations of motion are as follows: The interest in higher-order Lagrangians of type (7) is currently growing due to the dynamics of higher-order systems in various branches of theoretical physics, control theory, and applied mathematics being intensively studied. A brief survey is given in [4].
A similar development has also been noted in the theory of electrical circuits, stigmatized to some extent by the idea of the exclusivity of three fundamental elements, namely the inductor, capacitor, and resistor. A plethora of novel two-terminal devices with startling parameters and intriguing application potentials, such as the ovonic threshold switch, Josephson junction, or various diodes (Esaki, Gunn, IMPact ionization Avalanche Transit-Time (IMPATT)), appeared in the 1960s. Researchers in the field of nonlinear electronics were then faced with phenomena that were difficult to categorize within the current circuit theory, or they were classified as "anomalies", "exotic", or "paradox" [5]. It was a period when the intellectual basis of a new scheme began to be formed that would put all the above phenomena in order.
The axiomatic basis of the nonlinear circuit theory was laid down in 1969 by Leon Chua in the book [6]. The principles in [6] contributed to the prediction of a memristor as the fourth fundamental electrical element in 1971 [7]. In the concept of the so-called periodical table of fundamental elements from 1980 [8], the resistor, capacitor, inductor, and memristor are just special cases of more general (α,β) elements, also denoted as the HOEs (Higher-Order Elements). They preserve an unambiguous relationship between the variables v (α) and i (β) , where the indices α,β are the orders of derivatives/integrals (for positive/negative indices) of the terminal voltages v and currents i with respect to time. Until then, elements occupied the following positions in the new schematics: (0,0)-resistor (R), (0,−1)-capacitor (C), (−1,0)-inductor (L), and (−1,−1)-memristor (MR). After discovering the nanodevice with a memristive behavior in the Hewlett-Packard laboratories in 2008 [9], Chua appealed to search for the other missing elements from the table [10], namely the memcapacitor (MC) (−1,−2) and the meminductor (ML) (−2,−1). The electrical elements from the family of the so-called Frequency-Dependent Negative Resistors (FDNRs), known since the 1960s and gradually synthetized as active elements with the help of Bruton's work [11], behave as elements from Chua's table. The work [12] distinguishes between the FDNC (Frequency-Dependent Negative Conductor) and the FDNR (Frequency-Dependent Negative Resistor). Such resistors with negative resistances, providing quadratic dependence of the small-signal conductance or resistance on the frequency, are the (1,−1) or (−1,1) elements from Chua's table. In [13], Soliman introduces the FDPC (Frequency-Dependent Positive Conductance), whose conductance depends on the fourth power of frequency. In fact, it is the (2,−2) element. The mechanical inerter from 2002 [14] was recognized as the (1,0) element [15].
Some recent papers suggest modifications of the original Chua's table of fundamental elements [16][17][18]. The so-called storeyed structure of the HOEs [16] just transforms Chua's table into a different graphic form [16]. The storeyed structure and Chua's table, therefore, contain identical sets of the HOEs. Wang's triangular periodic table [17] contains only HOEs that occupy only three so-called ∆-diagonals of the original Chua's table (see Section 2) for ∆ = 1, 0, and 1. Since no other HOEs, the negative resistors being among them, are missing in Wang's table, this table cannot be utilized for our purposes. It is well documented via the impossibility of drawing a Lagrangian for circuits with dissipative elements [19].
The reason why Chua's and Wang's table contain different sets of elements consists in different conceptions of the fundamental element. While Chua's definition starts from an unambiguous axiomatic approach [8], Wang introduces additional limiting conditions, which must be fulfilled in order to denote the element as fundamental. These limitations are of postulating nature, namely: The fundamental electric quantities are the flux and the charge. A fundamental circuit element should link two electric quantities, at least one of which should be fundamental [17]. Such a tapered definition of the fundamental element results in the absence of quite a few of the HOEs in Wang's table. However, these missing elements are unsubstitutable in modeling the existing dynamic systems [20].
The work [18], starting from the intuitive approach, re-defines the fundamental element even more. For example, according to [18], the fundamental elements must be linear. Such a premise must inevitably exclude the memristor, whose constitutive relation is nonlinear, from the set of fundamental elements. Since the corresponding table of elements suggested in [18] contains only a subset of linear elements, it cannot be used for studying the variational principles in dynamic systems containing the HOEs.
For the above reasons, this work will, therefore, be based on Chua's approach. The Chua's table of fundamental elements became an important tool of the so-called predictive modeling [20]. In spite of some limits of this approach [21], it can be advantageous in the analysis and also synthesis, provided that the system be built exclusively from HOEs as the basic building blocks. One limit is due to the fact that only two-terminal HOEs appear in Chua's table. In order to model a certain class of systems, the above idea of two-terminal fundamental elements should, therefore, be generalized to multipoles or multiports [22]. The importance of predictive modeling is increasing when studying systems with higher-order dynamics. It can be documented in researching chaotic phenomena, where the higher-order dynamics is a necessary condition. According to [23], 427 important papers about the analysis or synthesis of chaotic circuits via memristors, memcapacitors, and meminductors were published within one decade after discovering the memristor in Hewlett-Packard labs in 2008 [9], 283 from them within the last three years of that decade.
In spite of such a precipitate development, two basic questions have still not been answered: Can the Lagrangian (7) of a system compounded from HOEs be constructed and under what conditions, and does Hamilton's variational principle hold for such systems. The first responses to these questions were the studies [24][25][26], which introduced Lagrange and Hamilton formalisms for circuits with memristors and memcapacitors. The follow-up paper [27] defines the state functions of general HOEs and derives the Euler-Lagrange equations from them. It also states in [27] that "The paper does not deal with the question of whether or not the extremal principle connected with the original idea of the Lagrangian is fulfilled." This paper aims to answer the above questions. Section 2 summarizes the current classification of the fundamental elements according to their positions on the ∆-diagonals of the table and the definitions of their state functions according to [27]. Recalled therein are some rules of the taxicab geometry, which hold for the table of elements, recently published in [28]. The following key Section 3 describes the derivation of Hamilton's variational principle for circuits containing elements from an arbitrary Σ-diagonal of the table. Section 4 illustrates that the principle does not hold for circuits containing elements from various Σ-diagonals. Section 5 transforms the Lagrangian of the circuit with HOEs into a form that is compatible with the hitherto used forms of the classical potential functions. In Section 6, the results are applied to the analysis and synthesis of the Pais-Uhlenbeck oscillator via HOEs.

Table of Higher-Order Elements
The basic attributes of Chua's table of HOEs are well known from the original paper [8]. The elements are drawn in Figure 1 as circles with the integer-type (α,β) coordinates in the rectangular grid. The (α,β) coordinates are not the only way of defining the location of the element in the table. There may be other coordinates (Σ,Δ), where Σ = α + β and Δ = β − α. All elements with the same value Σ/Δ lie on one of the Σ/Δ-diagonals; see Figure 1. The (Σ,Δ) element is located at the intersection of both diagonals.
It is well known that the Δ-diagonals determine the small-signal behavior of the elements located on these diagonals [8]. Consider the voltage representation of the element according to (9), and analyze the small-signal deflections of voltages and currents dv and di around the operating point Q lying on the constitutive relation. It is obvious that where mQ is the slope of the constitutive relation at point Q. The Fourier transform of Equation (10) and a simple rearrangement yield the formula for the small-signal impedance of the element It follows from (11) that elements lying on a given Δ-diagonal have the same small-signal characters. For Δ repeatedly increasing by one, the element type will change with a period of 4 among the positive resistor + R, inductor L, negative resistor − R, and the capacitor C. For all these cases, the elements are generally frequency-dependent according to (11).
Considering small-signal modeling, it follows from Figure 1 that a given Σ-diagonal can be occupied either by resistors (regularly alternating positive and negative resistors for even values of Σ), or by reactive elements (regularly alternating inductors and capacitors for odd values of Σ). It is also known that the number Σ determines the physical dimension of the corresponding state function [27]. The dissipative function with the dimension of power is a state function of resistors from the diagonal Σ = 0. The state functions of the inductors and capacitors from the diagonal Σ = −1 are energies, and the state function of the memristor (Σ = −2) is called action [7], i.e., the time integral of the energy. All these functions are defined in the same manner, thus as integrals of the constitutive relation of the element with respect to the corresponding constitutive variable. It means that they can be illustrated via the areas below the graphs of given constitutive relations; see Figure 2. The pair of state functions originating by the integration with respect to one or the other constitutive variable have been denoted in physics and electrical engineering as a function and co-function. In 1951, Cherry Every element is one-port, which predicts an unambiguous relation between its port variables v (α) and i (β) , i.e., it preserves the constitutive relation F(v (α) ,i (β) ) = 0, where F() is generally a nonlinear function, in all circumstances and regardless of the behavior of the surrounding network. The elements with explicit forms of the constitutive relation.
are called current-or voltage-controlled elements. The (α,β) coordinates are not the only way of defining the location of the element in the table. There may be other coordinates (Σ,∆), where Σ = α + β and ∆ = β − α. All elements with the same value Σ/∆ lie on one of the Σ/∆-diagonals; see Figure 1. The (Σ,∆) element is located at the intersection of both diagonals.
It is well known that the ∆-diagonals determine the small-signal behavior of the elements located on these diagonals [8]. Consider the voltage representation of the element according to (9), and analyze the small-signal deflections of voltages and currents dv and di around the operating point Q lying on the constitutive relation. It is obvious that where m Q is the slope of the constitutive relation at point Q. The Fourier transform of Equation (10) and a simple rearrangement yield the formula for the small-signal impedance of the element It follows from (11) that elements lying on a given ∆-diagonal have the same small-signal characters. For ∆ repeatedly increasing by one, the element type will change with a period of 4 among the positive resistor + R, inductor L, negative resistor − R, and the capacitor C. For all these cases, the elements are generally frequency-dependent according to (11).
Considering small-signal modeling, it follows from Figure 1 that a given Σ-diagonal can be occupied either by resistors (regularly alternating positive and negative resistors for even values It is also known that the number Σ determines the physical dimension of the corresponding state function [27]. The dissipative function with the dimension of power is a state function of resistors from the diagonal Σ = 0. The state functions of the inductors and capacitors from the diagonal Σ = −1 are energies, and the state function of the memristor (Σ = −2) is called action [7], i.e., the time integral of the energy. All these functions are defined in the same manner, thus as integrals of the constitutive relation of the element with respect to the corresponding constitutive variable. It means that they can be illustrated via the areas below the graphs of given constitutive relations; see Figure 2. The pair of state functions originating by the integration with respect to one or the other constitutive variable have been denoted in physics and electrical engineering as a function and co-function. In 1951, Cherry and Miller introduced the notions of energy and coenergy [29] for the reactive elements and the power functions content and cocontent [30] for resistive elements. The functions and co-functions are assigned to individual representations of the state functions of the four fundamental electrical elements, according to Table 1 [31]. and Miller introduced the notions of energy and coenergy [29] for the reactive elements and the power functions content and cocontent [30] for resistive elements. The functions and co-functions are assigned to individual representations of the state functions of the four fundamental electrical elements, according to Table 1 [31]. The generalization of state functions and cofunctions to arbitrary (α,β) elements is introduced in [27]. The state functions are integrals of the corresponding constitutive relations for the current or voltage representation of the element. The mechanism for classifying the representation (12) into functions and cofunctions is obvious from Table 1.

Type of the Element
It follows from the definition formulae (12) that the physical dimension of the state function of The state functions of all the elements located on the common Σ-diagonal have, therefore, the same physical dimension. It is a necessary condition for their additivity.
Let the symbols T * and V or + D and -D * denote the sums of state functions T * α,β and Vα,β of type L and type C element, or the functions + Dα,β and − D * α,β of type + R and -R elements from one Σ-diagonal,

Type of the Element
The generalization of state functions and cofunctions to arbitrary (α,β) elements is introduced in [27]. The state functions are integrals of the corresponding constitutive relations for the current or voltage representation of the element. The mechanism for classifying the representation (12) into functions and cofunctions is obvious from Table 1.
It follows from the definition formulae (12) that the physical dimension of the state function of a general (α,β) element is [Volt·Amper·sec −Σ ] or [Joule·sec -(Σ+1) ]. The state functions of all the elements located on the common Σ-diagonal have, therefore, the same physical dimension. It is a necessary condition for their additivity.
Let the symbols T * and V or + D and -D * denote the sums of state functions T * α,β and V α,β of type L and type C element, or the functions + D α,β and − D * α,β of type + R and -R elements from one Σ-diagonal, i.e., The summations are done for all circuit elements of L and C or + R and -R types from a given Σ-diagonal. Seeing that we work with generalized coordinates, which can be of various physical dimensions, the quantities (13a) and (13b) can be regarded as generalized energies and contents, respectively.
The regular alternation of potential functions and cofunctions within the even ( + D and -D * for elements of + R and -R types) and odd (T * and V for elements of L and C types) Σ-diagonals is demonstrated in Figure 3. The summations are done for all circuit elements of L and C or + R and -R types from a given Σdiagonal. Seeing that we work with generalized coordinates, which can be of various physical dimensions, the quantities (13a) and (13b) can be regarded as generalized energies and contents, respectively.
The regular alternation of potential functions and cofunctions within the even ( + D and -D * for elements of + R and -R types) and odd (T * and V for elements of L and C types) Σ-diagonals is demonstrated in Figure 3.
For the current representation of state function (12), the voltage variable v (α) = dSα,β/di (β) is equal to the slope of the state function of the element at the operating point. Considering only the elements within one given Σ-diagonal, each of them will have a unique index β. The voltage quantity v (α) of a specific element will be equal to the partial derivative of the corresponding total state functions (13a) or (13b) of the circuit with respect to a specific current variable i (β) .
The rules of the Taxicab geometry hold in the table of fundamental elements. According to these rules, the distance between the elements is given by the shortest path between them in the frame of the rectangular net. In such a type of geometry, the circle drawn around a central point has the form of a square whose diagonals occupy the horizontal and vertical positions. According to [28], the order  of the differential equation describing the behavior of one-port consisting of serial and parallel combinations of HOEs is equal to the radius of the smallest quarter-circle that incepts all the elements of the one-port. The center of the quarter-circle is occupied by the so-called hidden element with the coordinates (αMAX,βMIN) or (αMIN,βMAX) for a series or a parallel connection of the elements. The subscripts MAX and MIN denote the maximum and minimum coordinates that are occupied by the elements in the table. An example of constructing the quarter-circles for series and parallel connections of HOEs is shown in Figure 4.  For the current representation of state function (12), the voltage variable v (α) = dS α,β /di (β) is equal to the slope of the state function of the element at the operating point. Considering only the elements within one given Σ-diagonal, each of them will have a unique index β. The voltage quantity v (α) of a specific element will be equal to the partial derivative of the corresponding total state functions (13a) or (13b) of the circuit with respect to a specific current variable i (β) .
The rules of the Taxicab geometry hold in the table of fundamental elements. According to these rules, the distance between the elements is given by the shortest path between them in the frame of the rectangular net. In such a type of geometry, the circle drawn around a central point has the form of a square whose diagonals occupy the horizontal and vertical positions. According to [28], the order O of the differential equation describing the behavior of one-port consisting of serial and parallel combinations of HOEs is equal to the radius of the smallest quarter-circle that incepts all the elements of the one-port. The center of the quarter-circle is occupied by the so-called hidden element with the coordinates (α MAX ,β MIN ) or (α MIN ,β MAX ) for a series or a parallel connection of the elements. The subscripts MAX and MIN denote the maximum and minimum coordinates that are occupied by the elements in the table. An example of constructing the quarter-circles for series and parallel connections of HOEs is shown in Figure 4.
where the index j = 0 corresponds to the hidden element, v is the terminal voltage, m is the number of elements increased by the hidden element if appropriate, fi() is the constitutive relation (9) of the i th element, i = 0, ..., m. Note that, according to Figure 4, the order βj of the derivative of current in differential equation (14) will take values from βMIN to (βMIN + ). Thus, it is natural to introduce a new variable x = i (βMIN) , and the differential equation drawn for such a generalized current will contain its derivatives of orders 0 to , which will correspond to the distances of the elements from the hidden element in vertical direction. Similarly, let us introduce the generalized voltage u = v (αMAX) . According to (14), the orders of derivatives of generalized voltages fj will also be between 0 and  since they will correspond to the distances of the elements from the hidden element in the horizontal direction.
The symbols x and u can also be universally used for generalized voltages and currents of elements connected in parallel. The situation for a series and parallel connection is summarized in Table 2. For a parallel connection, where KC (β) L is applied, the driving variable is u = i (βMAX) flowing in the common node. The individual terms of differential equations will signify the generalized currents uj = ij (βMAX) flowing through individual elements.

Law
Hidden hidden element (-7,-3) where the index j = 0 corresponds to the hidden element, v is the terminal voltage, m is the number of elements increased by the hidden element if appropriate, f i () is the constitutive relation (9) of the i th element, i = 0, ..., m. Note that, according to Figure 4, the order β j of the derivative of current in differential equation (14) will take values from β MIN to (β MIN + O). Thus, it is natural to introduce a new variable x = i (βMIN) , and the differential equation drawn for such a generalized current will contain its derivatives of orders 0 to O, which will correspond to the distances of the elements from the hidden element in vertical direction. Similarly, let us introduce the generalized voltage u = v (αMAX) . According to (14), the orders of derivatives of generalized voltages f j will also be between 0 and O since they will correspond to the distances of the elements from the hidden element in the horizontal direction.
The symbols x and u can also be universally used for generalized voltages and currents of elements connected in parallel. The situation for a series and parallel connection is summarized in Table 2. For a parallel connection, where KC (β) L is applied, the driving variable is u = i (βMAX) flowing in the common node. The individual terms of differential equations will signify the generalized currents u j = i j (βMAX) flowing through individual elements.

Law
Hidden Element The uniform notation enables choosing both generalized voltages and currents as the generalized coordinates x i in the Lagrangian (7) and in the corresponding equations of motion (8).
The circuits with HOEs, where Hamilton's principle holds are governed by the set of differential equations (8). Consider that each equation represents KV (α) L or KC (β) L for some loop or node. Then the variables x i , i = 1,..., n, will be loop variables of i (β) type or nodal variables of v (α) type, respectively.
The Lagrange formalism (7) and (8) is based on the generalized coordinates x i . On the other hand, predictive modeling starts from different variables, namely from generalized voltages and currents of individual elements joined by the constitutive relations (9). In order to apply the Lagrange formalism to circuits with HOEs, the transform between these two sets of variables must be used. The following procedure is a generalization of the method described in [19]. Consider the above transform as a linear combination of the variables x i where x ε is the generalized voltage or current of a specific element ε, ε a i are the coefficients of the linear combination for the element ε. For example, if the element coincides with the loops No. 1 and 3, the corresponding coefficients ε a 1 and ε a 3 will be +1 or −1 according to the agreement or disagreement of the reference directions of the element and of the loop. If the element does not coincide with the loop, the corresponding coefficient is zero. A similar consideration applies to the nodal variables-the element ε directed between nodes 2 and 4 provides the coefficients ε a 2 = +1 and ε a 4 = −1, and all the other coefficients are zero.
The symbol ε will be specified as follows: all the elements from the circuit will be classified according to their positions in Chua's table (more elements of the same type will be counted as one element) and numbered from 0 to m. Zero will be assigned to the hidden element. Each element from the circuit will belong to one type from the set ε k , k = 0, ..., m.

Sufficient Condition for Hamilton's Principle in Circuits with HOEs
In this Section, the following proposition will be proved: Proposition 1. If the circuit is composed exclusively of the elements from Chua's table that are located on the common Σ-diagonal, then it conforms to Hamilton's principle.
The proof will be given for the KV (α) L representation of the circuit, according to Table 2. The proof for the dual KC (β) L representation is analogous.
Consider a loop in the circuit containing the (α,β) elements that are located on a general Σ-diagonal in Figure 5. The elements are numbered in ascending order from 0 to m, starting with the hidden element (α MAX ,β MIN ). The constitutive relations and the state functions of the elements of εk type, k = 0,1,2,…, m, are If the notation for generalized voltages and currents for KV (α) L from Table 2 is taken into consideration, then Equation (17) follows from the fact that the distance of the k th element from the hidden element along the β coordinate is equal to the number k (Considering the position of the element on the diagonal, the same holds for the α coordinate.) The formal variable x in (17) can be substituted by the generalized current of a specific element ε, which is of ε k type with the constitutive relation f k ().
Let us analyze the impact of the trajectory variation δxi , i = 1,..., n, on the state functions Sk. It follows from the definition (16) of the state function that For a particular element ε, the variation δx (k) will depend on which loops the element will be involved in. Considering Equation (15), it holds The index i of the term i fk indicates that the constitutive relation is completed by a correct sign according to the reference direction of the element within the i th loop, or it makes this relation zero if the element is not present in the loop.
The variations of the coordinates, not their derivatives, appear in the variational task (4). In the next step, the variations δxi (k) will, therefore, be transferred to δxi. This well-known procedure via repetitive integration of (19) by parts yields the result The constitutive relations and the state functions of the elements of ε k type, If the notation for generalized voltages and currents for KV (α) L from Table 2 is taken into consideration, then Equation (17) follows from the fact that the distance of the k th element from the hidden element along the β coordinate is equal to the number k (Considering the position of the element on the diagonal, the same holds for the α coordinate.) The formal variable x in (17) can be substituted by the generalized current of a specific element ε, which is of ε k type with the constitutive relation f k ().
Let us analyze the impact of the trajectory variation δx i , i = 1,..., n, on the state functions S k . It follows from the definition (16) of the state function that For a particular element ε, the variation δx (k) will depend on which loops the element will be involved in. Considering Equation (15), it holds The index i of the term i f k indicates that the constitutive relation is completed by a correct sign according to the reference direction of the element within the i th loop, or it makes this relation zero if the element is not present in the loop.
The variations of the coordinates, not their derivatives, appear in the variational task (4). In the next step, the variations δx i (k) will, therefore, be transferred to δx i . This well-known procedure via repetitive integration of (19) by parts yields the result If the boundary points of the original trajectory and the first derivatives of the (k−1) th orders are preserved, it follows from (15) that the first right-side term of (20) must be zero. Summing the terms in (20) separately for odd and for even indices k and the calculating their difference will result in An arbitrary number of elements with the constitutive functions f k and state functions S k can be present in the circuit. The summation of Equation (21) for all the elements ε in the circuits leads to the formula The sum with the summation index k in (22) represents KV (α) L for the i th loop, or one of the Euler-Lagrange equations (8). The formula (22) is therefore equal to zero, and the left-side integrand is the Lagrange function of the entire circuit. Then, if the boundary conditions are fulfilled for every i = 1, . . . , n, then the existing trajectory is an extremal of the action potential, or The following differential equation holds for a series circuit with one loop in which just one element of each type occurs: Note that only the even-order derivatives of the independent variable are present in the differential equation of the linear circuit in which Hamilton's principle applies. The reason is that if the constitutive relations f k from (22) are linear, then it holds where m k is the slope of the constitutive relation at the operating point.
The boundary conditions (23) indicate that the trajectory variation is ruled by stricter conditions than for the classical variational principle. All virtual trajectories must copy the first (m−1) orders of the curvatures of the original trajectory at both boundary points. Therefore, the variational principle of systems containing HOEs retains fewer degrees of freedom of the variation of the trajectory, as indicated in Figure 6.

Necessary Condition for Hamilton's Principle in Circuits with HOEs
Let us prove the following proposition:

Proposition 2. If a circuit, made up exclusively from the elements from Chua's table, has to comply with
Hamilton's principle, then all its elements must occur on a common Σ-diagonal.
Consider the contrary case: one element, εadd, with the constitutive relation fadd() and state function Sadd, is located beyond the diagonal, as shown in Figure 7, case (a) or (b). Let its distance from the hidden element be dα + dβ, where the subscripts denote the components of this distance along the α and the β coordinates. The distance of the element from the Σ-diagonal is with the difference dα − dβ being positive/negative below/above the diagonal, respectively. The constitutive relation fadd() introduces the coupling condition

Necessary Condition for Hamilton's Principle in Circuits with HOEs
Let us prove the following proposition: Proposition 2. If a circuit, made up exclusively from the elements from Chua's table, has to comply with Hamilton's principle, then all its elements must occur on a common Σ-diagonal.
Consider the contrary case: one element, ε add , with the constitutive relation f add () and state function S add , is located beyond the diagonal, as shown in Figure 7, case (a) or (b). Let its distance from the hidden element be d α + d β , where the subscripts denote the components of this distance along the α and the β coordinates. The distance of the element from the Σ-diagonal is with the difference d α − d β being positive/negative below/above the diagonal, respectively. The constitutive relation f add () introduces the coupling condition KV (α) L for the i th loop will be in the form However, the variation can be identically equal to zero only for dα = dβ, i.e., for a zero distance of the element from the -diagonal. Due to its time-domain differentiation or integration, the state function Sadd loses the ability to be a scalar potential function, which can generate the monogenic quantity fadd [1]. The generalized voltages fadd (d) and fadd (-d) , by which the terms below and above the diagonal contribute by (32) to KV (α) L, are polygenic quantities [1] that cannot be included in the Lagrangian. That is why Hamilton's principle holds only for d = 0, i.e., when all the elements of the circuit are located on the same -diagonal.
The proof holds for an arbitrary position of the hidden elements on the -diagonal and for an arbitrary number of elements. These parameters can be modified in order to achieve that the point, lying off the diagonal, will be situated within the characteristic quarter-circle, as illustrated in Figure 7.

Compatibility with the Classical Variational Principle
The variation of the action is then However, the variation can be identically equal to zero only for d α = d β , i.e., for a zero distance of the element from the Σ-diagonal. Due to its time-domain differentiation or integration, the state function S add loses the ability to be a scalar potential function, which can generate the monogenic quantity f add [1]. The generalized voltages f add (d) and f add (-d) , by which the terms below and above the Σ-diagonal contribute by (32) to KV (α) L, are polygenic quantities [1] that cannot be included in the Lagrangian. That is why Hamilton's principle holds only for d = 0, i.e., when all the elements of the circuit are located on the same Σ-diagonal.
The proof holds for an arbitrary position of the hidden elements on the Σ-diagonal and for an arbitrary number of elements. These parameters can be modified in order to achieve that the point, lying off the diagonal, will be situated within the characteristic quarter-circle, as illustrated in Figure 7.
with the difference dα − dβ being positive/negative below/above the diagonal, respectively. The constitutive relation fadd() introduces the coupling condition Figure 7. The loop is formed via interconnecting m + 1 elements located on a common Σ-diagonal (green circles) and one element (red circle) (a) below or (b) above the diagonal. The elements numbered from 0 to m are characterized by constitutive relations f 0 to f m . The added element beyond the diagonal has the constitutive relation f add . The corresponding state functions are S 0 to S m and S add .

Compatibility with the Classical Variational Principle
Basically, the Lagrangian (25) of the variational principle (24) can be of two different forms according to whether the circuit elements are located on the odd or the even Σ-diagonal of the table (see Figure 3).
The L and C elements occupy the odd Σ-diagonals, and the state functions S 0 to S m from (25) are the potential functions (13a). Their definition formulae together with (25) provide the resulting form of the Lagrangian for a circuit built only from elements of the odd Σ-diagonal where T* and V are the potential cofunction and function, representing parts of the circuit built exclusively from the elements of L and C types. The Lagrangian can, therefore, be obtained by subtracting the sum of the generalized energies of all capacitors from the sum of the generalized coenergies of all inductors in the circuit. The even Σ-diagonals contain the elements of + R and − R types, and the state functions S 0 to S m from (25) are the potential functions (13b), + D for positive, and − D * for negative resistive elements. Their definition formulae together with (25) give the resulting Lagrangian for a circuit built only from the elements from even Σ-diagonal where + D and − D* are the potential function and cofunction, representing parts of the circuit built exclusively from the elements of + R and − R types. The Lagrangian can, therefore, be obtained by subtracting the sum of generalized dissipative cofunctions of all negative resistors from the sum of generalized dissipative functions of all positive resistors in the circuit. Formulaes (36) and (37), defining the Lagrangian, can be modified, considering that the Euler-Lagrange equations generated by the Lagrangian do not depend on its sign. The sign of the Lagrangian only exerts influence on the type of the extreme of the action potential, which is reached during the movement along the existing trajectory.
Equations (36) and (37) hold for the current representation of the circuit. The voltage representation works with dual quantities, which means that the functions are replaced by cofunctions and vice versa [27].

Application: Pais-Uhlenbeck Oscillator
The necessary and sufficient condition for the validity of Hamilton's principle, namely that all the circuit elements must occupy one common Σ-diagonal, strictly holds only for nonlinear elements, i.e., elements with nonlinear constitutive relations. It is well known that the character of the linear element is not changed during the element movement along the ∆-diagonal [8,32]. Such a movement changes the α and β indices by the same value k, thus both the original (F(v (α) ,i (β) ) = 0) and the new (F(v (α+k) ,i (β+k) ) = 0) relations hold simultaneously. The linear (α,β) element can therefore be considered as an arbitrary (α + k,β + k) element for an arbitrary integer k. Taking also into account the linear elements, one has a chance to escape from the strict rule of one Σ-diagonal. The procedure is illustrated on the example of the synthesis of the Pais-Uhlenbeck (PU) oscillator via higher-order elements. The PU oscillator [33] is a toy model of recent higher-derivative theories. Its Lagrangian and differential equation are adopted from [34].
The PU oscillator is described by the second-order Lagrangian x, ..
and the resulting trajectory x(t) is the extremal of the action functional, which is based on the Lagrangian (38). The Euler-Lagrange equation, corresponding to this Lagrangian, is Equation (39) is a linear differential equation with constant parameters with the solution where and C + , C -, S + , and Sare constants given by the initial conditions. The system, therefore, oscillates at two different frequencies. For = 0, it changes to a classical harmonic oscillator with the oscillation frequency ω.
Since the oscillator fulfills Hamilton's principle, it should be synthetized via higher-order elements. Comparing the differential equation (39) and the relation (26), the synthesis can be performed via three elements with the constitutive relations x, f 2 ..
All three elements will neighbor on an arbitrary common Σ-diagonal. The generalized coordinate x (current variable for the series connection) and the physical size of generalized forces (voltage variables of individual in-loop elements) will be given by the specific positioning of the hidden element in Chua's table (see Figure 8).
All three elements will neighbor on an arbitrary common Σ-diagonal. The generalized coordinate x (current variable for the series connection) and the physical size of generalized forces (voltage variables of individual in-loop elements) will be given by the specific positioning of the hidden element in Chua's table (see Figure 8).  The given characteristic quarter-circles for series and parallel connections confirm that the order of the circuits is always four [28].
Qa is the quarter-circle for the series connection of the inductor L, capacitor C, and the (1,−2) element, whose small-signal model represents the inductance decreasing with the fourth power of frequency. In this case, the variable x is the double time-domain integral of current or the integral of charge σ = qdt. For nonlinear elements, i.e., elements whose constitutive relations f0(), f1(), and f2() are nonlinear functions, the series circuit in Figure 9a models the nonlinear PU oscillator, governed by the nonlinear differential equation  Figure 8 illustrates the characteristic quarter-circles Q a , Q b , and Q c , which represent three different implementations of the oscillator in Figure 9a-c. The given characteristic quarter-circles for series and parallel connections confirm that the order of the circuits is always four [28].
Q a is the quarter-circle for the series connection of the inductor L, capacitor C, and the (1,−2) element, whose small-signal model represents the inductance decreasing with the fourth power of frequency. In this case, the variable x is the double time-domain integral of current or the integral of charge σ = qdt. For nonlinear elements, i.e., elements whose constitutive relations f 0 (), f 1 (), and f 2 () are nonlinear functions, the series circuit in Figure 9a models the nonlinear PU oscillator, governed by the nonlinear differential equation The Lagrangian is a nonlinear form of the Lagrangian (37) as a sum of the coenergies of the inductive elements L and the (1,−2) minus the energy of the capacitor: .. Figure 8 is the quarter-circle for the series connection of the (−4,−7) element, which is concurrently the hidden element, and two other nearest neighbors of C and L types on the corresponding Σ-diagonal. If these two elements are linear, their character does not change during their movement along the ∆-diagonal. That is why the linear capacitor C and the linear inductor L can be used according to Figure 9b. In this case, the quantity x is the sevenfold time-domain integral of current or the quantity i (−7) (t). For a nonlinear (−4,−7) element, the Lagrangian is in the form L x, . x, ..
.. where the right-side terms of (45) are the coenergy of the (−4,−7) element, energy of the capacitor, and coenergy of the inductor, respectively.
If the differential equation (43) represents the equilibrium of the generalized velocities and not the equilibrium of generalized forces as we hitherto supposed, it can be considered the generalized Kirchhoff law for the current variables i (β) , and this equation can be modeled by a parallel connection of (α,β) elements. Figure 8 shows the quarter-circle Q c for a parallel connection of the resistor R and two elements, F 1 and F 2 , namely the frequency-dependent negative resistors. F 1 is the classical FDNC [12], whose small-signal resistance decreases with the square of the frequency. F 2 is the FDNR with opposite dependence, i.e., its resistance increases with the square of the frequency [12]. Comparing these effects with the frequency-dependent terms of the differential equation (39) yields the explanation of why the given circuitry can generate the solution of the equation of just this type. Removing the F 2 element, which represents the first term of the differential equation (39) or (43), leads to the classical harmonic oscillator for = 0.
If all three elements are nonlinear, then the Lagrangian has a general form (44). The individual right-side terms of (44) mean the content of the FDNR, cocontent of the resistor, and content of the FDNC.
The Lagrangian is a nonlinear form of the Lagrangian (37) as a sum of the coenergies of the inductive elements L and the (1,−2) minus the energy of the capacitor: Qb in Figure 8 is the quarter-circle for the series connection of the (−4,−7) element, which is concurrently the hidden element, and two other nearest neighbors of C and L types on the corresponding Σ-diagonal. If these two elements are linear, their character does not change during their movement along the Δ-diagonal. That is why the linear capacitor C and the linear inductor L can be used according to Figure 9b. In this case, the quantity x is the sevenfold time-domain integral of current or the quantity i (−7) (t). For a nonlinear (−4,−7) element, the Lagrangian is in the form ( ) ( ) where the right-side terms of (45) are the coenergy of the (−4,−7) element, energy of the capacitor, and coenergy of the inductor, respectively. If the differential equation (43) represents the equilibrium of the generalized velocities and not the equilibrium of generalized forces as we hitherto supposed, it can be considered the generalized Kirchhoff law for the current variables i (β) , and this equation can be modeled by a parallel connection of (α,β) elements.

Discussion
It is proved that Hamilton's variational principle holds for circuits built from higher-order elements that occupy the common Σ-diagonal of Chua's table. Furthermore, it is shown that this principle is voided if any other element outside this diagonal appears in the circuit. For the case of elements located on the odd Σ-diagonal, the Lagrangian is the difference between the generalized energies of the elements of L and C types. For the even Σ-diagonal, it is the difference between the generalized contents of the elements of + R and − R types. Created in this way, the Lagrangian has the form (7),