General Strong Polarization

Similar documents
General Strong Polarization

General Strong Polarization

General Strong Polarization

Locality in Coding Theory

Uncertain Compression & Graph Coloring. Madhu Sudan Harvard

Property Testing and Affine Invariance Part I Madhu Sudan Harvard University

Low-Degree Polynomials

Variations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra

Algebraic Codes and Invariance

Local Decoding and Testing Polynomials over Grids

On The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time Delays

ECE 6540, Lecture 06 Sufficient Statistics & Complete Statistics Variations

Two Decades of Property Testing Madhu Sudan Harvard

Support Vector Machines. CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington

COMPRESSION FOR QUANTUM POPULATION CODING

Review for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa

Worksheets for GCSE Mathematics. Quadratics. mr-mathematics.com Maths Resources for Teachers. Algebra

CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I

1. The graph of a function f is given above. Answer the question: a. Find the value(s) of x where f is not differentiable. Ans: x = 4, x = 3, x = 2,

Estimate by the L 2 Norm of a Parameter Poisson Intensity Discontinuous

Module 7 (Lecture 27) RETAINING WALLS

Angular Momentum, Electromagnetic Waves

Work, Energy, and Power. Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition

Approximate Second Order Algorithms. Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo

Lesson 24: True and False Number Sentences

Lecture 6. Notes on Linear Algebra. Perceptron

Chapter 6 Reed-Solomon Codes. 6.1 Finite Field Algebra 6.2 Reed-Solomon Codes 6.3 Syndrome Based Decoding 6.4 Curve-Fitting Based Decoding

A A A A A A A A A A A A. a a a a a a a a a a a a a a a. Apples taste amazingly good.

Lecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher

Lecture 11. Kernel Methods

(1) Correspondence of the density matrix to traditional method

A linear operator of a new class of multivalent harmonic functions

Terms of Use. Copyright Embark on the Journey

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa

Module 7 (Lecture 25) RETAINING WALLS

Classical RSA algorithm

Rotational Motion. Chapter 10 of Essential University Physics, Richard Wolfson, 3 rd Edition

Charge carrier density in metals and semiconductors

Introduction to time series econometrics and VARs. Tom Holden PhD Macroeconomics, Semester 2

ME5286 Robotics Spring 2017 Quiz 2

(2) Orbital angular momentum

(1) Introduction: a new basis set

Lesson 7: Linear Transformations Applied to Cubes

CHAPTER 2 Special Theory of Relativity

10.4 The Cross Product

International Journal of Mathematical Archive-5(3), 2014, Available online through ISSN

SECTION 7: FAULT ANALYSIS. ESE 470 Energy Distribution Systems

Worksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra

7.3 The Jacobi and Gauss-Seidel Iterative Methods

Gradient expansion formalism for generic spin torques

Heteroskedasticity ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD

Lecture No. 1 Introduction to Method of Weighted Residuals. Solve the differential equation L (u) = p(x) in V where L is a differential operator

A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions

Lecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method

Charged-Particle Interactions in Matter

Quantum Mechanics. An essential theory to understand properties of matter and light. Chemical Electronic Magnetic Thermal Optical Etc.

Multiple Regression Analysis: Inference ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD

3. Mathematical Modelling

CHAPTER 4 Structure of the Atom

Data Preprocessing. Jilles Vreeken IRDM 15/ Oct 2015

TECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES

A Circular Cylinder. 1 Introduction. RITA CHOUDHURY 1 and BIBHASH DEKA 2

Physics 371 Spring 2017 Prof. Anlage Review

APPROXIMATING A COMMON SOLUTION OF A FINITE FAMILY OF GENERALIZED EQUILIBRIUM AND FIXED POINT PROBLEMS

Math, Stats, and Mathstats Review ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD

Grover s algorithm. We want to find aa. Search in an unordered database. QC oracle (as usual) Usual trick

Module 4 (Lecture 16) SHALLOW FOUNDATIONS: ALLOWABLE BEARING CAPACITY AND SETTLEMENT

Extreme value statistics: from one dimension to many. Lecture 1: one dimension Lecture 2: many dimensions

Hopfield Network for Associative Memory

Non-Parametric Weighted Tests for Change in Distribution Function

Math 171 Spring 2017 Final Exam. Problem Worth

A semi-device-independent framework based on natural physical assumptions

Inferring the origin of an epidemic with a dynamic message-passing algorithm

PHY103A: Lecture # 4

Topics. Module 3 Lecture 10 SHALLOW FOUNDATIONS: ULTIMATE BEARING CAPACITY NPTEL ADVANCED FOUNDATION ENGINEERING-I

Control of Mobile Robots

A Step Towards the Cognitive Radar: Target Detection under Nonstationary Clutter

PHY103A: Lecture # 9

Introduction to Density Estimation and Anomaly Detection. Tom Dietterich

" = Y(#,$) % R(r) = 1 4& % " = Y(#,$) % R(r) = Recitation Problems: Week 4. a. 5 B, b. 6. , Ne Mg + 15 P 2+ c. 23 V,

Last Name _Piatoles_ Given Name Americo ID Number

Quantum Chosen-Ciphertext Attacks against Feistel Ciphers

Statistical Learning with the Lasso, spring The Lasso

Independent Component Analysis and FastICA. Copyright Changwei Xiong June last update: July 7, 2016

The Elasticity of Quantum Spacetime Fabric. Viorel Laurentiu Cartas

Photon Interactions in Matter

9. Switched Capacitor Filters. Electronic Circuits. Prof. Dr. Qiuting Huang Integrated Systems Laboratory

Eureka Lessons for 6th Grade Unit FIVE ~ Equations & Inequalities

SECTION 7: STEADY-STATE ERROR. ESE 499 Feedback Control Systems

Secondary 3H Unit = 1 = 7. Lesson 3.3 Worksheet. Simplify: Lesson 3.6 Worksheet

Big Bang Planck Era. This theory: cosmological model of the universe that is best supported by several aspects of scientific evidence and observation

International Journal of Mathematical Archive-4(7), 2013, Available online through ISSN

Solving Fuzzy Nonlinear Equations by a General Iterative Method

Aeroacoustic optimisation by means of adjoint sensitivity maps

Expectation Propagation performs smooth gradient descent GUILLAUME DEHAENE

10.4 Controller synthesis using discrete-time model Example: comparison of various controllers

Quantum state measurement

Integrating Rational functions by the Method of Partial fraction Decomposition. Antony L. Foster

Chapter 22 : Electric potential

Lecture 11: Polar codes construction

Transcription:

General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Gurswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) May 1, 018 G.Tech: General Strong Polarization 1 of 1

This Talk: Martingales and Polarization [0,1]-Martingale: Sequence of r.v.s XX 0, XX 1, XX, XX tt 0,1, EE XX tt+1 XX 0,, XX tt = XX tt Well-studied in algorithms: [MotwaniRaghavan, 4.4], [MitzenmacherUpfal, 1] E.g., XX tt = Expected approx. factor of algorithm given first tt random choices. Usually study concentration E.g. whp XX final [XX 0 ± εε] Today: Polarization: whp XX final εε, 1 εε Why? Polar codes lim tt XX tt = XX 0 lim tt XX tt = Bern XX 0 May 1, 018 G.Tech: General Strong Polarization of 1

Main Result: Definition and Theorem Strong Polarization: (informally) Pr XX tt (ττ, 1 ττ) εε if tt max OO(log 1/εε), oo(log 1/ττ) formally γγ > 0 ββ < 1, cc ss. tt. tt Pr XX tt (γγ tt, 1 γγ tt ) cc ββ tt A non-polarizing example: XX 0 = 1 ; XX tt+1 = XX tt ± tt Local Polarization: Both definitions qualitative! Variance in the middle: XX tt (ττ, 1 ττ) ττ > 0 σσ > 0 ss. tt. tt, XX tt ττ, 1 ττ VVVVVV XX tt+1 XX tt Suction at the ends: XX tt ττ, 1 ττ θθ > 0, cc <, ττ > 0 ss. tt. XX tt < ττ Pr XX tt+1 < XX tt cc σσ θθ Theorem: Local Polarization low Strong end condition. Polarization. Similar condition for high end May 1, 018 G.Tech: General Strong Polarization 3 of 1

Rest of the talk Background/Motivation: The Shannon Capacity Challenge Resolution by [Arikan 08], [Guruswami-Xia 13], [Hassani,Alishahi,Urbanke 13]: Polar Codes! Key ideas: Polar Codes & Arikan-Martingale. Implications of Weak/Strong Polarization. Our Contributions: Simplification & Generalization Local Polarization Strong Polarization Arikan-Martingale Locally Polarizes (generally) May 1, 018 G.Tech: General Strong Polarization 4 of 1

Shannon and Channel Capacity Shannon [ 48]: For every (noisy) channel of communication, capacity CC s.t. communication at rate RR < CC is possible and RR > CC is not. Needs lot of formalization omitted. Example: XX FF BBBBBB(pp) Acts independently on bits XX w.p. 1 pp 1 XX w.p. pp Capacity = 1 HH(pp) ; HH(pp) = binary entropy! Rate? Encode:FF kk nn FF nn kk lim nn RR nn nn HH pp = pp log 1 pp + 1 pp log 1 1 pp May 1, 018 G.Tech: General Strong Polarization 5 of 1

Achieving Channel Capacity Commonly used phrase. Many mathematical interpretations. Some possibilities: How small can nn be? Get RR > CC εε with polytime algorithms? [Luby et al. 95] Make running time poly nn εε? Open till 008 Shannon 48 : nn = OO CC RR Forney 66 :time = pppppppp nn, 1 εε Arikan 08: Invented Polar Codes Final resolution: Guruswami+Xia 13, Hassani+Alishahi+Urbanke 13 Strong analysis May 1, 018 G.Tech: General Strong Polarization 6 of 1

Polar Codes and Polarization Arikan: Defined Polar Codes, one for every tt Associated Martingale XX 0,, XX tt, ttth code: Length nn = tt If Pr XX tt γγ, 1 δδ εε then ttth code is εε + δδ - close to capacity, and Pr DDDDDDDDDDDD EEEEEEEEEEEE mm + eeeeeeeeee mm nn γγ Need γγ = oo 1 nn (strong polarization γγ = nnnnnn(nn)) Need εε = 1/nn Ω(1) ( strong polarization) May 1, 018 G.Tech: General Strong Polarization 7 of 1

Context for today s talk [Arikan 08]: Martingale associated to general class of codes. Martingale polarizes (weakly) implies code achieves capacity (weakly) Martingale polarizes weakly for entire class. [GX 13, HAU 13]: Strong polarization for one specific code. Why so specific? Bottleneck: Strong Polarization Rest of talk: Polar codes, Martingale and Strong Polarization May 1, 018 G.Tech: General Strong Polarization 8 of 1

Lesson 0: Compression Coding Linear Compression Scheme: HH, DD for compression scheme for bern pp nn if Linear map HH: FF nn FF mm Pr ZZ bern pp nn DD HH ZZ ZZ = oo 1 Hopefully mm HH pp + εε nn Compression Coding Let HH be such that HH HH = 0; Encoder: XX HH XX Error-Corrector: YY = HH XX + ηη YY DD HH YY = ww.pp. 1 oo(1) HH XX May 1, 018 G.Tech: General Strong Polarization 9 of 1

Question: How to compress? Arikan s key idea: Start with Polarization Transform : UU, VV UU + VV, VV Invertible so does nothing? If UU, VV independent, then UU + VV more random than either VV UU + VV less random than either Iterate (ignoring conditioning) End with bits that are almost random, or almost determined (by others). Output totally random part to get compression! May 1, 018 G.Tech: General Strong Polarization 10 of 1

Information Theory Basics Entropy: Defn: for random variable XX Ω, HH XX = pp xx log 1 where pp pp xx = Pr XX = xx xx xx Ω Bounds: 0 HH XX log Ω Conditional Entropy Defn: for jointly distributed XX, YY HH XX YY = pp yy HH(XX YY = yy) yy Chain Rule: HH XX, YY = HH XX + HH(YY XX) Conditioning does not increase entropy: HH XX YY HH(XX) Equality Independence: HH XX YY = HH XX XX YY May 1, 018 G.Tech: General Strong Polarization 11 of 1

Iteration by Example: AA, BB, CC, DD AA + BB + CC + DD, BB + DD, CC + DD, DD EE, FF, GG, HH EE + FF + GG + HH, FF + HH, GG + HH, HH BB + DD HH BB + DD AA + BB + CC + DD) BB + DD + FF + HH HH BB + DD + FF + HH AA + BB + CC + DD, EE + FF + GG + HH FF + HH HH FF + HH EE + FF + GG + HH) HH FF + HH FF + HH AA + BB + CC + DD, EE + FF + GG + HH, BB + DD + FF + HH May 1, 018 G.Tech: General Strong Polarization 1 of 1

The Polarization Butterfly PP nn UU, VV = PPnn UU + VV, PPnn VV UU ZZ 1 ZZ ZZ 1 + +1 ZZ + + WW 1 WW Polarization SS nn HH WW nn SS WW SS 0 SS HH pp + εε nn + ZZ nn WWnn Compression EE ZZ = PP nn ZZ SS VV +1 + +1 + WWnn +1 WWnn + Encoding time = OO nn log nn (given SS) To be shown: ZZ nn ZZ nn WW nn 1. Decoding = OO nn log nn. Polarization! May 1, 018 G.Tech: General Strong Polarization 13 of 1

The Polarization Butterfly: Decoding ZZ 1 ZZ ZZ 1 + +1 ZZ + + WW 1 WW Decoding idea: Given SS, WW SS = PP nn UU, VV Compute UU + VV Bern pp pp Determine VV? VV UU + VV ii Bern(rr ii ) +1 + +1 + + ZZ nn WWnn WWnn +1 WWnn + Key idea: Stronger induction! - non-identical product distribution Decoding Algorithm: Given SS, WW SS = PP nn UU, VV, pp 1,, pp nn - Compute qq 1,.. qqnn pp 1 pp nn ZZ nn ZZ nn WW nn - Determine UU + VV = DD WW SS +, qq 1 qqnn - Compute rr 1 rrnn pp 1 pp nn ; UU + VV - Determine VV = DD WW SS, rr 1 rrnn May 1, 018 G.Tech: General Strong Polarization 14 of 1

Polarization? Martingale? ZZ 1 ZZ +1 + ZZ nn ZZ 1 + +1 ZZ + +1 + ZZ nn + ZZ nn + WW 1 WW ZZ ii,0 ZZ ii,1 Z ii,tt ZZ ii,log nn WWnn WWnn +1 WWnn + WW nn Idea: Follow XX tt = HH ZZ ii,tt ZZ <ii,tt for random ii Martingale? ZZ jj,tt ZZ kk,tt XX tt can t distinguish ZZ jj,tt+1 ZZ kk,tt+1 ii = jj from ii = kk HH ZZ jj,tt, ZZ kk,tt = HH ZZ jj,tt+1, ZZ kk,tt+1 Chain Rule EE XX tt+1 XX tt = XX tt May 1, 018 G.Tech: General Strong Polarization 15 of 1

Remaining Tasks: 1. Prove that XX tt HH ZZ ii,tt ZZ <ii,tt polarizes locally. Prove that Local polarization Strong Polarization. Recall Strong Polarization: (informally) Pr XX tt (ττ, 1 ττ) εε if tt max OO(log 1/εε), oo(log 1/ττ) formally γγ > 0 ββ < 1, cc ss. tt. tt Pr XX tt (γγ tt, 1 γγ tt ) cc ββ tt Local Polarization: Variance in the middle: XX tt (ττ, 1 ττ) ττ > 0 σσ > 0 ss. tt. tt, XX tt ττ, 1 ττ VVVVVV XX tt+1 XX tt σσ Suction at the ends: XX tt ττ, 1 ττ θθ > 0, cc <, ττ > 0 ss. tt. XX tt < ττ Pr XX tt+1 < XX tt cc θθ May 1, 018 G.Tech: General Strong Polarization 16 of 1

Local (Global) Strong Polarization Step 1: Medium Polarization: αα, ββ < 1 Pr XX tt ββ tt, 1 ββ tt αα tt Proof: Potential Φ tt = EE min XX tt, 1 XX tt Variance+Suction ηη < 1 s.t. Φ tt ηη φφ tt 1 (Aside: Why XX tt? Clearly XX tt won t work. And XX tt increases!) Φ tt ηη tt Pr XX tt ββ tt, 1 ββ tt ηη ββ tt May 1, 018 G.Tech: General Strong Polarization 17 of 1

Local (Global) Strong Polarization Step 1: Medium Polarization: αα, ββ < 1 Pr XX tt ββ tt, 1 ββ tt αα tt Step : Medium Polarization + tt-more steps Strong Polarization: Proof: Assume XX 0 αα tt ; Want XX tt γγ tt w.p. 1 OO αα 1 tt Pick cc large & ττ s.t. Pr XX ii < XX ii 1 cc XX ii 1 ττ θθ.1: Pr ii ss. tt. XX ii ττ ττ 1 αα tt (Doob s Inequality).: Let YY ii = log XX ii ; Assuming XX ii ττ for all ii YY ii YY ii 1 log cc w.p. θθ; & kk w.p. kk. EEEE YY tt tt log cc θθ + 1 tt log γγ Concentration! Pr YY tt tt log γγ αα tt QED! May 1, 018 G.Tech: General Strong Polarization 18 of 1

Polarization of Arikan Martingale Ignoring conditioning: XX tt = HH pp XX tt+1 = HH pp pp w.p. ½ = HH pp HH pp pp w.p. ½ Variance in the middle: Continuity of HH and separation of pp pp from pp Suction at high end: XX tt = 1 εε Suction at low end: XX tt+1 = 1 OO(εε ) w.p. ½ XX tt = HH pp XX tt+1 = HH pp log pp w.p. ½ May 1, 018 G.Tech: General Strong Polarization 19 of 1

General Strong Polarization? So far: PP nn = HH l for HH = 1 1 1 0 What about other matrices? E.g., PP nn = HH 3 l for HH 3 = Open till our work! 1 1 0 1 0 0 0 0 1 Necessary conditions: HH invertible, not-lowertriangular Our work: Necessary conditions suffice (for local polarization)!. May 1, 018 G.Tech: General Strong Polarization 0 of 1

Conclusion Simple General Analysis of Strong Polarization But main reason to study general polarization: Maybe some matrices polarize even faster! This analysis does not get us there. But hopefully following the (simpler) steps in specific cases can lead to improvements. May 1, 018 G.Tech: General Strong Polarization 1 of 1

Thank You! May 1, 018 G.Tech: General Strong Polarization of 1