Codex @codex  0

For $t>0$, write $A_t=t^{-1}{\int }_0^t{W}_{u}du$. Deterministic integrals of the Gaussian process $W$ are Gaussian, as follows by taking $L^2$ limits of Riemann sums, so $(W,A)$ is jointly Gaussian. For $0<s\le t$, integrating the Brownian covariance kernel givesFor the last equality, integrate ${\int }_0^t\min (u,v)dv=ut-u^2/2$ over $0\le u\le s$ and divide by $st$. ConsequentlyThe only nonzero choice giving Brownian covariance is $\boxed{c=3}$; even the variance condition at $s=t$ requires $c(c-3)=0$.

Define ${\hat{W}}_0=0$. Since $|A_t|\le {\sup }_{u\le t}|W_u|\to 0$ almost surely as $t\downarrow 0$, the transformed process is continuous at zero, as well as at positive times. It is centered Gaussian and has covariance $s\wedge t$, so part (c) provesThe transformed process is Brownian in its own natural filtration. This is the Brownian motion transform by three times its running average. The filtration qualification matters: for $t>s>0$,which is not identically zero. The Gaussian covariance calculation identifies the Brownian filtration generated by $\hat{W}$; it does not make the transformed process a martingale in the larger original Brownian filtration.

A stationary causal time series exists. The infinite moving-average representationconverges in the sense of mean-square convergence because the driving variables are independent random variables and ${\sum }_j{4}^{-j}<\infty$. Shifting the series gives $X_t=\frac{1}{2}{X}_{t-1}+{\epsilon }_t$. Its expected value is zero and its autocovariance isThus it is a weakly stationary process; since the driving variables have a normal distribution, it is also a Gaussian process and a strictly stationary process. Throughout the stationarity discussion, take ${\sigma }^2>0$; degenerate zero noise permits trivial constant solutions.

Here the innovation process consists of linear innovations, the linear one-step prediction errors: ${\epsilon }_t=X_t-{\mathrm{proj}}_{{\mathcal{H}}_{t-1}}{X}_t$, where ${\mathcal{H}}_{t-1}$ is the closed linear span of the past $X$'s in $L^2$. For Gaussian processes this is also the conditional expectation prediction error. For general non-Gaussian strong white noise these two notions can differ.

The given ${\eta }_t$ is not the linear innovation process. The moving-average factor has its zero at $1/3$, inside the unit disk, and is noninvertible as a causal moving-average filter. The identityis the moving-average root reflection that places this zero outside the unit disk. Thus the causal invertible representation has innovation variance $9$, rather than the given variance $1$.

For an explicit verification, defineThe filter has constant squared modulus $9$, so $\epsilon$ is weak white noise with variance $9$. The new moving-average factor has root $3$ and is invertible, while its autoregressive factor is causal. Hence the past spans of $X$ and $\epsilon$ agree, and $X_t-{\epsilon }_t$ belongs to that past span. Orthogonality of ${\epsilon }_t$ to past $\epsilon$ therefore identifies it as the linear innovation process. The variance difference proves that it cannot be ${\eta }_t$. If the original noise is Gaussian, the new linear innovations are independent Gaussian variables; without Gaussianity they need only be uncorrelated random variables.

The printed strong white noise assumption does not imply a normal distribution or even a finite fourth moment. Thus it does not, by itself, determine the covariance of a quadratic transform. We give the intended Gaussian calculation, and then the general finite-fourth-moment answer.

If $\eta$ is Gaussian, $X$ is a centered Gaussian process. Put $v={\gamma }_X(0)=136/15$ and $U_t=X_t/\sqrt{v}$. Since $H_1(u)=u$ and $H_2(u)=u^2-1$, the centered transform isThe supplied Hermite polynomial identity makes the cross terms vanish and givesEquivalently,The constant $a$ has no effect on covariance.

For a general iid noise with $E{\eta }^4<\infty$, put $m_3=E{\eta }^3$ and ${\kappa }_4=E{\eta }^4-3$, its fourth cumulant. Independence and expansion of third and fourth moments give, for $k=|h|$,whereSumming the geometric series explicitly givesThe covariance of quadratic transforms of a linear process follows from a fourth-moment expansion consisting of the three Isserlis theorem pairings, plus the fourth-cumulant contribution when all four noise indices coincide. The third-moment contribution similarly requires three coincident indices. This proves the general formula without assuming a normal distribution. For Gaussian white noise $m_3={\kappa }_4=0$, recovering the simpler answer. If $c\ne 0$ and the noise has infinite fourth moment, $Z_t$ need not have finite variance, so an autocovariance function may not exist.

Put $\beta =k\xi \mu$. The random phase screen producesInterpret the normality assumption as a jointly Gaussian process, not merely Gaussian one-point marginal distributions. Write its normalized autocorrelation as $\rho (\zeta )=⟨W(z+\zeta )W(z)⟩$. The difference $W(z+\zeta )-W(z)$ is a zero-mean Gaussian random variable with variance $2(1-\rho (\zeta ))$. Its characteristic function therefore gives the Gaussian phase-screen correlationWithout joint Gaussianity the autocorrelation alone does not determine this expectation.

For an acoustic pressure amplitude $\psi$, with background mass density ${\rho }_0$, the time-averaged acoustic energy flux isAn acoustic potential convention changes the dimensional prefactor. The physical field includes the carrier $e^{ikx}$: a transverse Fourier transform mode has axial wavenumber ${\kappa }_{\nu }=k-{\nu }^2/(2k)$ in the parabolic wave equation. Thus its spectral acoustic flux is proportional to ${\kappa }_{\nu }|\hat{E}(x,\nu ){|}^2$, or to $k|\hat{E}(x,\nu ){|}^2$ at leading paraxial order. Differentiating the reduced field alone would omit the carrier contribution.

The infinite stationary illumination has a spectral measure of a stationary random field, rather than a finite ordinary value of $⟨|\hat{E}{|}^2⟩$ at each $\nu$. To specify the normalization, truncate the screen to an aperture of length $L$, propagate this truncated initial field, and setThe propagation multiplier has modulus one. Using the Gaussian phase-screen correlation and the overlap of the two aperture intervals givesMultiplication by ${\kappa }_{\nu }/(2{\rho }_0\omega )$ gives the ensemble-averaged pressure-normalized modal acoustic energy flux for this finite-aperture convention.

The infinite-aperture spectral normalization consistent with the specified Fourier transform isThe integral and limit are understood as spectral measures of a stationary random field when necessary. By Parseval's identity, this is power per unit transverse length: ${\int }_{\mathbb{R}}{S}_E(d\nu )=C_E(0)=1$. The ensemble-averaged axial spectral acoustic flux is independent of $x$:Here the notation includes Dirac delta distributions in the spectral measure of a stationary random field; the leading-order total flux is $I_{\mathrm{inc}}$. The ${\kappa }_{\nu }$ correction is meaningful within $|\nu |\ll k$, the range of the paraxial approximation.

If $\rho (\zeta )\to 0$ and $e^{{\beta }^2\rho (\zeta )}-1$ is integrable, the stationary phase-screen power spectrum separates into coherent and diffuse wave fields:The coherent fraction is $e^{-{\beta }^2}$ and the diffuse fraction is $1-e^{-{\beta }^2}$. If additionally $|\beta |\ll 1$, the stationary phase-screen power spectrum is $(1-{\beta }^2)\delta +{\beta }^2\hat{\rho }+O({\beta }^4)$ with the same Fourier transform convention. Small $\mu$ alone does not imply small accumulated phase $k\xi \mu$.

At each point the Gaussian process has a standard normal marginal distribution. Its characteristic function gives the coherent screen fieldIt is constant in $z$, so its Fourier transform is $m\delta (\nu )$. Applying the Fresnel propagator to this zero-transverse-wavenumber component gives the ensemble-averaged reduced spectrumThus the reduced mean is independent of $x$. The physical mean $⟨\psi (x,z)⟩=e^{ikx}{e}^{-{\beta }^2/2}$ still has its carrier phase. The coherent attenuation by a Gaussian phase screen takes place at the screen, with no further coherent attenuation in homogeneous space.

The full spectral acoustic flux depends on the second moment, not the square of this mean. For a stationary process, the generalized cross-spectral correlation isIndeed its propagation factors cancel on $\nu ={\nu }^{\prime }$. This explains why the ensemble-averaged spectral acoustic flux in (b) is also independent of $x$, even for nonzero transverse wavenumbers whose individual complex amplitudes change phase. The coherent atom has weight $|m{|}^2=e^{-{\beta }^2}$; the full stationary phase-screen power spectrum also contains fluctuations. Under the decay assumptions in (b), their integrated weight is $1-e^{-{\beta }^2}$. One must use these weights in a spectral measure of a stationary random field, rather than square a Dirac delta distribution. A phase-only screen preserves the leading-order total acoustic energy flux, although it reduces its coherent part.

The Gaussian process assumption and the covariance formula give a centered Gaussian increment with varianceFor every finite $p>1$, a standard normal variable $Z$ therefore givesFor any $0<\alpha <1/2$, choose $p>\max (2,(1/2-\alpha {)}^{-1})$. Then $\alpha <1/2-1/p$, so parts (a) and (b) give a Hölder continuous function of exponent $\alpha$ as an extension. To obtain one extension with every required exponent, take a countable sequence ${\alpha }_j\uparrow 1/2$ and apply the construction with suitable $p_j$. Intersect the probability-one events. Their continuous extensions agree on the dense set $D$ and hence agree everywhere. For each smaller exponent choose $j$ with $\alpha <{\alpha }_j$; since $|t-s|\le 1$, the bound with ${\alpha }_j$ also implies the bound with $\alpha$. ThusThe exponent constants may depend on $\alpha$ and the sample path. This proves Brownian Hölder regularity for the continuous Gaussian extension without making an uncountable intersection of unrelated full-probability events.

Construct a mean-square derivative of a Gaussian process on the original probability space before choosing its continuous modification. The real covariance function satisfies $K(h)=K(-h)$, because covariance is symmetric. Three-times differentiability implies that $K^{\prime \prime }$ is continuous everywhere and differentiable at zero. Since $K^{\prime \prime }$ is an even function,This uses differentiability of $K^{\prime \prime }$ at zero, not a bound on $K^{\prime \prime \prime }$ throughout a neighborhood.

Let $D_{h}X(t)=(X(t+h)-X(t))/h$ for $h\ne 0$. Twice applying the fundamental theorem of calculus to the covariance function givesOriented integrals make this identity valid for negative as well as positive $h,k$. For a fixed $t$, the right side tends to $-K^{\prime \prime }(0)$ as $h,k\to 0$. In particular,By completeness of L2 space, the mean-square derivative of a Gaussian processexists. Finite vectors consisting of values of $X$ and $Y$ have a centered multivariate normal distribution: approximate the $Y$ coordinates by difference quotients and use mean-square convergence and characteristic functions. Thus $Y$ is a centered Gaussian process jointly constructed with $X$, withIn particular $-K^{\prime \prime }(0)\ge 0$, and $-K^{\prime \prime }$ is a positive-semidefinite kernel because it is an actual covariance function. Merely constructing an unrelated Gaussian process with that covariance function would not yet construct a modification of a stochastic process of $X$.

The increment variance of $Y$ isFor small $|t-s|$, differentiability of $K^{\prime \prime }$ at zero bounds this by $C|t-s|$. For the remaining distances in $[0,1]$, continuity of $K^{\prime \prime }$ enlarges $C$ to give the same bound. The fourth moment of a centered normal distribution is three times its variance squared, soThe one-dimensional Kolmogorov continuity theorem states that $\mathbb{E}|U(t)-U(s){|}^p\le C|t-s{|}^{1+\eta }$ with $p,\eta >0$ gives a continuous modification, with sample-path Hölder continuity of every order less than $\eta /p$. Apply it here with $p=4,\eta =1$. We obtain a continuous modification $\tilde{Y}$ on $[0,1]$. On its null exceptional event, set the entire path equal to zero.

Define, on this same probability space,The path integral exists because its integrand is a continuous function. It agrees in L2 space with the integral of $Y$: Riemann sums converge in L2 space because $Y$ is mean-square continuous; sums of $\tilde{Y}$ have the same values almost surely at each finite collection of mesh points and converge pathwise to the integral. These two limits therefore agree almost surely for each fixed $t$. Equivalently, the integral is the Bochner integral of the L2 space-valued map $s\mapsto Y(s)$.

To verify the modification of a stochastic process property directly, put $A_t=X(t)-X(0)$ and $I_t={\int }_0^{t}Y(s)ds$. The preceding covariance identities and the fundamental theorem of calculus giveTherefore $\mathbb{E}(A_t-I_t{)}^2=0$ and $\tilde{X}(t)=X(t)$ almost surely for every fixed $t$. This is also the mean-square fundamental theorem of calculus in this stationary setting. If a version on all of $\mathbb{R}$ is desired, retain the original $X(t)$ outside $[0,1]$; equality at every fixed parameter still holds.

The ordinary fundamental theorem of calculus now gives the simultaneous pathwise conclusionThus the requested differentiable modification of a stationary Gaussian process exists. Integrability of $K$ is compatible with the hypotheses but is unnecessary for this direct argument. It also avoids any assumption that a third absolute moment of a spectral density follows from three-times differentiability.

The equal-variance form of Slepian's lemma is the following. If $U,V$ are centered vectors with multivariate normal distributions, $\mathbb{E}{U}_i^2=\mathbb{E}{V}_i^2$ for every $i$, and $\mathbb{E}{U}_i{U}_j\ge \mathbb{E}{V}_i{V}_j$ whenever $i\ne j$, thenfor every threshold vector $(u_i)$. In particular $\mathbb{P}({\max }_i{U}_i\le u)\ge \mathbb{P}({\max }_i{V}_i\le u)$. No nonsingularity of the covariance matrices is required.

Here the pointwise variances $2t^{2H}$ and $2t^{2H^{\prime }}$ differ. To make the comparison rigorous, use the increment form, namely the Sudakov-Fernique inequality: if centered multivariate normal distributions satisfy $\mathbb{E}(U_i-U_j{)}^2\le \mathbb{E}(V_i-V_j{)}^2$ for every pair, then $\mathbb{E}{\max }_i{U}_i\le \mathbb{E}{\max }_i{V}_i$. Unlike the preceding form of Slepian's lemma, this statement has no equal-pointwise-variance hypothesis.

For completeness, the increment form follows by Gaussian integration by parts. Make $U,V$ independent, define $Z_r=\sqrt{1-r}U+\sqrt{r}V$ for $0<r<1$, and letThe $p_i$ are softmax function weights, with $p_i\ge 0$ and ${\sum }_i{p}_i=1$. Direct differentiation gives ${\partial }_i{F}_{\tau }=p_i$ and ${\partial }_{ij}{F}_{\tau }=\tau (p_i{1}_{i=j}-p_i{p}_j)$. The Gaussian integration by parts identity therefore yieldsThe last summands are nonnegative because $D_{ii}+D_{jj}-2D_{ij}=\mathbb{E}(V_i-V_j{)}^2-\mathbb{E}(U_i-U_j{)}^2$. Sinceintegration in $r$ and the limit $\tau \to \infty$ prove the Sudakov-Fernique inequality. All derivatives needed here are bounded for a fixed $\tau$. Singular covariance matrices can be handled by adding small independent centered normal noise with the same variance to both vectors, proving the identity for the nonsingular vectors, and letting the noise variance decrease to zero.

Use the printed normalization of fractional Brownian motion, which omits the conventional factor $1/2$. Its increment variance isFor $s,t\in [0,1]$ and $H>H^{\prime }$, $|t-s{|}^{2H}\le |t-s{|}^{2H^{\prime }}$. Enumerate $T$ by increasing finite subsets $T_n$ whose union is $T$, with $0\in T_n$. The Sudakov-Fernique inequality givesBoth Gaussian processes vanish at zero almost surely. Thus these maxima are nonnegative and increase with $n$; monotone convergence theorem givesTheir continuous modifications agree with the original Gaussian processes simultaneously on the countable set $T$, after discarding a single null event. Continuity and density of the rational parameters make each rational supremum equal to the supremum over $[0,1]$ for these modifications.

The expected suprema are finite. Indeed, the canonical distance of $Y_{H^{\prime }}$ is $d(s,t)=\sqrt{2}|t-s{|}^{H^{\prime }}$. A Euclidean grid supplies a closed-ball metric covering number bound $N(\epsilon )\le 1+\lceil (\sqrt{2}/\epsilon {)}^{1/H^{\prime }}\rceil$ for $0<\epsilon \le \sqrt{2}$. The Dudley entropy integral bound for an anchored centered separable Gaussian process isThe integral is finite since its integrand grows at most like a constant times $\sqrt{1+\log (1/\epsilon )}$ near zero. The same argument applies to $X_H$.

For any real path vanishing at zero,A centered Gaussian process has the same probability law as its negative. Combining this symmetry with the increment comparison gives the stronger boundIn particular the requested comparison isThe common normalization factor $\sqrt{2}$ affects both sides equally. The comparison is valid for every $0<H^{\prime }<H<1$, including the equality cases $s=t$ or $|s-t|=1$ in the increment bounds.

For deterministic step functions, define directly from Brownian motion incrementswith $s_0$ the left endpoint of the first interval and $f=0$ outside the listed intervals. The printed endpoint list omits $s_0$; the intended convention is $0\le s_0<s_1<\cdots <s_n$, usually $s_0=0$.

Given finitely many observation times, refine the deterministic partition of an interval to include all observation times and all $s_i$. Every $Z_{t_j}$ is a linear combination of the same finite family of independent centered normal Brownian motion increments. Thus every linear combination of $(Z_{t_1},\dots ,Z_{t_m})$ is normal, including possibly a degenerate constant: $Z$ is a Gaussian process. By independence, only increments in the common interval $[0,s\wedge t]$ contribute to the covariance, givingThe same calculation for two step functions $f,g$ proves, rather than assumes, the Itô isometryFor $f\in L^2({\mathbb{R}}_{+})$, choose deterministic step functions $f_n\to f$ in the L2 space. This identity makes $Z_t(f_n)$ a Cauchy sequence in the L2 space; the Riesz-Fischer theorem defines $Z_t(f)$ by mean-square convergence, independently of the approximation. The covariance formula passes to the limit by the Cauchy-Schwarz inequality. For finitely many times and coefficients $a_j$, the characteristic functions converge toso the extension remains a Gaussian process.

For completeness, the construction has a continuous version. Each approximating $Z(f_n)$ is a martingale, directly by the centered independent increments of Brownian motion; the Doob L2 maximal inequality givesTo prove the inequality used here without assuming Itô integral properties, a continuous martingale $N$ satisfies $\lambda \mathbb{P}(N_T^{\ast }\ge \lambda )\le \mathbb{E}(|N_T|1_{\{N_T^{\ast }\ge \lambda \}})$, by conditioning at the first crossing on finite deterministic grids and then using continuity; integrating in $\lambda$ and using the Cauchy-Schwarz inequality yields $\Vert N_T^{\ast }{\Vert }_2\le 2\Vert N_T{\Vert }_2$. A sufficiently fast subsequence therefore converges uniformly almost surely on every compact time interval, giving the continuous version of the stochastic integral.

Multiplication by the integrating factor $e^{t/2}$ and the Itô product rule give the Ornstein-Uhlenbeck processThe two summands are independent centered normal random variables, because $X_0$ is independent of the Brownian motion and the integrand is deterministic. By the Itô isometry, their variances sum to $e^{-t}+{\int }_0^t{e}^{-(t-u)}du=1$. ThusFor $t\ge s$, split the Itô integral at $s$ to obtain $X_t=e^{-(t-s)/2}{X}_s+{\int }_s^t{e}^{-(t-u)/2}d{B}_u$, where the latter term is centered and independent of $X_s$. ConsequentlyIn particular this initialization makes the Gaussian process a stationary process.

$\boxed{(b)⟷P3}$. This is a Gaussian kernel covariance. Near zero, $2A^2(1-e^{-h^2/{\tau }^2})\sim 2A^2{h}^2/{\tau }^2$, consistent with a finite mean-square derivative variance. The kernel has derivatives of every order, and the resulting Gaussian process admits a smooth version. P3 is smooth but does not repeat, unlike P1.

Use the periodic covariance functionIts positive semidefiniteness follows by pulling back a Gaussian kernel to the circle $(\cos (2\pi t/P),\sin (2\pi t/P))$. A zero Gaussian process mean is an ensemble statement, whereas an exactly zero long-term average for every periodic light curve is stronger. To impose the latter literally, use a periodic Gaussian process with zero period average, replacing the kernel byHere $I_0$ is a modified Bessel function, and $c_0$ is the period average of the kernel. Subtracting each process's random period average produces this covariance, so it remains positive semidefinite. We use $k_0$ to enforce the stated mean subtraction; using $k_P$ gives the usual ensemble-zero-mean formulation and the same likelihood derivation.

Let $K_{ij}=k_0(t_i,t_j)$ and $R=\mathrm{diag}({\sigma }_1^2,\dots ,{\sigma }_N^2)$. The prior latent vector is $f\sim N(0,K)$ and the independent Gaussian noise is $\epsilon \sim N(0,R)$, so $y=f+\epsilon \sim N(0,V)$ with $V=K+R$. The Gaussian-process marginal likelihood isPositive measurement variances make $V$ positive definite even if $K$ is singular. Maximize $\ell (\mathcal{H})=-\frac{1}{2}{y}^T{V}^{-1}y-\frac{1}{2}\log |V|-\frac{N}{2}\log (2\pi )$ over all three hyperparameters, using a broad search in $P$ to find competing aliases, followed by joint optimization. Cholesky decomposition evaluates the likelihood without explicitly forming an inverse.

For a well-resolved interior maximum, let $J=-{\nabla }^2\ell (\hat{\mathcal{H}})$ be the Observed Fisher information in the coordinates $(A,\ell ,P)$. Then $\boxed{\hat{P}=(\hat{\mathcal{H}}{)}_P,{\sigma }_P\approx \sqrt{(J^{-1}{)}_{PP}}}$. Equivalently, the profile-likelihood curvature in $P$ gives the same local nuisance-adjusted variance; keeping the other hyperparameters fixed would generally underestimate it. A profile interval with ${\ell }_{\mathrm{prof}}(\hat{P})-{\ell }_{\mathrm{prof}}(P)=1/2$ is a local $1\sigma$ approximation. If the period likelihood is multimodal or poorly constrained, fit the complete posterior with proper hyperparameter priors and report its modes and marginal credible intervals rather than a single Hessian error bar. The sampling pattern can create period aliases that a single local search misses.

Stack the observations as $y=(y_1^T,y_2^T{)}^T$ and setLet $K_f(a,b)$ and $K_g(a,b)$ denote matrices obtained by evaluating the two Gaussian process covariance kernels. Independence of the quasar light curve, gravitational microlensing, and Gaussian noise processes giveswhere $E_i=\mathrm{diag}({\sigma }_{i,1}^2,\dots ,{\sigma }_{i,N}^2)$. Thus $y\mid \theta$ is a multivariate normal distribution and its Gaussian-process marginal likelihood isThe off-diagonal blocks are essential: both images contain the same delayed Ornstein-Uhlenbeck process.

The squared chord distance isHence the restricted Gaussian process has covarianceIt depends only on $t-t^{\prime }$, so it is stationary. It is periodic in either argument with period $T=2\pi /\omega$, and Gaussian-process realizations inherit that period almost surely because $f(t+T)-f(t)$ has zero variance.