AN IMPROVED VARIATIONAL MODE DECOMPOSITION METHOD FOR INTERNAL WAVES SEPARATION
|
|
- Marjorie Ramsey
- 6 years ago
- Views:
Transcription
1 3rd Eropean Signal Processing Conference (EUSIPCO AN IMPROVED VARIATIONAL MODE DECOMPOSITION METHOD FOR INTERNAL WAVES SEPARATION Jeremy Schmitt, Ernesto Horne, Nelly Pstelni, Sylvain Jobad, Philippe Odier Laboratoire de Physiqe de l Ecole Normale Spériere de Lyon, UMR CNRS 567, F-6936 Lyon, France, firstname.lastname@ens-lyon.fr ABSTRACT This paper proposes to revisit the -D Variational Mode Decomposition (-D VMD in order to separate the incident and reflected waves in experimental images of internal waves velocity field. -D VMD aims at splitting an image into a seqence of oscillating components which are centered arond specific spatial freqencies. In this wor we develop a proximal algorithm with local convergence garantees, allowing more flexibility in order to deal with modes having different spectral properties and to add some optional constraints modeling prior informations. Or method is compared with the standard -D VMD and with a Hilbert based strategy sally employed for processing internal waves images. Index Terms Optimisation, mode decomposition, internal waves.. INTRODUCTION The aim of this stdy is to extract reflecting and incident waves from experimental images modeling the velocity field of oceanic internal wave reflection. Internal waves are waves propagating in the bl of a flid whose density varies with depth (ocean, atmosphere, interior of stars []. De to their intrinsic properties, these waves, as they reflect pon an inclined slope, can ndergo a strong focalization []. An experiment designed to stdy this phenomenon consists in observing the waves throgh a measrement of the flid velocity [3]. The main challenge in this stdy is to separate the incident and reflected waves, which have the same temporal freqency. The incident wave is localized in the spatial Forier spectrm, while the reflected wave is not. The classical procedre to perform this separation is based on a Hilbert transform [] bt it fails to achieve this separation properly by introdcing some bondary artefacts and ndesirable oscillations. This qestion can be formlated as an inverse problem that consists in extracting K oscillating components, denoted ( K with R N N, from the observed data z This wor is spported by the ANR-3-BS3- TRES grant and ANR--BS-6- ONLITUR R N N sch that z = K + ε, = where ε models an additive noise. The stdy of internal wave reflection corresponds to the specific case K =, where z models the velocity field, and each component is centered arond an nnown freqency ω = (ω,, ω, so that, for every location (l, l,..., N,..., N,,(l,l = a,(l,l cos ( v,,(l,l l + v,,(l,l l + ϕ, where a R N N models the amplitde changes and, for i,, the mean vale of v,i R N N is close to ω,i, and ϕ is a phase term. In the recent literatre of component extraction for image analysis, efficient methods are -D synchrosqeezing [5] or -D empirical mode decomposition [6]. The drawbac of the first one is to reqire a spervised domain selection in order to extract each component while the second one may sffer lac of robstness when noise is involved. Another strategy has been proposed in [7] based on monogenic analysis and in [8] for estimating jointly ( K and (ω K. However, all these for methods are not designed to extract components that can be spread ot on the freqency spectrm in one direction. In this wor, we propose to revisit the wor by Dragomirestiy and Zosso, initially formlated for analyzing -D signals [9] and extended to image analysis in [8], in order to incorporate constraints imposed by internal waves analysis. In section, we focs on the design of a new criterion adapted to the problem of extracting reflecting and incident waves from velocity field data. In Section 3, we derive an algorithm solving the associated minimization problem insring some local convergence garantees and in Section, experimental reslts demonstrate the efficiency of the proposed method compared to state-of-the-art strategies. A conclsion is drawn in Section 5. Notation A matrix in R N N is denoted in bold as and its components are ( (l,l L with L =,..., N,..., N. The Forier transform of the image is denoted û and the freqency indices are (ν, ν V = N /,..., N / N /,..., N / /5/$3. 5 IEEE 98
2 3rd Eropean Signal Processing Conference (EUSIPCO. CRITERION The soltion proposed in [8] aims to estimate jointly ( K and (ω K by solving z K min (,ω K + α K D ( = =,(l,l e j(ω,l +ω, l (l,l L, (,where the D analytic signal is defined in the Forier domain, for every freqency (ν, ν V, as û,(ν,ν = ( + sign(ω,ν + ω, ν û,(ν,ν. The D analytic signal is chosen to set to zero one half-plane of the freqency domain relatively to the freqency vector ω. D models the gradient operator and the coefficient α > denotes a reglarization parameter allowing to adjst the selectivity of the filter. The limitations related to this criterion are essentially twofold. First, incident and reflected wave have different spectral behaviors. In particlar, the spectrm of the reflected wave is very compact horizontally bt not vertically. We introdce parameters α,i, depending of the mode and the coordinate i,, in order to separately adjst the horizontal and vertical spectral compacity of each mode. Second, we wold lie to add some prior information lined to the physics of internal waves. In particlar, the reflected wave only propagates along the reflection slope and it is nll far away from the slope. This prior information can be introdced by means of a penalty term f (. For instance it can be chosen as an indicator fnction i C ( whose vale is if C = R N ( (l, l S (l,l = and + otherwise. For sch a choice of the penalty f, we impose the component to be zero in the set of indices S. Sch a choice of C denotes a non-empty, closed and convex sbset of R N N. Bt more generally, we assme that f is at least convex, lower semi-continos and proper from R N N to ], + ]. According to previos remars, the criterion we derive is: min f ( + λ z,ω + D ( α,,(l,l e jω,l (l,l L K D ( + α,, ( =,(l,l e jω,l (l,l L where D and D denote respectively the horizontal and vertical gradient operator. The parameters α, and α, allowing to adjst the selectivity are chosen positive. The parameter λ allows s to adjst the fidelity to data. Similarly to (, this criterion is also non-convex. In what follows, we propose an algorithmic soltion to estimate a local minimizer. 3. ALGORITHM A poplar approach for solving ( consists in alternating the minimization over ( K and (ω K. For every iterations m and each mode,..., K, we denote φ,m ( = f ( + λ + i + α, D ( (l,l e jω[m] [m] i z + α, D ( (l,l e jω[m] ψ,m (ω = α, D (,[m+],(l,l + α, D (,[m+],(l,l l, (l,l L, l (l,l L e jωl (l,l L e jωl (l,l L (3, (. Based on this fnctions, an alternating procedre is smmarized in Algorithm. Algorithm Gass-Seidel algorithm For m =,,... For =,..., K [m+] Argmin φ,m ( R N N ω [m+] Argmin ψ,m (ω ω R This algorithmic strategy is sed in [8] when α = α, = α, and f =. For this configration, sch a minimization strategy provides satisfactory reslts however it involves parameters tricy to adjst. Moreover, it is well nown that sch an alternating minimization procedre reqires restrictive conditions to garantee convergence to a local minimizer as discssed in [, ]. A simple soltion to overcome these difficlties is to replace each minimisation step by a proximity operator step as sggested in [] and described in Algorithm. Algorithm has attractive convergence properties. The convergence to a global minimm is ensred in the case where the criterion is coercive and semi-algebraic []. Algorithm Proximal alternating algorithm Set µ > and ρ > For m =,,... For =,..., K [m+] ω [m+] = arg min φ,m ( + ρ [m] R N N = arg min ψ,m (ω + µ ω ω[m] ω R In order to solve each step efficiently, we rewrite them in the Forier domain. We denote F the Forier transform and F the inverse Forier transform. De to the Parseval/Plancherel theorem, both sbproblems can be written in 98
3 3rd Eropean Signal Processing Conference (EUSIPCO the Forier domain: û [m+] = arg min û û [m] û + f (F û ρ + λ û + û [m] i ẑ i ( + α, j(ν ω [m], ( + sign(ν û (ν,ν ( +α, j(ν ω [m], (+sign(ν û (ν,ν ω [m+] = arg min ω R µ ω[m] ω + α, ( j(ν ω ( + sign(ν û [m+],(ν,ν ( + α, j(ν ω ( + sign(ν û [m+],(ν,ν (ν,ν V (ν,ν V, (ν,ν V (ν,ν V The estimation step of ω [m+] has an analytic expression that is, for every i,, ω [m] ω [m+],i + 8µα,i,i = + 8µα,i (ν,ν Ω [m] where Ω [m] Ω [m] (û[m+] (ν,ν Ω ν [m] i,(ν,ν (û[m+],(ν,ν, (5 denotes the half-plane in the Forier domain: = ν = (ν, ν V ν ω [m], + ν ω [m],. On the other hand, the estimation of [m+] is obtained throgh forward-bacward iterations generating a seqence (v [p] p bild as v [p+] ( [p] = prox γρf v + γb γf AF v [p], (6 where the proximity operator step allows s to deal with the non-smooth fnction f. It is defined as ( x R N N prox f x = arg min y f (y + x y, which is redced to the projection onto a closed convex set when f = ι C. In (6, A models a linear operator defined as A A: û (ν,ν û (ν,ν if (ν, ν Ω [m], otherwise, with A (ν,ν = +ρλ+ρ(α, (ν ω, +α, (ν ω, and b R N N is sch that b (l,l = [m],(l + ρλ(z,l (l,l [m] i,(l,l i. Algorithm 3 Prox-D-VMD algorithm Choose the parameters λ >, ρ >, µ > Set the parameters α, >, α, > Initialize ( [] K, ( ω [] K For every m =,,... For every,..., K Set v [] = [m] Set γ ] ; / A [ For p =,,... Update v [p] Set [m+] as detailled in (6 = lim p + v [p] Compte ω [m+], and ω [m+], sing (5. EXPERIMENTAL RESULTS Experiments are performed on a critical reflection velocity field image z of size N = 3 and N = 38 illstrated on Fig.(a. Or goal is to extract the incident wave and the reflected wave (ths, K = from z. Or method is compared with a method based on Hilbert-transform, which is sally sed to process the internal waves velocity field [] and with the -D VMD proposed in [8] where the data fidelity constraint is enforced strictly throgh an agmented Lagrangian. An illstration of each strategy is illstrated in Fig.. One can observe that the Hilbert based strategy leads to mode mixing, while the standard -D VMD provides ndesirables oscillations in the top of the -D field. The proposed method removes all these artefacts. In order to evalate the impact of the proposed constraints, we perform three different experiments. In the first experiment, there is no optional constraint (f = f = and no directionality (α,i α. In the second experiment, we add the directionality constraint for the second component by setting α, = α/, since the spectrm of the reflected wave is expected to be more spread ot horizontally than vertically. In the third experiment, we combine the directionality with the zero constraint by setting S to be the top first 7 rows of (α, = α/ and f = i C. The reslts are illstrated in Figre, which shows mean vertical profiles of the incident wave (solid line and the reflected wave (dashed line. Mean profiles are obtained by averaging the absolte vale on each row of the images. In each plot, we compare the proposed soltion with the ones obtained with the Hilbert filtering and -D VMD. The left colmn presents the reslts of the proposed method where α is fixed and λ varies. On the right colmn, λ is fixed and α varies. In the first experiment, we can observe that for larger λ, the soltion is closer to the standard -D VMD that imposes the exact eqality between the data and the extracted components. Here, varying λ allows s some flexibility, for instance when noise is involved. A similar behavior is observed regarding the parameter α. One can observe that the basic Hilbert filtering strongly attenates the incident wave The iterations are smmarized in Algorithm 3. 98
4 3rd Eropean Signal Processing Conference (EUSIPCO Data z Incident wave Reflected wave Hilbert (a (b (c -D VMD (d (e Proposed method (f (g Fig.. Reslts on internal waves velocity field with critical wave reflection. st row: reslts obtained throgh a temporal and spatial Hilbert filtering. nd row: -D VMD. 3rd row: proposed method with α, = α/ and zero vale constraint on the 7 first rows of mode (f = i C. on the bondaries and leads to ndesirable oscillations on the reflected wave. In the second experiment, the directionality constraint allows s to redce the oscillations of the first IMF on the first 3 pixels near the slope. Moreover, adding this flexibility in the directionality allows s to be less sensitive to the choices of λ and α, which is why the profiles for different vales of parameters are sperposed. In the third experiment, we illstrate that the zero-constraint on the pixels of the second component which are far from the slope enables to cancel the ndesired oscillations of the reflected wave and the bondary effects. In or experiments, we observe that iterates are enogh to insre convergence. Withot the zero constraint, the convergence of the algorithm is fast (abot 3 seconds while with a zero constraint, the convergence is close to mintes (since the comptation of the modes involves forward-bacward iterations. In or experiments, the [] are set to zero, while the ω [] are initialized on the nit circle. however, we observe that the initialization of the and ω have no incidence on the final estimate. 5. CONCLUSION This paper presents an efficient method for the separation of incident and reflected wave in internal waves experimental images. This method revisits the -D VMD proposed by Dragomiretsiy and Zosso by adding stronger convergence garantees, giving more degrees of freedom for the setting of parameters and the possibility to add some prior nowledge throgh additional constraints. The proposed method has been tested on internal wave experimental images with critical reflection, and otperforms the classic Hilbert filtering since it strongly redces the bondary effects and the ndesired oscillations in the estimated reflected wave. REFERENCES [] B. R. Stherland, Internal gravity waves, Cambridge University Press,. [] OM Phillips, The dynamics of the pper ocean, Cambridge University Press, 966. [3] G. Bordes, A. Venaille, S. Jobad, P. Odier, and T. Daxois, Experimental observation of a strong mean flow indced by internal gravity waves., Physics of Flids, vol., no. 8, pp. Article No. 866, Agst. [] M.J. Mercier, N.B. Garnier, and T. Daxois, Reflection and diffraction of internal waves analyzed with the hilbert transform, Phys. Flids, vol., no. 8, pp. 866, 8. [5] M. Clasel, T. Oberlin, and V. Perrier, The monogenic synchrosqeezed wavelet transform: A tool for the decomposition/demodlation of AM-FM images, Appl. Comp. Harm. Anal.,, to appear. [6] J. Schmitt, N. Pstelni, P. Borgnat, P. Flandrin, and L. Condat, A -D Prony-Hang transform: A new tool for -D spectral analysis, IEEE Trans. Image Process., vol., pp , Oct.. [7] M. Unser, D. Sage, and D. Van De Ville, Mltiresoltion monogenic signal analysis sing the Riesz-Laplace wavelet transform, IEEE Trans. Image Process., vol. 8, no., pp. 8, November 9. [8] K. Dragomiretsiy and D. Zosso, Two-Dimensional Variational Mode Decomposition, in Energy Minimization Methods in Compter Vision and Pattern Recognition, Hong Kong, China, Janary [9] K. Dragomiretsiy and D. Zosso, Variational Mode Decomposition, IEEE Trans. Signal Process., vol. 6, no. 3, pp. 53 5, Febrary. [] T. F. Chan and C. K. Wong, Convergence of the alternating minimization algorithm for blind deconvoltion, Linear Algebra Appl, vol. 36, pp ,. [] J. Bolte, P. L. Combettes, and J.-C Pesqet, Alternating proximal algorithm for blind image recovery, in Proc. Int. Conf. Image Process., Hong Kong, China, September, pp , hal
5 3rd Eropean Signal Processing Conference (EUSIPCO α =, λ =,, 3 α =, 5,, λ =.8.6. Data Hilbert filtering D VMD Proposed method α =, λ = Proposed method α =, λ = Proposed method α =, λ = Data Hilbert filtering D VMD Proposed method α =, λ = Proposed method α = 5, λ = Proposed method α =, λ = Experiment Experiment Experiment Fig.. Mean profile over the vertical axis of the estimates for different vales of parameters α and λ. Blac solid line is the mean profile of data z. Color solid lines represent the extracted incident wave profiles with the different methods, while color dashed lines represent the extracted reflected waves profiles: brown for Hilbert filtering, red for -D VMD (α =, pin/ble/cyan for the proposed method. In the first colmn, λ = and α =, 5,. In the second colmn, α = and λ =,, 3. First row shows the profiles estimated by or method in the nidirectional and nconstrained case. Second row presents the profiles estimated with the proposed method in the directional case (α, = α/. Third row shows the profiles estimated with the proposed method in the directional and constrained case (f = i C. 98
2D HILBERT-HUANG TRANSFORM. Jérémy Schmitt, Nelly Pustelnik, Pierre Borgnat, Patrick Flandrin
2D HILBERT-HUANG TRANSFORM Jérémy Schmitt, Nelly Pustelnik, Pierre Borgnat, Patrick Flandrin Laboratoire de Physique de l Ecole Normale Suprieure de Lyon, CNRS and Université de Lyon, France first.last@ens-lyon.fr
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More informationLecture 8: September 26
10-704: Information Processing and Learning Fall 2016 Lectrer: Aarti Singh Lectre 8: September 26 Note: These notes are based on scribed notes from Spring15 offering of this corse. LaTeX template cortesy
More informationNonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator
Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational
More informationSolving a System of Equations
Solving a System of Eqations Objectives Understand how to solve a system of eqations with: - Gass Elimination Method - LU Decomposition Method - Gass-Seidel Method - Jacobi Method A system of linear algebraic
More informationAffine Invariant Total Variation Models
Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,
More informationStep-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter
WCE 007 Jly - 4 007 London UK Step-Size onds Analysis of the Generalized Mltidelay Adaptive Filter Jnghsi Lee and Hs Chang Hang Abstract In this paper we analyze the bonds of the fixed common step-size
More informationChapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More informationMove Blocking Strategies in Receding Horizon Control
Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationBridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control
Bridging the Gap Between Mltigrid, Hierarchical, and Receding-Horizon Control Victor M. Zavala Mathematics and Compter Science Division Argonne National Laboratory Argonne, IL 60439 USA (e-mail: vzavala@mcs.anl.gov).
More informationA State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics
A State Space Based Implicit Integration Algorithm for Differential Algebraic Eqations of Mltibody Dynamics E. J. Hag, D. Negrt, M. Ianc Janary 28, 1997 To Appear Mechanics of Strctres and Machines Abstract.
More informationFREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS
7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,
More information1 Undiscounted Problem (Deterministic)
Lectre 9: Linear Qadratic Control Problems 1 Undisconted Problem (Deterministic) Choose ( t ) 0 to Minimize (x trx t + tq t ) t=0 sbject to x t+1 = Ax t + B t, x 0 given. x t is an n-vector state, t a
More informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
More informationSafe Manual Control of the Furuta Pendulum
Safe Manal Control of the Frta Pendlm Johan Åkesson, Karl Johan Åström Department of Atomatic Control, Lnd Institte of Technology (LTH) Box 8, Lnd, Sweden PSfrag {jakesson,kja}@control.lth.se replacements
More informationCHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy
CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK Wassim Joini and Christophe Moy SUPELEC, IETR, SCEE, Avene de la Bolaie, CS 47601, 5576 Cesson Sévigné, France. INSERM U96 - IFR140-
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationRobust Tracking and Regulation Control of Uncertain Piecewise Linear Hybrid Systems
ISIS Tech. Rept. - 2003-005 Robst Tracking and Reglation Control of Uncertain Piecewise Linear Hybrid Systems Hai Lin Panos J. Antsaklis Department of Electrical Engineering, University of Notre Dame,
More informationBumpless transfer for discrete-time switched systems
29 American Control Conference Hyatt Regency Riverfront, St. Lois, O, USA Jne 1-12, 29 WeB12.3 Bmpless transfer for discrete-time switched systems I. ALLOCI, L. HETEL, J. DAAFOUZ, ember, IEEE, C. IUNG,
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationModel Discrimination of Polynomial Systems via Stochastic Inputs
Model Discrimination of Polynomial Systems via Stochastic Inpts D. Georgiev and E. Klavins Abstract Systems biologists are often faced with competing models for a given experimental system. Unfortnately,
More informationA Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model
entropy Article A Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model Hyenkyn Woo School of Liberal Arts, Korea University of Technology and Edcation, Cheonan
More informationAn extremum seeking approach to parameterised loop-shaping control design
Preprints of the 19th World Congress The International Federation of Atomatic Control An extremm seeking approach to parameterised loop-shaping control design Chih Feng Lee Sei Zhen Khong Erik Frisk Mattias
More informationRight Trapezoid Cover for Triangles of Perimeter Two
Kasetsart J (Nat Sci) 45 : 75-7 (0) Right Trapezoid Cover for Triangles of Perimeter Two Banyat Sroysang ABSTRACT A convex region covers a family of arcs if it contains a congrent copy of every arc in
More informationModelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationEfficient quadratic penalization through the partial minimization technique
This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC272754474, IEEE Transactions
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationStability of Model Predictive Control using Markov Chain Monte Carlo Optimisation
Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability
More informationReducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation
JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of
More informationA fundamental inverse problem in geosciences
A fndamental inverse problem in geosciences Predict the vales of a spatial random field (SRF) sing a set of observed vales of the same and/or other SRFs. y i L i ( ) + v, i,..., n i ( P)? L i () : linear
More informationAssignment Fall 2014
Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions
More informationSystem identification of buildings equipped with closed-loop control devices
System identification of bildings eqipped with closed-loop control devices Akira Mita a, Masako Kamibayashi b a Keio University, 3-14-1 Hiyoshi, Kohok-k, Yokohama 223-8522, Japan b East Japan Railway Company
More informationNon-collinear upconversion of infrared light
Non-collinear pconversion of infrared light Christian Pedersen, Qi H, Lasse Høgstedt, Peter Tidemand-Lichtenberg, and Jeppe Seidelin Dam * DTU Fotonik, Frederiksborgvej 399, 4000 Roskilde, Denmark * jdam@fotonik.dt.dk
More informationEvaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method
17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG
More informationA Model-Free Adaptive Control of Pulsed GTAW
A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,
More informationL 1 -smoothing for the Ornstein-Uhlenbeck semigroup
L -smoothing for the Ornstein-Uhlenbeck semigrop K. Ball, F. Barthe, W. Bednorz, K. Oleszkiewicz and P. Wolff September, 00 Abstract Given a probability density, we estimate the rate of decay of the measre
More informationMultivariable Ripple-Free Deadbeat Control
Scholarly Jornal of Mathematics and Compter Science Vol. (), pp. 9-9, October Available online at http:// www.scholarly-jornals.com/sjmcs ISS 76-8947 Scholarly-Jornals Fll Length Research aper Mltivariable
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationUpper Bounds on the Spanning Ratio of Constrained Theta-Graphs
Upper Bonds on the Spanning Ratio of Constrained Theta-Graphs Prosenjit Bose and André van Renssen School of Compter Science, Carleton University, Ottaa, Canada. jit@scs.carleton.ca, andre@cg.scs.carleton.ca
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationTheoretical and Experimental Implementation of DC Motor Nonlinear Controllers
Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,
More information), σ is a parameter, is the Euclidean norm in R d.
WEIGHTED GRAPH LAPLACIAN AND IMAGE INPAINTING ZUOQIANG SHI, STANLEY OSHER, AND WEI ZHU Abstract. Inspired by the graph Laplacian and the point integral method, we introdce a novel weighted graph Laplacian
More informationExperimental Study of an Impinging Round Jet
Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a
More informationAN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS
th Eropean Conference on Comptational Mechanics (ECCM ) 7th Eropean Conference on Comptational Flid Dynamics (ECFD 7) 5 Jne 28, Glasgow, UK AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD
More informationDownloaded 07/06/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 4, No., pp. A4 A7 c 8 Society for Indstrial and Applied Mathematics Downloaded 7/6/8 to 8.83.63.. Redistribtion sbject to SIAM license or copyright; see http://www.siam.org/jornals/ojsa.php
More informationStudy of the diffusion operator by the SPH method
IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait
More informationFRÉCHET KERNELS AND THE ADJOINT METHOD
PART II FRÉCHET KERNES AND THE ADJOINT METHOD 1. Setp of the tomographic problem: Why gradients? 2. The adjoint method 3. Practical 4. Special topics (sorce imaging and time reversal) Setp of the tomographic
More informationCHAPTER 8 ROTORS MOUNTED ON FLEXIBLE BEARINGS
CHAPTER 8 ROTORS MOUNTED ON FLEXIBLE BEARINGS Bearings commonly sed in heavy rotating machine play a significant role in the dynamic ehavior of rotors. Of particlar interest are the hydrodynamic earings,
More informationMulti-Voltage Floorplan Design with Optimal Voltage Assignment
Mlti-Voltage Floorplan Design with Optimal Voltage Assignment ABSTRACT Qian Zaichen Department of CSE The Chinese University of Hong Kong Shatin,N.T., Hong Kong zcqian@cse.chk.ed.hk In this paper, we stdy
More informationEstimating models of inverse systems
Estimating models of inverse systems Ylva Jng and Martin Enqvist Linköping University Post Print N.B.: When citing this work, cite the original article. Original Pblication: Ylva Jng and Martin Enqvist,
More informationNetwork Coding for Multiple Unicasts: An Approach based on Linear Optimization
Network Coding for Mltiple Unicasts: An Approach based on Linear Optimization Danail Traskov, Niranjan Ratnakar, Desmond S. Ln, Ralf Koetter, and Mriel Médard Abstract In this paper we consider the application
More informationGradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How
1 Gradient Projection Anti-windp Scheme on Constrained Planar LTI Systems Jstin Teo and Jonathan P. How Technical Report ACL1 1 Aerospace Controls Laboratory Department of Aeronatics and Astronatics Massachsetts
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationDetermining of temperature field in a L-shaped domain
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 0, (:-8 Determining of temperatre field in a L-shaped domain Oigo M. Zongo, Sié Kam, Kalifa Palm, and Alione Oedraogo
More informationA FIRST COURSE IN THE FINITE ELEMENT METHOD
INSTRUCTOR'S SOLUTIONS MANUAL TO ACCOMANY A IRST COURS IN TH INIT LMNT MTHOD ITH DITION DARYL L. LOGAN Contents Chapter 1 1 Chapter 3 Chapter 3 3 Chapter 17 Chapter 5 183 Chapter 6 81 Chapter 7 319 Chapter
More informationClassify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.
Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports
More informationApplying Fuzzy Set Approach into Achieving Quality Improvement for Qualitative Quality Response
Proceedings of the 007 WSES International Conference on Compter Engineering and pplications, Gold Coast, stralia, Janary 17-19, 007 5 pplying Fzzy Set pproach into chieving Qality Improvement for Qalitative
More informationThe SISTEM method. LOS ascending
The SISTEM method Simltaneos and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measrements A new global approach to obtain three-dimensional displacement maps by integrating
More informationCDS 110b: Lecture 1-2 Introduction to Optimal Control
CDS 110b: Lectre 1-2 Introdction to Optimal Control Richard M. Mrray 4 Janary 2006 Goals: Introdce the problem of optimal control as method of trajectory generation State the maimm principle and give eamples
More informationCOMPARATIVE STUDY OF ROBUST CONTROL TECHNIQUES FOR OTTO-CYCLE MOTOR CONTROL
ACM Symposim Series in Mechatronics - Vol. - pp.76-85 Copyright c 28 by ACM COMPARATIVE STUDY OF ROUST CONTROL TECHNIQUES FOR OTTO-CYCLE MOTOR CONTROL Marcos Salazar Francisco, marcos.salazar@member.isa.org
More informationPerformance analysis of the MAP equalizer within an iterative receiver including a channel estimator
Performance analysis of the MAP eqalizer within an iterative receiver inclding a channel estimator Nora Sellami ISECS Rote Menzel Chaer m 0.5 B.P 868, 308 Sfax, Tnisia Aline Romy IRISA-INRIA Camps de Bealie
More informationBLIND IMAGE WATERMARKING BASED ON SAMPLE ROTATION WITH OPTIMAL DETECTOR
7th Eropean Signal Processing Conference (EUSIPCO 009 Glasgow, Scotland, Agst -8, 009 BLIND IMAGE WATERMARKING BASED ON SAMPLE ROTATION WITH OPTIMAL DETECTOR S.M.E. Sahraeian, M.A. Akhaee, and F. Marvasti
More informationIncompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System
International Jornal of Compter Applications (97 8887) Volme 79 No., October Incompressible Viscoelastic Flow of a Generalised Oldroed-B Flid throgh Poros Medim between Two Infinite Parallel Plates in
More information1 The space of linear transformations from R n to R m :
Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations
More informationA Survey of the Implementation of Numerical Schemes for Linear Advection Equation
Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear
More informationB-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Quadratic Optimization Problems in Continuous and Binary Variables
B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Qadratic Optimization Problems in Continos and Binary Variables Naohiko Arima, Snyong Kim and Masakaz Kojima October 2012,
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationStair Matrix and its Applications to Massive MIMO Uplink Data Detection
1 Stair Matrix and its Applications to Massive MIMO Uplink Data Detection Fan Jiang, Stdent Member, IEEE, Cheng Li, Senior Member, IEEE, Zijn Gong, Senior Member, IEEE, and Roy S, Stdent Member, IEEE Abstract
More informationConvex Hypersurfaces of Constant Curvature in Hyperbolic Space
Srvey in Geometric Analysis and Relativity ALM 0, pp. 41 57 c Higher Edcation Press and International Press Beijing-Boston Convex Hypersrfaces of Constant Crvatre in Hyperbolic Space Bo Gan and Joel Sprck
More informationSensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters
Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed
More informationarxiv: v1 [physics.flu-dyn] 11 Mar 2011
arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department
More informationFast Algorithms for Restoration of Color Wireless Capsule Endoscopy Images
Fast Algorithms for Restoration of Color Wireless Capsle Endoscopy Images Haiying Li, W.-S. L, and Max Q.-H. Meng School of Control Science and Engineering, Shangdon University, Jinan, China Dept. of Electrical
More informationThe Real Stabilizability Radius of the Multi-Link Inverted Pendulum
Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract
More informationCFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion
CFD-Simlation thermoakstischer Resonanzeffekte zr Bestimmng der Flammentransferfnktion Ator: Dennis Paschke Technische Universität Berlin Institt für Strömngsmechanik nd Technische Akstik FG Experimentelle
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationAn optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model
An optimizing redced order FDS for the tropical Pacific Ocean redced gravity model Zhendong Lo a,, Jing Chen b,, Jiang Zh c,, Riwen Wang c,, and I. M. Navon d,, a School of Science, Beijing Jiaotong University,
More informationCHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES. George P. Yanev
Pliska Std. Math. Blgar. 2 (211), 233 242 STUDIA MATHEMATICA BULGARICA CHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES George P. Yanev We prove that the exponential
More informationDirect linearization method for nonlinear PDE s and the related kernel RBFs
Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract
More informationDiscussion Papers Department of Economics University of Copenhagen
Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,
More informationA Contraction of the Lucas Polygon
Western Washington University Western CEDAR Mathematics College of Science and Engineering 4 A Contraction of the Lcas Polygon Branko Ćrgs Western Washington University, brankocrgs@wwed Follow this and
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) MAE 5 - inite Element Analysis Several slides from this set are adapted from B.S. Altan, Michigan Technological University EA Procedre for
More informationKrauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.
Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationEVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE
13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth
More informationCubic graphs have bounded slope parameter
Cbic graphs have bonded slope parameter B. Keszegh, J. Pach, D. Pálvölgyi, and G. Tóth Agst 25, 2009 Abstract We show that every finite connected graph G with maximm degree three and with at least one
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationOn the Optimization of Numerical Dispersion and Dissipation of Finite Difference Scheme for Linear Advection Equation
Applied Mathematical Sciences, Vol. 0, 206, no. 48, 238-2389 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ams.206.6463 On the Optimization of Nmerical Dispersion and Dissipation of Finite Difference
More informationA TIME-FREQUENCY FORMULA FOR LMMSE FILTERS FOR NONSTATIONARY UNDERSPREAD CONTINUOUS-TIME STOCHASTIC PROCESSES
16th Eropean Signal Processing Conference (EUSIPCO 2008), Lasanne, Switzerland, Agst 25-29, 2008, copyright by EUASIP A TIME-FEQUENCY FOMULA FO LMMSE FILTES FO NONSTATIONAY UNDESPEAD CONTINUOUS-TIME STOCHASTIC
More informationJoint Transfer of Energy and Information in a Two-hop Relay Channel
Joint Transfer of Energy and Information in a Two-hop Relay Channel Ali H. Abdollahi Bafghi, Mahtab Mirmohseni, and Mohammad Reza Aref Information Systems and Secrity Lab (ISSL Department of Electrical
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More information