Periodic Functions & Orthogonal Functions — Question 10

PDF ↗

Question 10

Let u(x)=sin⁡xu(x)=\sin x and v(x)=sin⁡5xv(x)=\sin 5x on [0,2π][0,2\pi]. Compare the continuous average inner product with its equally spaced sampling version: ⟨f,g⟩=12π∫02πfgdx,⟨f,g⟩N=1N∑j=0N−1f(2πj/N)g(2πj/N).\langle f,g\rangle=\frac 1{2\pi}\int_0^{2\pi}fg\,dx,\qquad \langle f,g\rangle_N=\frac 1N\sum_{j=0}^{N-1}f(2\pi j/N)g(2\pi j/N). The sampled form may fail to be positive definite on functions even though it is positive definite on the vectors of sampled values.

Tasks

  1. Compute the continuous Gram matrix of u,vu,v; a Gram matrix has entries consisting of all pairwise inner products.

  2. For N=4N=4, list both sample vectors and compute the sampled Gram matrix. Exhibit a nonzero continuous function in their span with zero sampled norm, and sketch the curves with the common sampled points.

  3. For integer mm, derive Sm=N−1∑j=0N−1cos⁡(2πmj/N)S_m=N^{-1}\sum_{j=0}^{N-1}\cos(2\pi mj/N) using a finite geometric sum. Express every entry of the sampled Gram matrix in terms of S2,S4,S6,S10S_2,S_4,S_6,S_{10}.

  4. Find the smallest integer N≥3N\geq 3 for which the sampled Gram matrix equals the continuous one. Check every smaller candidate and explain why exactness for this pair does not imply exactness for all continuous functions.

Original worksheet page 1: question and worked solution for 8-3-010
Show solutionHide solution

Question 10 – Solution

Strategy. Distinguish continuous orthogonality from orthogonality of a finite set of sampled values.

Step 1: Compute the continuous matrix. Product-to-sum gives ⟨u,v⟩=0\langle u,v\rangle=0, while each sine has average square 1/21/2. Hence G=12I2\boxed{G=\tfrac 12 I_2}.

Step 2: Identify the four-point alias. At 0,π/2,π,3π/20,\pi/2,\pi,3\pi/2, both functions give (0,1,0,−1)(0,1,0,-1). Thus G4=12(1111).\boxed{G_4=\frac 12\begin{pmatrix}1&1\\1&1\end{pmatrix}.} The nonzero function u−vu-v vanishes at every sample, so its sampled norm is zero, although its continuous squared norm is one. The samples have lost the distinction.

Step 3: Derive the exact sampling formulas. Put z=e2πim/Nz=e^{2\pi im/N}. If NN divides mm, every summand equals one. Otherwise z≠1z\ne 1 and ∑j=0N−1zj=(1−zN)/(1−z)=0\sum_{j=0}^{N-1}z^j=(1-z^N)/(1-z)=0. Taking real parts yields Sm=1S_m=1 when N∣mN\mid m, and Sm=0S_m=0 otherwise. Product-to-sum now gives GN=12(1−S2S4−S6S4−S61−S10).\boxed{G_N=\frac 12\begin{pmatrix}1-S_2&S_4-S_6\\S_4-S_6&1-S_{10}\end{pmatrix}.}

Step 4: Find the smallest exact sample count. For N=3N=3 the off-diagonal entry is −1/2-1/2; for N=4N=4 it is 1/21/2. For N=5N=5, the second diagonal entry is zero. For N=6N=6 the off-diagonal entry is −1/2-1/2. None matches GG. For N=7N=7, none of 2,4,6,102,4,6,10 is divisible by seven, so G7=G\boxed{G_7=G} and seven is the smallest allowed count. This is exact only for the tested products: for example, sin⁡7x\sin 7x vanishes at all seven sample points but has continuous average square 1/21/2.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 8-3-010

Original worksheet layout. Use Enlarge or open the PDF for a closer view.