Improved Adaptive Time Delay Estimation Algorithm Based on Fourth-order Cumulants

Similar documents
Design of Two-Channel Low-Delay FIR Filter Banks Using Constrained Optimization

Figure 1 Siemens PSSE Web Site

11.5 MAP Estimator MAP avoids this Computational Problem!

Lecture 10. Erbium-doped fiber amplifier (EDFA) Raman amplifiers Have replaced semiconductor optical amplifiers in the course

RADIATION THERMOMETRY OF METAL IN HIGH TEMPERATURE FURNACE

Optimization of Microphone Array Geometry with Evolutionary Algorithm

Analysis the Transient Process of Wind Power Resources when there are Voltage Sags in Distribution Grid

Hybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification

Estimating Conditional Mean and Difference Between Conditional Mean and Conditional Median

Efficiency Optimal of Inductive Power Transfer System Using the Genetic Algorithms Jikun Zhou *, Rong Zhang, Yi Zhang

High Precision Feedback Control Design for Dual-Actuator Systems

Operational transconductance amplifier based voltage-mode universal filter

Enhancement of Robust Tracking Performance via Switching Supervisory Adaptive Control

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Research on DC resistivity for an arbitrarily anisotropic earth using circular scanning measurement

Liang Peng, Qiwei Yao Estimating conditional means with heavy tails

ROOT LOCUS. Poles and Zeros

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS

Robust Expectation-Maximization Algorithm for Multiple Wide-Band Acoustic Source Localization in the Presence of Non-Uniform Noise Variances

Hybrid models for hardware-in-the-loop simulation of hydraulic systems Part 2: experiments

Systems Analysis. Prof. Cesar de Prada ISA-UVA

Simulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System

Bandwidth expansion of a pressure control system for pneumatic anti-vibration apparatuses in presence of dead time

Midterm 3 Review Solutions by CC

PERFORMANCE MEASURES FOR SPARSE SPIKE INVERSION VS. BASIS PURSUIT INVERSION

A Search Space Optimization Technique for Improving Ambiguity Resolution and Computational Efficiency

p and transverse momentum p BF

GNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase

Active Disturbance Rejection Control for an Electro-statically Actuated MEMS Device

Practical issues of reverse time migration: true-amplitude gathers, noise removal and harmonic-source encoding

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis

Feedback Control Systems (FCS)

Design By Emulation (Indirect Method)

Iterative Decoding of Trellis-Constrained Codes inspired by Amplitude Amplification (Preliminary Version)

Numerical Simulation of Triaxial Compression Stress Paths Tests for Unsaturated Soil

skipping section 6.6 / 5.6 (generating permutations and combinations) concludes basic counting in Chapter 6 / 5

USEFUL TECHNIQUES FOR FIELD ANALYSTS IN THE DESIGN AND OPTIMIZATION OF LINEAR INDUCTION MOTORS

Narrow-Scope Searching FrFT Method Applied in the Elimination of Sonar Ranging Error

Lecture 6. Erbium-doped fiber amplifier (EDFA) Raman amplifiers Have replaced semiconductor optical amplifiers in the course

Signaling over MIMO Multi-Base Systems: Combination of Multi-Access and Broadcast Schemes

Temperature Sensitive Paint Heat Flux Measurements in Hypersonic Tunnels

Generation of Polarization Entangled Photon Pairs in a Planar Waveguide Divya Bharadwaj 1,* and K Thyagarajan 1 1

Fluid Detection in Tight Gas Sand from the Seismic Data

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Chapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem

Lecture 8 - SISO Loop Design

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems

Fuzzy Model of Human s Performance for Guarding a Territory in an Air Combat

Determination of the local contrast of interference fringe patterns using continuous wavelet transform

Fast Lithography Image Simulation By Exploiting Symmetries in Lithography Systems

Discovery Mass Reach for Excited Quarks at Hadron Colliders

AN ADAPTIVE SIGNAL SEARCH ALGORITHM IN GPS RECEIVER

OVERSHOOT FREE PI CONTROLLER TUNING BASED ON POLE ASSIGNMENT

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis

ON THE UNIQUENESS OF MEROMORPHIC FUNCTIONS SHARING THREE WEIGHTED VALUES. Indrajit Lahiri and Gautam Kumar Ghosh

Internal Model Control

Chapter Landscape of an Optimization Problem. Local Search. Coping With NP-Hardness. Gradient Descent: Vertex Cover

696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of

Analysis of Feedback Control Systems

Lecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)

Lecture 10 Filtering: Applied Concepts

UNSUPERVISED LEARNING WITH IMBALANCED DATA

White Rose Research Online URL for this paper: Version: Accepted Version

LTV System Modelling

Some Approaches to the Analysis of a Group of Repeated Measurements Experiment on Mahogany Tree with Heteroscedustic Model

Digital Control System

What lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine?

USING NONLINEAR CONTROL ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS

An estimation approach for autotuning of event-based PI control systems

Massachusetts Institute of Technology Dynamics and Control II

5.5 Application of Frequency Response: Signal Filters

Reliability of Cu Pillar Bumps for Flip-Chip Packages with Ultra low-k Dielectrics

Assessment of Performance for Single Loop Control Systems

The Winding Path to RL

Microwave magnetoelectric fields: An analytical study of topological characteristics

Price Protection with Consumer s Policy Behavior Beibei LI, Huipo WANG and Yunfu HUO

HIGHER-ORDER SPATIAL FDTD SCHEMES FOR EM PROPAGATION IN DISPERSIVE MEDIA

A Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking

MODELING OF NEGATIVE INFLUENCES AT THE SIGNAL TRANSMISSION IN THE OPTICAL MEDIUM. Rastislav Róka, Filip Čertík

A NEW YIELD CRITERION FOR ORTHOTROPIC SHEET METALS UNDER PLANE-STRESS CONDITIONS

Jul 4, 2005 turbo_code_primer Revision 0.0. Turbo Code Primer

Using Maple to Evaluate the Partial Derivatives of Two-Variables Functions

Linearteam tech paper. The analysis of fourth-order state variable filter and it s application to Linkwitz- Riley filters

Research Article FPGA Realization of Sensorless PMSM Speed Controller Based on Extended Kalman Filter

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis

DESIGN AND ANALYSIS OF FLEXIBLE STRUCTURES WITH PARAMETRIC UNCERTAINTIES FOR ACTIVE VIBRATION CONTROL USING ANSYS

Do Dogs Know Bifurcations?

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Simulation of Wound Rotor Synchronous Machine under Voltage Sags

CHAPTER 5. The Operational Amplifier 1

Übung zu Globale Geophysik I 05 Answers

Depth Discrimination for Low-Frequency Sources Using a Horizontal Line Array of Acoustic Vector Sensors Based on Mode Extraction

Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material

Efficient Methods of Doppler Processing for Coexisting Land and Weather Clutter

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

Comparison of Meteorological Radar Signal Detectability with Noncoherent and Spectral-Based Processing

A New Smith Predictor and Fuzzy Immune Control for Hybrid Networked Control Systems

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

Multi-dimensional Fuzzy Euler Approximation

Transcription:

Available online www.jocr.com Journal of hemical and Pharmaceutical Reearch, 016, 8(5):889-894 Reearch Article ISSN : 0975-784 ODEN(USA) : JPR5 Imroved Adative Time Delay Etimation Algorithm Baed on Fourth-order umulant Ling Tang and Yumei hen ollege of Automation and Electronic Information, Sichuan Univerity of Science & Engineering, Zigong 64000, hina ABSTRAT Time-delay etimation i widely ued in onar and radar alication. In order to etimate time delay accurately in variable ignal to noie ratio(snr), an imroved algorithm of adative time delay etimation() baed on fourth-order cumulant ha been rooed in thi aer. The ingle adative filter in the conventional model i earated into two adative ubunit connected in cacade, one i to track the delay, the other to track the SNR, o a to obtain the accurate ignal time delay etimation in the low and variable SNR. The theoretical analyi and imulation reult confirm the validity and feaibility of the rooed algorithm. Keyword: fourth-order cumulant, time delay, adative filter, variable SNR INTRODUTION Time Delay Etimation() a an active reearch in ignal roceing i widely ued in the field of radar, onar, communication and o on. The common method are hae method, double ectrum method, correlation method, adative filter arameter model method and o on. With the continuou develoment of ignal roceing method, a variety of algorithm are introduced into time delay etimation, which imrove the reciion and convergence, reduce the comutation. Becaue of blind Gau roerty high order accumulation i widely ued in time delay etimation. Liu Ying et al.[1] rooed an adative time delay etimation algorithm baed on four-order cumulant(fo-lms), which can effectively ure the influence of correlated Gau noie, but in the cae of low SNR the etimated reult are not atifactory. Aiming at thi roblem, a new algorithm of time delay etimation algorithm FO-E i rooed in the literature[], which can ure the imact of related or not related to Gau noie, and obtain the accurate time delay etimation of non Gau ignal in lower SNR. However when SNR change, the erformance of the ytem i greatly reduced. In thi aer, an imroved adative time delay etimation algorithm i rooed. The ingle adative filter i earated into two adative ubunit connected in cacade, one i to track the delay, the other to track the SNR, o a to imrove the erformance in the variable SNR. MATHEMATIAL MODEL onider one ignal beamed from a ditant ource, with interference and noie in the channel. Thu the ignal arameter model received on two atially earated enor[] i: x ( ( + w1 ( y ( a( k D) + w ( = (1) = () 889

Ling Tang and Yumei hen J. hem. Pharm. Re., 016, 8(5):889-894 In the formula ( i non Gau and zero mean tationary random ignal, ( k D) i the time delay of (, a i the fading factor that atified a (0,1], uually a = 1. w1 ( and w ( are the correlative zero mean Gau noie which i indeendent with (. By convolution theory: n= k D) = in n D) ( k n) ( () Where in v) = in(π v), then relace with a large oitive integer ( > D πv truncation error, we can get the amling value of the function. ), and ignoring the Inert the formula () into (): ( = a in ( k i) + n ( k = a in [ k i) n1 ( k i) ] + n ( y ) (4) The uroe of time delay etimation i to ure the Gau noie effectively by uing the limited obervation ignal ) k and y (, o a to etimate the time delay D. ALGORITHM ANALYSIS For tationary random roce with zero mean ( defined a follow: ( τ,0,0) cum[ k + τ ), ] x and y (, the FO of the random variable can be = (5) Define cro four order cumulant(fo) of the random variable and y( ( τ,0,0) cum[ y( k + τ ), ] = (6) Inert the formula (4) into (6), we get ( τ,0,0) = cum a in [ k i + τ ) w1 ( k i + τ )] + w ( k + τ ), (7) According to the roertie that the cumulant can be added and the FO of Gau noie i zero, the equation(7) can be orted out: [ ( k i τ ), )] ( τ,0,0) = a in cum x + k a: (8) That i: ( τ,0,0) = ( τ,0,0) (9) (,0,0) = a in ( τ i,0,0) = a in ( τ i,0,0) τ (10) Aarently the Gau noie contained in ( τ,0,0), ( τ,0,0) i ureed. The formula (10) can be exreed a a vector: = a4 S (11) xy x 890

Ling Tang and Yumei hen J. hem. Pharm. Re., 016, 8(5):889-894 Where, a i the contant, and S = [in D),in + 1 D), L,in D)] xy [ (,0,0), ( + 1,0,0), L, (,0,0 ] T = ) 4x = (0,0,0) (1,0,0) M (,0,0) ( 1,0,0) (0,0,0) M ( 1,0,0) L L O L T (,0,0) ( + 1,0,0) M (0,0,0) In the ue of adative time delay etimation method to get the time delay, the otimal weight coefficient hould be [4]: a σ T [in, i =, + 1, L,0, L, 1, ] σ + σ w (1) Baed on the maximum value of the weight coefficient, i obtained. However thi method can get the maximum error in the cae of low noie, which lead to inaccurate etimation. In thi cae, the original adative filter coefficient can be decouled into two art of g and in c, with uing in filter, at the ame time adding the gain g for the iteration to track the factor ( a σ ) ( σ + σ ) w o a to obtain the otimal weight coefficient. In thi way the etimation can be obtained directly without the incorrect maximum value. The ytem model i howed in Figure 1. y( FO calculation FO calculation (τ,0,0) g (τ,0,0) FIR filter with inc amling xzxx (τ,0,0) + e( + - Figure 1: Sytem model (τ,0,0 In Figure 1 ignal ) (τ,0,0 ignal ) : xzxx xzxx after an adative ower factor g ) and an adative FIR filter into the outut ) ) ( τ,0,0) = g in ( τ i,0,0) Adative error function i defined a: If = [ ( τ,0,0) ( τ,0,0)] xzxx J (14) τ =, + 1, L, τ, the error function i (1) 891

Ling Tang and Yumei hen J. hem. Pharm. Re., 016, 8(5):889-894 ) ) J = g in ( τ i,0,0) ( τ,0,0) τ = ) ) = ( g S ) ( g S ) T 4 x xy 4 x xy (15) Make the gradient of the error function to the arameter D ) and g ) are zero, the otimal arameter vector can be olved. The iterative roce of the otimal arameter vector i realized by uing the teeet decent method in the adative roce[5]. Each time the value of the current vector i changed by a negative gradient vector: ) ) µ d J D( k + 1) = D( ) ) g( D (16) µ d and g tability and eed. ) ) J g k + 1) = g( µ g ) g ( (17) µ i reectively the convergence factor of D ) ( and g ) (, which i ued to adjut the adative We can ee that the initial value of delay time and gain in the iteration are needed, initial value (0) [6] mut D ˆ (0), g ) (0) atifie 1.45 D D + 1. 45 delay time [7]. hould be between 0 and 1, o a to be converged to the real SIMULATIONS AND ANALYSIS The ignal ource i a tationary ignal with zero mean value generated by a non Gau random ignal generator. The real time delay i D=, initial time delay etimation D=, fading factor a i 1. The oitive order of the filter i =10, the amling oint are 000, and the noie i the correlated Gau white noie. The erformance of FO-E, FO-LMS and the imroved algorithm are comared uing MATLAB in the imulation. Exeriment 1 Time delay etimation under different SNR Aume SNR i 5dB or -5dB, the reult of different algorithm are hown in Figure. In the high SNR (SNR=5dB), three algorithm can quickly converge to the true value of time delay. But with SNR decreaing, the detection ability of the algorithm have been reduced. FO-LMS can not accurately etimate the time delay, comared with FO-E algorithm the rooed method obtain more accurate etimation ˆD.5 0 500 1000 1500 000 loo.5.5 0 500 1000 1500 000 loo 0 500 1000 1500 000 loo (a) FO-LMS(SNR=5dB) (b) FO-E(SNR=5dB) (c) The rooed algorithm(snr=5db) 89

Ling Tang and Yumei hen J. hem. Pharm. Re., 016, 8(5):889-894.5 0 500 1000 1500 000 loo.5.5 0 500 1000 1500 000 loo 0 500 1000 1500 000 loo (d) FO-LMS(SNR=-5dB) (e) FO-E(SNR=-5dB) (f) The rooed algorithm(snr=-5db) Figure : Performance comarion under correlated noie Exeriment Time delay etimation invariable SNR When SNR reduced from 0dB to -10dB, comare the delay etimation. The firt 1000 data with SNR=0dB, the lat 1000 data with SNR=-10dB.The otimal olution of the gain factor in FO-ETDGE method i g * = 0. 09, when SNR=-10dB the otimal olution i g * = 0. 5. In FO-E method, g i fixed to 1. In Figure we comare their erformance in variable SNR. Obviouly, FO-LMS algorithm get incorrect, the effect of FO EE i not good, which i becaue that any change in the filter arameter i conidered to be a change in the time delay, o can not be tracked in time. In contrat, the effect of the imroved algorithm i the bet, even in the cae of a large change to SNR, it can alo be converged to the etimated value..5.5.5 0 500 1000 1500 000 0 500 1000 1500 000 loo loo (a) FO-LMS (b) FO-E (c) The rooed algorithm Figure : in variable SNR ONLUSION 0 500 1000 1500 000 loo In thi aer an imroved adative time delay etimation algorithm baed on four order cumulant i rooed. The ingle adative filter i earated into two adative ubunit connected in cacade, which making the adative roce of time delay and ignal to noie ratio to be earated, avoiding the oor erformance of delay etimation in variable SNR. Simulation reult demontrate the effectivene of the rooed algorithm. Acknowledgment Thi work wa uorted by roject of Sichuan Provincial Deartment of Education(1ZB018) and roject of Artificial Intelligence Key Laboratory Of Sichuan Province(01RYY0). 89

Ling Tang and Yumei hen J. hem. Pharm. Re., 016, 8(5):889-894 REFERENES [1] LIU Ying, WANG Shu-xun, WANG Ben-ing, Journal of Sytem Simulation, 00, 14(6), 700-70 []LI ong-ying, WANG Jian-ying, YIN Zhong-ke, ZHANG Jiang-li, Journal of the hina railway ociety, 006, 8(6), 55-58 [] ZHANG Duan-jin ZHANG Zhong-hua GUO Jian-jun ZHANG De-jing, Journal Of Zhengzhou Univerity(Engineeringcience), 010,1(1), 10-106 [4]H..So and P.. hing. A novel contrained time delay etimator. Proc. Of ISP 9. Beijing hina. 199, 188-19 [5]XIONG Qiu-ben, WANG Jiang, YANG Jing-hu, Electronic Warfare Technology, 011, 6(), 10-15 [6]Faramarz Fekri, Mina Sartiii and Reell M. Merereau. IEEE Tranaction on Signal Proceing. May 005, 5(5) [7]H..So and P..hing. Performance analyi of ETDGE an efficient and unbiaed TDOA etimator. IEEE Proc.-Radar Navig. December 1998, 145(6), 5-0 894