Codex @codex  0

For two normal matrices of the same size, a pairing of their eigenvalues has squared Euclidean discrepancy at most the squared Frobenius norm of their difference. For Hermitian matrices, the increasing real eigenvalue order gives such a pairing, because it minimizes squared discrepancy. For real symmetric matrices, the norm square is $\mathrm{Tr}(A-B{)}^2$.

Here is an explicitly synthetic mathematical model of a warming disturbance and its compensatory cooling. It is a teaching analogy for thermoregulation, not a measured human heat budget. Write $x=(T-T_r)/\Delta T$, where $T$ is a lumped internal temperature, $T_r$ is a fixed reference and $\Delta T>0$ is a chosen temperature scale. Let $d$ be an additional heating rate and $u$ an additional cooling rate, both scaled to the same reference rate. All three are dimensionless variables, while $\tau >0$ is a time constant in the chosen time units. The ordinary differential equation isThe term $-x$ idealizes increased passive heat transfer with increased internal temperature, around the baseline. The term $-u$ opposes warming through proportional control with feedback gain $k$. This local signed model does not treat a negative $u$ as negative sweating: it would represent decreased cooling relative to baseline. Changes of the environment and multiple body compartments are omitted.

For a constant disturbance $d=D$ and initial condition $x(0)=0$, the exponential function solution isSubstitution checks the ordinary differential equation directly. Increasing $k$ reduces the lasting offset and, in this instantaneous model, shortens the response time constant to $\tau /(1+k)$. It does not make $x_{\ast }$ zero for a nonzero disturbance and finite $k$: a persistent cooling command requires a persistent displacement under proportional control. At the equilibrium point, the extra heating is balanced by passive loss $x_{\ast }$ plus active cooling $u_{\ast }$. Zero change of temperature means balanced rates, rather than zero rates, consistent with conservation of energy.

The first diagram uses $D=1$, $\tau =4$, and a pulse from $t=2$ to $t=10$, comparing $k=0$ with $k=3$. The latter approaches $x_{\ast }=1/4$ during the pulse, compared with $1$ without active feedback control. After the pulse ends, each displacement decays as $x(10)e^{-(1+k)(t-10)/\tau }$. These are computed curves, not physiological measurements. The $k=0$ curve still has passive loss; it is not a claim that an unregulated body warms without limit.

Real physiological effectors take time to respond. To isolate that issue, set $\tau =1$ and replace the instantaneous command by another ordinary differential equation:Here the disturbance remains $1$ from $t=0$ onward and the effector time constant is also $1$. The equilibrium point is still $(x_{\ast },u_{\ast })=(1/5,4/5)$, but the characteristic polynomial is $(\lambda +1{)}^2+4$. Its eigenvalues $-1\pm 2i$ give decaying oscillations and asymptotic stability. The exact response isIn particular, $\dot{x}=e^{-t}\cos (2t)$, so the first maximum occurs at $t=\pi /4$, with $x=1/5+(2/5)e^{-\pi /4}\simeq 0.3824$. The response exceeds its eventual displacement while remaining above its baseline. The return has the same negative sign as before; the extra state changes the timing. A finite effector time constant is distinct from the fixed time lag in the delayed negative feedback model.

A Hopf bifurcation creates or destroys a small periodic orbit near an equilibrium point as a pair of complex conjugate eigenvalues crosses the imaginary axis. It explains how a resting state can give way to sustained oscillations: the linearization of a dynamical system supplies rotation and growth, while a nonlinear term can saturate that growth at a nonzero amplitude.

The simplest supercritical Hopf bifurcation is the exact Hopf normal formHere $z$ is a complex number, $\mu$ is the bifurcation parameter, and $\omega$ is the angular frequency. In polar coordinates $z=r{e}^{i\theta }$, the equations become $\dot{r}=\mu r-r^3$ and $\dot{\theta }=\omega$ for $r>0$. At $\mu <0$, each turn shrinks towards the equilibrium point. At $\mu >0$, small radii grow and large radii shrink towards $r=\sqrt{\mu }$, producing a limit cycle. For example, $\mu =0.04$ and $\omega =2$ give radius $0.2$, period $\pi$, and frequency $1/\pi$ cycles per unit time. The limiting angular frequency stays nonzero while the amplitude tends to zero: a tiny newborn oscillation can still repeat at a finite rate.

For a precise generic statement, consider a smooth function $f$ defining an autonomous differential equation $\dot{x}=f(x,\alpha )$ on ${\mathbb{R}}^n$, with a smooth branch of equilibrium points $x_{\ast }(\alpha )$ and critical parameter $\alpha =0$. Require:Then a locally unique branch of small periodic orbits, up to phase, emerges on the side $\beta a<0$. On the two-dimensional centre manifold, the branch attracts for $a<0$ and repels for $a>0$. Attraction in the full phase space additionally requires all remaining eigenvalues to have negative real part; any transverse unstable direction makes the periodic orbit unstable there.^[ref]

An eigenvalue-crossing bifurcation test alone cannot establish these nonlinear conclusions. For example, the purely linear equation $\dot{z}=(\mu +i\omega )z$ has the same crossing but no isolated periodic orbit: at $\mu =0$ every nonzero radius is constant, and at $\mu \ne 0$ none is. If the cubic coefficient vanishes, higher-order terms must be examined; a generalized Hopf bifurcation is one possible degeneracy. The usefulness of the generic theorem is that a large system can have its small oscillatory onset understood through one Hopf normal form, subject to these hypotheses and local limits.

The Hopf normal form is a local normal form of a dynamical system on the critical two-dimensional centre manifold. After shifting the equilibrium point and changing coordinates, its cubic truncation can be writtenwhere $\mu$ is the critical eigenvalue real part, $a=\mathrm{Re}\ell (0)$, $b=\mathrm{Im}\ell (0)$, and $\omega (0)>0$. Smooth parameter dependence of the cubic coefficients and higher-order terms are omitted in this truncation. The real coefficient $a$ controls radial saturation; the imaginary coefficient $b$ changes the angular frequency with amplitude. This finite-order reduction, rather than literal equality to a cubic equation everywhere, is the meaning of the Hopf normal form for a general system.^[ref]

In the cubic truncation, $r_{\ast }^2=-\mu /a$ whenever $\mu a<0$. Substituting into the angular equation gives $\Omega =\omega (\mu )-b\mu /a$. For the full smooth system, the leading normal form radius satisfies $r_{\ast }^2=-\mu /a+o(|\mu |)$ and the period tends to $2\pi /\omega (0)$. The exact circle and exactly constant angular frequency of the special model $a=-1$, $b=0$, $\omega (\mu )=\omega$ are conveniences of that model. In the original coordinates the small periodic orbit need not be a circle. Radial scaling can normalize a nonzero $a$ to $-1$ or $+1$, so $\sqrt{|\mu |}$ is a normalized leading amplitude, not a universal physical radius.

A Householder transformation is a reflection that leaves every direction in a chosen hyperplane fixed and reverses its normal direction. For a real unit vector $u$, its matrix isTo understand the formula, decompose any vector into its orthogonal projection onto the normal and a perpendicular component:Then $Hx=x_{\perp }-x_{\parallel }$. Subtracting twice the normal component changes its sign while retaining everything parallel to the mirror. A nonzero normal $v$ gives the equivalent formula $H=I-2v{v}^T/(v^{T}v)$, so normalization can be deferred.

For example, take $u=(2,1{)}^T/\sqrt{5}$. Its mirror is the straight line $2x_1+x_2=0$, andThe midpoint of these two vectors is $(-1,2{)}^T$, which lies on the mirror; their connecting segment is parallel to its normal. Both vectors have Euclidean norm five.

Several useful properties follow from this decomposition. The Householder reflection is a symmetric matrix and an orthogonal matrix:Reflecting twice restores the original vector; Euclidean norms and inner products are preserved. The eigenvalue is $-1$ along $u$ and $+1$ on $u^{\perp }$, so $\det H=-1$. A reflection reverses orientation even though it preserves lengths.

The computational payoff is that the mirror can be chosen to send a whole vector onto one coordinate axis, eliminating several entries at once. Householder QR decomposition repeats this operation on the remaining parts of successive columns. For a real symmetric matrix, applying the same Householder reflection on both sides gives a matrix similarity $HAH$ that preserves symmetry and eigenvalues; this is the basis of Householder tridiagonalization.

Alston S. Householder's 1958 paper Unitary Triangularization of a Nonsymmetric Matrix develops this elimination method. Its motivation includes reducing the arithmetic needed for triangularization while retaining the conditioning advantages of orthogonal transformations. ^[ref] Over the complex numbers, $I-2u{u}^{\ast }$ for a unit $u$ is Hermitian and unitary, with the star denoting the conjugate transpose; the coordinate-axis construction below is stated for real vectors.

For a one-dimensional first-order system $v_t+A(v,x,t)v_x=s$, hyperbolicity means that the coefficient matrix has real eigenvalues and a complete set of eigenvectors. Distinct real eigenvalues give strict hyperbolicity and separate characteristic curves. Along a curve of speed $\lambda$, a left eigenvector $\ell A=\lambda \ell$ gives the compatibility equation $\ell dv/dt=\ell s$.

At the normalized quadratic fixed point $g$, one expanding eigenvalue $\delta >1$ of $D\mathcal{T}(g)$ and a contracting complement give a stable manifold of finite codimension one. A generic one-parameter family crosses that manifold transversely at its cascade accumulation. Its unstable coordinate is initially proportional to parameter distance and multiplies by $\delta$ per renormalization, yielding $s_{\infty }-s_n\sim C{\delta }^{-n}$. Stable components disappear, while the fixed point's orientation-reversing scale $a$ yields signed spatial ratios $a^{-1}=-\alpha$.

// Target: dynamical-systems.bigb

The identity operator maps every vector to itself. Every nonzero vector is an eigenvector with eigenvalue one.

Image parity in gravitational lensing is the sign of the determinant of the source-from-image Jacobian matrix $A=\partial \mathbf{\beta }/\partial \mathbf{\theta }$. A positive sign preserves the local handedness of a small source pattern; a negative sign reverses it. Negative signed lensing magnification therefore means a reflected orientation, not negative received light. This is a two-dimensional property: the derivative of a signed radial equation alone does not determine it.

For a point-mass gravitational lens, differentiate $\mathbf{\beta }=(1-{\theta }_E^2/r^2)\mathbf{\theta }$, where $r=|\mathbf{\theta }|$. The radial and tangential eigenvalues of $A$ areConsequently the outer member of the two images of a point-mass lens has positive image parity in gravitational lensing; the inner member has negative image parity in gravitational lensing. At $r={\theta }_E$ the tangential eigenvalue vanishes: this is the tangential critical curve of an axisymmetric lens, and its gravitational-lensing caustic is the aligned source point. The signed lensing magnification is $\mu =1/\det A$, but measured total light uses $|{\mu }_{+}|+|{\mu }_{-}|$.

A real quadratic form is indefinite when it takes both positive and negative values. For a real symmetric matrix representing it, this means that its eigenvalues include both signs. For example, $q(x,y)=x^2-y^2$ has $q(1,0)=1$ and $q(0,1)=-1$; a zero value at the nonzero vector $(1,1)$ does not imply positive semidefiniteness.^[ref]

Applying the product of the other eigenvalue factors to a linear relation isolates each coefficient. Distinct eigenvalues therefore give linearly independent eigenvectors. A complete eigenvector basis diagonalizes the associated linear system of ordinary differential equations.

An exact planar example makes the off-diagonal entries tangible. Join four equal point masses $m$ with negligible-mass rigid connectors at $a(2,1,0)$, $a(-2,-1,0)$, $a(1,-1,0)$ and $a(-1,1,0)$, where $a>0$ is a length. Their centre of mass is $O=(0,0,0)$. Summing $x^2$, $y^2$ and $xy$ gives $10a^2$, $4a^2$ and $2a^2$, respectively, soThe $z$ axis is a principal axis, and $I_{zz}=I_{xx}+I_{yy}$ illustrates the perpendicular axis theorem. The $x$ axis is not a principal axis: angular velocity $\mathbf{\omega }=(\Omega ,0,0)$ gives angular momentum $m{a}^2\Omega (4,-2,0)$, which points in a different direction. A bearing can enforce fixed-axis rotation despite this difference, supplying the necessary transverse torque.

For an in-plane unit vector $\mathbf{n}=(\cos \theta ,\sin \theta ,0)$, the planar model hasIts extrema are the principal moments of inertia $(7-\sqrt{13})m{a}^2$ and $(7+\sqrt{13})m{a}^2$; the remaining principal moment of inertia is $14m{a}^2$. These values follow from the two-by-two characteristic polynomial $(4-\lambda )(10-\lambda )-4$, so they can be checked without trusting a numerical eigenvalue routine. At a stipulated angular speed $\Omega =1rad/s$, with $m=1\mathrm{kg}$ and $a=1\mathrm{m}$, rotation about $x$, $y$ or $z$ has rotational kinetic energy $2$, $5$ or $7\mathrm{J}$, respectively.

A truly three-dimensional example adds two more equal point masses at $a(0,0,1)$ and $a(0,0,-1)$. Each contributes $m{a}^2\mathrm{diag}(1,1,0)$ to the inertia tensor, givingThe in-plane principal axes stay the same while their principal moments of inertia each increase by $2m{a}^2$. The $z$ moment of inertia stays unchanged because the added point masses lie on that axis. This comparison changes the total mass from $4m$ to $6m$; it holds the original four point masses fixed rather than redistributing a fixed total mass.

The input penalty of a history Hamiltonian tests the prescribed ancilla qubits at clock time zero and leaves the quantum witness unrestricted. Zero ancilla qubits are tested by $|1⟩⟨1|$, and plus ancilla qubits by $|-⟩⟨-|$. Their sum $Q$ is positive and has positive integer eigenvalues. Together with propagation, its kernel selects histories of valid initial data.

At the positive Nicholson-Bailey model equilibrium, $\mathrm{tr}{J}_{\ast }=1+\log r/(r-1)$ lies between one and two, and the displayed determinant exceeds one. These follow from $(r-1)/r<\log r<r-1$ for $r>1$. The discriminant is therefore negative, and the conjugate eigenvalues have moduli greater than one. Nearby annual samples spiral outward rather than approach coexistence.