HYBRID FAST ALGORITHM FOR S TRANSFORMS
|
|
- Abigail Fisher
- 5 years ago
- Views:
Transcription
1 8th European Signal Processing Conerence (EUSIPCO-) Aalborg Denar August -7 HYBRID FAST ALGORITHM FOR S TRASFORMS Soo-Chang Pei Jian-Jiun Ding Pai-Wei Wag and Wen Fu Wang Departent o Electrical Engineering ational Taiwan University o. Sec. Roosevelt Rd. 67 Taipei Taiwan R.O.C TEL: Fax: Eail: pei@cc.ee.ntu.edu.tw djj@cc.ee.ntu.edu.tw d99@ntu.edu.tw r9696@ntu.edu.tw ABSTRACT The S transor is useul in tie-requency analysis. In this paper we propose a hybrid algorith to ipleent it adaptively. Since the window size o the S transor varies with it is reasonable to use dierent algorith or dierent requency to ipleent it. In this paper we use the sub IDFT algorith in the low requency region and the sectioned convolution algorith in the high requency region to ipleent the S transor. Fro siulation our algorith reduces % o the coputation tie and uch iproves the eiciency o the S transor.. ITRODUCTIO The S transor [][] is deined as: ( τ t) iπ t S( τ ) = x()exp[ t ] e dt. () π It is useul or tie requency analysis. It is a odiication o the short tie Fourier transor (STFT also naed as the Gabor transor) []: iπ t G( τ ) = x( t)exp[ ( τ t) / ] e dt. () π Copared with the STFT the Gaussian as o the S transor varies with. This aes the S transor have adaptive resolution (high resolution in the low requency region and low resolution in the high requency region). The S transor is useul in signal decoposition [] power quality analysis [] seisoy [6] bioedical engineering [7][8] atospheric science [9] and other applications related to tie-requency analysis. In [] an algorith or ipleenting the S transor in the Fourier doain was introduced. Since the S transor in () can be viewed as the ollowing convolution operation iπ t St = xte () g t ro along t-axis () where g () t = exp[ t / ] () π iπ t FT[ xte ] Xu FT g t = e () the S transor can be ipleented as ollows []: = + () π u Step : Peror the discrete Fourier transor (DFT) or the saples o x(t): ( ) [ ] π X Δ = x n e where ΔΔ t = / () n= x [ (( n )) ] = x( nδ t ) () is the odulus operation and is soe integer that should be larger than the nuber o saple points o x(t). ote that i the tie duration o x(t) is [T T ] and T < then it is possible that n < and we should use ((n)) to re-arrange the index beore applying the DFT. Step : For each copute the -point inverse discrete Fourier transor (IDFT) as ollows. Then the result will be the S transor o x(t) at τ = nδ t and = Δ. π π j / ( Δt Δ ) (( + ) Δ ) S n X e e. (6) =/ Suppose that has M possible values. Then the above algorith requires one DFT M IDFTs and M ties o ultiplications o exp(π / ). Thus the coplexity is: ( M + ) + M (7) and the average coplexity or each is: ( M + ) + +. (8) M In this paper we ind that the eiciency o the S transor can be urther iproved. First notice that when is large the tie duration o g (t) = exp[ (τ t) /] is very sall. It eans that () is a convolution o x(t)exp(jπt) with a very short unction. The convolution o a long unction and a short unction can be ipleented by the sectioned convolution. By contrast when is sall since FT[g (t)] = exp[π u / ] decays very ast the bandwidth o g (t) is very sall. In this condition the input o the IDFT in (6) is very short and we can use several shorter length IDFTs instead o the -point IDFT in (6) to ipleent the S transor. Thereore one can ipleent the S transor adaptively according to the value o. (i) When is large the sectioned convolution algorith can be applied. See Section. (ii) When is sall the algorith o sub IDFT can be applied. See Section. EURASIP ISS
2 (iii) Even in the iddle requency region using both the sectioned convolution and the sub IDFT algoriths are ore eicient than using the original algorith. See Figs. and. We also peror several siulations in Section to show that the proposed hybrid algorith can uch iprove the eiciency o the S transor.. SECTIOED COVOLUTIO I HIGH FREQUECY REGIO Suppose that the nuber o sapling points or x(t) is : sapling xt () xn ( Δ t ) where n = n n +.. n +. (9) In theory the unction g (t) in () has ininite duration. However one can give a threshold: threshold =.ax ( g ( t) ) () and ignore the case where g (t) < threshold. It aes g (t) becoe a tie-liited unction with the tie duration o: B < t < B where B = / =.769 / () and the nuber o sapling points in the duration is.769 L = + Δ t () where eans the rounding down operation. Thereore () becoes the linear convolution o an -length sequence x(nδ t ) and an L-length sequence g(nδ t ). However when is large the tie-duration o g (t) is very short and the value o L in () is also very sall. In this case instead o () one can divide x(nδ t ) into several shorter sequences and use the sectioned convolution to ipleent the convolution o x(nδ t ) and g(nδ t ). That is Step : Divide x(nδ t ) into B subsections. The length o each subsection is where = /B and eans the rounding up operation. That is π ( n+ p+ n) xp [ n] = x ( n + p + n) Δt e () where n = p =. B and exp(jπ(n +p +n)/) coes ro substituting t = (n +p +n)δ t and = Δ into exp(jπt) in (). Moreover as in () = /Δ t Δ should be satisied. Step : Then we use the -point DFT / IDFT pair to ipleent the convolution in () instead o the -point DFT / IDFT pair used in () and (6) where should satisy > + L ( /B) + L. () This step can be accoplished by the ollowing sub steps: Step -: [ ] [ ] π X p = xp n e n= or p =. B. () π π j p = p = Step -: [ ] [ ] T n X e e or p =. B. (6) Step -: Qp [ n] = Tp [ n] or n = +L L = (L)/ Qp [ n] = Tp [ n+ ] or n = L L (7) and Q p [n] = otherwise. Then Q p [n] is near to Qp [ n] xp [ n] Δtg ( nδ t) (8) where * eans the linear convolution and g (t) is deined in (). Steps - and - can be proved as ollows π j n Δ t DFT { Δ } exp[ ] tg nδ t = e Δt π n ( DFT eans the -point DFT) exp( t / ) e dt π ut π t nδt u Δu ΔtΔ u= / = exp( π u / ) (ro ()) = exp( π ) Δ Δt π (ro u = Δ u = /Δ t and = Δ ) = exp( / ) (ro Δ t Δ = /). (9) Thereore ro (6) Tp [ n] = IDFT DFT [ ] x p n exp π { } { { [ ]} { }} IDFT DFT x n DFT Δg nδ. () p t t Thus T p [n] is the circular convolution o x p [n] and Δ t g (nδ t ) and Q p [n] in (7) is their linear convolution. Step : Su the linear convolution results o each subsection. Then we obtain the result o the S transor (denoted by S(nΔ t Δ )): (a) S( ( n + p + n) Δt Δ ) Qp [ n] when (i) L n < L p =. B (ii) n < L p = or (iii) L n < p = B where L = (L)/ (b) S( ( n + p + n) Δt Δ ) Qp [ n] + Qp [ n+ ] when n < L p =. B (c) S( ( n + p + n) Δt Δ ) Qp [ n] + Qp+ [ n ] when L n < p =. B. () Step can be proven ro the act that π n B x[ ne ] = xp [ nn p] () π n ( t ) [ ] t ( t) S nδ Δ x n e Δ g nδ (ro ()) 78
3 Table Coplexities o the original algorith and the proposed algorith in Section or ipleenting the S transor. Here Δ t = Δ =. and t. requency (/) + (relect the coplexity o the original algorith) [B( /B+L)/] ( /B+L) + +BLB (relect the coplexity o the proposed algorith) B = B = B = 6 B = 9 = = = = B ( Δt Δ ) p [ ] Δt ( Δt) S n x n n p g n B Qp n n p [ ] (ro (8)). () Then we analyze the coplexity o the proposed algorith. Suppose that there are M possible values o. Then the ultiplication o exp(π /( )) and the -point IDFT in (6) should be coputed M B ties. Since [ ] [ ] π ( n+ p) X X e thus in () and () we only have to calculate x p [n] and X p [] or = /. Thereore there are / ultiplications in () where is the nuber o input sapling points (see (9)) and the -point DFT in () should be perored B/ ties. The product o exp(jπ(n +p )/ ) in () can be erged with the product o exp(π /( )) in (6). Thereore the total coplexity o the proposed algorith is B + M B MB + + B MB ( / B L ) L () + / + MB / B+ L p + / = p + () and the average coplexity or each o the proposed algorith is near to B ( / ) + MB B L L M /( M) + B / B+ L B / B+ L L BL B (6) where is the nuber o input sapling points B is the nuber o sections and L is deined in () when M is suicient large. In coparison or the original algorith in [] the average coplexity is as in (8). In the case where L << i we choose B and the constraint that < is satisied then /B + L << and (6) will be uch less than (8). For exaple i Δ t and Δ are ixed to. ro () = and the value o (/) + in (8) is near to 66. I the duration o the input signal is t then the value o is. In Table we show the value o (6)) or dierent B and (ote that aects the value o L (see ()) and hence the value o (6)). It shows that the proposed algorith indeed has less coplexity than the original algorith. Especially when is large using larger B (i.e. dividing the input into ore sections) will be even ore eicient.. SUB IDFT ALGORITHM I LOW FREQUECY REGIO In previous section we described that when is large the sectioned convolution algorith is helpul or iproving the eiciency o the S transor. In this section we show that when is sall the sub DFT algorith will be ore helpul or iproving the eiciency. ote that when is sall since = Δ the value o in (6) is also sall. It eans that the exponential ter exp(π / ) decays very ast. I we ignore the case where exp(π / ) < then (6) can be rewritten as: ( Δt Δ ) = (( + ) Δ ) = π π j S n X e e (7) where.96 = = π. (8) Δt (7) is the -point IDFT whose input has + points. I is very sall + is uch less than. For exaple when Δ t = Δ =. and = the value o = (/Δ t Δ ) is but the value o in () is only. Since + = 6 << in this case using the -point IDFT to ipleent (7) directly would be very ineicient. Instead we can use the sub IDFT algorith described as ollows to ipleent (7) and iprove the eiciency. First ind an integer L such that (a) L + and (b) B = / L is an integer. (9) Then (7) can be re-expressed as: π π j ( ab + c) = π π π j c j a L (( + ) Δt Δ ) = (( + ) Δ ) S ab c X e e = ( ) = { X ( + ) Δ e e } e () where c = ((n)) B a = (nc)/b and n = ab + c. () Then ro the act that the coplexity o B ties o the (/B )-point IDFTs is less than that o one -point IDFT we can use the ollowing way to ipleent the -point IDFT instead o using (6) directly to iprove the eiciency. (Step in () is unchanged.) (Modiication or Step ): First calculate = ( + Δ ) π j c Xc X ( ) e e π or c = B. () Then peror the L -point IDFT or each o X c () (c = B ): 79
4 . x -axis coplexity. Line Line Line Line High requency region: (Using the sectioned convolution with larger B) Middle requency region: (Using the sectioned convolution with saller B) Fig. The coplexities o using three dierent algoriths or ipleenting the S transor. Line : the original algorith Lines and : using the sectioned convolution algorith in Section (B = or line and B = or Line ) Line : using the sub IDFT algorith in Section. (( ) t ) L c S ab + c Δ Δ = IDFT X. π j a L Xc e = = () ote that i there are M possible values or then in () the L -point IDFT should be perored B M ties. Thereore the coplexity o the odiied algorith is: (ro Step ) + L MB ( + ) (ro ()) + MB L (ro ()) + M + M L. () Here we use the act that B ( +) B L = (ro (9)). Thus the average coplexity or each o the proposed sub IDFT algorith is: + L+ L+ () M i M is suicient large. Copared with (8) since L = /B << it is obvious that (/) L is uch less than (/). Thereore using the proposed sub IDFT algorith can indeed iprove the eiciency o the S transor especially in the condition where is sall (ote that since L = + and is proportional to (see (8)) i is sall the value o L in () is also sall). In Fig. we copare the eiciencies o the three algoriths or ipleenting the S transor. (i.e. original algorith the sectioned convolution algorith in Section and the sub IDFT algorith in Section ) Here Δ t = Δ =. (i.e. = (/ Δ t Δ ) = ) (6) and t. Line is (/) + which is the coplexity o the original algorith. It is invariant with. Lines and are [B( /B+L)/] ( /B+L) + +BLB when B = and respectively. They are the coplexities when using the sectioned convolution algorith in Section. Line is (/) L + which is the coplexity o the sub IDFT algorith described in this section. Low requency region: Using the sub IDFT algorith Fig. Using dierent algoriths in dierent regions to ipleent the S transor where B eans the nuber o sections when using the sectioned convolution algorith. Fro Fig. it is obvious that when is sall using the sub IDFT algorith proposed in this Section will be ore eicient or ipleenting the S transor. When is in the iddle region it is proper to use the sectioned convolution algorith in Section with saller B to ipleent the S transor. In the high requency region it is proper to use the sectioned convolution algorith with larger B as in Fig.. Thereore we can use the hybrid algorith to ipleent the S transor. For dierent the algorith or ipleenting the S transor is also dierent.. SIMULATIOS t-axis We peror several siulations to copare the eiciencies o the proposed hybrid algorith and the original algorith to ipleent the S transor. There are two input signals: Fig. (a): x(t) = cos(πt) or t x(t) = cos(6πt) or t x(t) = cos(πt) or t (7) Fig. (a): x(t) = cos(πt) or t x(t) = cos(πt) or t x(t) = cos(πt) or t. (8) Their S transors (coputed by the original algorith in (6) and (7)) are shown in Figs. (b) and (b). Then we use the proposed hybrid algorith to copute the S transors. In Fig. (b) we use the sub IDFT algorith or < and use the sectioned convolution with or. We plot the results in Figs. (c) and (c). The results are all the sae as those o using the original algorith. It proos that the proposed hybrid algorith is valid. Then we show the coputation tie in Table. The results show that the proposed hybrid algorith saves over % o the coputation tie and is uch ore eicient than the original algorith. 7
5 - (b) (a) S transor - (c) Fast algorith o the S transor Table Coparing the coputation ties o the original and the proposed hybrid algoriths or ipleenting the S transor Coputation tie Signal in Fig. Signal in Fig. Original algorith Hybrid algorith Iproveent. sec.7 sec.% 8. sec.7 sec.8% REFERECES requency - - requency - - [] R. G. Stocwell L. Mansinha and R. P. Lowe Localization o the coplex spectru: the S transor IEEE Tran. Signal Processing vol. no. Apr Fig. (a) The input signal x(t) (deined in (7)) (b) the S transor o x(t) coputed by the original algorith and (c) the S transor o x(t) coputed by the proposed hybrid algorith. [] S. Ventosa C. Sion M. Schiel J. J. Dañobeitia and A. Mànuel The S-Transor Fro a Wavelet Point o View IEEE Trans. Signal Process vol. 6 no. 7 pp July 8. (a) [] M. J. Bastiaans Gabor s expansion o a signal into Gaussian eleentary signals Proc. IEEE vol. 68 pp [] L. J. Stanovic T. Thayaparan and M. Daovic Signal Decoposition by Using the S-Method EUSIPCO Sept.. (b) The S transor - (c) Fast algorith o the S transor [] P. K. Dash B. K. Panigrahi and G. Panda Power quality analysis using S-transor IEEE Trans. Power Deliv. vol. 8 pp. 6- Apr.. requency (Hz) requency (Hz) Fig. (a) The input signal x(t) (deined in (8)) (b) the S transor o x(t) coputed by the original algorith and (c) the S transor o x(t) coputed by the proposed hybrid algorith.. COCLUSIOS In this paper we propose a hybrid algorith to iprove the ipleentation eiciency o the S transor. We ind that when is sall it is proper to use the sub IDFT algorith instead o the original ethod to ipleent the S transor. When is large it is proper to use the sectioned convolution algorith and the nuber o sections is increased when grows larger. The proposed algorith uch reduces the coputation tie o the S transor and is very helpul or tie-requency analysis. [6] J. H. Gao W. C. Chen Y. M. Li and F. Tian Generalized S-transor and seisic response analysis o thin interbeds Chinese J. Geophysics Chinese Ed. vol. 6 pp. 6- July. [7] G. Livanos. Ranganathan and J. Jiang Heart sound analysis using the S transor IEEE Coputers in Cardioy pp [8] S. Assous A. Hueau M. Tartas P. Abraha and J. P. L Huillier S-Transor Applied to Laser Doppler Flowetry Reactive Hypereia Signals IEEE Trans. Bioed. Eng. vol. pp.-7 June 6. [9] Y. Portnyagin E. Merzlyaov C. Jacobi. Mitchell H. Muller A. Manson W. Singer P. Hoann and A. Fachrutdinova Soe results o S-transor analysis o the transient planetary-scale wind oscillations in the lower therosphere Earth Planets and Space vol. pp
N-Point. DFTs of Two Length-N Real Sequences
Coputation of the DFT of In ost practical applications, sequences of interest are real In such cases, the syetry properties of the DFT given in Table 5. can be exploited to ake the DFT coputations ore
More informationA method to determine relative stroke detection efficiencies from multiplicity distributions
A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,
More informationA practical approach to real-time application of speaker recognition using wavelets and linear algebra
A practical approach to real-tie application o speaker recognition using wavelets and linear algebra Duc Son Pha, Michael C. Orr, Brian Lithgow and Robert Mahony Departent o Electrical and Coputer Systes
More informationFourier Series Summary (From Salivahanan et al, 2002)
Fourier Series Suary (Fro Salivahanan et al, ) A periodic continuous signal f(t), - < t
More informationQ ESTIMATION WITHIN A FORMATION PROGRAM q_estimation
Foration Attributes Progra q_estiation Q ESTIMATION WITHIN A FOMATION POGAM q_estiation Estiating Q between stratal slices Progra q_estiation estiate seisic attenuation (1/Q) on coplex stratal slices using
More informationBEE604 Digital Signal Processing
BEE64 Digital Signal Processing Copiled by, Mrs.S.Sherine Assistant Professor Departent of EEE BIHER. COTETS Sapling Discrete Tie Fourier Transfor Properties of DTFT Discrete Fourier Transfor Inverse Discrete
More informationFast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials
Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter
More informationSPECTRUM sensing is a core concept of cognitive radio
World Acadey of Science, Engineering and Technology International Journal of Electronics and Counication Engineering Vol:6, o:2, 202 Efficient Detection Using Sequential Probability Ratio Test in Mobile
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationStructured Illumination Super-Resolution Imaging Achieved by Two Steps based on the Modulation of Background Light Field
017 nd International Seinar on Applied Physics, Optoelectronics and Photonics (APOP 017) ISBN: 978-1-60595-5-3 Structured Illuination Super-Resolution Iaging Achieved by Two Steps based on the Modulation
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationEfficient Filter Banks And Interpolators
Efficient Filter Banks And Interpolators A. G. DEMPSTER AND N. P. MURPHY Departent of Electronic Systes University of Westinster 115 New Cavendish St, London W1M 8JS United Kingdo Abstract: - Graphical
More informationSexually Transmitted Diseases VMED 5180 September 27, 2016
Sexually Transitted Diseases VMED 518 Septeber 27, 216 Introduction Two sexually-transitted disease (STD) odels are presented below. The irst is a susceptibleinectious-susceptible (SIS) odel (Figure 1)
More informationArithmetic Unit for Complex Number Processing
Abstract Arithetic Unit or Coplex Nuber Processing Dr. Soloon Khelnik, Dr. Sergey Selyutin, Alexandr Viduetsky, Inna Doubson, Seion Khelnik This paper presents developent o a coplex nuber arithetic unit
More informationDECONVOLUTION VERSUS CONVOLUTION A COMPARISON FOR MATERIALS WITH CONCENTRATION GRADIENT
Materials Structure, vol. 7, nuber () 43 DECONVOLUTION VERSUS CONVOLUTION A COMPARISON FOR MATERIALS WITH CONCENTRATION GRADIENT David Raaa Departent o Electronic Structures, Faculty o Matheatics and Physics,
More informationA New Algorithm for Reactive Electric Power Measurement
A. Abiyev, GAU J. Soc. & Appl. Sci., 2(4), 7-25, 27 A ew Algorith for Reactive Electric Power Measureent Adalet Abiyev Girne Aerican University, Departernt of Electrical Electronics Engineering, Mersin,
More informationFingerprint Verification Based on Invariant Moment Features and Nonlinear BPNN
800 International Journal o Control, Autoation, and Systes, vol. 6, no. 6, pp. 800-808, Deceber 008 Fingerprint Veriication Based on Invariant Moent Features and Nonlinear BPNN Ju Cheng Yang and Dong Sun
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationComplexity reduction in low-delay Farrowstructure-based. filters utilizing linear-phase subfilters
Coplexity reduction in low-delay Farrowstructure-based variable fractional delay FIR filters utilizing linear-phase subfilters Air Eghbali and Håkan Johansson Linköping University Post Print N.B.: When
More informationAcoustic Source Localization and Discrimination in Urban Environments
Acoustic Source Localization and Discriination in Urban Environents Manish Kushwaha, Xenoon Koutsouos, Peter Volgyesi and Aos Ledeczi Institute or Sotware Integrated Systes (ISIS) Departent o Electrical
More informationCHAPTER 8 : A PARALLEL TWO DIMENSIONAL WAVELET PACKET TRANSFORM AND ITS APPLICATION TO MATRIX-VECTOR MULTIPLICATION.
HAPTR 8 : A PARALLL TWO DIMNSIONAL WAVLT PAKT TRANSFORM AND ITS APPLIATION TO MATRIX-VTOR MULTIPLIATION. ric GOIRAND Mladen Victor WIKRHAUSR Washington University Departent o Matheatics Brookings Drive,
More informationA remark on a success rate model for DPA and CPA
A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance
More informationFunctions: Review of Algebra and Trigonometry
Sec. and. Functions: Review o Algebra and Trigonoetry A. Functions and Relations DEFN Relation: A set o ordered pairs. (,y) (doain, range) DEFN Function: A correspondence ro one set (the doain) to anther
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More informationTime-frequency plane behavioural studies of harmonic and chirp functions with fractional Fourier transform (FRFT)
Maejo International Journal of Science and Technology Full Paper ISSN 95-7873 Available online at www.ijst.ju.ac.th Tie-frequency plane behavioural studies of haronic and chirp functions with fractional
More informationA MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION
A eshsize boosting algorith in kernel density estiation A MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION C.C. Ishiekwene, S.M. Ogbonwan and J.E. Osewenkhae Departent of Matheatics, University
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationFigure 1: Equivalent electric (RC) circuit of a neurons membrane
Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of
More informationP032 3D Seismic Diffraction Modeling in Multilayered Media in Terms of Surface Integrals
P032 3D Seisic Diffraction Modeling in Multilayered Media in Ters of Surface Integrals A.M. Aizenberg (Institute of Geophysics SB RAS, M. Ayzenberg* (Norwegian University of Science & Technology, H.B.
More informationOptical spectroscopic techniques
What happens i a saple is illuinated by light? (absorbed light) Optical spectroscopic techniques illuinating light saple transitted light Absorptionspectroscopies: UV-VIS, IR László Seller eitted light
More informationOn Constant Power Water-filling
On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives
More informationThis model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.
CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationOn the Mixed Discretization of the Time Domain Magnetic Field Integral Equation
On the Mixed Discretization of the Tie Doain Magnetic Field Integral Equation H. A. Ülkü 1 I. Bogaert K. Cools 3 F. P. Andriulli 4 H. Bağ 1 Abstract Tie doain agnetic field integral equation (MFIE) is
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationAn Improved Particle Filter with Applications in Ballistic Target Tracking
Sensors & ransducers Vol. 72 Issue 6 June 204 pp. 96-20 Sensors & ransducers 204 by IFSA Publishing S. L. http://www.sensorsportal.co An Iproved Particle Filter with Applications in Ballistic arget racing
More informationWarning System of Dangerous Chemical Gas in Factory Based on Wireless Sensor Network
565 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 59, 07 Guest Editors: Zhuo Yang, Junie Ba, Jing Pan Copyright 07, AIDIC Servizi S.r.l. ISBN 978-88-95608-49-5; ISSN 83-96 The Italian Association
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationc 2013 Society for Industrial and Applied Mathematics
SIAM J. MATRIX ANAL. APPL. Vol. 34, No. 3, pp. 1213 123 c 213 Society or Industrial and Applied Matheatics χ 2 TESTS FOR THE CHOICE OF THE REGULARIZATION PARAMETER IN NONLINEAR INVERSE PROBLEMS J. L. MEAD
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationECS455: Chapter 5 OFDM
ECS455: Chapter 5 OFDM 5.4 Cyclic Prefix (CP) 1 Dr.Prapun Suksopong prapun.co/ecs455 Office Hours: BKD 3601-7 Tuesday 9:30-10:30 Friday 14:00-16:00 Three steps towards odern OFDM 1. Mitigate Multipath
More informationPULSE-TRAIN BASED TIME-DELAY ESTIMATION IMPROVES RESILIENCY TO NOISE
PULSE-TRAIN BASED TIME-DELAY ESTIMATION IMPROVES RESILIENCY TO NOISE 1 Nicola Neretti, 1 Nathan Intrator and 1,2 Leon N Cooper 1 Institute for Brain and Neural Systes, Brown University, Providence RI 02912.
More information7. Discrete Fourier Transform (DFT)
7 Discrete ourier Trasor (DT) 7 Deiitio ad soe properties Discrete ourier series ivolves to seueces o ubers aely the aliased coeiciets ĉ ad the saples (T) It relates the aliased coeiciets to the saples
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More informationMultiscale Entropy Analysis: A New Method to Detect Determinism in a Time. Series. A. Sarkar and P. Barat. Variable Energy Cyclotron Centre
Multiscale Entropy Analysis: A New Method to Detect Deterinis in a Tie Series A. Sarkar and P. Barat Variable Energy Cyclotron Centre /AF Bidhan Nagar, Kolkata 700064, India PACS nubers: 05.45.Tp, 89.75.-k,
More informationSymmetric Discrete Orthonormal Stockwell Transform
Symmetric Discrete Orthonormal Stockwell Transform Yanwei Wang and Jeff Orchard Department of Applied Mathematics; David R. Cheriton School of Computer Science, University Avenue West, University of Waterloo,
More informationIntroduction to Multiresolution Analysis of Wavelets
Introuction to Multiresolution Analysis o Wavelets Aso Ray Proessor o Mechanical Engineering The Pennsylvania State University University Par PA 68 Tel: (84) 865-6377 Eail: axr@psu.eu Organization o the
More informationSHORT TIME FOURIER TRANSFORM PROBABILITY DISTRIBUTION FOR TIME-FREQUENCY SEGMENTATION
SHORT TIME FOURIER TRANSFORM PROBABILITY DISTRIBUTION FOR TIME-FREQUENCY SEGMENTATION Fabien Millioz, Julien Huillery, Nadine Martin To cite this version: Fabien Millioz, Julien Huillery, Nadine Martin.
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More information( )( )( )( )( ) University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2010
Hoework Assignent 4: Due at 5 p.. 7/9/ Text Probles: 6.4, 6.9, 6, 6.3 University o Washington Departent o Cheistry Cheistry 45/456 Suer Quarter P6.9) Jaes Watt once observed that a hard-working horse can
More informationDepartment of Electronic and Optical Engineering, Ordnance Engineering College, Shijiazhuang, , China
6th International Conference on Machinery, Materials, Environent, Biotechnology and Coputer (MMEBC 06) Solving Multi-Sensor Multi-Target Assignent Proble Based on Copositive Cobat Efficiency and QPSO Algorith
More informationDETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION
DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION Masaki WAKUI 1 and Jun IYAMA and Tsuyoshi KOYAMA 3 ABSTRACT This paper shows a criteria to detect
More informationA Graphical Model Approach to Source Localization in Wireless Sensor Networks
A Graphical Model Approach to Source Localization in Wireless Sensor Networks Manish Kushwaha, Xenoon Koutsoukos Institute or Sotware Integrated Systes (ISIS) Departent o Electrical Engineering and Coputer
More informationlecture 36: Linear Multistep Mehods: Zero Stability
95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,
More informationDonald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More informationSmith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtime
Sith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtie Oscar Caacho, a, * Francisco De la Cruz b a Postgrado en Autoatización e Instruentación. Grupo en Nuevas Estrategias
More informationMobile Communications TCS 455
Mobile Counication TCS 455 Dr. Prapun Sukopong prapun@iit.tu.ac.th Lecture 24 1 Office Hour: BKD 3601-7 Tueday 14:00-16:00 Thurday 9:30-11:30 Announceent Read Chapter 9: 9.1 9.5 Section 1.2 fro [Bahai,
More informationHomework 3 Solutions CSE 101 Summer 2017
Hoework 3 Solutions CSE 0 Suer 207. Scheduling algoriths The following n = 2 jobs with given processing ties have to be scheduled on = 3 parallel and identical processors with the objective of iniizing
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationAn analytical relation between relaxation time spectrum and molecular weight distribution
An analytical relation between relaxation tie spectru and olecular weight distribution Wolfgang Thi, Christian Friedrich, a) Michael Marth, and Josef Honerkap b) Freiburger Materialforschungszentru, Stefan-Meier-Straße
More informationPOD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION EQUATION IN NEURO-MUSCULAR SYSTEM
6th European Conference on Coputational Mechanics (ECCM 6) 7th European Conference on Coputational Fluid Dynaics (ECFD 7) 1115 June 2018, Glasgow, UK POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION
More informationStructural optimization of an automobile transmission case to minimize radiation noise using the model reduction technique
Journal o Mechanical Science and Technology 25 (5) (2011) 1247~1255 www.springerlink.co/content/1738-494x DOI 10.1007/s12206-011-0135-3 Structural optiization o an autoobile transission case to iniize
More informationA Division Algorithm Using Bisection Method in Residue Number System
International Journal of Coputer, Consuer and Control IJ3C), Vol., No. 03) 59 A Division Algorith Using Bisection Method in Residue Nuber Syste * Chin-Chen Chang and Jen-Ho Yang Abstract. Introduction
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 3 Signals & Systems Pro. Mark Fowler Discussion #9 Illustrating the Errors in DFT Processing DFT or Sonar Processing Example # Illustrating The Errors in DFT Processing Illustrating the Errors in
More information&& && F( u)! "{ f (x)} = f (x)e # j 2$ u x. f (x)! " #1. F(u,v) = f (x, y) e. f (x, y) = 2D Fourier Transform. Fourier Transform - review.
2D Fourier Transfor 2-D DFT & Properties 2D Fourier Transfor 1 Fourier Transfor - review 1-D: 2-D: F( u)! "{ f (x)} = f (x)e # j 2$ u x % & #% dx f (x)! " #1 { F(u) } = F(u)e j 2$ u x du F(u,v) = f (x,
More informationCooperative spectrum sensing with secondary user selection for cognitive radio networks over Nakagami-m fading channels
IET Counications Research Article Cooperative spectru sensing with secondary user selection or cognitive radio networks over Nakagai- ading channels ISSN 1751-8628 Received on 14th May 2015 Revised on
More informationThe Weierstrass Approximation Theorem
36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined
More informationAssessment of wind-induced structural fatigue based on joint probability density function of wind speed and direction
The 1 World Congress on Advances in Civil, Environental, and Materials Research (ACEM 1) eoul, Korea, August 6-3, 1 Assessent of wind-induced structural fatigue based on oint probability density function
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More informationIntroduction to Analog And Digital Communications
Introduction to Analog And Digital Communications Second Edition Simon Haykin, Michael Moher Chapter Fourier Representation o Signals and Systems.1 The Fourier Transorm. Properties o the Fourier Transorm.3
More informationTracking using CONDENSATION: Conditional Density Propagation
Tracking using CONDENSATION: Conditional Density Propagation Goal Model-based visual tracking in dense clutter at near video frae rates M. Isard and A. Blake, CONDENSATION Conditional density propagation
More informationDEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS
ISSN 1440-771X AUSTRALIA DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS An Iproved Method for Bandwidth Selection When Estiating ROC Curves Peter G Hall and Rob J Hyndan Working Paper 11/00 An iproved
More informationTransfer from Arbitrary Pure State to Target Mixed State for Quantum Systems
reprints o the 8th FAC World Congress ilano (taly August 8 - epteber, ranser ro Arbitrary ure tate to arget ixed tate or Quantu ystes J. WEN,. CONG Departent o Autoation, University o cience and echnology
More informationAn improved self-adaptive harmony search algorithm for joint replenishment problems
An iproved self-adaptive harony search algorith for joint replenishent probles Lin Wang School of Manageent, Huazhong University of Science & Technology zhoulearner@gail.co Xiaojian Zhou School of Manageent,
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationThe case for cycling the representer method with nonlinear models
The case for cycling the representer ethod with nonlinear odels Hans Ngodock Scott Sith Gregg Jacobs Naval Research Lab Stennis Space Center 1 Outline The representer ethod (long story short) The EnKF/S
More informationRandomized Recovery for Boolean Compressed Sensing
Randoized Recovery for Boolean Copressed Sensing Mitra Fatei and Martin Vetterli Laboratory of Audiovisual Counication École Polytechnique Fédéral de Lausanne (EPFL) Eail: {itra.fatei, artin.vetterli}@epfl.ch
More informationFinite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields
Finite fields I talked in class about the field with two eleents F 2 = {, } and we ve used it in various eaples and hoework probles. In these notes I will introduce ore finite fields F p = {,,...,p } for
More informationRational Filter Wavelets*
Ž Journal of Matheatical Analysis and Applications 39, 744 1999 Article ID jaa19996550, available online at http:wwwidealibraryco on Rational Filter Wavelets* Kuang Zheng and Cui Minggen Departent of Matheatics,
More informationKONINKL. NEDERL. AKADEMIE VAN WETENSCHAPPEN AMSTERDAM Reprinted from Proceedings, Series A, 61, No. 1 and Indag. Math., 20, No.
KONINKL. NEDERL. AKADEMIE VAN WETENSCHAPPEN AMSTERDAM Reprinted fro Proceedings, Series A, 6, No. and Indag. Math., 20, No., 95 8 MATHEMATIC S ON SEQUENCES OF INTEGERS GENERATED BY A SIEVIN G PROCES S
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationMiddlesex University Research Repository
Middlesex University Research Repository An open access repository o Middlesex University research http://eprints.dx.ac.uk Vien, Quoc-Tuan and Stewart, Brian G. and Tianield, Huaglory and Nguyen, Huan
More informationGeneral Properties of Radiation Detectors Supplements
Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent
More informationIAENG International Journal of Computer Science, 42:2, IJCS_42_2_06. Approximation Capabilities of Interpretable Fuzzy Inference Systems
IAENG International Journal of Coputer Science, 4:, IJCS_4 6 Approxiation Capabilities of Interpretable Fuzzy Inference Systes Hirofui Miyajia, Noritaka Shigei, and Hiroi Miyajia 3 Abstract Many studies
More informationModified Systematic Sampling in the Presence of Linear Trend
Modified Systeatic Sapling in the Presence of Linear Trend Zaheen Khan, and Javid Shabbir Keywords: Abstract A new systeatic sapling design called Modified Systeatic Sapling (MSS, proposed by ] is ore
More informationDamage Detection using Stochastic Subspace Identification
Daage Detection using Stochastic Subspace Identiication S. H. Si 1 and B. F. Spencer, Jr. 2 1 Departent o Civil and Environental Engineering, University o Illinois at Urbana-Chapaign, Urbana, IL 6181,
More informationWhen Short Runs Beat Long Runs
When Short Runs Beat Long Runs Sean Luke George Mason University http://www.cs.gu.edu/ sean/ Abstract What will yield the best results: doing one run n generations long or doing runs n/ generations long
More informationApplication of Differential Geometry in Magnetic Levitation Control
Appl. Math. In. Sci. 8, No. 5, 2437-2442 (2014) 2437 Applied Matheatics & Inoration Sciences An International Journal http:dx.doi.org10.12785ais080539 Application o Dierential Geoetry in Magnetic Levitation
More informationHIGH RESOLUTION NEAR-FIELD MULTIPLE TARGET DETECTION AND LOCALIZATION USING SUPPORT VECTOR MACHINES
ICONIC 2007 St. Louis, O, USA June 27-29, 2007 HIGH RESOLUTION NEAR-FIELD ULTIPLE TARGET DETECTION AND LOCALIZATION USING SUPPORT VECTOR ACHINES A. Randazzo,. A. Abou-Khousa 2,.Pastorino, and R. Zoughi
More informationNonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy
Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo
More informationA VARIABLE STEP-SIZE FREQUENCY-DOMAIN ADAPTIVE FILTERING ALGORITHM FOR STEREOPHONIC ACOUSTIC ECHO CANCELLATION
18th European Signal Processing Conference EUSIPCO-2010 Aalborg, Denark, August 23-27, 2010 A VARIABLE STEP-SIZE FREQUENCY-DOMAIN ADAPTIVE FILTERING ALGORITHM FOR STEREOPHONIC ACOUSTIC ECHO CANCELLATION
More informationA new type of lower bound for the largest eigenvalue of a symmetric matrix
Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box
More informationPHY 171. Lecture 14. (February 16, 2012)
PHY 171 Lecture 14 (February 16, 212) In the last lecture, we looked at a quantitative connection between acroscopic and icroscopic quantities by deriving an expression for pressure based on the assuptions
More informationOptimization of Flywheel Weight using Genetic Algorithm
IN: 78 7798 International Journal o cience, Engineering and Technology esearch (IJET) Volue, Issue, March 0 Optiization o Flywheel Weight using Genetic Algorith A. G. Appala Naidu*, B. T.N. Charyulu, C..C.V.aanaurthy
More information