Misplaced Pages

Rudin–Shapiro sequence

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Rudin-Shapiro sequence)

In mathematics, the Rudin–Shapiro sequence, also known as the Golay–Rudin–Shapiro sequence, is an infinite 2-automatic sequence named after Marcel Golay, Harold S. Shapiro, and Walter Rudin who investigated its properties.

Definition

Each term of the Rudin–Shapiro sequence is either 1 {\displaystyle 1} or − 1 {\displaystyle -1} . If the binary expansion of n {\displaystyle n} is given by

n = ∑ k ≥ 0 ϵ k ( n ) 2 k , {\displaystyle n=\sum _{k\geq 0}\epsilon _{k}(n)2^{k},}

then let

u n = ∑ k ≥ 0 ϵ k ( n ) ϵ k + 1 ( n ) . {\displaystyle u_{n}=\sum _{k\geq 0}\epsilon _{k}(n)\epsilon _{k+1}(n).}

(So u n {\displaystyle u_{n}} is the number of times the block 11 appears in the binary expansion of n {\displaystyle n} .)

The Rudin–Shapiro sequence ( r n ) n ≥ 0 {\displaystyle (r_{n})_{n\geq 0}} is then defined by

r n = ( − 1 ) u n . {\displaystyle r_{n}=(-1)^{u_{n}}.}

Thus r n = 1 {\displaystyle r_{n}=1} if u n {\displaystyle u_{n}} is even and r n = − 1 {\displaystyle r_{n}=-1} if u n {\displaystyle u_{n}} is odd.

The sequence u n {\displaystyle u_{n}} is known as the complete Rudin–Shapiro sequence, and starting at n = 0 {\displaystyle n=0} , its first few terms are:

0, 0, 0, 1, 0, 0, 1, 2, 0, 0, 0, 1, 1, 1, 2, 3, ... (sequence A014081 in the OEIS)

and the corresponding terms r n {\displaystyle r_{n}} of the Rudin–Shapiro sequence are:

+1, +1, +1, −1, +1, +1, −1, +1, +1, +1, +1, −1, −1, −1, +1, −1, ... (sequence A020985 in the OEIS)

For example, u 6 = 1 {\displaystyle u_{6}=1} and r 6 = − 1 {\displaystyle r_{6}=-1} because the binary representation of 6 is 110, which contains one occurrence of 11; whereas u 7 = 2 {\displaystyle u_{7}=2} and r 7 = 1 {\displaystyle r_{7}=1} because the binary representation of 7 is 111, which contains two (overlapping) occurrences of 11.

Historical motivation

The Rudin–Shapiro sequence was introduced independently by Golay, Rudin, and Shapiro. The following is a description of Rudin's motivation. In Fourier analysis, one is often concerned with the L 2 {\displaystyle L^{2}} norm of a measurable function f : [ 0 , 2 π ) → [ 0 , 2 π ) {\displaystyle f\colon [0,2\pi )\to [0,2\pi )} . This norm is defined by

| | f | | 2 = ( 1 2 π ∫ 0 2 π | f ( t ) | 2 d t ) 1 / 2 . {\displaystyle ||f||_{2}=\left({\frac {1}{2\pi }}\int _{0}^{2\pi }|f(t)|^{2}\,\mathrm {d} t\right)^{1/2}.}

One can prove that for any sequence ( a n ) n ≥ 0 {\displaystyle (a_{n})_{n\geq 0}} with each a n {\displaystyle a_{n}} in { 1 , − 1 } {\displaystyle \{1,-1\}} ,

sup x ∈ R | ∑ 0 ≤ n < N a n e i n x | ≥ | | ∑ 0 ≤ n < N a n e i n x | | 2 = N . {\displaystyle \sup _{x\in \mathbb {R} }\left|\sum _{0\leq n<N}a_{n}e^{inx}\right|\geq \left|\left|\sum _{0\leq n<N}a_{n}e^{inx}\right|\right|_{2}={\sqrt {N}}.}

Moreover, for almost every sequence ( a n ) n ≥ 0 {\displaystyle (a_{n})_{n\geq 0}} with each a n {\displaystyle a_{n}} is in { − 1 , 1 } {\displaystyle \{-1,1\}} ,

sup x ∈ R | ∑ 0 ≤ n < N a n e i n x | = O ( N log ⁡ N ) . {\displaystyle \sup _{x\in \mathbb {R} }\left|\sum _{0\leq n<N}a_{n}e^{inx}\right|=O({\sqrt {N\log N}}).}

However, the Rudin–Shapiro sequence ( r n ) n ≥ 0 {\displaystyle (r_{n})_{n\geq 0}} satisfies a tighter bound: there exists a constant C > 0 {\displaystyle C>0} such that

sup x ∈ R | ∑ 0 ≤ n < N r n e i n x | ≤ C N . {\displaystyle \sup _{x\in \mathbb {R} }\left|\sum _{0\leq n<N}r_{n}e^{inx}\right|\leq C{\sqrt {N}}.}

It is conjectured that one can take C = 6 {\displaystyle C={\sqrt {6}}} , but while it is known that C ≥ 6 {\displaystyle C\geq {\sqrt {6}}} , the best published upper bound is currently C ≤ ( 2 + 2 ) 3 / 5 {\displaystyle C\leq (2+{\sqrt {2}}){\sqrt {3/5}}} . Let P n {\displaystyle P_{n}} be the n-th Shapiro polynomial. Then, when N = 2 n − 1 {\displaystyle N=2^{n}-1} , the above inequality gives a bound on sup x ∈ R | P n ( e i x ) | {\displaystyle \sup _{x\in \mathbb {R} }|P_{n}(e^{ix})|} . More recently, bounds have also been given for the magnitude of the coefficients of | P n ( z ) | 2 {\displaystyle |P_{n}(z)|^{2}} where | z | = 1 {\displaystyle |z|=1} .

Shapiro arrived at the sequence because the polynomials

P n ( z ) = ∑ i = 0 2 n − 1 r i z i {\displaystyle P_{n}(z)=\sum _{i=0}^{2^{n}-1}r_{i}z^{i}}

where ( r i ) i ≥ 0 {\displaystyle (r_{i})_{i\geq 0}} is the Rudin–Shapiro sequence, have absolute value bounded on the complex unit circle by 2 n + 1 2 {\displaystyle 2^{\frac {n+1}{2}}} . This is discussed in more detail in the article on Shapiro polynomials. Golay's motivation was similar, although he was concerned with applications to spectroscopy and published in an optics journal.

Properties

The Rudin–Shapiro sequence can be generated by a 4-state automaton accepting binary representations of non-negative integers as input. The sequence is therefore 2-automatic, so by Cobham's little theorem there exists a 2-uniform morphism φ {\displaystyle \varphi } with fixed point w {\displaystyle w} and a coding τ {\displaystyle \tau } such that r = τ ( w ) {\displaystyle r=\tau (w)} , where r {\displaystyle r} is the Rudin–Shapiro sequence. However, the Rudin–Shapiro sequence cannot be expressed as the fixed point of some uniform morphism alone.

There is a recursive definition

{ r 2 n = r n r 2 n + 1 = ( − 1 ) n r n {\displaystyle {\begin{cases}r_{2n}&=r_{n}\\r_{2n+1}&=(-1)^{n}r_{n}\end{cases}}}

The values of the terms rn and un in the Rudin–Shapiro sequence can be found recursively as follows. If n = m·2 where m is odd then

u n = { u ( m − 1 ) / 4 if  m ≡ 1 ( mod 4 ) u ( m − 1 ) / 2 + 1 if  m ≡ 3 ( mod 4 ) {\displaystyle u_{n}={\begin{cases}u_{(m-1)/4}&{\text{if }}m\equiv 1{\pmod {4}}\\u_{(m-1)/2}+1&{\text{if }}m\equiv 3{\pmod {4}}\end{cases}}}
r n = { r ( m − 1 ) / 4 if  m ≡ 1 ( mod 4 ) − r ( m − 1 ) / 2 if  m ≡ 3 ( mod 4 ) {\displaystyle r_{n}={\begin{cases}r_{(m-1)/4}&{\text{if }}m\equiv 1{\pmod {4}}\\-r_{(m-1)/2}&{\text{if }}m\equiv 3{\pmod {4}}\end{cases}}}

Thus u108 = u13 + 1 = u3 + 1 = u1 + 2 = u0 + 2 = 2, which can be verified by observing that the binary representation of 108, which is 1101100, contains two sub-strings 11. And so r108 = (−1) = +1.

A 2-uniform morphism φ {\displaystyle \varphi } that requires a coding τ {\displaystyle \tau } to generate the Rudin-Shapiro sequence is the following: φ : a → a b b → a c c → d b d → d c {\displaystyle {\begin{aligned}\varphi :a&\to ab\\b&\to ac\\c&\to db\\d&\to dc\end{aligned}}} τ : a → 1 b → 1 c → − 1 d → − 1 {\displaystyle {\begin{aligned}\tau :a&\to 1\\b&\to 1\\c&\to -1\\d&\to -1\end{aligned}}}

The Rudin–Shapiro word +1 +1 +1 −1 +1 +1 −1 +1 +1 +1 +1 −1 −1 −1 +1 −1 ..., which is created by concatenating the terms of the Rudin–Shapiro sequence, is a fixed point of the morphism or string substitution rules

+1 +1 → +1 +1 +1 −1
+1 −1 → +1 +1 −1 +1
−1 +1 → −1 −1 +1 −1
−1 −1 → −1 −1 −1 +1

as follows:

+1 +1 → +1 +1 +1 −1 → +1 +1 +1 −1 +1 +1 −1 +1 → +1 +1 +1 −1 +1 +1 −1 +1 +1 +1 +1 −1 −1 −1 +1 −1 ...

It can be seen from the morphism rules that the Rudin–Shapiro string contains at most four consecutive +1s and at most four consecutive −1s.

The sequence of partial sums of the Rudin–Shapiro sequence, defined by

s n = ∑ k = 0 n r k , {\displaystyle s_{n}=\sum _{k=0}^{n}r_{k}\,,}

with values

1, 2, 3, 2, 3, 4, 3, 4, 5, 6, 7, 6, 5, 4, 5, 4, ... (sequence A020986 in the OEIS)

can be shown to satisfy the inequality

3 5 n < s n < 6 n  for  n ≥ 1 . {\displaystyle {\sqrt {{\frac {3}{5}}n}}<s_{n}<{\sqrt {6n}}{\text{ for }}n\geq 1\,.}

Let ( s n ) n ≥ 0 {\displaystyle (s_{n})_{n\geq 0}} denote the Rudin–Shapiro sequence on { 0 , 1 } {\displaystyle \{0,1\}} , in which case s n {\displaystyle s_{n}} is the number, modulo 2, of occurrences (possibly overlapping) of the block 11 {\displaystyle 11} in the base-2 expansion of n {\displaystyle n} . Then the generating function

S ( X ) = ∑ n ≥ 0 s n X n {\displaystyle S(X)=\sum _{n\geq 0}s_{n}X^{n}}

satisfies

( 1 + X ) 5 S ( X ) 2 + ( 1 + X ) 4 S ( X ) + X 3 = 0 , {\displaystyle (1+X)^{5}S(X)^{2}+(1+X)^{4}S(X)+X^{3}=0,}

making it algebraic as a formal power series over F 2 ( X ) {\displaystyle \mathbb {F} _{2}(X)} . The algebraicity of S ( X ) {\displaystyle S(X)} over F 2 ( X ) {\displaystyle \mathbb {F} _{2}(X)} follows from the 2-automaticity of ( s n ) n ≥ 0 {\displaystyle (s_{n})_{n\geq 0}} by Christol's theorem.

The Rudin–Shapiro sequence along squares ( r n 2 ) n ≥ 0 {\displaystyle (r_{n^{2}})_{n\geq 0}} is normal.

The complete Rudin–Shapiro sequence satisfies the following uniform distribution result. If x ∈ R ∖ Z {\displaystyle x\in \mathbb {R} \setminus \mathbb {Z} } , then there exists α = α ( x ) ∈ ( 0 , 1 ) {\displaystyle \alpha =\alpha (x)\in (0,1)} such that

∑ n < N exp ⁡ ( 2 π i x u n ) = O ( N α ) {\displaystyle \sum _{n<N}\exp(2\pi ixu_{n})=O(N^{\alpha })}

which implies that ( x u n ) n ≥ 0 {\displaystyle (xu_{n})_{n\geq 0}} is uniformly distributed modulo 1 {\displaystyle 1} for all irrationals x {\displaystyle x} .

Relationship with one-dimensional Ising model

Let the binary expansion of n be given by

n = ∑ k ≥ 0 ϵ k ( n ) 2 k {\displaystyle n=\sum _{k\geq 0}\epsilon _{k}(n)2^{k}}

where ϵ k ( n ) ∈ { 0 , 1 } {\displaystyle \epsilon _{k}(n)\in \{0,1\}} . Recall that the complete Rudin–Shapiro sequence is defined by

u ( n ) = ∑ k ≥ 0 ϵ k ( n ) ϵ k + 1 ( n ) . {\displaystyle u(n)=\sum _{k\geq 0}\epsilon _{k}(n)\epsilon _{k+1}(n).}

Let

ϵ ~ k ( n ) = { ϵ k ( n ) if  k ≤ N − 1 , ϵ 0 ( n ) if  k = N . {\displaystyle {\tilde {\epsilon }}_{k}(n)={\begin{cases}\epsilon _{k}(n)&{\text{if }}k\leq N-1,\\\epsilon _{0}(n)&{\text{if }}k=N.\end{cases}}}

Then let

u ( n , N ) = ∑ 0 ≤ k < N ϵ ~ k ( n ) ϵ ~ k + 1 ( n ) . {\displaystyle u(n,N)=\sum _{0\leq k<N}{\tilde {\epsilon }}_{k}(n){\tilde {\epsilon }}_{k+1}(n).}

Finally, let

S ( N , x ) = ∑ 0 ≤ n < 2 N exp ⁡ ( 2 π i x u ( n , N ) ) . {\displaystyle S(N,x)=\sum _{0\leq n<2^{N}}\exp(2\pi ixu(n,N)).}

Recall that the partition function of the one-dimensional Ising model can be defined as follows. Fix N ≥ 1 {\displaystyle N\geq 1} representing the number of sites, and fix constants J > 0 {\displaystyle J>0} and H > 0 {\displaystyle H>0} representing the coupling constant and external field strength, respectively. Choose a sequence of weights η = ( η 0 , … , η N − 1 ) {\displaystyle \eta =(\eta _{0},\dots ,\eta _{N-1})} with each η i ∈ { − 1 , 1 } {\displaystyle \eta _{i}\in \{-1,1\}} . For any sequence of spins σ = ( σ 0 , … , σ N − 1 ) {\displaystyle \sigma =(\sigma _{0},\dots ,\sigma _{N-1})} with each σ i ∈ { − 1 , 1 } {\displaystyle \sigma _{i}\in \{-1,1\}} , define its Hamiltonian by

H η ( σ ) = − J ∑ 0 ≤ k < N η k σ k σ k + 1 − H ∑ 0 ≤ k < N σ k . {\displaystyle H_{\eta }(\sigma )=-J\sum _{0\leq k<N}\eta _{k}\sigma _{k}\sigma _{k+1}-H\sum _{0\leq k<N}\sigma _{k}.}

Let T {\displaystyle T} be a constant representing the temperature, which is allowed to be an arbitrary non-zero complex number, and set β = 1 / ( k T ) {\displaystyle \beta =1/(kT)} where k {\displaystyle k} is the Boltzmann constant. The partition function is defined by

Z N ( η , J , H , β ) = ∑ σ ∈ { − 1 , 1 } N exp ⁡ ( − β H η ( σ ) ) . {\displaystyle Z_{N}(\eta ,J,H,\beta )=\sum _{\sigma \in \{-1,1\}^{N}}\exp(-\beta H_{\eta }(\sigma )).}

Then we have

S ( N , x ) = exp ⁡ ( π i N x 2 ) Z N ( 1 , 1 2 , − 1 , π i x ) {\displaystyle S(N,x)=\exp \left({\frac {\pi iNx}{2}}\right)Z_{N}\left(1,{\frac {1}{2}},-1,\pi ix\right)}

where the weight sequence η = ( η 0 , … , η N − 1 ) {\displaystyle \eta =(\eta _{0},\dots ,\eta _{N-1})} satisfies η i = 1 {\displaystyle \eta _{i}=1} for all i {\displaystyle i} .

See also

Notes

  1. ^ John Brillhart and Patrick Morton, winners of a 1997 Lester R. Ford Award (1996). "A Case Study in Mathematical Research: The Golay–Rudin–Shapiro Sequence". Amer. Math. Monthly. 103 (10): 854–869. doi:10.2307/2974610. JSTOR 2974610.{{cite journal}}: CS1 maint: numeric names: authors list (link)
  2. Weisstein, Eric W. "Rudin–Shapiro Sequence". MathWorld.
  3. ^ Pytheas Fogg (2002) p.42
  4. Everest et al (2003) p.234
  5. Golay, M.J.E. (1949). "Multi-slit spectrometry". Journal of the Optical Society of America. 39 (437–444): 437–444. doi:10.1364/JOSA.39.000437. PMID 18152021.
  6. Golay, M.J.E. (1951). "Static multislit spectrometry and its application to the panoramic display of infrared spectra". Journal of the Optical Society of America. 41 (7): 468–472. doi:10.1364/JOSA.41.000468. PMID 14851129.
  7. Rudin, W. (1959). "Some theorems on Fourier coefficients". Proceedings of the American Mathematical Society. 10 (6): 855–859. doi:10.1090/S0002-9939-1959-0116184-5.
  8. Shapiro, H.S. (1952). "Extremal problems for polynomials and power series". Master's Thesis, MIT.
  9. Salem, R.; Zygmund, A. (1954). "Some properties of trigonometric series whose terms have random signs". Acta Mathematica. 91: 245–301. doi:10.1007/BF02393433. S2CID 122999383.
  10. Allouche and Shallit (2003) p. 78–79
  11. Allouche and Shallit (2003) p. 122
  12. Brillhart, J.; Morton, P. (1978). "Über Summen von Rudin–Shapiroschen Koeffizienten". Illinois Journal of Mathematics. 22: 126–148. doi:10.1215/ijm/1256048841.
  13. Saffari, B. (1986). "Une fonction extrémale liée à la suite de Rudin–Shapiro". C. R. Acad. Sci. Paris. 303: 97–100.
  14. Allouche, J.-P.; Choi, S.; Denise, A.; Erdélyi, T.; Saffari, B. (2019). "Bounds on Autocorrelation Coefficients of Rudin-Shapiro Polynomials". Analysis Mathematica. 45 (4): 705–726. arXiv:1901.06832. doi:10.1007/s10476-019-0003-4. S2CID 119168430.
  15. Finite automata and arithmetic, Jean-Paul Allouche
  16. Allouche and Shallit (2003) p. 192
  17. Allouche and Shallit (2003) p. 352
  18. Müllner, C. (2018). "The Rudin–Shapiro sequence and similar sequences are normal along squares". Canadian Journal of Mathematics. 70 (5): 1096–1129. arXiv:1704.06472. doi:10.4153/CJM-2017-053-1. S2CID 125493369.
  19. Allouche and Shallit p. 462–464
  20. Allouche and Shallit (2003) p. 457–461

References

Category: