According to Kennedy's Robust op-amp realization of Chua's circuit(1992), the differential equations satisfied by several physical quantities in Chua's circuit are
C1dvC1dt=G(vC2−vC1)−g(vC1)C2dvC2dt=G(vC1−vC2)+iLLdiLdt=−vC2
Where G=1/R. The circuit diagram is given at the bottom, where the quantity represented by each letter is indicated clearly.
When Chua invented his circuit, he deliberately designed it in such a way that the three equilibrium points were unstable. But I don't know how he did this, although I can analyze the stability using eigenvalues.
For small x, f(x)=−mx for some positive m. So the equations above are linear for samll x. In order to analyze the stability of the zero solution vC1=vC2=iL=0, I write x=[vC1,vC2,iL]T and M=[−ab0c−cd0−f0] where all letters are positive constants, and a=1/RC1,b=1/RC1+m/C1 and so on. So the differential equations above can be written as ˙x=Mx.
The characteristics equation det is \lambda^3+(a+c)\lambda^2+(ac-bc+fd)\lambda+fda=0.
If the trivial solution \mathbf x=0 is a stable solution, then according to a theorem(see here), the real parts of all eigenvalues of M have to be negative. Thus it is not hard to deduce that the coefficients of the characteristic equation must be all positive. Therefore ac-bc+fd>0. Substituting the expressions for a,b,\ldots,f in terms of R,C_1,... and rearranging gives m<\frac{RC_1}{L}.
Questions:
- So does it mean that when m\geq\frac{RC_1}{L} the zero solution is always unstable? Does this result help us to choose the values of parameters?
- Intuitively, for the trivial solution to be unstable, we want the Chua's diode to develop power at a faster rate than the resistor consuming power. If we ignore the capacitors for a moment, then this means that m<1/R. However, the R in m\geq\frac{RC_1}{L} is on the numerator rather than the denominator. why are the two inequalities (m\geq\frac{RC_1}{L} and m<1/R) so different? this seems quite counter-intuitive. Also, intuitively, for the trivial solution to be unstable, we want the Chua's diode to develop power at a faster rate than the resistor consuming power. If we ignore the capacitors for a moment, then this means that m<1/R. However, the R in m\geq\frac{RC_1}{L} is on the numerator rather than the denominator. why are the two inequalities (m\geq\frac{RC_1}{L} and m<1/R) so different? this seems quite counter-intuitive.
Also, can I improve my results? This is just a sufficient condition for the equilibrium to be unstable. Can I somehow make it a necessary condition as well?
Or are there any better methods I can use instead of evaluating eigenvalues?
EDIT: It appears that the condition m\geq\frac{RC_1}{L} should be also sufficient if all eigenvalues are real.
No comments:
Post a Comment