arxiv: v1 [cs.it] 5 Aug 2008

Transcription

1 Generalized Cross-correlation Properties of Chu Sequences arxiv: v1 [cs.it] 5 Aug 2008 Jae Won Kang, Younghoon Whang, Hyo Yol Park and Kwang Soon Kim, Senior Member, IEEE Department of Electrical and Electronic Engineering, Yonsei University, 134 Shinchon-dong, Seodaemun-gu, Seoul , Korea Tel: Fax: ks.kim@yonsei.ac.kr Corresponding Author Abstract In this paper, we analyze the cross-correlation properties for Chu sequences, which provide information on the distribution of the maximum magnitudes of the cross-correlation function. Furthermore, we can obtain the number of available sequences for a given maximum magnitude of the cross-correlation function and the sequence length. Chu sequences, cross-correlation function. Index Terms

2 I. INTRODUCTION IN general, it is desired to design a set of sequences with an impulsive autocorrelation function and a zero cross-correlation function for many practical applications. However, according to the Welch bound, the Sarwate bound, the Sidelnikov bound, the Massey bound and other bounds [1] [4], it was shown to be impossible to construct such an ideal set of sequences. Therefore, searching large families of sequences with good auto-correlation function and cross-correlation function properties has been one of the most interesting topics in sequence design. For evaluating the correlation properties, one good choice is to use the maximum sidelobe magnitude of the autocorrelation function and the maximum magnitude of the cross-correlation function, which are respectively denoted as ˆθ a and ˆθ c in this paper. Here, the following questions arise naturally: how many pairs of sequences are available for a given maximum values of ˆθ a and ˆθ c and what is the distribution of the magnitude of the cross-correlation function? Among well known good sequences are Kasami [5], Gold [5], Chu [7] [9] and complex four-phase [10] sequences. For Kasami and Gold sequences, it was shown that there are N + 1 sequences satisfying ˆθ a = 1 and ˆθ c = 1 + 2/N [5][6][11][12], where N is the sequence length. For four-phase sequences, the number of sequences satisfying ˆθ a = 1 + N and ˆθ c = 1 + N is N + 2 [10][13][14]. On the other hand, the autocorrelation function of Chu sequences is known to be zero except at the lag of an integer multiple of the sequence length [7] [9][15] [17]. A set of Chu sequences with length N is defined as C = {a r 0 < r < N, gcdn, r) = 1}, where the kth element of a r, a r k), is defined as ) jπ rk2, N even, N a r k) = ) rkk + 1) jπ, N odd. N The periodic autocorrelation function with lag τ, θ r τ), of the sequence a r is defined as 1) θ r τ) = N τ 1 N τ 1 a r k)a r k + τ)+ k=n τ a r k)a r k + τ N). 2)

3 In [8], it was shown that the periodic autocorrelation function of Chu sequences satisfies N, τ mod N = 0, θ r τ) = 0, τ mod N 0. 3) Let a r and a s be any two Chu sequences with length N. Then, the cross-correlation function θ r,s τ) of a r and a s with lag τ is defined as θ r,s τ) = = N τ 1 N 1 a r k)a sk + τ) + N τ 1 k=n τ a r k)a sk + τ N) a r k)a s k + τ), 4) where the last equality comes from the fact that a r k + d) = a r k + d + N) for an arbitrary integer d [8]. In [1] and [18], it was shown that the maximum magnitude of the cross-correlation function ˆθ c can be lower-bounded as a function of the sequence length and the maximum magnitude of the autocorrelation function, ˆθ a. By using this lower-bound, the optimum correlation properties of a set of sequences can be defined and it follows that the lower bound of ˆθ c is equal to N when ˆθ a equals zero. Certain pairs of Chu sequences, a r and a s, meet this lower-bound when gcdr s, N) = 1. However, in order to obtain more Chu sequences with relatively low cross-correlation values, we need to investigate more general cross-correlation properties. In this paper, we derive general properties for cross-correlation function of Chu sequences. Using the derived properties, we can obtain the magnitude distribution of the cross-correlation function. Here, the maximum magnitude denotes the maximum magnitude value of the cross-correlation function of two given Chu sequences among all possible lags and its distribution is taken over all possible pairs of Chu sequences. In addition, the number of available sequences can be obtained for a given value of ˆθ c and the given sequence length. The remaining of this paper is organized as follows. In Section II, the magnitude of cross-correlation function of Chu sequences are described. In Section III, the distribution of the maximum magnitude of the cross-correlation function and the number of available Chu sequences for given maximum crosscorrelation value and the sequence length are investigated. Finally, Section IV concludes this paper.

4 II. CHARACTERISTIC OF THE CROSS-CORRELATION FUNCTION OF CHU SEQUENCES In order to investigate the cross-correlation function of Chu sequences in detail, we need to find what are the possible values that the cross-correlation function of Chu sequences can take, which are given in the following theorem. Definition 1: Let r and s be positive integers satisfying 0 < r, s < N, gcdn, r) = 1 and gcdn, s) = 1. Also define g r,s = gcdn, r s), u r,s = N/g r,s and v r,s = r s)/g r,s. Then u r,s is relatively prime with v r,s. Also, for a given lag τ, we can rewrite it as τ = i τ g r,s +d τ, where i τ = τ/g r,s and d τ = τ i τ g r,s. Theorem 1: The magnitude of the cross-correlation function θ r,s τ), θ r,s τ), is given as Ng r,s δ K d τ ), N and u r,s v r,s even, or N odd, θ r,s τ) = Ng r,s δ K d τ g ) r,s, N even and u r,s v r,s odd, 2 0, otherwise, where δ K ) is the Kroneker delta function. To prove Theorem 1, the following lemmas are useful. Lemma 1 [7][8]: The hth primitive root of unity ξ h can be defined as ξ h = j2π u h), where u is any integer relatively prime to h. Then, for any integer v, 0 < v h 1, h 1 ξ ±vk h = 0, ξ h 1.

5 Lemma 2: The squared magnitude of the cross-correlation function is given as θ r,s τ) 2 = g r,s 1 u r,s g r,s g r,s 1 u r,s g r,s The proof of Lemma 2 is given in Appendix A. 1) ur,svr,sm2 j2π smd ), N even, g r,s 1) vrmur,s+1) j2π smd ), N odd. g r,s Now, the proof of Theorem 1 is given as follows. Proof of Theorem 1: First, consider the case when N and u r,s v r,s are even. Then, from Lemma 2, we obtain g r,s 1 θ r,s τ) 2 = u r,s g r,s j2π smd τ g r,s ). 5) If d τ = 0, θ r,s i τ u r,s ) 2 = u r,s g 2 r,s = Ng r,s. When, d τ 0, since s is relatively prime with g r,s, θ r,s i τ g r,s + d τ ) 2 = 0. Now, consider the case when N is even and u r,s v r,s is odd. When d τ = g r,s /2, we obtain from Lemma 2 that θ r,s i τ g r,s + g r,s 2 ) 2 = ur,sg r,s g r,s 1 { j2π mu r,sv r,s m + s) 2 )}. 6) We know that s is odd because s is relatively prime with N. If m is odd, u r,s v r,s m is odd and u r,s v r,s m+s is even. On the other hand, if m is even, u r,s v r,s m is also even. Thus, mu r,s v r,s m + s) is always even and it shows that θ r,s i τ g r,s + g r,s /2) 2 = u r,s g 2 r,s = Ng r,s. When d τ g r,s /2, from Lemma 2, we can rewrite θ r,s τ) 2 as θ r,si τ g r,s + g r,s 2 + d τ ) 2 g r,s 1 = ur,s g r,s g r,s 1 = u r,s g r,s { j2π j2π smd τ g r,s mu r,sv r,s m + s) 2 ) + smd τ g r,s )} 7) = 0,

6 where d τ = d τ g r,s /2 and the last equality comes from the fact that s is relatively prime with g r,s. Finally, consider the case where N is odd. Then g r,s and u r,s should be odd. Then, from Lemma 2, we obtain which concludes the proof. g r,s 1 θ r,s i τ g r,s + d τ ) 2 = u r,s g r,s = Ng r,s δ k d τ ), j2π smd τ g r,s ) 8) III. DISTRIBUTION OF THE MAXIMUM MAGNITUDES OF THE CROSS-CORRELATION FUNCTION A. The uniform property Theorem 1 tells us that the characteristic of the cross-correlation function of two Chu sequences, a r and a s, depends only on g r,s = gcdr s, N). For example, when g r,s = 1, ˆθ c meets the lower bound of N. On the other hand, when g r,s = N, ˆθ c becomes the largest value of N. However, it has not yet been investigated how many sequences are available for a given values of ˆθ c and the sequence length. To answer the question, it is required to investigate the distribution of the maximum magnitude values of the cross-correlation function. Definition 2: Any given integer N can be represented as N = k i=1 p c i i, where p i denotes the ith smallest prime factor of N. Let us define µ N = {n 0 < n < N, gcdn, N) = 1} as the index set of Chu sequences of length N. Also, for a given integer c, define the following sets and function as follows. U N,c = {n c 0 n < N} R N,c = {n c n µ N } D N,c = {n c 0 n < N and n / µ N } P m N,c = {np m c 0 n < N/p m } G N,x S) = {n n S and gcdn, N) = x} for a given integer set S. A : The cardinality of a set A.

7 From Theorem 1, we can see that the maximum magnitude of the cross-correlation function between a r and a s is ˆθ r,s = max θ r,s τ) = g r,s N. Thus, for given N, s µ N and x, it is easily seen that τ G N,x R N,s ) = {r s g r,s = gcdr s, N) = x and r µ N } is the set of differences between s and all Chu sequence indices whose maximum squared magnitude of the cross-correlation function with a s is equal to x. Then, G N,x R N,s ) is the number of available Chu sequences satisfying ˆθ r,s 2 /N = x. Then, the main result of this subsection is given in the following theorem. Theorem 2: Let 1 s s N be two different integers relatively prime with N. Then, G N,x R N,s ) = G N,x R N,s ). Theorem 2 indicates that the distribution of the maximum magnitudes of the cross-correlation function for a given Chu sequence set can be obtained by fixing one sequence arbitrarily and examining the cross-correlation functions with the other sequences. The following Lemmas 3 5 are useful the proof of Theorem 2. Lemma 3: For any two different integers c and c, G N,x U N,c ) = G N,x U N,c ). Proof: It has been proved that gcdc + mn, N) = gcdc, N) [19][20]. Then, it is easily seen that { gcd c + 1, N), gcd c + 2, N),, gcd c + N, N) } = { gcd1, N), gcd2, N),, gcdn, N) } for any integer c. Therefore, G N,x U N,c ) = G N,x U N,c ). Lemma 4: Let a and b be positive integers satisfying gcda, b) = 1. Also, for an arbitrary positive integer m, define C = {na c k n < k }, where k is an arbitrary integer and k = k + mb. Then, C contains exactly m integer multiples of b.

8 Proof: For the ith element c i = k i a c of C, we can represent it as c i = qk i )b + ek i ), where qk i ) = c i /b and eb i ) = c i mod b. Note that such a pair of qk) and ek) is unique for a given c i [20]. Let d ij = c i c j. Then d ij = k i k j )a = {qk i ) qk j )}b + ek i ) ek j ). Thus, ek i ) = ek j ) implies that k i k j ) is an integer multiple of b and vice versa because a is relatively prime with b. Now, consider the partition {C r, 0 r < m}, where C r = {k r,i a + c k + rb k r,i < k + r + 1)b}. Then, each C r contain exactly one element that is an integer multiple of b since ek r,i ), k +rb k r,i < k +r +1)b, are all distinct and 0 ek r,i ) < b, which concludes the proof. Lemma 5: Let 1 s s < N be two different integers relatively prime with N. Then, if x is a divisor of N, G N,x D N,s ) = G N,x D N,s ). Proof: For a given N = k i=1 p c i i, since D N,s = k i=1 PN,s i from Definition 2, G N,x D N,s ) = ki=1 G N,x PN,s) i. Then, GN,x D N,s ) can be rewritten as k ) k 1 G N,x D N,s ) = G N,x P i k N,s G N,x P i 1 N,s P i 2 N,s) + i=1 i 1 =1 i 2 =i k k 1 + 1) k 2 G N,x PN,s) im + 1)k 1 G k 9) N,x PN,s) i. i 1 =1 i 2 =i 1 +1 i k 1 =i k 2 +1 m=1 i=1 If gcdx, p m ) = 1, it is seen easily that there always exist N/p m x) integer multiples of x among the elements in P m N,s from Lemma 4. If gcdx, p m ) 1, p m should be a divisor of x since p m is a prime number. Thus, there is no integer multiple of x among the elements in P m N,s since x is relatively prime with s. Thus, G N,x PN,s m ) does not depend on s as long as s is relatively prime with N, i.e., G N,x P m N,s ) = G N,x P m N,s ). Let M be an arbitrary subset of {1,, k}. Also, let mi denote the ith element of M. Then, for a given index set M, define L M N,s = M i=1 P m i N,s = {nl m s 0 n < N/l m }, where l m = M i=1 p mi. Similarly, if gcdx, l m ) = 1, there always exist N/ l m x) integer multiples of x among elements in L M N,s from Lemma 4. If gcdx, l m ) 1, x should be an integer multiple of p mi for some m i M. Since s is relatively prime with all p mi, m i M, there is no integer multiple of x among the elements in L M N,s. Thus, G N,x L M N,s) does not depend on s as long as s is relatively prime with N.

9 Thus, from 9), G N,x D N,s ) = G N,x D N,s ). The proof of Theorem 2 is now given as follows. Proof of Theorem 2: When x is not a divisor of N, G N,x R N,s ) = regardless of s. When x is a divisor of N, from Lemmas 3 and 5, we have already seen that G N,x U N,s ) and G N,x D N,s ) does not depend on s as long as s is relatively prime with N and x is a divisor of N. Also, from Definition 2, since R N,s D N,s =, G N,x R N,s ) = G N,x U N,s ) G N,x D N,s ). Thus, G N,x R N,s ) does not depend on s as long as s is relatively prime with N. B. The distribution In this subsection, the distribution of the maximum magnitudes of the cross-correlation function is investigated. The main result of this subsection is given as follows. Theorem 3: For N = k i=1 p c i i, 0 < x N, any s relatively prime with N, G N,x R N,s ) is given as G N,x R N,s ) = i M x φ x i) j M k x Φ x j), 10) where M x, M k x, φ xi) and Φ x i) are defined as follows. When x is a divisor of N, M x denotes the set of indices of the prime factors of x and M k x = {1,, k} M x. Then, x can be represented as x = i M x p n ix) i, where 0 < n i x) c i. Also, φ x i) and Φ x i) are defined as φ x i) p c i n i x) i Φ x j) p c j j 2p c j 1 j. 1 δ K c i n i x))) p c i n i x) 1 i, When x is not a divisor of N, M x, φ x i) = 0 and Φ x i) = 0, so that G N,x R N,s ) = 0. The following Lemmas 6 10 are useful for the proof of Theorem 3.

10 Lemma 6 Euler function [20]): Let ϕn) be the number of positive integers that are relatively prime with N = k i=1 p c i i among {n 0 < n < N}. Then, ϕn) is given as ϕn) = N k i=1 ) 1 p 1 i. 11) Lemma 7: If x is not a divisor of N, G N,x R N,s ) = 0 since gcdr s, N) cannot be equal to x. Lemma 8: When N is a prime number, for s µ N and 0 < x N, G N,x R N,s ) is given as 1, if x = N, G N,x R N,s ) = N 2, if x = 1, 0, otherwise. 12) Proof: Since N is a prime number, gcdr s, N) = 1 when s r and gcdr s, N) = N when s = r. Thus, G N,x R N,s ) = 0 for 1 < x < N. Since µ N = R N,s = N 1 from Definition 2, G N,1 R N,s ) = N 2 and G N,N R N,s ) = 1. Lemma 9: For N = k i=1 p i, its divisor x and gcds, N) = 1, we can denote M x be the set of indices of the prime factors of x so that x = i M x p i and M k x = {1,, k} M x. Then, G N,x R N,s ) is given as G N,x R N,s ) = 1, if x = N, p i 2), i Mx k if x = p i N, i M x 0, otherwise. 13) Proof: When k = 1, 13) holds from Lemma 8. Suppose that 13) holds for any divisor x of N when k = K and let N = Np K+1. Then, the set of all divisors of N is given as {y = x or xp K+1 x {all divisors of N}}. Note that µ N = { n 0 < n < N, gcdn, N ) = 1} = { n + mn n µ N, 0

11 m < p K+1, gcdn + mn, p K+1 ) = 1}. Thus, R N,s = { n s n µ N } = p K+1 1 R m N,s, where R m N,s = { n + mn s n µ N, gcdn + mn, p K+1 ) = 1 }. When y = x for a divisor x of N, G N,yR m N,s) = { n + mn n G N,x R N,s ), gcdn + mn, p K+1 ) = 1, gcdn s + mn, p K+1 ) = 1 } since gcdn s + mn, N ) = x implies gcdn s, N) = x and gcdn s + mn, N ) = x if gcdn s, N) = x and gcdn s + mn, p K+1 ) = 1. Similarly, when y=xp K+1 for a divisor x of N, G N,yR m N,s ) = { n + mn n G N,xR N,s ), gcdn + mn, p K+1 ) = 1, gcdn s + mn, p K+1 ) = p K+1 }. Now, define B N n) = { n + mn 0 m < p K+1 }. Then, it is easily seen that p K+1 1 G N,yR N,s) = G N,yRN m,s) = n G N,x R N,s ) B N n) α s N,N x, y), 14) where α s N,N x, y) = { ζ ζ n G N,x R N,s ) B N n), gcdζ +s, p K+1 ) 1 or gcdζ, p K+1 ) y/x }. From Lemma 4, it is easily seen that α s N,N x, y) B Nn) = 2 when y = x since there is exactly one element in B N n) for each of the two conditions : gcdζ+s, p K+1 ) 1 and gcdζ, p K+1 ) 1 and the elements are different due to the fact that gcdζ +s, p K+1 ) 1 and gcds, N) = 1 implies gcdζ, p K+1 ) = 1. Similarly, p K+1 1 elements are not integer multiples of p K+1 in B N n) when y = xp K+1 since the first condition gcdζ + s, p K+1 ) 1 and gcds, N ) = 1 implies the second condition gcdζ, p K+1 ) = 1 p K+1. Thus, from 14), G N,yR N,s) = p K+1 G N,x R N,s ) 2 G N,x R N,s ) = p K+1 2) i M K x p i 2) = i My K+1 p i 2) when y = x and G N,yR N,s) = p K+1 G N,x R N,s ) p K+1 1) G N,x R N,s ) = i M p x K i 2) = i My K+1 p i 2) when y = xp K+1, which concludes the proof. Lemma 10: For N = p c 1 1, its divisor x = p n 1x) 1, 0 < n 1 x) c 1 and s µ N, G N,x R N,s ) is given as φ x 1), if x = p n 1x) 1, G N,x R N,s ) = Φ x 1), ifx = 1, 15) 0, otherwise.

12 Proof: When n 1 x) = c 1, G N,N R N,s ) = { r s gcdr s, N) = N, r µ N } = { 0 }, which implies G N,N R N,s ) = 1. When 0 < n 1 x) < c 1, G n N,p 1 x)r N,s ) = { r s gcdr s, N) = p n 1x) 1, r 1 µ N } = { r s r = mp n 1x) 1 + s, 0 m < p c 1 n 1 x) 1, gcdm, p 1 ) = 1 }. Since there are p c 1 n 1 x) 1 1 integer multiples of p 1 among { m 0 m < p c 1 n 1 x) 1 }, G n N,p 1 x)r N,s ) = p c 1 n 1 x) 1 p c 1 n 1 x) 1 1. When x = 1, G N,1 R N,s ) = { r s gcdr s, N) = 1, r µ N }. Note that G N,1 R N,s ) = R N,s N x=2 G N,x R N,s ) = µ N c 1 i=1 G N,p i 1 R N,s ). Since µ N = p c p 1 1) from Lemma 6 and G N,N R N,s ) = 1, G N,1 R N,s ) = p c 1 1 p c c 1 1 i=1 p c 1 i 1 p c 1 i 1 1 ) 1 = p c 1 1 2p c Now, the proof of Theorem 3 is given as follows. Proof of Theorem 3: When k = 1, 10) holds from Lemma 10. Suppose that 10) holds for any divisor x of N when k = K and let N = Np c K+1 K+1. Then, the set of all divisors of N is given as a { xp l K+1 x {all divisors of N}, 0 l c K+1 }. Note that µ N = { n 0 < n < N, gcdn, N ) = 1} = { n + mn n µ N, 0 m < p c K+1, gcdn + mn, p K+1 ) = 1}. Thus, R N,s = { n s n µ N } = m {0,,p c K+1 K+1 1} Rm N,s, where Rm N,s = { n + mn s n µ N, gcdn + mn, p K+1 ) = 1 }. When y = x for a divisor x of N, G N,yR m N,s ) = { n + mn n G N,xR N,s ), gcdn + mn, p K+1 ) = 1, gcdn s + mn, p K+1 ) = 1 } since gcdn s + mn, N ) = x implies gcdn s, N) = x and gcdn s+mn, N ) = x if gcdn s, N) = x and gcdn s+mn, p K+1 ) = 1. Also, when y = xp l K+1 for a divisor x of N and 1 l c K+1, G N,yR m N,s) = { n+mn n G N,x R N,s ), gcdn+mn, p K+1 ) = 1, gcdn s + mn, p l K+1 ) = pl K+1 }. Now, define B N n) = { n + mn 0 m < p c K+1 K+1 }. Then, it is easily seen that G N,yR N,s) = = p ck+1 K+1 1 n G N,x R N,s ) G N,yR m N,s ) B N n) α s N,N x, y), 16) where α s N,N x, y) = { ζ ζ n G N,x R N,s ) B N n), gcdζ + s, p K+1 ) 1 or gcdζ, p c K+1) y/x }. Similarly to the proof of Lemma 9, α s N,N x, y) B Nn) = 2p c K+1 1 K+1 when y = x and G N,yR N,s)

13 then is given as G N,yR N,s) =p c K+1 K+1 φ x i) i M x = i M y φ y i) j M K+1 y j M K x Φ x j) 2p c K+1 1 K+1 φ x i) Φ x j) i M x j Mx K Φ y j), since M y = M x, My K+1 = Mx K {K + 1}, φ y i) = φ x i) for i M x, Φ y j) = Φ x j) for j Mx K and Φ y K +1) = p c K+1 K+1 2p c K+1 1 K+1 so that i M Φ y K+1 y j) = p c K+1 K+1 2p c ) K+1 1 K+1 i M x Φ x i). Also, similarly to the proof of Lemma 9, α s N,N x, y) B Nn) = p c K+1 K+1 1 when y = xp c K+1 K+1 and G N,yR N,s) is then given as G N,yR N,s) = p c K+1 K+1 φ x i) Φ x j) p c K+1 K+1 1) φ x i) Φ x j) i M x j Mx K i M x j Mx K = since M y = M x {K + 1}, M K+1 y i M x φ x i) j M K x Φ x j) = i M y φ y i) j M K+1 y Φ y j), 17) 18) = M K x, φ y i) = φ x i) for i M x and φ y K + 1) = 1 so that i M y φ y i) = i M x φ x i) and i My K+1 Φ y j) = i M Φ x xj) for j M K K x = My K+1. When y = xp l K+1 for 0 < l < c K+1, G N,yR N,s) = pc K+1 1 K+1 G N,yR m N,s ) = { ζ ζ n G N,x R N,s ) B N,yn), ζ = mp l K+1 + s, 0 m < p c K+1 l K+1, gcdm, p K+1 ) = 1}. Similar to the proof of Lemma 10, G N,yR N,s) is given as since M y = M x {K +1}, M K+1 y G N,yR N,s) = p c K+1 l K+1 p c ) K+1 l 1 K+1 GN,x R N,s ) = p c K+1 l K+1 p c ) K+1 l 1 K+1 φ x i) i M x j Mx K = i M y φ y i) j M K+1 y Φ y j), Φ x j) 19) = Mx K, φ yi) = φ x i) for i M x and φ y K +1) = p c K+1 l K+1 p c K+1 l 1 K+1 so that i M y φ y i) = p c K+1 l K+1 p c ) K+1 l 1 K+1 i M x φ x i) and Φ y j) = Φ x j) for j Mx K = My K+1, which concludes the proof. C. Number of Available Chu sequences In this section, the number of available Chu sequences satisfying a given maximum magnitude of the cross-correlation is investigated. The main result of this subsection is given in the following theorems.

14 Theorem 4: Let a partial Chu sequence set C A be defined as C A = {a r r A } for a given partial index set A µ N and ˆθ A c be the maximum magnitude of the cross-correlation among sequences in C A. Then, ˆθ A c θ if and only if any two elements in C A, a r and a s, satisfy gcdr s, N) θ 2 /N. Proof: The proof is apparent since ˆθ r,s = N gcdr s, N) from Theorem 1. Theorem 5: For a given N θ N, let X N,θ = { all divisors of N = k i=1 p c i i greater than θ 2 /N }, x min = min X N,θ, M xmin = {n 1 n < x min, gcdn, N) = 1}, and x ϕmin = arg min x XN,θ ϕ x). Then, the largest cardinality among partial Chu sequence sets satisfying ˆθ A c θ is lower bounded by M xmin and upper bounded by ϕx ϕmin ). Proof: The lower bound is apparent from Theorem 4 since the difference of any two elements in M xmin is smaller than x min. Now, consider the upper bound part. Define R xϕmin = n Nxϕmin R n x ϕmin, 20) where R n x ϕmin = {mx ϕmin + n 0 m < N/x ϕmin } and N xϕmin = {n 1 n < x ϕmin, gcdn, x ϕmin ) = 1}. Then, it is apparent that µ n R xϕmin and N xϕmin = ϕx ϕmin ). Let A be any partial Chu sequence set satisfying ˆθ A c θ. Then, A Rn x ϕmin 1 because the difference of any two distinct elements in R n x ϕmin is at least x min. Thus, we can pick at most one element in R n x ϕmin for each n N xϕmin, which proves the upper bound. Lemma 11: R n x ϕmin µ n for any n N xϕmin. Proof: Since gcdn, x ϕmin ) = 1, gcdmx ϕmin +n, N) = gcdn, N/x ϕmin m)x ϕmin ) = gcdn, N/x ϕmin m). If gcdn, N) = 1, gcdmx ϕmin + n, N) = 1 when m = 0. If gcdn, N) = g > 1, gcdmx ϕmin + n, N) = 1 when 1 m = Ng 1)/x ϕmin g) < N/x ϕmin, which concludes the proof.

15 Note that Lemma 11 implies that the upper bound in Theorem 5 is tight. Also, although not shown licitly, exhaustive search showed that we can find a partial Chu sequence set with cardinality equal to the upper bound for the sequence length N up to Thus, although not proved, we may conjecture that the largest cardinality among partial Chu sequence sets satisfying ˆθ c A θ is ϕx ϕmin ). Example 1 N = 143, θ = 1430) : In this example, x min = x ϕmin = 11, R xϕmin = {11m + n 0 m < 13 and 1 n < 11} and ϕ x ϕmin ) = 10. It is easily found that C A, A = {1, 2,, 10}, satisfies ˆθ c A θ and A = ϕ x ϕmin ). Example 2 N = 154, θ = 1540) : In this example, x min = 11, x ϕmin = 14, R xϕmin = {14m + n 0 m < 11, 1 n < 14, gcdn, 14) = 1} and ϕ x ϕmin ) = 6. Note that {n 1 n < 14, gcdn, 14) = 1} = {1, 3, 5, 9, 11, 13} and 11 is not relatively prime with N so that it cannot be included in a partial Chu sequence set. Instead of 11, pick = 39 to construct A = {1, 3, 5, 9, 13, 39}. Then, it is easily found that C A satisfies ˆθ c A θ and A = ϕ x ϕ min ). The three subfigures in Fig. 1 shows the number of available Chu sequences satisfying ˆθ A /N 2 < θn 2 for θn 2 = 0.01, 0.05, 0.1, when N is around 512, 1024, and 2048, respectively. Also, Table I gives the prime factors of N used in Fig. 1. It is well known that N is preferred to be a prime number, which is confirmed from the results. Also, the results indicate that, in cases we must choose N among nonprime numbers by some reason, the number of available Chu sequences tends to increase as the smallest difference between prime factors increases. Thus, it is preferred to choose N composed of two prime factors with relatively large difference e.g., 508, 514, 515, 2045, 2049, 2051 as shown in Table I. IV. CONCLUSIONS In this paper, we analyzed generalized cross-correlation properties for Chu sequences. From the analysis, it was obtained that i) the magnitude of the cross-correlation function between any two Chu sequences for all possible lags, ii) the distribution of the maximum magnitude of the cross-correlation among a given

16 Chu sequence set and iii) the number of available Chu sequences satisfying a given cross-correlation constraint. APPENDIX A. Proof of Lemma 1 When N is an even number, we can rewrite θ r,s τ) as ) ) N 1 θ r,s τ) = jπ rk2 s k + τ)2 jπ N N ) ur,sg r,s 1 { = π sτ2 vr,s k 2 j2π sk i )} τg r,s + d τ ). u r,s g r,s Then, the squared magnitude, θ r,s τ) 2, is given as 21) θ r,s τ) 2 = u r,sg r,s 1 u r,sg r,s 1 l=0 { j2π vr,s k 2 2u r,s sk i τg r,s + d τ ) u r,s g r,s )} { j2π sl iτ g r,s + d τ ) v r,sl 2 )}. u r,s g r,s 2u r,s 22) The last term in 22) is periodic with period u r,s g r,s because { sl + ur,s g r,s )i τ g r,s + d τ ) j2π v r,s l + u r,s g r,s ) 2 )} u r,s g r,s 2u r,s { sliτ g r,s + d τ ) = j2π v r,sl 2 )} jπv r,s u r,s gr,s) 2 {j2πiτ g r,s + d τ + lg r,s v r,s )} u r,s g r,s 2u r,s { sliτ g r,s + d τ ) = j2π v r,sl 2 )}, u r,s g r,s 2v r,s 23) where the last equality comes from the fact that g r,s is always even so that jπv r,s u r,s g 2 r,s) = 1 because r s should be even when N is even. For a periodic function fx) with period N, it is easily seen that N 1 l=0 fl) = N+k 1 e=k fe). 24)

17 Then, from 23) and 24), 22) can be rewritten as u r,sg r,s 1 { θ r,s τ) 2 vr,s k 2 = j2π sk i )} τg r,s + d τ ) u r,sg r,s+k 1 { se iτ g r,s + d τ ) j2π v r,se 2 )} e=k u r,s g r,s 2u r,s u r,sg r,s 1 { vr,s k 2 = j2π sk i )} τg r,s + d τ ) u r,sg r,s 1 { se + k) iτ g r,s + d τ ) j2π v r,se + k) 2 )} e=0 u r,s g r,s 2u r,s u r,sg r,s 1 { vr,s e 2 = j2π se i )} ur,sg τg r,s + d τ ) r,s 1 j2π v ) r,ske. u r,s e=0 The last term can be divided in two terms in 25), when e = mu r,s and e mu r,s for 0 m < g r,s. Therefore, it can be ressed as θ r,s τ) 2 = { vr,s e 2 j2π se i τg r,s + d τ ) e=mu r,s + { vr,s e 2 j2π se i τg r,s + d τ ) e mu r,s )} ur,sg r,s 1 )} ur,sg r,s 1 j2π v ) r,ske u r,s j2π v r,ske u r,s When e mu r,s for 0 m < g r,s, the last term is equal to 0 from Lemma 1 in 26) because u r,s is relatively prime with v r,s. Accordingly, we can rewrite 26) as follows θ r,s τ) 2 = { vr,s e 2 j2π se i τg r,s + d τ ) e=mu r,s g r,s 1 { vr,s u r,s m 2 = j2π sm i τg r,s + d τ ) 2 g r,s g r,s 1 ) smdτ = u r,s g r,s 1) vr,sur,sm2. g r,s )} ur,sg r,s 1 )} ur,sg r,s 1 e=0 j2π v r,ske u r,s ) ) j2πv r,s me) When N is an odd number, we can derive θ r,s τ) 2 in similar to case when N is even. Accordingly, we. 25) 26) 27) can rewrite θ r,s τ) 2 as ) ) N 1 rk k + 1) s k + τ) k τ) θ r,s τ) = jπ jπ N N = jπ s ) ur,sg τ2 + τ) r,s 1 { vr,s k 2 + k) j2π sk i )} τg r,s + d τ ). u r,s g r,s 28)

18 Then, the squared magnitude, θ r,s τ) 2, is given as θ r,s τ) 2 = u r,sg r,s 1 u r,sg r,s 1 { j2π l=0 vr,s k 2 + k) sk i )} { τg r,s + d τ ) sl iτ g r,s + d τ ) j2π v r,s l 2 )} + l). u r,s g r,s 2u r,s The last term in 29) is periodic with period u r,s g r,s because { j2π s l + u } r,sg r,s ) i τ g r,s + d τ ) u r,s g r,s jπv r,s l + ur,s g r,s ) 2 + l + u r,s g r,s ) ) u r,s { sl iτ g r,s + d τ ) = j2π v r,s l 2 )} + l) jπv r,s g r,s u r,s g r,s +1))j2πi τ g r,s +d τ +lg r,s v r,s )) u r,s g r,s 2u r,s { sl iτ g r,s + d τ ) = j2π v r,s l 2 )} + l), u r,s g r,s 2u r,s where the last equality comes from the fact that g r,s +1 is always even so that jπv r,s g r,s u r,s g r,s + 1)) = 1 because r s should be odd when N is odd. Then, from 24) and 30), 29) can be rewritten as u r,sg r,s 1 { θ r,s τ) 2 vr,s k 2 + k) = j2π sk i )} τg r,s + d τ ) u r,sg r,s+k 1 { se iτ g r,s + d τ ) j2π v r,s e 2 )} + e) e=k u r,s g r,s 2u r,s u r,sg r,s 1 { vr,s e 2 + e) = j2π se i )} ur,sg τg r,s + d τ ) r,s 1 j2π v ) r,ske. u r,s e=0 The last term can be divided in two terms in 31), when e = mu r,s and e mu r,s for 0 m < g r,s. 29) 30) 31) Therefore, it can be ressed as θ r,s τ) 2 = { vr,s e 2 + e) j2π e=mu r,s + e mu r,s { j2π 2u r,s se i τg r,s + d τ ) u r,s g r,s )} ur,sg r,s 1 vr,s e 2 + e) se i )} ur,sg τg r,s + d τ ) r,s 1 j2π v ) r,ske u r,s j2π v r,ske u r,s ). 32)

19 When e mu r,s for 0 m < g r,s, the last term is equal to 0 from Lemma 1 in 32) because u r,s is relatively prime with v r,s and we can rewrite 32) as follows θ r,s τ) 2 = { vr,s e 2 + e) j2π se i τg r,s + d τ ) = e=mu r,s g r,s 1 { vr,s mu r,s + 1) j2π sm i τg r,s + d τ ) 2 g r,s ) 1) vr,smur,s+1) smdτ. g r,s 1 = u r,s g r,s On that way we have proved Lemma 1. g r,s )} ur,sg r,s 1 )} ur,sg r,s 1 e=0 j2π v r,ske u r,s ) j2πv r,s me) 33) ACKNOWLEDGMENT This research was supported by the MKE Ministry of Knowledge Economy), Korea, under the ITRC Information Technology Research Center) support program supervised by the IITA Institute of Information Technology Assessment) IITA-2008-C )) REFERENCES [1] D. V. Sarwate, Bounds on correlation and autocorrelation of sequences,, IEEE Trans. Inform. Theory, vol. IT-25, pp Nov [2] L. R. Welch, Lower bounds on the maximum cross correlation of signals,, IEEE Trans. Inform. Theory, vol. IT-20, No. 3, pp , May [3] V. M. Sidelnikov, On mutual correlation of sequences,, Soviet Math Doklady, vol. IT-12, No. 1, pp , [4] J. L. Massey, On Welch s bound for the crosscorrelation of a sequence set,, IEEE ISIT 90, pp. 385, San Diego, CA, USA, Sep [5] T. Kasami, Weight Distribution Formula for Same Class of Cyclic Codes,, Coordinated Science Lab., Univ. Illinois at Urbana- Champaign, Urbana, II., Tech. Rep. R285, [6] R. Gold, Maximal recursive sequences with 3-valued recursive crosscorrelation functions,, IEEE Trans. Inform. Theory, vol. IT-14, pp , Jan [7] R.L. Frank and S. A. Zadoff, Phase shift pulse codes with good periodic correlation properties,, IRE Trans. Inform. Theory, vol. IT-8, pp , Oct [8] D.C. Chu, Polyphase codes with good periodic correlation properties,, IEEE Trans. Inform. Theory, vol. IT-18, pp , July [9] R. L. Frank, Comments on Polyphase codes with good correlation properties,, IEEE Trans. Inform. Theory, vol. IT-19, pp. 244, Mar. 1973

20 [10] S. Bozta, R. Hammons and P. V. Kummer, 4-phase sequences with near optimum correlation properties,, IEEE Trans. Inform. Theory, vol. 38, pp , May [11] D. V. Sarwate and M. B. Pursley, Crosscorrelation properties of pseudo random and related sequencs,, Proc. IEEE, vol. 68, pp , May [12] S. I. Park, S. R. Park, I. Song, and N. Suehiro, Multiple acces interference reduction for QS-CDMA system with a novel clss of polyphase sequences,, IEEE Trans. Inform. Theory, vol. 46, pp , July 2000 [13] P. Fan and M. Darnell, Sequence Design for Communications Applications,, Wiley and RSP, New York, 1996) ISBN x [14] S. R. Park, I. Song, S. Yoon and J. Lee, A New polyphase sequence with perfect even and good odd cross-correlation functions for DS/CDMA systems,,, IEEE Trans. Vehicular Tech., vol. 51, No. 5, pp , Sep [15] B. M. Popovic, Generalized chirp-like polyphase sequences with optimum correlation properties,, IEEE Trans. Inform. Theory, vol. IT-38, pp , July A.B Nov [16] Y. Liu and P. Fan, Modified Chu sequences with smaller alphabet sized,, Electron. Lett.,vol. 40 No. 10, pp A.B, May [17] L. Lu and V. K. Dubey, Extended orthogonal polyphase codes for multicarrier CDMA System,,, IEEE Comm. Lett., vol. 8, No. 12, pp , Dec [18] D.Y. Peng, and P.Z. Fan, Generalized Sarwate bounds on the periodic autocorrelation and crosscorrelations of binary sequences,, Electron. Lett.,vol. 3824), pp. [19] R. Kwiatek and G. Zwara, The divisibility of intergers and integer relatively primes,, Formalized Mathematics., vol. 1, No. 5, pp , Nov [20] D. M. Burton, Elementary number theory,, Allyn and Bacon, 1980) ISBN

21 TABLE I THE PRIME FACTORS OF ADJACENT NUMBERS FOR N = 512, 1024 AND 2048 N Prime factors 3, , ,3, 5,17 7, ,19 2,257 5, ,3, 43 11, 47 N Prime factors ,3,5, , 7,73 3,11, ,41 2, 3 3,19 13, ,257 3, 7 3 N Prime factors 3 2, ,7,73 5, 409 2,3, 11,31 23, ,683 2, 5 2,41 7, ,3 3,

22 Number of available sequences θ N 2 =0.01 θ N 2 =0.05 θ N 2 = Length of sequences Number of available sequences θ N 2 =0.01 θ N 2 =0.05 θ N 2 = Length of sequences Number of available sequences θ N 2 =0.01 θ N 2 =0.05 θ N 2 = Length of sequences Fig. 1. Number of available sequences for a given maximum magnitude bound of the cross-correlation function