Shaping Low-Density Lattice Codes Using Voronoi Integers

Size: px
Start display at page:

Download "Shaping Low-Density Lattice Codes Using Voronoi Integers"

Transcription

1 Shaping Low-Density Lattice Codes Using Voronoi Integers Nuwan S. Ferdinand Brian M. Kurkoski Behnaam Aazhang Matti Latva-aho University of Oulu, Finland Japan Advanced Institute of Science and Technology Rice University, USA University of Oulu, Finland Information Theory Workshop Hobart, Tasmania, Australia November 2014

2 Nested Lattice Codes Achieve Capacity Lattice codes can achieve the capacity of AWGN channel [Erez and Zamir 04] Nested lattice codes: Want which is simultaneously good for coding and shaping Other information theoretic results using lattices: Lattices for relay channel e.g. [Song-Devroye 13] Two-way (Bidirectional) relay channel e.g. [Wilson et al.] Compute-forward relaying [Nazer-Gastpar 11] How to move from information theory to practical lattice codes? 2/24

3 Capacity-Approaching Lattice Constructions Recent high-dimension lattice constructions approach capacity Construction A with LDPC codes Construction D with turbo codes, spatially coupled LDPC Lattices based on polar codes Low-Density Lattice Codes [Sommer et al. 2008] Common claim: within few tenth of db of unconstrained capacity: No assumption about the channel power constraint. 3/24

4 Satisfy Power Constraint with Nested Lattices Good for correcting errors Good for quantization (satisfy the power constraint) high complexity 4/24

5 1.53 db Shaping Gain of Sphere over Cube Separate lattice and shaping region B contribution to signal power: Average Power B x 2 dx Depends only on shape of B (normalized second moment) nv (B) 2 n +1 G(B) M n V ( ) Depends only on coding lattice Λ Shaping Gain B B lim n G(cube) G(n-sphere) = e 6 = 1.53 db lim n G(n-sphere) = 1 2 e G(cube) = /24

6 for Well-Known Lattices Z 1 A 2D D4 D E 7 E Λ 16 Λ lattice dimension n

7 Satisfy Power Constraint with Nested Lattices n = 1, 2 n = 8 n = 23 n = 1000 n = 10 4 n = 10 5 E8 Leech small n Well-known lattices Weak coding gain efficient shaping algorithms Good shaping gain (0.65~1.0 db) Large n BP-based lattices Strong coding gain Inefficient shaping algorithms Uncertain coding gains: two cases: 0.4 db shaping gain 7/24

8 It Would Be Great If Find a construction that: Has the capacity-approaching coding gain high-dimension lattices Has the shaping gains and implementation complexity of a well-known lattice like E8. Must overcome the problem of mismatch in dimensions 8/24

9 Outline Key result: a lattice construction technique for shaping LDLC lattices Elements of the technique: 1. Voronoi Integers Shape integers using small-dimension lattices 2. Systematic lattice encoding: lattice point is nearby corresponding integer Results Full 0.65 db shaping gain of the E8 lattice. (2.1 db from 1/2 log(snr+1) ) Competing nested LDLCs obtained only 0.4 db, using higher complexity First, review 1.53 shaping gain result and LDLC lattices 9/24

10 Low-Density Lattice Codes LDLC lattices introduced by Sommer, Shalvi and Feder [IT 2008] LDLC have a sparse inverse generator matrix H Gaussian Belief-propagation decoding High dimension, n = 100, 1000, 10000, Come within 0.6 db of unconstrained capacity LDLCs for the power-constrained channel [Sommer et al ITW 2009] H matrix in triangular form, use M algorithm for quantization Obtained 0.4 db gain over hypercube (out of 1.53 db) 10/24

11 LDLC Latin Square Construction Inverse generator H = G 1 has constant row and column weight d. Latin square: each row/column {h 1, h 2,,h d } with random ±, h 1 h 2 h d Choose h 1 = 1 Example: {1, 1/2, 1/3} (forces determinant to be 1) Random sign changes d = 7 gives good performance BP convergence condition: 11/24

12 Nested Lattice Codes With LDLCs Sommer [ITW 2009]: Triangular construction Construct dimension n = 10,000 Put 1's on main diagonal, make triangular 90% Latin square weight d: {h 1, h 2,,h d } Quantization/shaping using M-Algorithm Such a lower triangular matrix can be generated with simi Complexity is O(ndM), but M is large Shaping gain of 0.4 db over hypercube Modest shaping gain for high complexity 12/24

13 Proposed Construction u 1 c 1 u P Z n Splitter u 2 Z m { s c 2 Combine integers c P Z n Systematic shaping x H 1 pc kq Substract offset a x 1 x a u n{m c n{m Voronoi Integers Systematic lattice encoding Offset to reduce average power 13/24

14 Voronoi Integers Z m / shape 5 Under systematic shaping, if the integers are 4 shaped, then lattice code will be shaped. 3 0,3 2 3,3 0,2 1,1 L is a small-dimensional lattice 1 2,3 3,2 0,1 1,0 2,7 quantization, i.e. shaping, is easy c 2 0 1,3 2,2 3,1 0,0 1,7 2,6 3,5 1 1,2 2,1 3,0 0,7 1,6 2,5 3,4 Define Voronoi Integers Z m / shape 2 2,0 3,7 0,6 1,5 2,4 set of integers inside fundamental region 3 3,6 0,5 1,4 4 0, c 1 14/24

15 13 12 Voronoi Integers Example /24

16 13 12 Voronoi Integers Example /24

17 13 12 Voronoi Integers Example /24

18 Systematic Lattice Encoding Requirement Triangular H with 1 s on diagonal 18/24

19 Systematic Lattice Encoding Example using 4 c = (2,4) (3,4) (4,4) Recall: c = round(x) Note Voronoi volume det(h) = 1 3 (3,3) (4,3) and the integer grid also has vol. 1 (2,3) No gaps (4,2) 2 (3,2) (2,2) /24

20 Systematic Lattice Encoding Example using 4 c = (2,4) (3,4) (4,4) Recall: c = round(x) Note Voronoi volume det(h) = 1 3 (3,3) (4,3) and the integer grid also has vol. 1 (2,3) No gaps (4,2) 2 (3,2) (2,2) /24

21 20/24

22 Average Transmit Power 0.7 E8 lattice alone No LDLC Transmit power, Gain over hypercube Gain over hypercube shaping a=a opt Z8/LE8 with LDLC systematic a=0 c (Z8/LE8) x (Z8/LE8 with LDLC) 0.2 E8 shaping gain bound log 2 (L) 21/24

23 Power-Constrained AWGN Channel AWGN channel with average power constraint 5 bits/dimension coding: LDLC lattice dimension n = 10,000 shaping: E8 lattice with m = 8 Compare with M-Algorithm LDLC shaping of Sommer et al 22/24

24 0.65 db Gain Over Hypercube Shaping Hypercube shaping [7] Nested lattice shaping 7] Proposed shaping Uniform input capacity AWGN capacity 0.15 db better than M algorithm, SER 10 3 and much lower complexity 0.4 db Average SNR in db 0.65dB 23/24

25 Conclusion Lattices are an alternative to finite-field codes for AWGN Shaping techniques to obtain 1.53 db are accessible Coset codes/nested lattice codes, high complexity We proposed: Voronoi integers using low-dimension lattices Systematic lattice shaping for LDLCs High coding gain of LDLCs, good shaping gain of E8 lattice 24/24

Lattice Coding and its Applications in Communications

Lattice Coding and its Applications in Communications Lattice Coding and its Applications in Communications Alister Burr University of York alister.burr@york.ac.uk Introduction to lattices Definition; Sphere packings; Basis vectors; Matrix description Codes

More information

Rewriting Codes for Flash Memories Based Upon Lattices, and an Example Using the E8 Lattice

Rewriting Codes for Flash Memories Based Upon Lattices, and an Example Using the E8 Lattice Rewriting Codes for Flash Memories Based Upon Lattices, and an Example Using the E Lattice Brian M. Kurkoski kurkoski@ice.uec.ac.jp University of Electro-Communications Tokyo, Japan Workshop on Application

More information

Leech Constellations. of Construction-A Lattices

Leech Constellations. of Construction-A Lattices Leech Constellations of Construction-A Lattices Nicola di Pietro and Joseph J. Boutros, Senior Member, IEEE arxiv:1611.04417v3 [cs.it] 2 Aug 2017 Abstract The problem of communicating over the additive

More information

Iterative Encoding with Gauss-Seidel Method for Spatially-Coupled Low-Density Lattice Codes

Iterative Encoding with Gauss-Seidel Method for Spatially-Coupled Low-Density Lattice Codes 202 IEEE International Symposium on Information Theory Proceeings Iterative Encoing with Gauss-Seiel Metho for Spatially-Couple Low-Density Lattice Coes Hironori Uchikawa, Brian M. Kurkoski, Kenta Kasai

More information

Integer Low-Density Lattices based on Construction A

Integer Low-Density Lattices based on Construction A Integer Low-Density Lattices based on Construction A Nicola di Pietro, Joseph J. Boutros, Gilles Zémor, Loïc Brunel Mitsubishi Electric R&D Centre Europe, Rennes, France Email: {n.dipietro, l.brunel}@fr.merce.mee.com

More information

Chapter 2 Uncertainty Analysis and Sampling Techniques

Chapter 2 Uncertainty Analysis and Sampling Techniques Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying

More information

Cross-Packing Lattices for the Rician Fading Channel

Cross-Packing Lattices for the Rician Fading Channel Cross-Packing Lattices for the Rician Fading Channel Amin Sakzad, Anna-Lena Trautmann, and Emanuele Viterbo Department of Electrical and Computer Systems Engineering, Monash University. Abstract We introduce

More information

Lecture outline. Monte Carlo Methods for Uncertainty Quantification. Importance Sampling. Importance Sampling

Lecture outline. Monte Carlo Methods for Uncertainty Quantification. Importance Sampling. Importance Sampling Lecture outline Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford KU Leuven Summer School on Uncertainty Quantification Lecture 2: Variance reduction

More information

Finite-length analysis of the TEP decoder for LDPC ensembles over the BEC

Finite-length analysis of the TEP decoder for LDPC ensembles over the BEC Finite-length analysis of the TEP decoder for LDPC ensembles over the BEC Pablo M. Olmos, Fernando Pérez-Cruz Departamento de Teoría de la Señal y Comunicaciones. Universidad Carlos III in Madrid. email:

More information

Monte Carlo Methods for Uncertainty Quantification

Monte Carlo Methods for Uncertainty Quantification Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)

More information

1102 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 51, NO. 3, MARCH Genyuan Wang and Xiang-Gen Xia, Senior Member, IEEE

1102 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 51, NO. 3, MARCH Genyuan Wang and Xiang-Gen Xia, Senior Member, IEEE 1102 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 51, NO 3, MARCH 2005 On Optimal Multilayer Cyclotomic Space Time Code Designs Genyuan Wang Xiang-Gen Xia, Senior Member, IEEE Abstract High rate large

More information

Contents Critique 26. portfolio optimization 32

Contents Critique 26. portfolio optimization 32 Contents Preface vii 1 Financial problems and numerical methods 3 1.1 MATLAB environment 4 1.1.1 Why MATLAB? 5 1.2 Fixed-income securities: analysis and portfolio immunization 6 1.2.1 Basic valuation of

More information

An advanced method for preserving skewness in single-variate, multivariate, and disaggregation models in stochastic hydrology

An advanced method for preserving skewness in single-variate, multivariate, and disaggregation models in stochastic hydrology XXIV General Assembly of European Geophysical Society The Hague, 9-3 April 999 HSA9.0 Open session on statistical methods in hydrology An advanced method for preserving skewness in single-variate, multivariate,

More information

Monte Carlo Methods in Finance

Monte Carlo Methods in Finance Monte Carlo Methods in Finance Peter Jackel JOHN WILEY & SONS, LTD Preface Acknowledgements Mathematical Notation xi xiii xv 1 Introduction 1 2 The Mathematics Behind Monte Carlo Methods 5 2.1 A Few Basic

More information

Uncertainty Analysis with UNICORN

Uncertainty Analysis with UNICORN Uncertainty Analysis with UNICORN D.A.Ababei D.Kurowicka R.M.Cooke D.A.Ababei@ewi.tudelft.nl D.Kurowicka@ewi.tudelft.nl R.M.Cooke@ewi.tudelft.nl Delft Institute for Applied Mathematics Delft University

More information

Lattices from equiangular tight frames with applications to lattice sparse recovery

Lattices from equiangular tight frames with applications to lattice sparse recovery Lattices from equiangular tight frames with applications to lattice sparse recovery Deanna Needell Dept of Mathematics, UCLA May 2017 Supported by NSF CAREER #1348721 and Alfred P. Sloan Fdn The compressed

More information

Square Grid Benchmarks for Source-Terminal Network Reliability Estimation

Square Grid Benchmarks for Source-Terminal Network Reliability Estimation Square Grid Benchmarks for Source-Terminal Network Reliability Estimation Roger Paredes Leonardo Duenas-Osorio Rice University, Houston TX, USA. 03/2018 This document describes a synthetic benchmark data

More information

Monte Carlo Methods for Uncertainty Quantification

Monte Carlo Methods for Uncertainty Quantification Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)

More information

High Dimensional Edgeworth Expansion. Applications to Bootstrap and Its Variants

High Dimensional Edgeworth Expansion. Applications to Bootstrap and Its Variants With Applications to Bootstrap and Its Variants Department of Statistics, UC Berkeley Stanford-Berkeley Colloquium, 2016 Francis Ysidro Edgeworth (1845-1926) Peter Gavin Hall (1951-2016) Table of Contents

More information

Pathloss and Link Budget From Physical Propagation to Multi-Path Fading Statistical Characterization of Channels. P r = P t Gr G t L P

Pathloss and Link Budget From Physical Propagation to Multi-Path Fading Statistical Characterization of Channels. P r = P t Gr G t L P Path Loss I Path loss L P relates the received signal power P r to the transmitted signal power P t : P r = P t Gr G t L P, where G t and G r are antenna gains. I Path loss is very important for cell and

More information

Monte Carlo Methods in Financial Engineering

Monte Carlo Methods in Financial Engineering Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures

More information

Square-Root Measurement for Ternary Coherent State Signal

Square-Root Measurement for Ternary Coherent State Signal ISSN 86-657 Square-Root Measurement for Ternary Coherent State Signal Kentaro Kato Quantum ICT Research Institute, Tamagawa University 6-- Tamagawa-gakuen, Machida, Tokyo 9-86, Japan Tamagawa University

More information

ELEMENTS OF MONTE CARLO SIMULATION

ELEMENTS OF MONTE CARLO SIMULATION APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the

More information

Show that the column rank and the row rank of A are both equal to 3.

Show that the column rank and the row rank of A are both equal to 3. hapter Vectors and matrices.. Exercises. Let A 2 5 4 3 2 4 2 2 3 5 4 2 4 3 Show that the column rank and the row rank of A are both equal to 3. 2. Let x and y be column vectors of size n, andleti be the

More information

Table of Contents. Kocaeli University Computer Engineering Department 2011 Spring Mustafa KIYAR Optimization Theory

Table of Contents. Kocaeli University Computer Engineering Department 2011 Spring Mustafa KIYAR Optimization Theory 1 Table of Contents Estimating Path Loss Exponent and Application with Log Normal Shadowing...2 Abstract...3 1Path Loss Models...4 1.1Free Space Path Loss Model...4 1.1.1Free Space Path Loss Equation:...4

More information

ROM Simulation with Exact Means, Covariances, and Multivariate Skewness

ROM Simulation with Exact Means, Covariances, and Multivariate Skewness ROM Simulation with Exact Means, Covariances, and Multivariate Skewness Michael Hanke 1 Spiridon Penev 2 Wolfgang Schief 2 Alex Weissensteiner 3 1 Institute for Finance, University of Liechtenstein 2 School

More information

CPSC 540: Machine Learning

CPSC 540: Machine Learning CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x

More information

Modeling Volatility Risk in Equity Options: a Cross-sectional approach

Modeling Volatility Risk in Equity Options: a Cross-sectional approach ICBI Global Derivatives, Amsterdam, 2014 Modeling Volatility Risk in Equity Options: a Cross-sectional approach Marco Avellaneda NYU & Finance Concepts Doris Dobi* NYU * This work is part of Doris Dobi

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days

Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days 1. Introduction Richard D. Christie Department of Electrical Engineering Box 35500 University of Washington Seattle, WA 98195-500 christie@ee.washington.edu

More information

GENERATION OF STANDARD NORMAL RANDOM NUMBERS. Naveen Kumar Boiroju and M. Krishna Reddy

GENERATION OF STANDARD NORMAL RANDOM NUMBERS. Naveen Kumar Boiroju and M. Krishna Reddy GENERATION OF STANDARD NORMAL RANDOM NUMBERS Naveen Kumar Boiroju and M. Krishna Reddy Department of Statistics, Osmania University, Hyderabad- 500 007, INDIA Email: nanibyrozu@gmail.com, reddymk54@gmail.com

More information

Financial Risk Management

Financial Risk Management Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #3 1 Maximum likelihood of the exponential distribution 1. We assume

More information

EC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods

EC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions

More information

Thursday, March 3

Thursday, March 3 5.53 Thursday, March 3 -person -sum (or constant sum) game theory -dimensional multi-dimensional Comments on first midterm: practice test will be on line coverage: every lecture prior to game theory quiz

More information

Volatility Spillovers and Causality of Carbon Emissions, Oil and Coal Spot and Futures for the EU and USA

Volatility Spillovers and Causality of Carbon Emissions, Oil and Coal Spot and Futures for the EU and USA 22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Volatility Spillovers and Causality of Carbon Emissions, Oil and Coal

More information

NOISE VARIANCE ESTIMATION IN DS-CDMA AND ITS EFFECTS ON THE INDIVIDUALLY OPTIMUM RECEIVER

NOISE VARIANCE ESTIMATION IN DS-CDMA AND ITS EFFECTS ON THE INDIVIDUALLY OPTIMUM RECEIVER NOISE VRINCE ESTIMTION IN DS-CDM ND ITS EFFECTS ON THE INDIVIDULLY OPTIMUM RECEIVER R. Gaudel, F. Bonnet, J.B. Domelevo-Entfellner ENS Cachan Campus de Ker Lann 357 Bruz, France. Roumy IRIS-INRI Campus

More information

Sphere Packings, Lattices and Groups

Sphere Packings, Lattices and Groups J.H. Conway N.J.A. Sloane Sphere Packings, Lattices and Groups Third Edition With Additional Contributions by E. Bannai, R.E. Borcherds, J. Leech, S.P. Norton, A.M. Odlyzko, R.A. Parker, L. Queen and B.B.

More information

Alternative VaR Models

Alternative VaR Models Alternative VaR Models Neil Roeth, Senior Risk Developer, TFG Financial Systems. 15 th July 2015 Abstract We describe a variety of VaR models in terms of their key attributes and differences, e.g., parametric

More information

BCJR Algorithm. Veterbi Algorithm (revisted) Consider covolutional encoder with. And information sequences of length h = 5

BCJR Algorithm. Veterbi Algorithm (revisted) Consider covolutional encoder with. And information sequences of length h = 5 Chapter 2 BCJR Algorithm Ammar Abh-Hhdrohss Islamic University -Gaza ١ Veterbi Algorithm (revisted) Consider covolutional encoder with And information sequences of length h = 5 The trellis diagram has

More information

Write legibly. Unreadable answers are worthless.

Write legibly. Unreadable answers are worthless. MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are

More information

Gamma Distribution Fitting

Gamma Distribution Fitting Chapter 552 Gamma Distribution Fitting Introduction This module fits the gamma probability distributions to a complete or censored set of individual or grouped data values. It outputs various statistics

More information

Quasi-Monte Carlo for Finance

Quasi-Monte Carlo for Finance Quasi-Monte Carlo for Finance Peter Kritzer Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Linz, Austria NCTS, Taipei, November 2016 Peter Kritzer

More information

Zhen Sun, Milind Dawande, Ganesh Janakiraman, and Vijay Mookerjee

Zhen Sun, Milind Dawande, Ganesh Janakiraman, and Vijay Mookerjee RESEARCH ARTICLE THE MAKING OF A GOOD IMPRESSION: INFORMATION HIDING IN AD ECHANGES Zhen Sun, Milind Dawande, Ganesh Janakiraman, and Vijay Mookerjee Naveen Jindal School of Management, The University

More information

Fast Convergence of Regress-later Series Estimators

Fast Convergence of Regress-later Series Estimators Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser

More information

ANALYSIS OF THE BINOMIAL METHOD

ANALYSIS OF THE BINOMIAL METHOD ANALYSIS OF THE BINOMIAL METHOD School of Mathematics 2013 OUTLINE 1 CONVERGENCE AND ERRORS OUTLINE 1 CONVERGENCE AND ERRORS 2 EXOTIC OPTIONS American Options Computational Effort OUTLINE 1 CONVERGENCE

More information

Computational Finance Binomial Trees Analysis

Computational Finance Binomial Trees Analysis Computational Finance Binomial Trees Analysis School of Mathematics 2018 Review - Binomial Trees Developed a multistep binomial lattice which will approximate the value of a European option Extended the

More information

EENG473 Mobile Communications Module 3 : Week # (11) Mobile Radio Propagation: Large-Scale Path Loss

EENG473 Mobile Communications Module 3 : Week # (11) Mobile Radio Propagation: Large-Scale Path Loss EENG473 Mobile Communications Module 3 : Week # (11) Mobile Radio Propagation: Large-Scale Path Loss Practical Link Budget Design using Path Loss Models Most radio propagation models are derived using

More information

Non-Data-Aided Parameter Estimation in an Additive White Gaussian Noise Channel

Non-Data-Aided Parameter Estimation in an Additive White Gaussian Noise Channel on-data-aided Parameter Estimation in an Additive White Gaussian oise Channel Fredrik Brännström Department of Signals and Systems Chalmers University of Technology SE-4 96 Göteborg, Sweden Email: fredrikb@s.chalmers.se

More information

FIT5124 Advanced Topics in Security. Lecture 1: Lattice-Based Crypto. I

FIT5124 Advanced Topics in Security. Lecture 1: Lattice-Based Crypto. I FIT5124 Advanced Topics in Security Lecture 1: Lattice-Based Crypto. I Ron Steinfeld Clayton School of IT Monash University March 2016 Acknowledgements: Some figures sourced from Oded Regev s Lecture Notes

More information

Information Processing and Limited Liability

Information Processing and Limited Liability Information Processing and Limited Liability Bartosz Maćkowiak European Central Bank and CEPR Mirko Wiederholt Northwestern University January 2012 Abstract Decision-makers often face limited liability

More information

Financial Computing with Python

Financial Computing with Python Introduction to Financial Computing with Python Matthieu Mariapragassam Why coding seems so easy? But is actually not Sprezzatura : «It s an art that doesn t seem to be an art» - The Book of the Courtier

More information

Solving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function?

Solving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function? DOI 0.007/s064-006-9073-z ORIGINAL PAPER Solving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function? Jules H. van Binsbergen Michael W. Brandt Received:

More information

Implementing Models in Quantitative Finance: Methods and Cases

Implementing Models in Quantitative Finance: Methods and Cases Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1

More information

CPSC 540: Machine Learning

CPSC 540: Machine Learning CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}

More information

Lecture 8 : The dual lattice and reducing SVP to MVP

Lecture 8 : The dual lattice and reducing SVP to MVP CSE 206A: Lattice Algorithms and Applications Spring 2007 Lecture 8 : The dual lattice and reducing SVP to MVP Lecturer: Daniele Micciancio Scribe: Scott Yilek 1 Overview In the last lecture we explored

More information

Path Loss Models and Link Budget

Path Loss Models and Link Budget Path Loss Models and Link Budget A universal path loss model P r dbm = P t dbm + db Gains db Losses Gains: the antenna gains compared to isotropic antennas Transmitter antenna gain Receiver antenna gain

More information

I. More Fundamental Concepts and Definitions from Mathematics

I. More Fundamental Concepts and Definitions from Mathematics An Introduction to Optimization The core of modern economics is the notion that individuals optimize. That is to say, individuals use the resources available to them to advance their own personal objectives

More information

Solving the Stochastic Steady-State Diffusion Problem Using Multigrid

Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Program Advisor: Howard Elman Department of Computer Science May 5, 2016 Tengfei

More information

THE ENERGY EFFICIENCY OF THE ERGODIC FADING RELAY CHANNEL

THE ENERGY EFFICIENCY OF THE ERGODIC FADING RELAY CHANNEL 7th European Signal Processing Conference (EUSIPCO 009) Glasgow, Scotland, August 4-8, 009 THE ENERGY EFFICIENCY OF THE ERGODIC FADING RELAY CHANNEL Jesús Gómez-Vilardebó Centre Tecnològic de Telecomunicacions

More information

Approximate methods for dynamic portfolio allocation under transaction costs

Approximate methods for dynamic portfolio allocation under transaction costs Western University Scholarship@Western Electronic Thesis and Dissertation Repository November 2012 Approximate methods for dynamic portfolio allocation under transaction costs Nabeel Butt The University

More information

GENERATION OF APPROXIMATE GAMMA SAMPLES BY PARTIAL REJECTION

GENERATION OF APPROXIMATE GAMMA SAMPLES BY PARTIAL REJECTION IASC8: December 5-8, 8, Yokohama, Japan GEERATIO OF APPROXIMATE GAMMA SAMPLES BY PARTIAL REJECTIO S.H. Ong 1 Wen Jau Lee 1 Institute of Mathematical Sciences, University of Malaya, 563 Kuala Lumpur, MALAYSIA

More information

Lindner, Szimayer: A Limit Theorem for Copulas

Lindner, Szimayer: A Limit Theorem for Copulas Lindner, Szimayer: A Limit Theorem for Copulas Sonderforschungsbereich 386, Paper 433 (2005) Online unter: http://epub.ub.uni-muenchen.de/ Projektpartner A Limit Theorem for Copulas Alexander Lindner Alexander

More information

Asymptotic results discrete time martingales and stochastic algorithms

Asymptotic results discrete time martingales and stochastic algorithms Asymptotic results discrete time martingales and stochastic algorithms Bernard Bercu Bordeaux University, France IFCAM Summer School Bangalore, India, July 2015 Bernard Bercu Asymptotic results for discrete

More information

Is Greedy Coordinate Descent a Terrible Algorithm?

Is Greedy Coordinate Descent a Terrible Algorithm? Is Greedy Coordinate Descent a Terrible Algorithm? Julie Nutini, Mark Schmidt, Issam Laradji, Michael Friedlander, Hoyt Koepke University of British Columbia Optimization and Big Data, 2015 Context: Random

More information

Constrained Sequential Resource Allocation and Guessing Games

Constrained Sequential Resource Allocation and Guessing Games 4946 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 54, NO. 11, NOVEMBER 2008 Constrained Sequential Resource Allocation and Guessing Games Nicholas B. Chang and Mingyan Liu, Member, IEEE Abstract In this

More information

Computational Methods in Finance

Computational Methods in Finance Chapman & Hall/CRC FINANCIAL MATHEMATICS SERIES Computational Methods in Finance AM Hirsa Ltfi) CRC Press VV^ J Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor &

More information

Model Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16

Model Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16 Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint

More information

A Systematic Global Macro Fund

A Systematic Global Macro Fund A Systematic Global Macro Fund Correlation and Portfolio Construction January 2013 Working Paper Lawson McWhorter, CMT, CFA Head of Research Abstract Trading strategies are usually evaluated primarily

More information

Package multiassetoptions

Package multiassetoptions Package multiassetoptions February 20, 2015 Type Package Title Finite Difference Method for Multi-Asset Option Valuation Version 0.1-1 Date 2015-01-31 Author Maintainer Michael Eichenberger

More information

Course information FN3142 Quantitative finance

Course information FN3142 Quantitative finance Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken

More information

ARCH Models and Financial Applications

ARCH Models and Financial Applications Christian Gourieroux ARCH Models and Financial Applications With 26 Figures Springer Contents 1 Introduction 1 1.1 The Development of ARCH Models 1 1.2 Book Content 4 2 Linear and Nonlinear Processes 5

More information

Online Appendix: Extensions

Online Appendix: Extensions B Online Appendix: Extensions In this online appendix we demonstrate that many important variations of the exact cost-basis LUL framework remain tractable. In particular, dual problem instances corresponding

More information

Prentice Hall Connected Mathematics, Grade 7 Unit 2004 Correlated to: Maine Learning Results for Mathematics (Grades 5-8)

Prentice Hall Connected Mathematics, Grade 7 Unit 2004 Correlated to: Maine Learning Results for Mathematics (Grades 5-8) : Maine Learning Results for Mathematics (Grades 5-8) A. NUMBERS AND NUMBER SENSE Students will understand and demonstrate a sense of what numbers mean and how they are used. Students will be able to:

More information

Outline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.

Outline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0. Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization

More information

Introduction to Sequential Monte Carlo Methods

Introduction to Sequential Monte Carlo Methods Introduction to Sequential Monte Carlo Methods Arnaud Doucet NCSU, October 2008 Arnaud Doucet () Introduction to SMC NCSU, October 2008 1 / 36 Preliminary Remarks Sequential Monte Carlo (SMC) are a set

More information

Financial Risk Management

Financial Risk Management Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given

More information

Omitted Variables Bias in Regime-Switching Models with Slope-Constrained Estimators: Evidence from Monte Carlo Simulations

Omitted Variables Bias in Regime-Switching Models with Slope-Constrained Estimators: Evidence from Monte Carlo Simulations Journal of Statistical and Econometric Methods, vol. 2, no.3, 2013, 49-55 ISSN: 2051-5057 (print version), 2051-5065(online) Scienpress Ltd, 2013 Omitted Variables Bias in Regime-Switching Models with

More information

International Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN

International Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18,   ISSN International Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL

More information

Multilevel quasi-monte Carlo path simulation

Multilevel quasi-monte Carlo path simulation Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein

More information

Portfolio Management and Optimal Execution via Convex Optimization

Portfolio Management and Optimal Execution via Convex Optimization Portfolio Management and Optimal Execution via Convex Optimization Enzo Busseti Stanford University April 9th, 2018 Problems portfolio management choose trades with optimization minimize risk, maximize

More information

Crystal Lattices. Daniele Toffoli December 7, / 42

Crystal Lattices. Daniele Toffoli December 7, / 42 Crystal Lattices Daniele Toffoli December 7, 2016 1 / 42 Outline 1 Bravais Lattice: definitions and examples 2 Unit cell: primitive, conventional and Wigner-Seitz 3 Crystal structure: lattices with basis

More information

A First Course in Probability

A First Course in Probability A First Course in Probability Seventh Edition Sheldon Ross University of Southern California PEARSON Prentice Hall Upper Saddle River, New Jersey 07458 Preface 1 Combinatorial Analysis 1 1.1 Introduction

More information

Monte-Carlo Planning: Introduction and Bandit Basics. Alan Fern

Monte-Carlo Planning: Introduction and Bandit Basics. Alan Fern Monte-Carlo Planning: Introduction and Bandit Basics Alan Fern 1 Large Worlds We have considered basic model-based planning algorithms Model-based planning: assumes MDP model is available Methods we learned

More information

Workshop Session 1: Theoretical understanding of flow-based method, intuitiveness, hybrid-coupling

Workshop Session 1: Theoretical understanding of flow-based method, intuitiveness, hybrid-coupling Workshop Session 1: Theoretical understanding of flow-based method, intuitiveness, hybrid-coupling by Joel HOEKSEMA (APX), Adrien ATAYI (RTE) and Pieter SCHAVEMAKER (e-bridge) 1 Workshop 1 Agenda I. Flow-based

More information

Approximating a multifactor di usion on a tree.

Approximating a multifactor di usion on a tree. Approximating a multifactor di usion on a tree. September 2004 Abstract A new method of approximating a multifactor Brownian di usion on a tree is presented. The method is based on local coupling of the

More information

A Comparison of Universal and Mean-Variance Efficient Portfolios p. 1/28

A Comparison of Universal and Mean-Variance Efficient Portfolios p. 1/28 A Comparison of Universal and Mean-Variance Efficient Portfolios Shane M. Haas Research Laboratory of Electronics, and Laboratory for Information and Decision Systems Massachusetts Institute of Technology

More information

Computational Finance Improving Monte Carlo

Computational Finance Improving Monte Carlo Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal

More information

16 MAKING SIMPLE DECISIONS

16 MAKING SIMPLE DECISIONS 247 16 MAKING SIMPLE DECISIONS Let us associate each state S with a numeric utility U(S), which expresses the desirability of the state A nondeterministic action A will have possible outcome states Result

More information

Chapter 10 Inventory Theory

Chapter 10 Inventory Theory Chapter 10 Inventory Theory 10.1. (a) Find the smallest n such that g(n) 0. g(1) = 3 g(2) =2 n = 2 (b) Find the smallest n such that g(n) 0. g(1) = 1 25 1 64 g(2) = 1 4 1 25 g(3) =1 1 4 g(4) = 1 16 1

More information

Strategies for Improving the Efficiency of Monte-Carlo Methods

Strategies for Improving the Efficiency of Monte-Carlo Methods Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful

More information

LECTURE 3: FREE CENTRAL LIMIT THEOREM AND FREE CUMULANTS

LECTURE 3: FREE CENTRAL LIMIT THEOREM AND FREE CUMULANTS LECTURE 3: FREE CENTRAL LIMIT THEOREM AND FREE CUMULANTS Recall from Lecture 2 that if (A, φ) is a non-commutative probability space and A 1,..., A n are subalgebras of A which are free with respect to

More information

Relevant parameter changes in structural break models

Relevant parameter changes in structural break models Relevant parameter changes in structural break models A. Dufays J. Rombouts Forecasting from Complexity April 27 th, 2018 1 Outline Sparse Change-Point models 1. Motivation 2. Model specification Shrinkage

More information

Introduction to the Maximum Likelihood Estimation Technique. September 24, 2015

Introduction to the Maximum Likelihood Estimation Technique. September 24, 2015 Introduction to the Maximum Likelihood Estimation Technique September 24, 2015 So far our Dependent Variable is Continuous That is, our outcome variable Y is assumed to follow a normal distribution having

More information

Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014

Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014 NTNU Page 1 of 5 Institutt for fysikk Contact during the exam: Professor Ingve Simonsen Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014 Allowed help: Alternativ D All written material This

More information

Trust Region Methods for Unconstrained Optimisation

Trust Region Methods for Unconstrained Optimisation Trust Region Methods for Unconstrained Optimisation Lecture 9, Numerical Linear Algebra and Optimisation Oxford University Computing Laboratory, MT 2007 Dr Raphael Hauser (hauser@comlab.ox.ac.uk) The Trust

More information

Probability distributions relevant to radiowave propagation modelling

Probability distributions relevant to radiowave propagation modelling Rec. ITU-R P.57 RECOMMENDATION ITU-R P.57 PROBABILITY DISTRIBUTIONS RELEVANT TO RADIOWAVE PROPAGATION MODELLING (994) Rec. ITU-R P.57 The ITU Radiocommunication Assembly, considering a) that the propagation

More information

A model reduction approach to numerical inversion for parabolic partial differential equations

A model reduction approach to numerical inversion for parabolic partial differential equations A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 3, Mikhail Zaslavsky 3 University of Michigan, Ann

More information

M.S. in Quantitative Finance & Risk Analytics (QFRA) Fall 2017 & Spring 2018

M.S. in Quantitative Finance & Risk Analytics (QFRA) Fall 2017 & Spring 2018 M.S. in Quantitative Finance & Risk Analytics (QFRA) Fall 2017 & Spring 2018 2 - Required Professional Development &Career Workshops MGMT 7770 Prof. Development Workshop 1/Career Workshops (Fall) Wed.

More information

Geostatistical Inference under Preferential Sampling

Geostatistical Inference under Preferential Sampling Geostatistical Inference under Preferential Sampling Marie Ozanne and Justin Strait Diggle, Menezes, and Su, 2010 October 12, 2015 Marie Ozanne and Justin Strait Preferential Sampling October 12, 2015

More information

Research Statement. Dapeng Zhan

Research Statement. Dapeng Zhan Research Statement Dapeng Zhan The Schramm-Loewner evolution (SLE), first introduced by Oded Schramm ([12]), is a oneparameter (κ (0, )) family of random non-self-crossing curves, which has received a

More information