Implementation of Orthogonal Frequency Coded SAW Devices Using Apodized Reflectors

Similar documents
Weighted SAW Reflector Gratings for Orthogonal Frequency Coded SAW ID Tags and Sensors

696 ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 55, no. 3, march 2008

Design, Analysis And Implementation Of Orthogonal Frequency Coding In Saw Devices Used For Spread Spectrum Tags And Sensors

Weighted SAW Reflector Gratings for Orthogonal Frequency Coded SAW ID Tags and Sensors DIV

DESIGN OF GLOBAL SAW RFID TAG DEVICES C. S. Hartmann, P. Brown, and J. Bellamy RF SAW, Inc., 900 Alpha Drive Ste 400, Richardson, TX, U.S.A.

Ultra-wideband Orthogonal Frequency Coded Saw Correlators

Wireless Passive SAW Sensors using Coded Spread Spectrum Techniques

Implementation of On-Off Passive Wireless Surface Acoustic Wave Sensor Using Coding and Switching Techniques DIV

Design And Simulation For Encoded Pn-ofc Saw Sensor Systems

MAKING TRANSIENT ANTENNA MEASUREMENTS

200-GHz 8-µs LFM Optical Waveform Generation for High- Resolution Coherent Imaging

ACOUSTIC MICRO IMAGING ANALYSIS METHODS FOR 3D PACKAGES

High-Speed Optical Modulators and Photonic Sideband Management

Paper VI. Non-synchronous resonators on leaky substrates. J. Meltaus, V. P. Plessky, and S. S. Hong. Copyright 2005 IEEE.

Interference of Chirp Sequence Radars by OFDM Radars at 77 GHz

UWB SAW sensors and tags

PWST SAW - Sensor System

Bragg and fiber gratings. Mikko Saarinen

Modern radio techniques

Antenna Measurements using Modulated Signals

ME scope Application Note 01 The FFT, Leakage, and Windowing

Spread Spectrum (SS) is a means of transmission in which the signal occupies a

Direct-Conversion I-Q Modulator Simulation by Andy Howard, Applications Engineer Agilent EEsof EDA

A SHEAR WAVE TRANSDUCER ARRAY FOR REAL-TIME IMAGING. R.L. Baer and G.S. Kino. Edward L. Ginzton Laboratory Stanford University Stanford, CA 94305

Time and Frequency Domain Windowing of LFM Pulses Mark A. Richards

Keysight Technologies Pulsed Antenna Measurements Using PNA Network Analyzers

SIDELOBES REDUCTION USING SIMPLE TWO AND TRI-STAGES NON LINEAR FREQUENCY MODULA- TION (NLFM)

Passive Polymer. Figure 1 (a) and (b). Diagram of a 1-3 composite (left) and a 2-2 composite (right).

Chapter 1. Overview. 1.1 Introduction

Ultrasound Beamforming and Image Formation. Jeremy J. Dahl

Time reversal: a flexible approach for identifying and measuring surface acoustic wave delay lines acting as wireless, passive sensors

Ultrasonic Linear Array Medical Imaging System

Radar-Verfahren und -Signalverarbeitung

Ginzton Laboratory, W. W. Hansen Laboratories of Physics Stanford University, Stanford, CA 94305

Liquid sensor probe using reflecting SH-SAW delay line

A STUDY ON NON-CONTACT ULTRASONIC TECHNIQUE FOR ON-LINE INSPECTION OF CFRP

Experiment 2 Effects of Filtering

A NOVEL DIGITAL BEAMFORMER WITH LOW ANGLE RESOLUTION FOR VEHICLE TRACKING RADAR

Low Power LFM Pulse Compression RADAR with Sidelobe suppression

NEW LASER ULTRASONIC INTERFEROMETER FOR INDUSTRIAL APPLICATIONS B.Pouet and S.Breugnot Bossa Nova Technologies; Venice, CA, USA

QUESTION BANK SUBJECT: DIGITAL COMMUNICATION (15EC61)

Performance of Wideband Mobile Channel with Perfect Synchronism BPSK vs QPSK DS-CDMA

SEPTEMBER VOL. 38, NO. 9 ELECTRONIC DEFENSE SIMULTANEOUS SIGNAL ERRORS IN WIDEBAND IFM RECEIVERS WIDE, WIDER, WIDEST SYNTHETIC APERTURE ANTENNAS

Implementing Orthogonal Binary Overlay on a Pulse Train using Frequency Modulation

Lecture 3: Wireless Physical Layer: Modulation Techniques. Mythili Vutukuru CS 653 Spring 2014 Jan 13, Monday

Lecture 9: Spread Spectrum Modulation Techniques

Advanced 3G & 4G Wireless Communication Prof. Aditya K. Jaganathan Department of Electrical Engineering Indian Institute of Technology, Kanpur

Ultra-wideband Spread Spectrum Communications using Software Defined Radio and Surface Acoustic Wave Correlators

Impulse Response as a Measurement of the Quality of Chirp Radar Pulses

Characterization of L5 Receiver Performance Using Digital Pulse Blanking

International Distinguished Lecturer Program

DIVERSE RADAR PULSE-TRAIN WITH FAVOURABLE AUTOCORRELATION AND AMBIGUITY FUNCTIONS

SAW Reflective Transducers and Antennas for Orthogonal Frequency Coded SAW Sensors

CHAPTER 5 FINE-TUNING OF AN ECDL WITH AN INTRACAVITY LIQUID CRYSTAL ELEMENT

The Effect of Notch Filter on RFI Suppression

Application Note 106 IP2 Measurements of Wideband Amplifiers v1.0

Introduction to SAW Filter Theory & Design Techniques

Effect of coupling conditions on ultrasonic echo parameters

Part A: Spread Spectrum Systems

High-overtone Bulk Acoustic Resonator (HBAR) as passive sensor: towards microwave wireless interrogation

DIGITAL BEAM-FORMING ANTENNA OPTIMIZATION FOR REFLECTOR BASED SPACE DEBRIS RADAR SYSTEM

Study on Multi-tone Signals for Design and Testing of Linear Circuits and Systems

Simulating and Testing of Signal Processing Methods for Frequency Stepped Chirp Radar

EWGAE 2010 Vienna, 8th to 10th September

Termination Insensitive Mixers By Howard Hausman President/CEO, MITEQ, Inc. 100 Davids Drive Hauppauge, NY

A 1-W GaAs Class-E Power Amplifier with an FBAR Filter Embedded in the Output Network

Mechanical Spectrum Analyzer in Silicon using Micromachined Accelerometers with Time-Varying Electrostatic Feedback

Spread Spectrum Techniques

Signals and Systems Lecture 9 Communication Systems Frequency-Division Multiplexing and Frequency Modulation (FM)

UNIT- 7. Frequencies above 30Mhz tend to travel in straight lines they are limited in their propagation by the curvature of the earth.

Photonic Signal Processing(PSP) of Microwave Signals

Plane wave excitation by taper array for optical leaky waveguide antenna

Amplitude independent RF instantaneous frequency measurement system using photonic Hilbert transform

Radar Signatures and Relations to Radar Cross Section. Mr P E R Galloway. Roke Manor Research Ltd, Romsey, Hampshire, United Kingdom

Extending Acoustic Microscopy for Comprehensive Failure Analysis Applications

Know how Pulsed Doppler radar works and how it s able to determine target velocity. Know how the Moving Target Indicator (MTI) determines target

Waveform Multiplexing using Chirp Rate Diversity for Chirp-Sequence based MIMO Radar Systems

Passive, Wireless SAW OFC Strain Sensor and Software Defined Radio Interrogator

Multi-Path Fading Channel

Electronically Steerable planer Phased Array Antenna

Comparison of FMCW-LiDAR system with optical- and electricaldomain swept light sources toward self-driving mobility application

Part A: Spread Spectrum Systems

Switch-less Dual-frequency Reconfigurable CMOS Oscillator using One Single Piezoelectric AlN MEMS Resonator with Co-existing S0 and S1 Lamb-wave Modes

Application of Ultrasonic Guided Waves for Characterization of Defects in Pipeline of Nuclear Power Plants. Younho Cho

Communication using Synchronization of Chaos in Semiconductor Lasers with optoelectronic feedback

THE PROPAGATION OF PARTIAL DISCHARGE PULSES IN A HIGH VOLTAGE CABLE

Evaluation of RF power degradation in microwave photonic systems employing uniform period fibre Bragg gratings

Characterization of Silicon-based Ultrasonic Nozzles

Improving CDM Measurements With Frequency Domain Specifications

Narrow- and wideband channels

Mobile Radio Propagation Channel Models

Agilent PNA Microwave Network Analyzers

MEASUREMENT OF RAYLEIGH WAVE ATTENUATION IN GRANITE USING

EE390 Final Exam Fall Term 2002 Friday, December 13, 2002

Standard Guide for Evaluating Characteristics of Ultrasonic Search Units 1

Chapter 4. Part 2(a) Digital Modulation Techniques

Constant-Envelope Variations of OFDM and OFDM-CDMA

Understanding How Frequency, Beam Patterns of Transducers, and Reflection Characteristics of Targets Affect the Performance of Ultrasonic Sensors

LWA Beamforming Design Concept

A STUDY OF AM AND FM SIGNAL RECEPTION OF TIME MODULATED LINEAR ANTENNA ARRAYS

Transcription:

Implementation of Orthogonal Frequency Coded SAW Devices Using Apodized Reflectors Derek Puccio, Don Malocha, Nancy Saldanha Department of Electrical and Computer Engineering University of Central Florida Orlando, FL 3286-245 rick.puccio@ieee.org Abstract Recently, orthogonal frequency coding (OFC) has been presented as a novel method for coding SAW tags and sensors []. Orthogonal frequency coding is a spread spectrum technique and has been shown to provide enhanced processing gain and reduced time ambiguity resulting in greater range and increased sensitivity when compared with single carrier frequency devices. The sensor works both as a tag and a sensor with the ability to send back tagged sensor information in a multi-sensor environment. The tag information is provided by a series of reflectors which map into a known chip sequence. The time-chip- sequence is coded by differing OFC and PN sequences. Therefore, the implementation of an OFC sensor requires reflectors having differing local carrier frequencies. In the case of narrow fractional bandwidths or high reflectivity (such as on LiNbO3), it is desirable to adjust the reflectivity per electrode in the various chips. For varying system requirements, the use of weighted reflectors is an option; both apodization and variable weighted reflectors are investigated. Experimental results on cosine weighted apodized reflectors will be compared to COM predictions on YZ LiNbO 3. This paper presents several OFC SAW device embodiments that are employed using wideband input transducers and multiple weighted reflector gratings. The devices operate in the differential mode using gratings on either side of the transducer. The advantages of using the weighted reflector gratings in OFC SAW devices are discussed. OFC temperature sensor device experimental results are presented using weighted reflectors to obtain the proper time coded response and compared to predictions. I. INTRODUCTION Orthogonal frequency coding (OFC) is a spread spectrum coding technique that has been successfully applied to SAW tags and sensors [2]. OFC offers several advantages over single frequency SAW tags and sensors including enhanced processing gain, increased range using chirp interrogation, and improved security. In addition, OFC relies on the use of several frequencies making simultaneous sensing and tagging possible in multiple sensor environments. OFC SAW devices are implemented using reflectors with periodicities chosen to match the chip frequency of interest. As a result, an OFC SAW device contains a series of reflectors whose center frequencies correspond to the OFC signal of interest. In the case of high reflectivity or narrow bandwidth materials, it is desirable to control the reflection and transmission characteristics of each reflector. The goal of this research is to study weighted reflectors as a method of controlling the reflected and transmitted SAW energy. In addition, arbitrary pulse shapes can be achieved using reflector weighting. Furthermore, pulse shaping can be used in phase shift keying, such as in phase and quadrature channels in minimum shift keying (MSK). Several steps must be taken in order to properly implement weighted reflectors. Both withdrawal weighting and apodization result in structures that contain various metallized and free surface regions. As a result, an accurate description of SAW velocity within these regions is critical to successful reflector design. Experimental velocity results as a function of normalized metal thickness and reflector duty cycle are given, and a discussion of the results is provided. Additionally, a two dimensional COM model has been created for analysis and design of apodized reflectors, and experimental results are compared with predictions for cosine weighted reflectors. Finally, experimental and simulated responses of an OFC SAW sensor using cosine weighted reflectors are given. The responses are applied to a simulated transceiver and results are discussed. II. REVIEW OF OFC TAGS AND SENSORS Orthogonal frequency coding (OFC) makes use of orthogonal frequencies in a manner similar to M-ary FSK. The orthogonality integral requires a relationship between local chip frequencies and their bandwidths. An example of a signal utilizing orthogonal frequencies is the linear stepped chirp shown in Figure. Each of the seven chips has a unique carrier center frequency as indicated by the individual This work is supported by the Florida Solar Energy Center (FSEC) and NASA. D. Puccio is supported by the NASA Graduate Student Researcher s Program (GSRP).

.8 Magnitude (Linear).6.4 f 3 Piezoelectric Substrate f 5 f f 6 f 2 f 4 f τ τ 2 f f 4 f 2 f 6 f f 5 f 3.2 2.2.4.6.8.2.4.6.8 Normalized Frequency.5 Magnitude (db) 3 4 5 6 Experimental COM Simulated.6.8.2.4.6.8 2 2.2 2.4.5.5.5 2 3 4 5 6 7 Normalized Time (Chip Lengths) 2 3 4 5 6 7 Normalized Time (Chip Lengths) Figure. Example of OFC uniform chip frequency responses (upper), stepped chirp response (middle), orthogonal frequency coded waveform (lower) chip frequency responses also shown in Figure. The frequency responses indicate that the center frequency of each chip is coincident with nulls of all other chip frequency responses. A level of coding is achieved by shuffling the chips so that the local carrier frequencies are no longer contiguous in time as shown in Figure. The order in which the orthogonal frequencies are used determines the code. The coded signal and the linear stepped chirp occupy the same bandwidth and all coding information is contained within the signal phase. OFC can be applied to SAW tags and sensors using shorted periodic reflector gratings as shown in Figure 2. Each reflector s periodicity is chosen so that the reflector center frequency corresponds to the OFC chip of interest. The device schematic shown is that of a temperature sensor which uses identical reflector banks on either side of a wideband transducer. However, different free space delays are employed on each side of the transducer. The impulse response of this device contains two identical approximations of the OFC signal of interest as shown in Figure 2. During matched filtering, this device produces two compressed pulses. The separation between the pulses is proportional to device temperature. III. SAW VELOCITY EXTRACTION Weighted reflectors can be employed using withdrawal weighting or apodization. Both techniques rely on eliminating portions of the metallized surface used in a uniform reflector. As a result, weighted reflector designs must account for velocity perturbations in order to keep wave fronts from different sections of the reflector phase coherent. For this reason, reflector velocity data were extracted for YZ LiNbO 3. Reflector velocities were extracted using the following equation derived from a transmission line model analogy [3]. h a a v B a vg, = vfs + sin () λ p p vfs π Y p where, v = free surface velocity fs a p Time (us) Figure 2. Orthogonal frequency coded SAW temperature sensor schematic (upper) and sensor impulse response (lower). = reflector grating duty cycle h = normalized metal thickness λ This equation is a function of two independent variables, the reflector grating duty cycle and normalized metal thickness. The difference between metallized and free surface velocities, v, is proportional to the normalized metal thickness. The transmission line model uses a shunt susceptance term, B /Y, to account for stored energy at the front and back of each reflector electrode, and this term is proportional to the square of the normalized metal thickness. The unknown terms in () were extracted using a test mask with several test structures. The test devices used a simple delay line with a reflector grating situated nearby. This embodiment is used in order to isolate just the reflector response using signal processing [3]. The test mask contained thirty five of these devices each having different combinations of reflector grating duty cycle and center frequency wavelength. Seven wavelengths were used between 8 and 6 µm, and five reflector grating duty cycles between 4% and 6% were employed. The mask was fabricated on YZ LiNbO 3, and signal processing techniques were applied to the swept frequency responses to extract reflector grating velocity data. A three dimensional optimization was performed using these data in order to obtain the unknown terms in (). The results for several metal thicknesses are shown in Figure 3, and there is good

Figure 3. YZ LiNbO 3 velocity vs. reflector duty cycle for several normalized metal thickensses. Curvature increases with metal thickness due to stored energy at edges of electrodes. agreement between measured and predicted results. Note that the curvature of each line is proportional to the normalized metal thickness. This is a direct result of the stored energy effect which has a large influence for thicker films. IV. WEIGHTED REFLECTORS A. Reflector Weighting Techniques Reflector weighting can be accomplished using withdrawal weighting or apodization. Previous studies of withdrawal weighting focused on eliminating unwanted modes in resonant cavities by shaping reflector frequency responses [4,5,6]. For OFC tag and sensor applications, the time domain response is of primary concern, and is better approximated through apodization. Therefore, the goals of this research are simulation and implementation of apodized reflectors. B. Implementation of Apodized Reflectors In an apodized reflector, each electrode covers a portion of the beam aperture defined by a sampled time window. As an example, Figure 4 shows a schematic of a sixteen period cosine weighted reflector. Note that the first eight electrodes are connected to the bottom bus bar, and the rest are connected at the top. This is done in order to reduce group delay variations over the reflector aperture. In order to simulate such a device, a two dimensional COM model was developed. The model is designed to segment the beam aperture in to uniform tracks as shown. The sum of each track s swept frequency response yields the overall reflector response. When designing apodized reflectors, special care must be taken to ensure phase coherence of the reflected waves from different tracks. As an example, consider an incident wave on the left side of the reflector shown in Figure 4. A portion of the wave will be reflected from track A, and after a short delay, another portion is reflected from Figure 4. Schematic of cosine weighted apodized reflector. Reflector is simulated by creating uniform tracks and summing responses from each. track B. This delay is implemented using a free surface which has a slightly higher velocity than within the electrodes. As a result, if no correction is made, the reflected waves from tracks A and B will have a small phase difference. Therefore, the initial free space delays must be extended slightly in order to keep the reflected waves coherent. The resulting electrodes will be made up of several smaller rectangles that are staggered across the aperture as shown. C. Cosine Weighted Reflector Results Using the apodization technique described, a weighted reflector was designed and fabricated using a cosine window for operation on YZ LiNbO 3 at 25 MHz. The reflector contained 24 periods, and was placed near a simple delay line as in [3]. The device s swept frequency response was obtained and signal processing applied in order to isolate the reflector response. The device was also simulated using the 2-D COM model, and the simulated and experimental reflector responses are plotted in Figure 5. Overall, there is good agreement between the simulated and predicted responses. However, the experimental response is nonsymmetrical, and has wider bandwidth than the prediction. To gain a better understanding of the differences between the simulated and measured reflector responses, an FFT was applied to both. In Figure 6, the time domain responses reveal that the experimental reflector response is shorter than the simulated which Figure 5. 2-D COM simulated and experimental responses of 24 period cosine weighted apodized reflector on YZ LiNbO 3.

.4.3 Experimental response COM simulated response Experimental window COM simulated window.2 Grating Reflection. -. -.2 -.3 -.4.4.42.44.46.48.5.52.54 Time (µs) Figure 6. 2-D COM simulated and experimental cosine weighted reflector time domain responses. Experimental response has shorter duration due to velocity errors. is expected due to the difference in their bandwidths. It is believed that the errors are caused by phase differences of the reflected waves from each track of the reflector. This result underlines the importance of accurate velocity information when designing weighted reflectors. V. OFC TEMPERATURE SENSOR USING COSINE REFLECTORS Implementation of an OFC SAW temperature sensor was accomplished using cosine weighted reflectors. Previous studies have defined a set of orthogonal frequencies for uniform weighting [,2]. Using a similar approach, orthogonal frequencies can be defined for several window types. Orthogonal frequencies for cosine weighting are defined as f n = 2 n τ -. Using this definition, a cosine weighted OFC temperature sensor was designed for operation on YZ LiNbO 3. The sensor was implemented using three cosine weighted reflectors as shown in Figure 7, and occupied a 24.5% fractional bandwidth centered at 25 MHz. The device was simulated using the 2-D COM model and verified by experiment using devices fabricated on YZ LiNbO 3. The simulated and experimental responses have been transformed into the time domain using an FFT, and are shown in Figure 8. The plot shows the cosine weighted reflector responses from either side of the transducer, and there is good agreement between measured and predicted responses. The experimental response was then applied to a simulated transceiver, and a compressed pulse response from one reflector bank is shown in Figure 9. The autocorrelation of the matched filter is also shown for comparison. Note the two pulses located one half chip length away on either side of the compressed pulse. These pulses are undesired, and subsequent calculations have shown that the level of these responses can be significantly reduced by implementing the device using in phase and quadrature channels along with PN coding. Lastly, the OFC sensor was tested over temperature Figure 7. Schematic drawing of three chip cosine weighted OFC sensor. Identical weighted reflector banks are placed on either side of wideband transducer. Free space delays are not drawn to scale. Magnitude (db) 2 25 3 35 4 Experimental COM Simulation 45.8.2.4.6.8 2 2.2 Time (us) Figure 8. Simulated and experimental time domain responses of three chip OFC SAW temperature sensor. Sensor implemented using cosine weighted apodized reflectors. between ºC and ºC, and the resulting compressed pulse responses are shown in Figure. Using an adaptive matched filter, the pulse amplitude remains constant as temperature is varied. VI. CONCLUSION This paper demonstrated the use of weighted reflectors in orthogonal frequency coded SAW tags and sensors. A review of OFC was given, and its application to SAW tags and sensors was demonstrated. The importance of accurate SAW velocity information was discussed, and velocity profiles were extracted for YZ LiNbO 3. A description of the extraction technique was given, and a discussion of the results was provided. Reflector withdrawal weighting and apodization were both considered as weighting techniques, and apodization was chosen as the optimal method for weighting OFC SAW devices. A description was given of apodized reflector design criteria including correcting errors due to velocity perturbations. The final apodized reflector design was optimized to reduce group delay variations across the beam aperture and ensure phase coherence of the reflected waves from each uniform section of the reflectors. In addition, a two dimensional COM model was developed in order to properly simulate the apodized reflector designs. The apodization design criteria were then applied to cosine weighted reflectors, and experimental devices

The 25 MHz OFC SAW sensor was fabricated on YZ LiNbO 3, and its swept frequency response agreed well with 2-D COM model predictions. The sensor was then tested over temperature and the responses were applied to a simulated transceiver which uses an adaptive matched filter. Using the adaptive matched filter, results were given showing that the compressed pulse amplitudes remained constant as temperature varied. The use of weighted reflectors for OFC sensors has been introduced, but additional work is warranted. Phase matching across the beam is difficult and more details need investigation to compensate exactly. Other weighting windows and embodiments may also show advantages for differing sensor applications. Figure 9. Compressed pulse response. ACKNOWLEDGMENT The authors wish to thank NASA for supporting Derek Puccio through the Graduate Student Researcher s Program (GSRP). We also wish to thank the Florida Solar Energy Center (FSEC) and NASA for supporting this work. Our NASA point of contact is Dr. Robert Youngquist. Figure. Cosine weighted OFC SAW sensor compressed pulses. Temperature varied between ºC and ºC. Adaptive matched filter yields uniform pulse amplitude as temperature varies. fabricated on YZ LiNbO 3 were compared with simulated responses generated by the 2-D COM model. The simulated and experimental reflector responses were very similar; however, subsequent analysis revealed small errors due to phase differences of the reflected waves from each uniform track of the reflector. These errors revealed the importance of accurate material velocity information when designing weighted reflectors. Lastly, an OFC SAW temperature sensor was implemented using cosine weighted reflectors. REFERENCES [] D.C. Malocha, D. Puccio, D. Gallagher, Orthogonal frequency coding for SAW device applications, Proc. IEEE International Ultrasonics Symposium, 24, pp. 82-85. [2] D. Puccio, D.C. Malocha, D. Gallagher, J. Hines, SAW sensors using orthogonal frequency coding, Proc. IEEE International Frequency Control Symposium, 24, pp. 37-3. [3] P.V. Wright, Modeling and experimental measurements of the reflection properties of SAW metallic gratings, Proc. IEEE International Ultrasonics Symposium, 984, pp. 54-63. [4] T. Omori, J. Akasaka, M. Arai, K. Hashimoto, M. Yamaguchi, Optimisation of weighted SAW grating reflectors with minimised time delay deviation, Proc. IEEE International Frequency Control Symposium and PDA Exhibition, 2, pp. 666-67. [5] W.J. Tanski, SAW resonators utilizing withdrawal weighted reflectors, IEEE Transactions on Sonics and Ultrasonics, Vol. 26, No. 6, pp. 44-4, November 979. [6] P.D. White, R.F. Mitchell, R. Stevens, P. Moore, M. Redwood, Synthesis and design of weighted reflector banks for SAW resonators, Proc. IEEE International Ultrasonics Symposium, 978, pp. 634-638.