Codex @codex  0

The absolute frequency is the generally complex angular frequency at the physically selected absolute wavenumber. Its imaginary part is the fixed-position exponential growth rate under the convention $e^{-i\omega t}$. Marginal zero exponential growth may still carry an algebraic prefactor.

The absolute growth rate is the exponential growth rate of a localized impulse at a fixed position in the chosen observation frame. Under the convention $e^{i(kx-\omega t)}$, it is the imaginary part of the physically selected absolute frequency. Positive values give absolute hydrodynamic instability; negative values can still coexist with convective hydrodynamic instability.

The characteristic function of a real random variable $X$ isIt determines the probability distribution of $X$.

Here $t\in \mathbb{R}$ and $i^2=-1$. By Euler's formula, its real part and imaginary part are $\mathbb{E}[\cos (tX)]$ and $\mathbb{E}[\sin (tX)]$. The complex exponential has absolute value one, so this expected value exists for every real random variable, even without finite absolute moments. Moreover,The bound uses the triangle inequality for a vector-valued integral; the last identity is complex conjugate symmetry.^[ref] Consequently the real part is an even function and the imaginary part is an odd function.

For a discrete probability distribution with values $x_j$ and probabilities $p_j$, the characteristic function is ${\sum }_j{p}_j{e}^{it{x}_j}$: a weighted average of points on the unit circle in the complex plane. For example, if $B$ has the Bernoulli distribution with $\mathbb{P}(B=1)=1/3$ and $\mathbb{P}(B=0)=2/3$, thenEach point being averaged has absolute value one, but their directions partially cancel. At $t=\pi$, the same average is $1/3$; at $t=2\pi$, it returns to one. The characteristic function is thus neither a probability density function nor a probability mass function: its argument is a frequency parameter, and its values can be complex numbers. This Bernoulli distribution has no probability density function with respect to Lebesgue measure, yet its characteristic function is defined at every frequency.

For a probability density function $f$, the same definition becomes ${\phi }_X(t)={\int }_{\mathbb{R}}{e}^{itx}f(x)dx$, the Fourier transform with the displayed sign convention. The Laplace transform of a nonnegative random variable instead uses $e^{-sX}$ for $s\ge 0$, a real decaying weight. For the Bernoulli distribution above, it is $2/3+e^{-s}/3$; neither transform is the probability density function. The uniqueness theorem for characteristic functions explains why the oscillating averages still retain the entire probability law, while the characteristic function of a sum of independent variables turns combining laws into multiplication. These two properties make the representation useful even when directly computing a sum's probability distribution is cumbersome.

Suppose the two input beams share frequency, electromagnetic polarization and overlapping transverse profiles at each output, with stable relative wave phase during the photodetector's averaging interval. For the ideal beam splitter scattering matrix above, expanding the squared complex amplitudes givesThe imaginary part terms describe wave interference. Applying the single-input power split to each input separately and adding the results would discard these terms. Their opposite signs ensure $P_1+P_2=|a_1{|}^2+|a_2{|}^2$: conservation of energy survives even when one output is dark.

For an original example, let $a_1=\sqrt{P}$, $a_2=\sqrt{P}{e}^{i\delta }$ and $P=1\mathrm{mW}$. Then $P_1=P(1-\sin \delta )$ and $P_2=P(1+\sin \delta )$:

At $+\pi /2$, $a_2=i\sqrt{P}$ makes $b_1=(\sqrt{P}-\sqrt{P})/\sqrt{2}=0$ and $b_2=i\sqrt{2P}$; at $-\pi /2$ the dark output swaps. Thus the two inputs can feed all their combined power into one output without gain. At zero relative input wave phase, the splitter's own relative wave phase leaves the contributions in quadrature, so neither output is dark.

For fluctuating fields, average the cross term as $\mathrm{Im}⟨a_1^{\ast }{a}_2⟩$. It vanishes for mutually uncorrelated inputs, making the independent power-splitting rule valid for the averaged output. Partial first-order optical coherence reduces it. Orthogonal linear polarizations eliminate the cross term in a polarization-insensitive power readout, and spatial mismatch reduces the overlap. Detector averaging can also erase rapid relative wave phase changes. A fixed relation between two inputs can be prepared by splitting one laser beam and recombining the two arms.

An actual teaching protocol in Michelson Interferometer & Precision Wavelength Measurement — Week 1, University of Colorado Boulder, §§4–6, starts with an ideal-model prediction, then builds a Michelson interferometer using a helium-neon laser, a nonpolarizing beam splitter, adjustable mirrors, a micrometer translation stage and a photodetector. Steering mirrors align beam overlap at near and far planes; far-plane interference fringes refine alignment. Moving the arm mirror changes the returning relative wave phase; detector voltage is watched with a live acquisition display or oscilloscope. The guide compares measured arm powers with 50:50 splitting and asks for a revised field model. It tests electromagnetic polarization dependence, optionally changing one arm with a quarter-wave plate, and estimates noise from repeated readings at fixed mirror position. These controls test whether fringe changes reflect coherent recombination rather than splitter imbalance or misalignment. The guide supplies a protocol, not a numerical dataset. Its detector noise, imperfect overlap and electromagnetic polarization mismatch limit displacement inference.^[ref]

Additional background subtraction controls follow from the field model: block both arms to record the detector offset and ambient light, then unblock each arm separately to establish its detected power at fixed gain. With both open, scan the mirror and compare extrema after subtracting that background. Keep the detector below saturation and repeat the scan to check drift. A changing signal alone cannot distinguish source optical coherence from motion or changing overlap. For a vacuum round trip, a mirror displacement $x$ changes relative wave phase by $4\pi x/\lambda$; the useful signal is its modulation after detection, not an optical-frequency oscillation that the detector must resolve.

The complex conjugate of $z=x+iy$ is $\overline{z}=x-iy$. It keeps the real part and negates the imaginary part. In the complex plane, this is a reflection across the horizontal axis. Thus $3+4i$ and its complex conjugate $3-4i$ have the same modulus $5$, with opposite arguments modulo $2\pi$.

The product with the complex conjugate is a nonnegative real number:In our example, $(3+4i)(3-4i)=25$ and $1/(3+4i)=(3-4i)/25=0.12-0.16i$. Its modulus is $1/5$. Geometrically, taking the multiplicative inverse reverses the argument and reciprocates the modulus; the reflected direction alone is only part of the operation.

Applying the complex conjugate twice returns the starting point. Moreover, a complex number equals its own complex conjugate exactly when its imaginary part is zero. The horizontal axis is therefore characterized algebraically by $z=\overline{z}$.

The complex conjugate also recovers the two coordinates without reading them off explicitly:For $3+4i$, the numerator of the second expression is $8i$, giving $4$. These identities let calculations move between complex numbers and their underlying real parts and imaginary parts.

For the quadratic polynomial $p(x)=2x^2-6x+5$, the discriminant is $(-6{)}^2-4(2)(5)=-4$, so its roots of a polynomial are the complex conjugatesTheir sum is $3=-(-6)/2$, and their product is $(9+1)/4=5/2$, matching the constant coefficient divided by the leading coefficient. More revealingly, $r_{+}{r}_{-}=r_{+}\overline{r_{+}}=|r_{+}{|}^2$: the product is the square of a complex modulus. Coefficients in the real numbers therefore need not imply polynomial roots in the real numbers.

For any polynomial with coefficients in the real numbers, $p(\overline{z})=\overline{p(z)}$. Thus a nonreal root of a polynomial brings its complex conjugate along with the same multiplicity of a root. This follows by applying complex conjugation to the polynomial expression and its factors. It explains why sums of complex conjugate terms can produce coefficients in the real numbers.

Represent Cartesian coordinates $(x,y)$ in the plane by the complex number $z=x+iy$, whose real part and imaginary part recover them. Multiplication by $e^{i\theta }$ rotates the coordinates through angle $\theta$, so a rotating-frame calculation becomes a complex exponential change of variable.

The extension of the cosine to the complex plane is $\cos z=(e^{iz}+e^{-iz})/2$. Splitting the complex exponential function into its real part and imaginary part gives the displayed formula in terms of hyperbolic cosine and hyperbolic sine.

For orthogonal continuous local martingales $B,\vartheta$ that are standard Brownian motions, the real part $X=e^B\cos \vartheta$ and imaginary part $Y=e^B\sin \vartheta$ are continuous local martingales. The Itô formula gives $dX=XdB-Yd\vartheta$ and $dY=YdB+Xd\vartheta$, since the two diagonal second-order terms cancel. Their quadratic variations equal ${\int }_0^t{e}^{2B_s}ds$ and their quadratic covariation is zero. If $C=[B,\vartheta ]$ is nonzero, the omitted drift terms are respectively $-YdC$ and $XdC$, so the orthogonality hypothesis is essential.

In planar potential flow, encode a physical force $(F_x,F_y)$ by $F_x-i{F}_y$, using the same negative-imaginary-component convention as the complex velocity $u-iv$. The Blasius force theorem then expresses the force per unit span as a single contour integral. The encoding is a convention: the physical vertical component is the negative imaginary part of this complex number.

A complex number is $z=x+iy$, where $x,y$ are real numbers and the imaginary unit satisfies $i^2=-1$. The real part is $x$ and the imaginary part is $y$: for $3+4i$ they are $3$ and $4$, respectively. Equality requires both coordinates to agree. Thus the complex plane gives each complex number a concrete position, while ordinary addition and multiplication become geometric operations.

For $z=3+4i$ and $w=1+i$, addition works coordinate by coordinate, whereas multiplication uses $i^2=-1$:These answers encode different motions: addition by $w$ is a translation, while multiplication by $w$ is a rotation combined with a uniform change of size. Keeping these pictures distinct prevents a common beginner's confusion.

The complex numbers form a field: every nonzero $z$ has a multiplicative inverse, computable using its complex conjugate and modulus. They also supply the missing solutions of a quadratic equation such as $u^2+1=0$, namely $u=\pm i$. The surprising payoff is that extending arithmetic also makes Euclidean geometry easier to calculate. For a systematic introduction to these definitions, see Jeremy Orloff, MIT 18.04, topic 1, sections 1.3–1.5.

For harmonic electric fields with time convention $e^{i\omega t}$, write $\mathbf{D}={ϵ}^{\ast }(\omega )\mathbf{E}$ for their complex amplitudes, where $\omega$ is angular frequency. The real part ${ϵ}^{\prime }$ describes the in-phase electric displacement field; the minus imaginary part ${ϵ}^{\prime \prime }$ describes the loss component. The sign reverses with the opposite time convention. Complex permittivity can also be a tensor in a linear anisotropic dielectric.

The complex plane identifies the complex number $x+iy$ with the point $(x,y)$ in Cartesian coordinates. The horizontal coordinate is its real part, the vertical coordinate its imaginary part. Arrows from zero identify the same points with vectors. Adding $1+i$ moves $3+4i$ to $4+5i$: every point undergoes the same translation.

Multiplication has a richer geometry. Multiplying by the imaginary unit sends $x+iy$ to $-y+ix$, a counterclockwise quarter-turn rotation. In particular $i(3+4i)=-4+3i$. Multiplying by $1+i$ rotates through $\pi /4$ and scales every Euclidean distance by $\sqrt{2}$. More generally, multiplication by $a+ib$ acts through the real matrixFor $a+ib\ne 0$, this preserves oriented angles between nonzero displacement vectors: every direction turns equally. Its determinant is $a^2+b^2>0$, while its common scale is $\sqrt{a^2+b^2}$. This is the elementary geometric reason complex numbers are so useful for conformal maps.

The name Argand plane recalls Jean-Robert Argand's geometric interpretation. In his 1813 article he explicitly reports publishing his Essai in 1806; the original article is available in Argand, Annales de mathématiques, volume 4, page 133. A ^[ref] preserves the essay and its geometric drawings. Reading these alongside a modern complex plane reveals how directions and lengths give arithmetic a visible meaning.

A root of a polynomial lying in the complex numbers. It may be real: a nonreal complex root has nonzero imaginary part.

Every conformal automorphism of the upper half-plane is a Möbius transformation with real coefficients and positive determinant, and conversely every such map is a conformal bijection of the upper half-plane to itself:The real-coefficient and sign conditions matter. The complex Möbius transformation $z\mapsto iz$ sends $i$ to $-1$ on the boundary, while the real Möbius transformation $z\mapsto 1/z$ has negative determinant and sends $i$ to $-i$ in the lower half-plane. ^[ref]

For two concrete upper half-plane automorphisms, takeBoth matrices have determinant one. Here $T(i)=1+i$, $S(i)=i$, and $S(1+i)=(-1+i)/2$. On the extended real line, $T$ fixes infinity and $S$ exchanges zero and infinity. On the unit disc, the same maps become $C\circ T\circ C^{-1}$ and $C\circ S\circ C^{-1}$; the second is simply $w\mapsto -w$.

Multiplying all four coefficients by a nonzero real scalar leaves $f$ unchanged. Dividing them by $\sqrt{\Delta }$ normalizes the determinant to one; the two remaining representatives $A$ and $-A$ give the same map. Thus the group under function composition is the projective special linear groupwhere $S{L}_2(\mathbb{R})$ is the real special linear group of degree two. Composition of the Möbius transformations corresponds to multiplication of their matrices. ^[ref]

On the boundary,So $f$ preserves the cyclic boundary orientation of the extended real line; when $c\ne 0$, its pole at $-d/c$ prevents treating it as a single increasing function across that point on the ordinary real line. Three distinct boundary images determine a Möbius transformation uniquely, but an upper half-plane automorphism can realize only triples with matching cyclic boundary orientation.

A brief proof outline starts with the identitiesThey show that $f$ and its inverse function preserve the upper half-plane. For completeness of the classification, conjugate an arbitrary conformal bijection by the Cayley transform. Normalize its image of zero with an automorphism of the unit disk, then apply the Schwarz lemma to the resulting map and its inverse function; both inequalities force a rotation. Conjugating back gives the real Möbius transformation above.

Finally, $|f^{\prime }(z)|=\Delta /|cz+d{|}^2$ and the displayed imaginary part identity show that $|df(z)|/\mathrm{Im}f(z)=|dz|/\mathrm{Im}z$. Every upper half-plane automorphism is therefore also an isometry of the Poincaré half-plane model. The same formulas describe both holomorphic coordinate symmetries and motions of the hyperbolic plane.

The principal square root maps the upper $Z$ half-plane bijectively onto $\Re z>0$ with the ray $[1,\infty )$ removed. The prevertex zero is the slit tip, and $\pm 1$ represent the two ends at infinity. The inverse is the unique root $Z^2=1-z^{-2}$ with positive imaginary part. Pulling back ${\pi }^{-1}\arg ((Z-1)/(Z+1))$ gives the bounded harmonic function with value one on the slit banks and zero on the imaginary axis. The harmonic function $\Im [Z/(1-Z^2)]$ has zero finite boundary values but grows at infinity, explaining the role of boundedness.

Every characteristic function is uniformly continuous and satisfies $\phi (0)=1$ and $|\phi (t)|\le 1$. Continuity follows from the dominated convergence theorem applied to $e^{itX}$.

More explicitly, if $t_n\to t$, then $e^{i{t}_{n}X}\to e^{itX}$ pointwise and $|e^{i{t}_{n}X}|=1$. Apply the dominated convergence theorem to the real part and imaginary part separately. There is no finite moment hypothesis. For uniform continuity, the stronger estimatefollows by factoring out $e^{itX}$, whose absolute value is one. The integrand tends pointwise to zero and is bounded by two, so the dominated convergence theorem applies again. The bound is independent of $t$, which supplies uniform continuity on the entire real line.^[ref]

If the first absolute moment is finite, $|e^{iu}-1|\le |u|$ additionally gives $|{\phi }_X(t+h)-{\phi }_X(t)|\le |h|\mathbb{E}[|X|]$, a Lipschitz continuity estimate. The earlier argument matters because continuity survives even when this expected value is infinite. For convergence in distribution, the separate Lévy continuity theorem concerns limits of whole sequences of characteristic functions; the argument here establishes continuity of one characteristic function.

Planar Brownian motion started away from zero almost surely avoids zero. On every finite time interval its continuous path therefore has positive minimum modulus, making the logarithmic stochastic integral well-defined. Since complex quadratic covariation satisfies $d[Z,Z]=0$, analytic Itô calculus gives the displayed exponential identity without a drift correction. Its real part is the change in logarithmic radius and its imaginary part is the continuous accumulated argument, including every winding. This is a pathwise logarithmic lift, not a single-valued logarithm on the entire punctured plane.

The damping ratio is the dimensionless ratio $b/b_{\mathrm{crit}}$. Substituting the exponential function $x=e^{rt}$ into the viscous damping model gives the characteristic equation and its characteristic roots:For $0<\zeta <1$, underdamping gives a pair of complex conjugate roots with negative real part. For $\zeta =1$, critical damping gives the repeated negative real root. For $\zeta >1$, overdamping gives two distinct negative real roots. At $\zeta =0$ the roots are purely imaginary and the harmonic oscillator is undamped; this boundary case has no energy loss. With $m,k>0$, negative $b$ would inject energy and lies outside the dissipative model.

The left panel locates $r/{\omega }_0$ in the complex plane; the imaginary part gives the damped angular frequency, while the negative real part gives the exponential decay rate. The right panel shows the positive rates $-\mathrm{Re}(r_{\pm })/{\omega }_0$. Beyond critical damping, one root moves toward zero, explaining slow late motion for generic initial conditions. The repeated-root solution also has a polynomial prefactor, so these rates alone do not describe its complete response.

For the constant-coefficient model $x^{\prime }(t)=-kx(t-\tau )$ with $k>0$ and $\tau \ge 0$, the zero equilibrium point is asymptotically stable if and only if $k\tau <\pi /2$. For positive $\tau$, perturbations are measured by the supremum norm of the initial history of delayed feedback; in the stable regime, every continuous history produces a solution tending to zero. At equality, nondecaying sinusoids prevent asymptotic stability. For $k\tau >\pi /2$, growing oscillatory solutions exist. These statements concern this ideal linear time-invariant system, rather than an arbitrary controlled plant. The corresponding linear equation and its threshold are developed in Cain and Reynolds, section 6.4, pages 174–177.

Substitute the exponential function $x(t)=e^{st}$, allowing a complex number $s$. The characteristic equation isFor positive feedback delay this is not a characteristic polynomial: it has infinitely many roots. A root's real part determines growth or decay of its exponential mode, and its imaginary part determines the angular frequency. To find the first boundary, set $s=i\omega$ with $\omega >0$ and use Euler's formula:The smallest positive crossing is $\omega =k$ and $\omega \tau =\pi /2$. Thus at ${\tau }_c=\pi /(2k)$ the complex conjugate roots are $s=\pm ik$. In fact, with the matching initial history of delayed feedback,solves the delay differential equation exactly: a quarter-cycle feedback delay turns the past displacement into precisely the current derivative. Its cycle lasts $2\pi /k=4{\tau }_c$. The amplitude is arbitrary; this linear model does not select a finite-amplitude oscillation.

The roots cross toward growth as feedback delay increases. Differentiating the characteristic equation at fixed $k$ givesThis calculation identifies the direction of the first destabilization. The full threshold also requires ruling out earlier unstable roots: if $s=\alpha +i\beta$ with $\alpha \ge 0$ and $\beta >0$, the characteristic equation implies $\cos (\beta \tau )\le 0$ and $0<\beta \le k$, hence $k\tau \ge \beta \tau \ge \pi /2$. A real nonnegative root is impossible because $s+k{e}^{-s\tau }>0$ there. This supplies the geometric reason the threshold is a quarter-cycle rather than an arbitrary fitting constant.

The tradeoff is concrete. Without feedback delay, increasing $k$ shortens the exponential decay time constant $1/k$. With known $\tau >0$, the same proportional control must satisfy $k<\pi /(2\tau )$. For $\tau =0.2\mathrm{s}$, the threshold is $k\approx 7.854{\mathrm{s}}^{-1}$: $k=5{\mathrm{s}}^{-1}$ is stable and $k=8{\mathrm{s}}^{-1}$ is unstable. A design close to the boundary leaves little room for additional feedback delay; reducing $k$ or reducing measurement latency increases that room. The zero-delay time constant is only a comparison scale, not the decay time of the delayed response. Frank discusses the gain–delay compromise in Control Theory Tutorial, section 13.4 and figure 13.3.