FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection
|
|
- Winifred Wilcox
- 5 years ago
- Views:
Transcription
1 FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection
2 Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection
3 Math Review Symmetry: Given a unitary representation of a group G on a vector space V, we say that a vector v V is invariant under the action of G if for a g G: ρ g v = v The set of G-invariant vector V G is a vector space.
4 Math Review Symmetry: The inear map π G is a projection onto V G, if: π G v V G for a v V π G v = v for a v V G v, w π G w = 0 for a v V G, w V. The map π G is the map sending a vector v to the cosest G-invariant vector.
5 Math Review Symmetry: The measure of symmetry of a vector v with respect to the group G is the size of its projection onto the space of G-invariant vectors: Sym 2 v, G = π G v 2
6 Math Review Convoution: Given two functions f(p) and g(p), we define the convoution of the two functions to be: f g q = f p g q p dp
7 Math Review Convoution: If we hod the function g fixed we can define a map from the space of functions back into itsef: C g f = f g Caim: The map C g is a inear operator.
8 Math Review Convoution: If we hod the function g fixed we can define a map from the space of functions back into itsef: C g f = f g Caim: Given functions f and h and scaars α and β: C g αf + βh q = α f p + β h p g q p dp = α f p g q p dp + β h p g q p dp = α C g f + β C g h
9 Math Review Convoution: Assume that the function g is rea-vaued and radia, i.e. the vaue of g at a point p is competey determined by the distance of p from the origin: g p = g( p ) Exampe: The function g is a Gaussian f g f g =
10 Math Review Convoution: Assume that the function g is rea-vaued and radia, i.e. the vaue of g at a point p is competey determined by the distance of p from the origin: g p = g( p ) Caim: In this case, C g is sef-adjoint (i.e. symmetric).
11 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Expanding the eft side, we get: C g f, h = C g f p h(p) dp
12 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Writing out the operator C g, we get: C g f, h = C g f p h(p) dp C g f, h = f g p h p dp
13 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Expressing the convoution as an integra gives: C g f, h = f g p h p dp C g f, h = f q g p q dq h p dp
14 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Changing the order of integration, we get: C g f, h = f q g p q dq h p dp C g f, h = f q g p q h p dp dq
15 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Using the fact that g is rea-vaued and radia: C g f, h = f q g p q h p dp dq C g f, h = f q h p g(q p) dp dq
16 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Using the equation for convoution, we get: C g f, h = f q h p g(q p) dp dq C g f, h = f q h g q dq
17 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h Using the equation for C g, we get: C g f, h = f q h g q dq C g f, h = f q C g h q dq
18 Math Review Convoution: Proof: We need to show that for any functions f and h: C g f, h = f, C g h And finay, using the equation for the dot-product: C g f, h = f q C g h q dq C g f, h = f, C g h
19 Outine Math Review Spherica Convoution Axia Symmetry Detection
20 Spherica Convoution/Correation In the case of the circe we used convoution / correation for two different tasks:
21 Spherica Convoution/Correation In the case of the circe we used convoution / correation for two different tasks: 1. We used convoution for operations ike smoothing circuar functions
22 Spherica Convoution/Correation In the case of the circe we used convoution / correation for two different tasks: 1. We used convoution for operations ike smoothing circuar functions 2. We used correation for operations ike aignment and symmetry detection k-fod Symmetry
23 Spherica Convoution/Correation Up to now, we thought of these two operations as essentiay the same. The situation changes as we move to functions on a sphere.
24 Spherica Convoution/Correation When we perform an operation ike smoothing, the input is: A function on the circe defining the signa, and A function on the circe defining the smoothing fiter The output of the operation is: A function on the circe = Signa Fiter Smoothed Signa
25 Spherica Convoution/Correation When we perform an operation ike aignment, the input is: Two functions on a circe The output is: A function on the space of 2D rotations
26 Spherica Convoution/Correation In the case of a circe, the situation is simper because the space of rotations is itsef a circe: There is a one-to-one mapping from points on a circe to rotations, with a point on a circe with ange θ corresponding to a rotation by an ange of θ. In the case of the sphere, the situation becomes more compicated: The sphere is a 2D space whie the rotations are a 3D space, so there can t be a one-to-one mapping.
27 Spherica Convoution In the case of a circe, we compute the vaue of the smoothed function at p by rotating the fiter so that (1,0) maps to p and then we compute the inner product of the signa with the rotated fiter. = *
28 Spherica Convoution In the case of a circe, we compute the vaue of the smoothed function at p by rotating the fiter so that (1,0) maps to p and then we compute the inner product of the signa with the rotated fiter. p =,
29 Spherica Convoution We can try an appy the same type of approach to the case of spherica functions. (0,1,0) Signa Fiter
30 Spherica Convoution We woud ike to define a new function on the sphere whose vaue at the point p is obtained by: Finding a rotation R that maps the North poe to p and then compute the inner product of the signa with the rotated fiter. p, R(0,1,0) Signa Rotated Fiter
31 Spherica Convoution The probem is that there are many different rotations that send the North poe to the point p, so this does not ead to a we-defined notion of smoothing. p, S(0,1,0) Signa Rotated Fiter
32 Spherica Convoution Reca: If we have two rotations R and S mapping the North poe to the point p, the rotations must differ by an initia rotation about the y-axis: S = R R y ψ (0,1,0) S R 1 p
33 Spherica Convoution Reca: Thus, we can make the notion of smoothing wedefined by ensuring that the initia rotation about the y-axis does not change the fiter.
34 Spherica Convoution Reca: This means that we can extend the circuar notion of smoothing to the sphere if we ensure that the fiter is symmetric about the y-axis: (0,1,0) Signa Fiter
35 Spherica Convoution Reca: This means that we can extend the circuar notion of smoothing to the sphere if we ensure that the fiter is symmetric about the y-axis: (0,1,0) If R and S are rotations mapping the North poe to p, then the rotation of the fiter by either R or S wi give Signa the same spherica Fiter function!
36 Spherica Convoution Convoution: Using the Euer ange representation, we know that the rotation taking the North poe to the point p = Φ(θ, φ) is the rotation: R θ, φ = R y θ R z φ y z x p θ φ
37 Spherica Convoution Convoution: Thus, given A spherica function f θ, φ A spherica fiter g(θ, φ) that is rotationaysymmetric about the y-axis The convoution of f with g at p = Φ θ, φ can be expressed by rotating g so the North poe gets mapped to p and computing the inner product: f g θ, φ = f, ρ R θ,φ (g)
38 Spherica Convoution Convoution: Expressing the spherica functions f and g in terms of the spherica harmonic basis, we get: f θ, φ = f, m Y m (θ, φ) g θ, φ = m= m= g, m Y m (θ, φ)
39 Spherica Convoution Convoution: Reca that the spherica harmonics can be expressed as a compex exponentia in θ times a poynomia in cos φ: Y m θ, φ = P m cos φ e imθ So a rotation by α degrees about the y-axis acts on the (, m)-th spherica harmonics by: ρ Ry α Y m = e imα Y m
40 Spherica Convoution Convoution: Thus, if the fiter g is rotationay symmetric about the y-axis, any rotation about the y-axis must not change g. That is, for a α we must have: ρ Ry α g = g Or in terms of the spherica harmonics: m= g, m Y m (θ, φ) = m= g, m = g, m e imα g, m e imα Y m θ, φ
41 Spherica Convoution Convoution: g, m = g, m e imα Thus, either: e imα = 1 for a α m = 0, or g, m = 0 Thus, in terms of the spherica harmonics, we get: g θ, φ = g, 0 Y 0 θ, φ
42 Spherica Convoution Y 0 0 θ, φ Im Y 1 1 θ, φ Y 1 0 θ, φ Re Y 1 1 θ, φ Im Y 2 2 θ, φ Im Y 2 1 θ, φ Y 2 0 θ, φ Re Y 2 1 θ, φ Re Y 2 2 θ, φ Im Y 3 3 θ, φ Im Y 3 2 θ, φ Im Y 3 1 θ, φ Y 3 0 θ, φ Re Y 2 1 θ, φ Re Y 3 2 θ, φ Re Y 3 3 θ, φ
43 Spherica Convoution g θ, φ = g, 0 Y 0 (θ, φ) Y 0 0 θ, φ Im Y 1 1 θ, φ Y 1 0 θ, φ Re Y 1 1 θ, φ Im Y 2 2 θ, φ Im Y 2 1 θ, φ Y 2 0 θ, φ Re Y 2 1 θ, φ Re Y 2 2 θ, φ Im Y 3 3 θ, φ Im Y 3 2 θ, φ Im Y 3 1 θ, φ Y 3 0 θ, φ Re Y 2 1 θ, φ Re Y 3 2 θ, φ Re Y 3 3 θ, φ
44 Spherica Convoution Convoution: Thus, the expression for the functions in terms of their spherica harmonic decomposition becomes: f θ, φ = f, m Y m (θ, φ) m= g θ, φ = g, 0 Y 0 (θ, φ) and we get an expression for the convoution: f g θ, φ = f, m Y m 0, ρ R θ,φ g, 0 Y m=
45 Spherica Convoution Convoution: By everaging the conjugate-inearity of the inner product and using the fact that the transformation ρ R is inear, we get: f g θ, φ = f g θ, φ = m=, m= f, m Y m, ρ R θ,φ g, 0 Y 0 f, m g(, 0) Y m, ρ R θ,φ Y 0
46 Spherica Convoution Convoution: Additionay, we know that: A rotation of an -th frequency function wi sti be an -th frequency function The space of -th frequency functions is orthogona to the space of -th frequency functions (if ) Thus, for a, we have: Y m, ρ R Y m = 0
47 Spherica Convoution Convoution: This ets us simpify the expression for the convoution: f g θ, φ = f g θ, φ =, m= m= f, m g(, 0) Y m, ρ R θ,φ Y 0 f, m g(, 0) Y m, ρ R θ,φ Y 0
48 Spherica Convoution Convoution: f g θ, φ = m= f, m g(, 0) Y m, ρ R θ,φ Y 0 To compute the convoution, we need to be abe to evauate the inner product: Y m, ρ R θ,φ Y 0
49 Spherica Convoution Convoution: What is the meaning of the function: Y m, ρ R θ,φ Y 0 This is a function on the sphere whose vaue at the point p = Φ(θ, φ) is the inner product of Y m with the rotation of Y 0, where the rotation takes the North poe to p.
50 Spherica Convoution Convoution: We woud ike to show that this function acts very simpy on the spherica harmonics: Y m, ρ R θ,φ Y 0 = λ Y m θ, φ Im Y 2 2 θ, φ Im Y 2 1 θ, φ Y 2 0 θ, φ Re Y 2 1 θ, φ Re Y 2 2 θ, φ
51 Spherica Convoution Convoution: Let s consider the operator C that maps spherica functions to spherica functions, defined by: C f θ, φ = f, ρ R θ,φ Y 0 As before, it turns out this map is a symmetric inear operator on the space of functions. Thus, there exists an orthonorma basis with respect to which C is diagona.
52 Spherica Convoution Convoution: Let s consider the operator C that maps spherica functions to spherica functions, defined by: C f θ, φ = f, ρ R θ,φ Y 0 This operator aso has the property that it commutes with rotations: Rotating a spherica function and then convoving with Y 0 is the same as first convoving with Y 0 and then rotating.
53 Spherica Convoution Convoution: So, as with the Lapacian, we have a case in which we are given a symmetric operator which commutes with rotations.
54 Spherica Convoution L: a symmetric operator R SO(3): a rotation V λ : the space of e. functions of L with e.vaue λ R L f = L R f λ R f = L R f f V λ R f V λ
55 Spherica Convoution Convoution: Thus, the subspace of -th frequency functions is a space of functions that are eigenvectors of C, a with the same eigenvaue: C Y m = λ, Y m Thus we have: C Y m = Y m 0 if λ otherwise
56 Spherica Convoution Convoution: Putting this a together, we get: Y m, ρ R θ,φ Y 0 = λ Y m (θ, φ) Thus, the equation for the convoution becomes: f g θ, φ = f g θ, φ = m= m= f, m g, 0 Y m, ρ R θ,φ Y 0 f, m g, 0 λ Y m (θ, φ)
57 Spherica Convoution Convoution: f g θ, φ = m= f, m g, 0 λ Y m (θ, φ) Thus, the convoution of f with g can be obtained by mutipying the (, m)-th spherica harmonic coefficients of f by λ g(, 0). As in the case of functions on a circe, this means that convoution in the spatia domain amounts to mutipication in the frequency domain.
58 Spherica Convoution Convoution: In order to be abe to use the convoution theorem for spherica functions, we need to know what the eigenvaues λ are.
59 Spherica Convoution Convoution: In order to be abe to use the convoution theorem for spherica functions, we need to know what the eigenvaues λ are. It turns out that these are: λ = 4π 2 + 1
60 Spherica Convoution Convoution: Which gives us the equation: f g, m = 4π f, m g(, 0) =
61 Outine Math Review Spherica Convoution Axia Symmetry Detection
62 Axia Symmetry Detection Given a spherica function f, we woud ike to compute the measure of the axia symmetry of f with respect to every axis through the origin. f(θ, φ) AxiaSym f (θ, φ)
63 Axia Symmetry Detection What is the measure of the axia symmetry of f about the y-axis? f(θ, φ)
64 Axia Symmetry Detection What is the measure of the axia symmetry of f about the y-axis? We know that f is axiay symmetric about the yaxis if it can be expressed as the sum of the Y 0 : Y 0 0 θ, φ Im Y 1 1 θ, φ Y 1 0 θ, φ Re Y 1 1 θ, φ Im Y 2 2 θ, φ Im Y 2 1 θ, φ Y 2 0 θ, φ Re Y 2 1 θ, φ Re Y 2 2 θ, φ
65 Axia Symmetry Detection What is the measure of the axia symmetry of f about the y-axis? We know that f is axiay symmetric about the yaxis if it can be expressed as the sum of the Y 0. We aso know that for m 0: Y 0, Y m = 0 So the projection onto the space of functions that are axiay symmetric about the y-axis is: π y k m= f, m Y m = f, 0 Y 0
66 Axia Symmetry Detection What is the measure of the axia symmetry of f about the y-axis? Thus, the measure of the axia symmetry of f about the y-axis is defined as: YAxiaSym 2 f = = f, 0 Y 0 f, 0 2 2
67 Axia Symmetry Detection More generay, we woud ike to be abe to compute the measure of the axia symmetry of f with respect to any axis. To compute the symmetry measure about the ine through p = Φ(θ, φ) we: Rotate so that p goes to the North poe, and Compute the symmetry measure about the y-axis. R(p) R p p
68 Axia Symmetry Detection More generay, we woud ike to be abe to compute the measure of the axia symmetry of f with respect Since to the any rotation axis. R(θ, φ) maps the To compute North the poe symmetry to p, the measure rotation we about are the ine through interested p = Φ(θ, in φ) is we: the inverse, R 1 (θ, φ). Rotate so that p goes to the North poe, and Compute the symmetry measure about the y-axis. R(p) R p p
69 Axia Symmetry Detection Using the fact that the spherica harmonics form an orthonorma basis, we know that the (, m)-th harmonic coefficient of f is defined by: f, m = f, Y m Thus, to compute the measure of axia symmetry about the axis through p we need to compute: AxiaSym f 2 θ, φ = ρ R 1 θ,φ f, Y 0 2
70 Axia Symmetry Detection Using the fact that ρ is a unitary representation we can re-write this equation as: AxiaSym f 2 θ, φ = AxiaSym f 2 θ, φ = ρ R 1 θ,φ f, Y 0 2 f, ρ R θ,φ Y 0 2
71 Axia Symmetry Detection Expressing f in terms of its spherica harmonic decomposition, we get: AxiaSym f 2 θ, φ = AxiaSym f 2 θ, φ = m= f, ρ R θ,φ Y 0 2 f(, m) Y m, ρ R θ,φ Y 0 2
72 Axia Symmetry Detection Appying the identity: Y m, ρ R θ,φ Y 0 = 4π Y m θ, φ we get an expression for the symmetry measure: AxiaSym f 2 θ, φ = AxiaSym f 2 θ, φ = 4π m= f(, m) Y m, ρ R θ,φ Y 0 m= f, m Y m θ, φ 2 2
73 Axia Symmetry Detection AxiaSym f 2 θ, φ = 4π m= f, m Y m θ, φ 2 Thus, the measure of axia symmetry can be computed by taking the weighted sum of the squares of the frequency components of f.
74 Axia Symmetry Detection AxiaSym f 2 θ, φ = 4π m= f, m Y m θ, φ 2 Initia Function Frequency Decomposition = = Axia Symmetry Descriptor Weighted Square Norms
FFTs in Graphics and Vision. Fast Alignment of Spherical Functions
FFTs in Graphics and Vision Fast Alignment of Spherical Functions Outline Math Review Fast Rotational Alignment Review Recall 1: We can represent any rotation R in terms of the triplet of Euler angles
More informationLECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationDavid Eigen. MA112 Final Paper. May 10, 2002
David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More information2M2. Fourier Series Prof Bill Lionheart
M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationLECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationLegendre Polynomials - Lecture 8
Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha
More informationAnouncements. Assignment 3 has been posted!
Anouncements Assignment 3 has been posted! FFTs in Graphics and Vision Correlation of Spherical Functions Outline Math Review Spherical Correlation Review Dimensionality: Given a complex n-dimensional
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationFourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form
Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationPhysics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions
Physics 116C Hemhotz s an Lapace s Equations in Spherica Poar Coorinates: Spherica Harmonics an Spherica Besse Functions Peter Young Date: October 28, 2013) I. HELMHOLTZ S EQUATION As iscusse in cass,
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationSTATISTICS ON HILBERT S SIXTEENTH PROBLEM
STATISTICS ON HILBERT S SIXTEENTH PROBLEM ANTONIO LERARIO AND ERIK LUNDBERG Abstract. We study the statistics of the number of connected components and the voume of a random rea agebraic hypersurface in
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationPhysics 505 Fall Homework Assignment #4 Solutions
Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments
More informationPhysics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27
Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationRadiation Fields. Lecture 12
Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence
More informationThe state vectors j, m transform in rotations like D(R) j, m = m j, m j, m D(R) j, m. m m (R) = j, m exp. where. d (j) m m (β) j, m exp ij )
Anguar momentum agebra It is easy to see that the operat J J x J x + J y J y + J z J z commutes with the operats J x, J y and J z, [J, J i ] 0 We choose the component J z and denote the common eigenstate
More informationWeek 6 Lectures, Math 6451, Tanveer
Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n
More informationMidterm 2 Review. Drew Rollins
Midterm 2 Review Drew Roins 1 Centra Potentias and Spherica Coordinates 1.1 separation of variabes Soving centra force probems in physics (physica systems described by two objects with a force between
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationForces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment
Forces of Friction When an object is in motion on a surface or through a viscous medium, there wi be a resistance to the motion This is due to the interactions between the object and its environment This
More informationTRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS
Vo. 39 (008) ACTA PHYSICA POLONICA B No 8 TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS Zbigniew Romanowski Interdiscipinary Centre for Materias Modeing Pawińskiego 5a, 0-106 Warsaw, Poand
More informationSampling with Bessel Functions
amping with Besse Functions K.I. Kou, T. Qian and F. ommen Abstract. The paper deas with samping of σ-bandimited functions in R m with Cifford-vaued, where bandimitedness means that the spectrum is contained
More informationMath 220B - Summer 2003 Homework 1 Solutions
Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary
More informationComputing Spherical Transform and Convolution on the 2-Sphere
Computing Spherica Transform and Convoution on the 2-Sphere Boon Thye Thomas Yeo ythomas@mit.edu May, 25 Abstract We propose a simpe extension to the Least-Squares method of projecting sampes of an unknown
More informationThe Hessian of Axially Symmetric Functions on SE(3) and Application in 3D Image Analysis
The Hessian of Axiay Symmetric Functions on SE(3) and Appication in 3D Image Anaysis M.H.J. Janssen 1, T.C.J. Dea Haije 1, F.C. Martin 1, E.J. Bekkers 1, R. Duits 1 1 Eindhoven University of Technoogy,
More information11.1 One-dimensional Helmholtz Equation
Chapter Green s Functions. One-dimensiona Hemhotz Equation Suppose we have a string driven by an externa force, periodic with frequency ω. The differentia equation here f is some prescribed function) 2
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationLecture 8 February 18, 2010
Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More information3D and 6D Fast Rotational Matching
3D and 6D Fast Rotationa Matching Juio Kovacs, Ph.D. Department of Moecuar Bioogy The Scripps Research Institute 10550 N. Torrey Pines Road, Mai TPC6 La Joa, Caifornia, 9037 Situs Modeing Workshop, San
More informationMAT 167: Advanced Linear Algebra
< Proem 1 (15 pts) MAT 167: Advanced Linear Agera Fina Exam Soutions (a) (5 pts) State the definition of a unitary matrix and expain the difference etween an orthogona matrix and an unitary matrix. Soution:
More informationFFTs in Graphics and Vision. The Laplace Operator
FFTs in Graphics and Vision The Laplace Operator 1 Outline Math Stuff Symmetric/Hermitian Matrices Lagrange Multipliers Diagonalizing Symmetric Matrices The Laplacian Operator 2 Linear Operators Definition:
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationFFTs in Graphics and Vision. Rotational and Reflective Symmetry Detection
FFTs in Graphics and Vision Rotational and Reflective Symmetry Detection Outline Representation Theory Symmetry Detection Rotational Symmetry Reflective Symmetry Representation Theory Recall: A group is
More informationCluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects.
Stear Dynamics & Structure of Gaaxies handout # Custer modeing Just as a sef-gravitating coection of objects. Coisions Do we have to worry about coisions? Gobuar custers ook densest, so obtain a rough
More informationON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS
SIAM J. NUMER. ANAL. c 1992 Society for Industria Appied Mathematics Vo. 6, No. 6, pp. 1716-1740, December 1992 011 ON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS G. BEYLKIN
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationDipartimento di Matematica
Dipartimento di Matematica V. CASARINO TWO-PARAMETER ESTIMATES FOR JOINT SPECTRAL PROJECTIONS ON S n Rapporto interno N. 3, febbraio 7 Poitecnico di Torino Corso Duca degi Abruzzi, 4-9 Torino-Itaia TWO-PARAMETER
More informationIntroduction to LMTO method
1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More information1.2 Partial Wave Analysis
February, 205 Lecture X.2 Partia Wave Anaysis We have described scattering in terms of an incoming pane wave, a momentum eigenet, and and outgoing spherica wave, aso with definite momentum. We now consider
More informationRestricted weak type on maximal linear and multilinear integral maps.
Restricted weak type on maxima inear and mutiinear integra maps. Oscar Basco Abstract It is shown that mutiinear operators of the form T (f 1,..., f k )(x) = R K(x, y n 1,..., y k )f 1 (y 1 )...f k (y
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationGaussian Curvature in a p-orbital, Hydrogen-like Atoms
Advanced Studies in Theoretica Physics Vo. 9, 015, no. 6, 81-85 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/astp.015.5115 Gaussian Curvature in a p-orbita, Hydrogen-ike Atoms Sandro-Jose Berrio-Guzman
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f( q ) = exp ( βv( q )), which is obtained by summing over
More informationORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION
J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationParallel-Axis Theorem
Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationLecture 17 - The Secrets we have Swept Under the Rug
Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics
More informationTHE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES
THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia
More information(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].
PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform
More informationSrednicki Chapter 51
Srednici Chapter 51 QFT Probems & Soutions A. George September 7, 13 Srednici 51.1. Derive the fermion-oop correction to the scaar proagator by woring through equation 5., and show that it has an extra
More informationA nodal collocation approximation for the multidimensional P L equations. 3D applications.
XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.
More informationT.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA
ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network
More informationDiscrete Techniques. Chapter Introduction
Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationDiscrete Techniques. Chapter Introduction
Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various
More informationANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP
ANALOG OF HEAT EQUATION FOR GAUSSIAN MEASURE OF A BALL IN HILBERT SPACE GYULA PAP ABSTRACT. If µ is a Gaussian measure on a Hibert space with mean a and covariance operator T, and r is a} fixed positive
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f(q ) = exp ( βv( q )) 1, which is obtained by summing over
More information18. Atmospheric scattering details
8. Atmospheric scattering detais See Chandrasekhar for copious detais and aso Goody & Yung Chapters 7 (Mie scattering) and 8. Legendre poynomias are often convenient in scattering probems to expand the
More informationIntroduction to Riemann Solvers
CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid
More informationIn Coulomb gauge, the vector potential is then given by
Physics 505 Fa 007 Homework Assignment #8 Soutions Textbook probems: Ch. 5: 5.13, 5.14, 5.15, 5.16 5.13 A sphere of raius a carries a uniform surface-charge istribution σ. The sphere is rotate about a
More informationDeterminantal point process models on the sphere
Determinanta point process modes on the sphere Jesper Møer 1, Morten Niesen 1, Emiio Porcu 2 and Ege Rubak 1 1 Department of Mathematica Sciences, Aaborg University, Denmark jm@math.aau.dk. mniesen@math.aau.dk,
More informationPolar Snakes: a fast and robust parametric active contour model
Poar Snakes: a fast and robust parametric active contour mode Christophe Coewet To cite this version: Christophe Coewet. Poar Snakes: a fast and robust parametric active contour mode. IEEE Int. Conf. on
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationRelated Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage
Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is
More informationJoel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS.
Joe Broida UCSD Fa 009 Phys 30B QM II Homework Set DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. You may need to use one or more of these: Y 0 0 = 4π Y 0 = 3 4π cos Y
More informationMat 1501 lecture notes, penultimate installment
Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the
More information1 Heat Equation Dirichlet Boundary Conditions
Chapter 3 Heat Exampes in Rectanges Heat Equation Dirichet Boundary Conditions u t (x, t) = ku xx (x, t), < x (.) u(, t) =, u(, t) = u(x, ) = f(x). Separate Variabes Look for simpe soutions in the
More informationHILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is
More informationAbsolute Value Preconditioning for Symmetric Indefinite Linear Systems
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce
More informationTerm Test AER301F. Dynamics. 5 November The value of each question is indicated in the table opposite.
U N I V E R S I T Y O F T O R O N T O Facuty of Appied Science and Engineering Term Test AER31F Dynamics 5 November 212 Student Name: Last Name First Names Student Number: Instructions: 1. Attempt a questions.
More informationFRIEZE GROUPS IN R 2
FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the
More informationarxiv: v3 [math.ca] 8 Nov 2018
RESTRICTIONS OF HIGHER DERIVATIVES OF THE FOURIER TRANSFORM MICHAEL GOLDBERG AND DMITRIY STOLYAROV arxiv:1809.04159v3 [math.ca] 8 Nov 018 Abstract. We consider severa probems reated to the restriction
More information