Use the part-preserving bipartite four-cycle density convention
Here is the indicator function of adjacency, and repeated labels are allowed. This is the homomorphism density convention appropriate to the four-cycle norm. We derive the needed counting and operator norm facts directly.
Let be the -by- biadjacency matrix, let be the all-ones vector, let , and put . The degree hypothesis says
so the centered matrix annihilates on both sides. In particular, , and
The two terms have zero products in either order. Expanding the matrix trace directly counts four adjacency factors:
Since and , we obtain
The given upper bound therefore implies
This computation proves the centered four-cycle identity for a biregular graph, rather than assuming a four-cycle norm identity.
The matrix is a positive semidefinite matrix. By the spectral theorem for real symmetric matrices it has nonnegative eigenvalues , and
Moreover, for the Euclidean norm,
Indeed , whose maximum on unit vectors is ; the Cauchy-Schwarz inequality characterization also shows that and have the same operator norm. Taking fourth roots thus gives
where the error parameter is nonnegative, as in the claimed bound.
Finally use the indicator vectors , . Their Euclidean norms are and , so the Cauchy-Schwarz inequality and the operator norm estimate yield
This is the required discrepancy bound. It also holds when either subset is empty or . In fact, because annihilates constants, replacing by gives the stronger bound . This is the spectral discrepancy bound for a biregular graph.
For the reality of the eigenvalues, one may initially allow a nonzero complex eigenvector with . Since is real symmetric, it is a Hermitian matrix, and
is real: the numerator equals its own complex conjugate and the denominator is positive. Once is real, a nonzero real or imaginary part of supplies a real eigenvector. We work with the real unit eigenvectors of the problem.
For two such eigenvectors, symmetry gives
Distinct eigenvalues therefore give zero dot product. Together with the prescribed unit lengths,
These vectors form an orthonormal basis of .
Expand a real unit vector as . Then , and its Rayleigh quotient is
Moreover,
Every coefficient on the right is strictly positive, so equality forces for . Equality occurs precisely along the top eigendirection:
To obtain a unit vector in using just , let and denote their first coordinates. If , take
Then and . Hence has first coordinate zero and unit length. If , simply take . In either case the Rayleigh quotient of the constructed vector is
The continuous quadratic form attains its maximum on the compact space , so
For , deletion of the first row and column gives , while . The spectral theorem for real symmetric matrices applied to identifies its largest eigenvalue with this maximum:
Combining the upper bound for every unit vector with the constructed lower bound proves
This is largest-eigenvalue interlacing for a principal submatrix: restricting a quadratic form to a coordinate hyperplane cannot exceed its old maximum, but leaves some direction in the span of the top two eigenvectors.
Use the spectral theorem for real symmetric matrices to write , with orthogonal and , . Set and generate
Then . The characteristic function is
which identifies the multivariate normal distribution . Equivalently, every linear combination is normal, with the required mean and variance. A Cholesky decomposition can be used for a positive-definite matrix. The spectral construction also works for singular , producing a degenerate Gaussian distribution on an affine subspace; in that case there need not be a density on all of .
The quadratic form has symmetric matrix
whose characteristic polynomial is . All its eigenvalues are positive, so is positive definite and is an ellipsoid. The spectral theorem for real symmetric matrices gives semi-axis lengths
The nearest points lie in the eigendirection of the largest eigenvalue. Since is an -eigenvector, they are
Use the spectral theorem for real symmetric matrices to write , with orthogonal and diagonal and nonnegative. For any sign vector, , so
This is an exact identity for every sign choice. Put . When , the Rademacher rounding for a semidefinite relaxation gives
Hence is feasible. In the finite spanning case from part (a), makes , and the spanning condition then gives . If a nonzero instead had , it would be an unbounded feasible direction rather than a vector to divide by zero.
Let be the tridiagonal matrix with diagonal two and adjacent entries minus one, incorporating zero Dirichlet boundary conditions. The actual PDF scheme is
An orthonormal basis of eigenvectors of has components , and substitution gives
They form an orthonormal basis by the spectral theorem for real symmetric matrices. Hence is invertible for every , and the amplification eigenvalues are . For the discrete L2 norm ,
This bound has constant one independent of time step and mesh, which is the required stability of a numerical method. Therefore every positive Courant number is stable. If a forcing or local-defect sequence is added, the discrete variation-of-constants formula gives a bound by the initial norm plus the sum of those perturbation norms; this also justifies the global error conclusion in part (a).
Under an orthogonal change of coordinates represented by , the component matrix changes by
Thus is similar to , so its characteristic polynomial, eigenvalues, and their algebraic multiplicities are unchanged.
If is symmetric, the spectral theorem for real symmetric matrices makes it diagonalizable, so algebraic and geometric multiplicities agree. Consequently the dimensions of its eigenspaces are also independent of the coordinate frame.
For the stated tensor, choose a unit eigenvector belonging to the simple eigenvalue . Its perpendicular plane is the eigenspace of , so
Hence
Let be the adjacency matrix of a graph. Since is a regular graph of degree , every row of sums to . Therefore, for ,
so is a graph eigenvalue.
If , then the quadratic form of the graph Laplacian gives
Every summand is nonnegative, so along every edge. The graph is connected, hence all coordinates of are equal and is a scalar multiple of . Thus the -eigenspace is one-dimensional. Since is a symmetric matrix, the spectral theorem for real symmetric matrices makes it diagonalizable, so the algebraic multiplicity is also one. This proves the constant eigenvector of a regular graph result.
Suppose the strongly regular graph has vertices, degree , common neighbours for every adjacent pair and common neighbours for every nonadjacent pair. By the walk count from powers of an adjacency matrix, counts common neighbours of and , while . Hence the adjacency-matrix relation for a strongly regular graph is
where is the all-ones matrix.
The spectral theorem for real symmetric matrices gives an orthonormal eigenbasis for . Part (a) says that the constant direction is the one-dimensional -eigenspace; every other eigenvector belongs to , so . If , the displayed identity gives
There are at most two possible roots in addition to . Thus has at most three distinct eigenvalues.
Because every column of has sample mean zero, every column of also has sample mean zero. The spectral theorem for real symmetric matrices gives
Hence for the sample covariance is , while the sample variance of is
Thus the are pairwise uncorrelated sample principal components, ordered by decreasing sample variance.
By the spectral theorem for real symmetric matrices, write with orthogonal and . Part a and the tensor-product property of the Fourier transform give
Thus
Applying the pullback rule from part c to the orthogonal change of variables, for which , replaces by and by . Since determinant and signature are invariant under orthogonal conjugation,
Expanding the determinant gives . Solving each eigenvalue equation yields
Their pairwise inner products are , and . Their squared Euclidean norms are , respectively, so division by gives an orthonormal basis. This also follows from the spectral theorem for real symmetric matrices: symmetry makes eigenvectors for distinct eigenvalues orthogonal.
For either sign, the dot product with is
The remaining pair also satisfies
Hence the three eigenvectors are mutually orthogonal. This also follows from the spectral theorem for real symmetric matrices, since and its eigenvalues are distinct, but the displayed products verify the actual chosen vectors directly.
Choose two perpendicular unit directions in the eigenspace:
Their pairwise dot products vanish and their Euclidean norms are one. They form an orthonormal basis of eigenvectors, with eigenvalues , as guaranteed for a real symmetric matrix by the spectral theorem for real symmetric matrices.
The pairwise dot products of the normalized vectors are
Their squared norms are respectively , , and , all one. Hence they form an orthonormal basis. This also follows from the spectral theorem for real symmetric matrices, because is symmetric and its three eigenvalues are distinct.
By the spectral theorem for real symmetric matrices, choose the normalized eigenvectors to form an orthonormal basis. With , orthonormality gives and . For every vector , expand and apply :
Since this equality holds for every , the spectral decomposition is
Each is the rank-one orthogonal projection onto its eigenvector direction. In a degenerate eigenspace, normalized eigenvectors must be chosen mutually orthogonal; normalization alone does not guarantee this resolution of the identity.
The spectral theorem for real symmetric matrices supplies an orthonormal basis of real eigenvectors. Put these in the columns of , so and . In the question's coordinate convention choose and . Since ,
Reorder the eigenvalues with positive, negative and zero. For nonzero eigenvalues set , and retain the remaining coordinates as . This is an invertible real transformation , where is the diagonal rescaling; the combined matrix gives . It gives
with the zero coordinates absent. The general allowed range is . The printed restriction is not valid for every real symmetric matrix: a negative-definite or zero form has no positive square. The invariance of these counts is Sylvester's law of inertia.
For the specified form, a particularly simple change is
Subtracting the first two squared coordinates verifies the identity directly, and proves invertibility. Thus , and there is one zero direction. The final rescaling need not remain orthogonal.
Choose real, as permitted by the spectral theorem for real symmetric matrices. Matrix powers commute, so
The proposed partner is nonzero because
Also the real scalar equals its own transpose , so it vanishes. Therefore and are orthogonal eigenvectors with the same eigenvalue. This is the basic pairing in the paired eigenspaces of the square of a real skew-symmetric matrix.