Least squares and motion. Nuno Vasconcelos ECE Department, UCSD

Similar documents
Least squares and motion. Nuno Vasconcelos ECE Department, UCSD

Lecture 12: Introduction to nonlinear optics II.

MECE 3320 Measurements & Instrumentation. Static and Dynamic Characteristics of Signals

Control Systems (Lecture note #6)

Dr. Junchao Xia Center of Biophysics and Computational Biology. Fall /21/2016 1/23

Mathematical Statistics. Chapter VIII Sampling Distributions and the Central Limit Theorem

CHAPTER: 3 INVERSE EXPONENTIAL DISTRIBUTION: DIFFERENT METHOD OF ESTIMATIONS

Anouncements. Conjugate Gradients. Steepest Descent. Outline. Steepest Descent. Steepest Descent

Chap 2: Reliability and Availability Models

Almost unbiased exponential estimator for the finite population mean

Improved Exponential Estimator for Population Variance Using Two Auxiliary Variables

Gauge Theories. Elementary Particle Physics Strong Interaction Fenomenology. Diego Bettoni Academic year

Improved Exponential Estimator for Population Variance Using Two Auxiliary Variables

COMPLEX NUMBERS AND ELEMENTARY FUNCTIONS OF COMPLEX VARIABLES

3.4 Properties of the Stress Tensor

8. Queueing systems. Contents. Simple teletraffic model. Pure queueing system

LECTURE 6 TRANSFORMATION OF RANDOM VARIABLES

An N-Component Series Repairable System with Repairman Doing Other Work and Priority in Repair

The Variance-Covariance Matrix

Two-Dimensional Quantum Harmonic Oscillator

Almost Unbiased Exponential Estimator for the Finite Population Mean

k of the incident wave) will be greater t is too small to satisfy the required kinematics boundary condition, (19)

Algorithms to Solve Singularly Perturbed Volterra Integral Equations

Major: All Engineering Majors. Authors: Autar Kaw, Luke Snyder

Frequency Response. Response of an LTI System to Eigenfunction

The Linear Regression Of Weighted Segments

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

Introduction to Laplace Transforms October 25, 2017

1973 AP Calculus BC: Section I

Introduction to logistic regression

Continous system: differential equations

Survival Analysis for Randomized Clinical Trials II Cox Regression. Ziad Taib Biostatistics AstraZeneca February 26, 2008

Consider a system of 2 simultaneous first order linear equations

State Observer Design

Let's revisit conditional probability, where the event M is expressed in terms of the random variable. P Ax x x = =

Phys Nov. 3, 2017 Today s Topics. Continue Chapter 2: Electromagnetic Theory, Photons, and Light Reading for Next Time

Supplementary Figure 1. Experiment and simulation with finite qudit. anharmonicity. (a), Experimental data taken after a 60 ns three-tone pulse.

Series of New Information Divergences, Properties and Corresponding Series of Metric Spaces

Boyce/DiPrima 9 th ed, Ch 7.6: Complex Eigenvalues

9. Simple Rules for Monetary Policy

Reliability analysis of time - dependent stress - strength system when the number of cycles follows binomial distribution

Delay-Dependent State Estimation for Time Delay Systems

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

Chapter 4. Continuous Time Markov Chains. Babita Goyal

Variable Satellite Usage in GPS Receivers

Unbalanced Panel Data Models

Lecture Y4: Computational Optics I

Summary: Solving a Homogeneous System of Two Linear First Order Equations in Two Unknowns

Estimation Theory. Chapter 4

Midterm Examination (100 pts)

Solution set Stat 471/Spring 06. Homework 2

Entropy Equation for a Control Volume

Fourier Series: main points

For the plane motion of a rigid body, an additional equation is needed to specify the state of rotation of the body.

innovations shocks white noise

(1) Cov(, ) E[( E( ))( E( ))]

The far field calculation: Approximate and exact solutions. Persa Kyritsi November 10th, 2005 B2-109

Density estimation III. Linear regression.

( A) ( B) ( C) ( D) ( E)

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

CHAPTER 7. X and 2 = X

Lecture 4: Laplace Transforms

Fundamentals of Speech Recognition Suggested Project The Hidden Markov Model

The Method of Steepest Descent for Feedforward Artificial Neural Networks

Inference on Curved Poisson Distribution Using its Statistical Curvature

Numerical Method: Finite difference scheme

Course 10 Shading. 1. Basic Concepts: Radiance: the light energy. Light Source:

Part B: Transform Methods. Professor E. Ambikairajah UNSW, Australia

Convergence tests for the cluster DFT calculations

Handout on. Crystal Symmetries and Energy Bands

Mellin Transform Method for the Valuation of the American Power Put Option with Non-Dividend and Dividend Yields

The rise of neural networks. Deep networks. Why many layers? Why many layers? Why many layers? 24/03/2017

Midterm exam 2, April 7, 2009 (solutions)

Counting the compositions of a positive integer n using Generating Functions Start with, 1. x = 3 ), the number of compositions of 4.

Option Pricing in a Fractional Brownian Motion Environment

The algebraic immunity of a class of correlation immune H Boolean functions

Classification. Linear Classification. What is a Linear Disciminant? Representing Classes. Decision Boundaries. What can be expressed?

EMPIRICAL STUDY IN FINITE CORRELATION COEFFICIENT IN TWO PHASE ESTIMATION

Notation for Mixed Models for Finite Populations

Response of LTI Systems to Complex Exponentials

Existence of Nonoscillatory Solutions for a Class of N-order Neutral Differential Systems

Linear System Review. Linear System Review. Descriptions of Linear Systems: 2008 Spring ME854 - GGZ Page 1

ELEN E4830 Digital Image Processing

A NOVEL DIFFERENCE EQUATION REPRESENTATION FOR AUTOREGRESSIVE TIME SERIES

STRUCTURAL FAULT DETECTION OF BRIDGES BASED ON LINEAR SYSTEM PARAMETER AND MTS METHOD

Asymptotic Behavior of Solutions of Nonlinear Delay Differential Equations With Impulse

EE 232 Lightwave Devices. Photodiodes

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

Final Exam : Solutions

Overview. Splay trees. Balanced binary search trees. Inge Li Gørtz. Self-adjusting BST (Sleator-Tarjan 1983).

QR factorization. Let P 1, P 2, P n-1, be matrices such that Pn 1Pn 2... PPA

Part I- Wave Reflection and Transmission at Normal Incident. Part II- Wave Reflection and Transmission at Oblique Incident

Midterm Exam. Tuesday, September hour, 15 minutes

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

Homework: Introduction to Motion

Comparisons of the Variance of Predictors with PPS sampling (update of c04ed26.doc) Ed Stanek

MODEL QUESTION. Statistics (Theory) (New Syllabus) dx OR, If M is the mode of a discrete probability distribution with mass function f

EEE 303: Signals and Linear Systems

Mathematical Formulation

Reliability of time dependent stress-strength system for various distributions

Transcription:

Las squars ad moo uo Vascoclos ECE Dparm UCSD

Pla for oda oda w wll dscuss moo smao hs s rsg wo was moo s vr usful as a cu for rcogo sgmao comprsso c. s a gra ampl of las squars problm w wll also wrap up dscusso o las squars roduc wo ps of moo smao bloc machg dffral mhods wll al abou moo ambgus ad local vs global moo

Las squars a las squars problm s o whr w hav wo varabls XY rlad b a uow fuco Y gx a rag s D { } a modl Y f;φ whr Φ s a vcor of paramrs h goal s: o fd h modl paramrs ha lad o h bs appromao o h obsrvd daa.. o drm h caocal ampl s h problm of fg a l o a s of pos hr * ε m Φ [ f Φ ] Φ a b ad f ; a b a + b 3

Two ma cass o-lar las squars fφ o lar o Φ.g. lar las squars fφ lar o Φ.g. f ; Φ s f ; Φ s o: all ha mars s lar o Φ boh olar o ohr lar modls: polomals spls ural wors Fourr dcomposos c. 4

5 o-lar las squars mos dffcul cas opmal soluo f ad ol f: grad of ε s zro Hssa of ε gav df gral hs has o closd form umrcal soluo.g. grad dsc pc al sma Φ ra < Φ z z z T ε Φ ε ε ε Φ ε ε ε ε ε L L Φ Φ Φ Φ + ε α εφ Φ ε

6 Lar las squars closd form soluo wr soluo s gv b ormal quaos.g. for a l f; + Φ Γ Φ Φ f f L T T Γ Γ Γ Φ Γ L

7 Vr powrful Q: wha s h bs lar appromao of a po squc b DFT sl poals? o g las squars soluo w d Γ 4 4 L L L

8 Bs Fourr appromao hs mas ha hs s orhoormal.. Γ T Γ ad.. h bs appromao ar h DFT coffcs assocad wh h poals Γ 4 4 L L L T T T Γ Γ Γ Γ Φ...

Sgal appromao Q: wha s h bad-pass flr h whos oupu bs appromas a sgal h frquc rag Ω? w hav s ha mus hav DFT Y X hc opmal flr has DFT Ω ohrws Ω H ohrws.. s h dal bad-pass flr of bad Ω uv: dal bs appromao LS ss! 9

oo smao s a mpora praccal ampl of LS problms ma applcaos: rcogo: ma vs ar characrzd b h p of moo.g. walg vs rug srog clus abou sc srucur.g. wh w roa a 3D obc moo of a pl drmd b how far h 3D po s from camra sgmao hgs ha mov oghr blog o h sam obc algm oc w ow h moo w ca alg mags a squc.g. h ASA paoramas comprsso sma moo alg mags rasm ol rror c

oo smao cosdr h followg wo mags m

oo smao cosdr h followg wo mags m +

oo smao goal: gv mags ad + for ach pl fd uv whch mmzs dffrc [ u v ] D + problm: mpossbl o solv from o pl alo wo uows uv o quao uv -u-v m 3

Fudamal law moo ca ol b solvd ovr a ghborhood d a las wo pls mas ss o cosdr mor ad mmz h avrag rror hs s las squars ε [ u v + ] uv -u-v m 4

Bloc machg fac s a o-lar las squars problm sc -u-v s a o-lar fuco of uv soluo : bloc machg for ach bloc + do a hausv sarch for h closs mach vr commo comprsso.g. PEG 5

Bloc machg s compuaoall sv d o compu h squard rror bw h bloc ad a collco of blocs h prvous mag dos o alwas produc good moo smas.g. ma machs ca b quall good hs s a problm for all moo smao mhods: moo ca b ambguous wh masurd locall.g. b machg wdows? 6

oo ambgus clarl w cao drm h moo of a fla ghborhood for a dg ghborhood w ca ol drm o of h wo compos h wo compos ar uqul dfd ol wh h ghborhood coas D mag srucur hs s calld h aprur problm?????? 7

8 Dffral mhods w ca a las lma h compl problm b loog for a closd-form soluo o problm: hs s a o-lar fuco of uv soluo: clarl h problm s du o hs quao ca b mad lar o uv b a Talor srs appromao v u + [ ] + d d v u * m ε v u v u

9 Dffral mhods whch lads o o: w ow how o compu hs rms A s h dffrc bw coscuv frams B s.. a fuco of h mag grad 4444 4 3 4444 4 444 4 3 444 4 B A v u + + + T A + T v u B T

Dffral mhods w hus hav ad h las squars problm s o: sc s cosa w om hs s ow lar las squars w ca us us our formula rcall ha + u v * ε + [ + u + v ]

Lar las squars f h h LS soluo s: wr soluo s gv b ormal quaos [ ] * m Φ Φ f ε Φ Γ Φ Φ f f L T T Γ Γ Γ Φ

Las squars soluo for moo sad of w hav ad wr [ ] * Φ f ε [ ] * + + v u ε Φ Φ v u f f

3 Las squars soluo h ormal quaos ar ladg o h soluo v u L L L L v u

4 Las squars soluo wh s hs wll dfd? o ha s h mar ha w usd o dc corrs s vrbl ol wh s wo gvalus ar o-zro H T wdow

prsg Wdows rcall: cosa wdow small gvalus dg wdow o mdum o small flow wdow o larg o small corr wdow wo larg gvalus H H H H hs cofrms wha w had alrad s: moo ca ol b compud uambguousl wh h ghborhood coas D formao.g. corrs 5

6 summar [UV] lsm w compu grads - for ach pl l wdow compu ma U u V v rur UV v u { } w w w w + +

7 Problms rcall w usd h Talor srs appromao hs s a good appromao ol for small uv o avod hs problm w d o us pramds v u v u +

Hrarchcal smao algorhm: do moo smao usg ad o oba u v warp wh u v : wpd -u -v up-sampl b o g wpd warp moo smao u v do moo smao usg ad o oba u v warp wh u v warp moo smao u v c. wpd... 8

Hrarchcal smao ach sag mprovs h mach u v soluo: upsampl all u v o full rsoluo add o oba uv u v + o ha small dsplacms a low rsoluo ar larg dsplacms a full rsoluo combs lar wh abl o sma larg dsplacms... 9

oo modls so far w hav dal local moo ach pl movs b slf raslao + u v local moo s h mos grc.g. r lavs blowg h wd o mpora alrav cas s ha of global moo moo of all pls sasfs o commo quao usuall du o camra moo: pag roao zoomg zoom roao 3

3 mpora cass po a m warpd o po a m + mpora global moos ar raslao b uv roao b θ scalg b s s θ uv s + v u ' ' θ θ θ θ cos s s cos ' ' s s ' '

3 Aff rasformaos hs ar all spcal cass of h aff rasformao moo of r mag dscrbd b Φ abcdf T ca accou for raslao roao scalg ad shar + f d c b a ' ' raslao roao uform scal ouform scal sharg

33