GB2447972A - Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain - Google Patents

Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain Download PDF

Info

Publication number
GB2447972A
GB2447972A GB0706263A GB0706263A GB2447972A GB 2447972 A GB2447972 A GB 2447972A GB 0706263 A GB0706263 A GB 0706263A GB 0706263 A GB0706263 A GB 0706263A GB 2447972 A GB2447972 A GB 2447972A
Authority
GB
United Kingdom
Prior art keywords
symbol
frequency offset
correlation
pilots
correlating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
GB0706263A
Other versions
GB0706263D0 (en
Inventor
Zhaohui Liu
Mike Page
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Holdings Corp
Original Assignee
Matsushita Electric Industrial Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Matsushita Electric Industrial Co Ltd filed Critical Matsushita Electric Industrial Co Ltd
Priority to GB0706263A priority Critical patent/GB2447972A/en
Publication of GB0706263D0 publication Critical patent/GB0706263D0/en
Publication of GB2447972A publication Critical patent/GB2447972A/en
Withdrawn legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2662Symbol synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2662Symbol synchronisation
    • H04L27/2663Coarse synchronisation, e.g. by correlation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2662Symbol synchronisation
    • H04L27/2665Fine synchronisation, e.g. by positioning the FFT window
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2675Pilot or known symbols
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2676Blind, i.e. without using known symbols
    • H04L27/2678Blind, i.e. without using known symbols using cyclostationarities, e.g. cyclic prefix or postfix
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2602Signal structure
    • H04L27/261Details of reference signals
    • H04L27/2613Structure of the reference signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

The invention synchronises the Symbol Timing in OFDM Systems using scattered pilots. Coarse symbol timing is achieved by autocorrelating <B>2</B> with a cyclic prefix and detecting a correlation peak <B>3</B>. At the same time, a frequency offset can be estimated <B>4</B>. Received samples are frequency offset compensated <B>5</B>, and a symbol index is determined <B>6</B> which allows a pattern of transmitted pilots to be identified. These frequency domain pilot sub-carrier symbols are then transformed to the time domain via an IFFT (Inverse Fast Fourier Transform) <B>7</B>. Given a search window around the coarse timing point, the frequency compensated samples are cross correlated with the time domain pilot sub-carriers waveform <B>8</B> to generate an estimate of the channel impulse response (CIR). Then the absolute value of the correlation is accumulated within a sliding window <B>9</B> and the position of the maximum <B>10</B> is regarded as the optimum timing point.

Description

Method and Appatus for Symbol Timing in OFDM Systems The present
invention relates to symbol timing and particularly symbol timing in orthogonal frequency division multiplexing (OFDM) systems with scattered pilots.
It is known in OFDM systems that in order to avoid inter-symbol interference (ISI) in a multipath channel, a guard interval (cyclic prefix) is included in the transmitted OFDM signal. The symbol timing is selected to minimize the inter-symbol interference.
Examples of OFDM communication systems based on scattered pilots are the Digital Video Broadcasting (DVB) system and the Digital Radio Mondiale (DRM) system. In general, there is no training symbol (such as a NULL symbol) in such an OFDM system. Instead, scattered pilots are provided for the purpose of channel estimation, so that a coherent receiver can be realised.
There are several approaches for symbol timing of an OFDM receiver: 1) Symbol timing based on correlation with the cyclic prefix.
2) Symbol timing based on correlation with training symbols.
3) Symbol timing based on maximizing the energy of the channel impulse response * 20 (CIR) estimate within a sliding window.
: It should be noted that there is usually no training symbol for a scattered pilot OFDM *:. system, so approach 2 may not be used.
Approach 1 can only give an approximate timing position for a multi-path channel.
Approach 2 cannot give the optimum timing position when the first path is not the strongest.
For approach 3, the CIR is usually calculated via interpolation and inverse fast Fourier transformation (IFFT) which is a frequency domain method and so may suffer from path "aliasing" where the incorrect signal is picked up for use in the lEFT. EP-A-1276291 relates to a frequency domain method which may suffer from such aliasing.
An object of the present invention is to achieve optimum symbol timing in a multi-path environment for a scattered pilots based OFDM system, even when the first path is not the strongest path and avoiding the problem of path aliasing.
The present invention provides a method of determining the symbol timing in an OFDM system, the method having a coarse timing and fine timing stage which are carried out with the aid of scattered pilots present in the symbol.
In particular, firstly, coarse symbol timing can be achieved by carrying out a correlation with a cyclic prefix. At the same time, a frequency offset can be estimated. Secondly, from the frequency offset compensated samples, the symbol index can be estimated and thus the known pilot sub-carriers can be transformed to the time domain via an IFFT. Thirdly, given a search window around the coarse timing point, the frequency compensated input samples can be correlated with the IFFT of the pilot sub-carriers. Finally, the absolute value of the correlation is accumulated within a sliding window and the position of the maximum is regarded as the optimum timing point. This approach has no path "alias" and is always able to synchronise to the first arriving path even if it is not the strongest path.
The present invention also provides an apparatus comprising means for performing the method.
In order that the present invention be more readily understood, an embodiment will be described with reference to the accompanying figures in which: Fig. 1 shows a schematic diagram according to a preferred embodiment of the present invention; Fig. 2 shows an OFDM symbol with cyclic prefix; Fig. 3 shows an example of the arrangement of scattered pilots in OFDM symbols; Fig. 4 shows a graph where windows are used to search for the optimum timing point; A preferred embodiment of symbol timing approach is shown Fig. 1. The figure shows the steps after analogue to digital conversion of the incoming signal in a receiver of digital broadcasting system which utilises OFDM with scattered pilots. Such a system is DVB and DRM but the invention is not limited to use in such systems.
After A/D conversion, the input samples are provided to memory 1 and auto-correlation 2 is carried out within a sliding window. The coarse symbol timing and frequency offset estimation 4 can be found from the position of maximum correlation 3. After frequency offset compensation 5, the symbol index 6 can be estimated from the location of the scattered pilots. Then the waveform of the pilot sub-carriers in the time domain can be constructed by carrying out IFFT 7. Applying cross-correlation 8 between the input samples from the frequency offset compensation 5 and the constructed pilot waveform, the amplitude of channel impulse response is estimated. Accumulation 9 within a sliding window will ensure synchronisation to the first arnval path even if it is not the strongest path. This is the overall method carried out in accordance with the preferred embodiment. The steps will be * considered in detail below. *. ** * * Se..
:* Fig. 2 shows the transmitted OFDM signal with its cyclic prefix. By calculating the auto- *:. 20 correlation of thc input samples, the coarse symbol timing and frequency offset estimate can be found. The details of such steps can be found in many references, e.g., reference [1]. a.. a
Figure 3 shows the pilot position within each OFDM symbol. In general, the power of each pilot is boosted and their phases and positions are changed depending on the symbol index.
After obtaining the coarse symbol timing, the scattered pilots can be used to obtain the optimal timing. To use this information, the symbol index must first be detected, since the positions and phases of the scattered pilots are different for each symbol. This can be achieved by the following steps: 1) The input samples are frequency offset compensated, generating a sample vector x[n] 2) An FFT is applied to x[n], generating the frequency domain output Y, where s is the symbol index and k is the sub-carrier index 3) The absolute value of YSk is accumulated at each presumed position, the maximum corresponding to the pilot pattern (i.e. the pilot sub-carrier positions in the symbol).
This is because the pilot sub-carriers are power boosted.
4) Once the pilot pattern has been detected, the symbol index can be searched within the set that comprises of all symbols with the detected pilot pattern. (Note that as shown in figure 3, several symbols may share the same pattern of pilot positions, this pattern repeating after a pilot pattern cycle). Denote the set as F and the symbol index. can be found from the following: for each s E P . exp(-j( - ) end = arg max(z) S. Where c is the pilot pattern cycle, is the phase of the pilot at sub-carrier k and 0SS* * . *5** symbol s. S...
The estimated symbol index corresponds to the maximum of z,.
Once the symbol index is known, then the waveform of pilots in the time domain may be * generated using an IFFT. Given a search window around the coarse symbol timing point, the : recreated pilot waveform may be cross-correlated with the frequency offset compensated input samples, and the optimum timing can be found by maximizing the energy of the CW within a sliding window. As shown in Figure 4, assuming the coarse symbol timing point is, n0, the length of the accumulation window is L, the search window is W and let p[n] be the IFFT of the pilot sub-carriers, then the optimal timing n can be found as follows: (1) Calculate the cross-correlation for n=OtoW-1 y[n] = p [k]x[k + n0 + n -w/2] end (2) Calculate the energy of the CIR within a sliding window for n=OtoW-L z[nJ = y[k + n]j end (3) flop1 =argmaxz[nJ As the CuR is estimated in time domain (rather than the frequency domain), any path "alias" is avoided.
In general, the length of accumulation window L should be equal to the channel size, which is not known. However, as the transmission mode is related to the channel characteristics, L should be no more than the length of the cyclic prefix. On the other hand, because the :. scattered pilots are regularly spaced, its time domain waveform is cyclic. So if the
. accumulation window is too large, peaks introduced by periodicity will also be involved. S....DTD: This can be solved by including irregularly spaced pilot sub-carriers to break up the periodic S...
sequence. In general, more pilot sub-carriers are available for the first symbol for the * 15 purpose of frame synchronization. For example, in DRM, gain pilots are regularly spaced * : and extra, irregularly spaced time reference pilots are available for the first symbol. *5S
Further approaches can be taken to mitigate the effects of noise. For instance, the output of the cross-correlation can be compared with a threshold; if the output is less than the threshold, it can be regarded as mainly caused by noise and set to zero.
References 1) M. Speth and etc., "Optimum receiver design for wireless broad-band systems using OFDM -part II," IEEE Trans. Comm., vol. 49, No. 4, pp. 571 -578, 2001 2) Wei Chee Lim and etc., "Joint channel estimation and OFDM synchronization in multipath fading, "in Proc. ICC'04, June, 2004, pp. 983 -987 * * * *** S... * . S...
S I... S.. S. S
S S S * S.
S S..
S

Claims (5)

1. A method of determining symbol timing in a scattered pilot based OFDM system comprising: receiving an OFDM symbol having a cyclic prefix; calculating auto-correlation of the received symbol in order to determine coarse symbol timing and a frequency offset estimation; performing frequency offset compensation of the received OFDM symbol on the basis of the estimated frequency offset; detecting the symbol index from the location of scattered pilots in the symbol; generating a waveform of the pilots in the time domain using an inverse Fast Fourier Transform (IFFT); estimating an amplitude of a channel impulse response (CIR) by cross correlating the generated waveform and frequency offset compensated OFDM symbol accumulating an absolute value of the cross-correlation within a sliding window, and determining the maximum value wherein the maximum value is the optimum symbol timing point.
2. The method of claim 1 wherein the length of the sliding window is related to the length of the cyclic prefix.
3. The method of claim 1 or 2 wherein the OFDM symbol includes irregularly spaced scattered pilots.
4. The method of any preceding claim wherein the output of the cross-correlation is compared with a threshold, and if the output is less than the threshold, it is set to zero.
5. An apparatus for determining symbol timing in a scattered pilot based OFDM system, comprising: means for receiving an OFDM symbol having a cyclic prefix; means for auto-correlating the received symbol in order to determine coarse symbol timing and a frequency offset estimation; means for performing frequency offset compensation of the received OFDM symbol on the basis of the estimated frequency offset; means for detecting the symbol index from the location of scattered pilots in the symbol; means for generating a waveform of the pilots in the time domain using an inverse Fast Fourier Transform (1FF!'); means for estimating an amplitude of a channel impulse response (CIR) by cross correlating the generated waveform and frequency offset compensated OFDM symbol; means for accumulating an absolute value of the cross-correlation within a sliding window and determining the maximum value, wherein the maximum value is the optimum symbol timing point.
GB0706263A 2007-03-30 2007-03-30 Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain Withdrawn GB2447972A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
GB0706263A GB2447972A (en) 2007-03-30 2007-03-30 Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
GB0706263A GB2447972A (en) 2007-03-30 2007-03-30 Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain

Publications (2)

Publication Number Publication Date
GB0706263D0 GB0706263D0 (en) 2007-05-09
GB2447972A true GB2447972A (en) 2008-10-01

Family

ID=38050565

Family Applications (1)

Application Number Title Priority Date Filing Date
GB0706263A Withdrawn GB2447972A (en) 2007-03-30 2007-03-30 Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain

Country Status (1)

Country Link
GB (1) GB2447972A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010131138A1 (en) 2009-05-14 2010-11-18 Koninklijke Philips Electronics, N.V. Robust sensing of dvb-t/h transmissions
US8194799B2 (en) 2009-03-30 2012-06-05 King Fahd University of Pertroleum & Minerals Cyclic prefix-based enhanced data recovery method
CN101924725B (en) * 2009-06-17 2013-01-30 国民技术股份有限公司 Frame synchronization method and device for OFDM system
US8867634B2 (en) 2009-03-12 2014-10-21 Thomson Licensing Method and appratus for spectrum sensing for OFDM systems employing pilot tones
WO2017036193A1 (en) * 2015-09-01 2017-03-09 中兴通讯股份有限公司 Method and device for estimating frequency offset
EP2437450A4 (en) * 2009-06-30 2017-06-07 ZTE Corporation Device and method for estimating time offset in orthogonal frequency division multiplexing (ofdm) system
CN111884979A (en) * 2020-07-30 2020-11-03 电子科技大学 OFDM smart grid impulse noise resistant symbol synchronization method
CN111884978A (en) * 2020-07-30 2020-11-03 电子科技大学 OFDM (orthogonal frequency division multiplexing) anti-impulse noise symbol synchronization method
US20220408285A1 (en) * 2018-10-24 2022-12-22 Lg Electronics Inc. Method and device for sidelink terminal to detect sidelink signal in wireless communication system

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11606240B1 (en) * 2021-09-30 2023-03-14 Silicon Laboratories Inc. Using preamble portion having irregular carrier spacing for frequency synchronization

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
IEEE Transactions on Communications, Volume 50, Number 4, April 2002, IEEE, Landstrom D et al, "Symbol time offset estimation in coherent OFDM systems", pages 545-549, available from http://ieeexplore.ieee.org/iel5/26/21498/00996067.pdf *
Proc. IEEE 63rd Vehicular Technology Conference (VTC) 2006, 7-10 May 2006, Jongkyung Kim et al, "Synchronization and channel estimation in cyclic postfix based OFDM system", Vol. 4 Pages 2028-32, available from http://ieeexplore.ieee.org/iel5/11096/35444/01683202.pdf?tp=&isnumber=&arnumber=1683202 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8867634B2 (en) 2009-03-12 2014-10-21 Thomson Licensing Method and appratus for spectrum sensing for OFDM systems employing pilot tones
US8194799B2 (en) 2009-03-30 2012-06-05 King Fahd University of Pertroleum & Minerals Cyclic prefix-based enhanced data recovery method
US8942336B2 (en) 2009-05-14 2015-01-27 Koninklijke Philips N.V. Robust sensing of DVB-T/H transmissions
KR101633628B1 (en) * 2009-05-14 2016-06-27 코닌클리케 필립스 엔.브이. Robust sensing of dvb-t/h transmissions
JP2012527143A (en) * 2009-05-14 2012-11-01 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Robust detection of DVB-T / H communication
CN102422634A (en) * 2009-05-14 2012-04-18 皇家飞利浦电子股份有限公司 Robust sensing of DVB-T/H transmissions
KR20120036839A (en) * 2009-05-14 2012-04-18 코닌클리케 필립스 일렉트로닉스 엔.브이. Robust sensing of dvb-t/h transmissions
CN102422634B (en) * 2009-05-14 2014-12-03 皇家飞利浦电子股份有限公司 Robust sensing of DVB-T/H transmissions
WO2010131138A1 (en) 2009-05-14 2010-11-18 Koninklijke Philips Electronics, N.V. Robust sensing of dvb-t/h transmissions
CN101924725B (en) * 2009-06-17 2013-01-30 国民技术股份有限公司 Frame synchronization method and device for OFDM system
EP2437450A4 (en) * 2009-06-30 2017-06-07 ZTE Corporation Device and method for estimating time offset in orthogonal frequency division multiplexing (ofdm) system
WO2017036193A1 (en) * 2015-09-01 2017-03-09 中兴通讯股份有限公司 Method and device for estimating frequency offset
US10148463B2 (en) 2015-09-01 2018-12-04 Xi'an Zhongxing New Software Co., Ltd. Method and device for estimating frequency offset
US20220408285A1 (en) * 2018-10-24 2022-12-22 Lg Electronics Inc. Method and device for sidelink terminal to detect sidelink signal in wireless communication system
CN111884979A (en) * 2020-07-30 2020-11-03 电子科技大学 OFDM smart grid impulse noise resistant symbol synchronization method
CN111884978A (en) * 2020-07-30 2020-11-03 电子科技大学 OFDM (orthogonal frequency division multiplexing) anti-impulse noise symbol synchronization method

Also Published As

Publication number Publication date
GB0706263D0 (en) 2007-05-09

Similar Documents

Publication Publication Date Title
GB2447972A (en) Synchronising OFDM Symbols by auto-correlating cyclic prefix then cross correlating scattered pilots in the time domain
CN108040028A (en) OFDM system anti-interference signal detection and synchronization method based on local sequence cross-correlation detection
EP1729467A2 (en) Apparatus and method for transmitting/receiving preamble signal in a wireless communication system
US8520783B2 (en) Receiving device
CN107231176A (en) A kind of OFDM MFSK underwater sound communications broadband Doppler shift method based on subcarrier energy
JP2006504359A (en) Method and apparatus for synchronizing training sequences
KR20010002477A (en) Symbol timing and frequency synchronizing device for OFDM signals and method thereof
JP5144687B2 (en) OFDM Code Timing Recovery Method and System
CN101836411A (en) Use time division multiplexed pilots in the broadcasting ofdm system synchronously
WO2008148078A2 (en) Ofdm-based device and method for performing synchronization in the presence of interference signals
KR101467252B1 (en) Method for estimating ofdm integer frequency offset, ofdm integer frequency offset estimator and ofdm receiver system
KR101110312B1 (en) Apparatus and method for frequency offset estimation
CN109495414B (en) Frequency offset estimation method, device and equipment and computer readable storage medium
JP3905541B2 (en) Delay profile estimation apparatus and correlator
CN101843061B (en) A method to discriminate a real echo peak from an aliased echo peak
CN101836410A (en) Use time division multiplexed pilots in the broadcasting ofdm system synchronously
KR101128287B1 (en) Ofdm receiver being capable of estimating timing error in multi-path fading channel, ofdm system having the same and timing error estimation method thereof
Nasraoui et al. Robust brute force and reduced complexity approaches for timing synchronization in IEEE 802.11 a/g WLANs
JPH08223132A (en) Receiver for digital transmission signal and digital transmission system
KR20020089793A (en) Synchronous signal detection apparatus for ofdm receiver
CN104618294B (en) OFDM integer frequency offsets Method of fast estimating based on training sequence and system
KR20090131543A (en) Method and apparatus for estimating integer frequency offset in orthogonal frequency division multiplexing system using differential combination
CN101997792B (en) Cyclic prefix-based blind channel length estimation method
Xiao et al. Time synchronization algorithms for IEEE 802.11 OFDM systems
Laourine et al. Joint carrier and sampling frequency offset estimation based on harmonic retrieval

Legal Events

Date Code Title Description
WAP Application withdrawn, taken to be withdrawn or refused ** after publication under section 16(1)