Matrix Transformation and Solutions of Wave Equation of Free Electromagnetic Field

Size: px
Start display at page:

Download "Matrix Transformation and Solutions of Wave Equation of Free Electromagnetic Field"

Transcription

1 Matrix Transformation and Solutions of Wave Equation of Free Electromagnetic Field Xianzhao Zhong Meteorological College of Yunnan Province, Kunming, 6508, China Abstract In this paper, the generalized differential wave equation for free electromagnetic field is transformed and formulated by means of matrixes. Then Maxwell wave equation and the second form of wave equation are deduced by matrix transformation. The solutions of the wave equations are discussed. Finally, two differential equations of vibration are established and their solutions are discussed. Key words: matrix transformation, wave equation. PACS: De 1 Introduction By means of matrix establishment and transformation [1], the generalized wave equation for free electromagnetic field [] is transformed and formulated in terms of a diagonal matrix. Then both wave equations, Maxwell wave address: zhxzh46@163.com 1

2 equation [3] and the second form of wave equation are deduced from the the generalized wave equation in terms of matrixes. Solutions of the two kinds of wave equations are discussed. For Maxwell wave equation, the wave amplitude F 0 is independent of wave number k and wave frequency ω. For the generalized wave equation, it is a function of wave number k and wave frequency ω. Solving the second form of wave equation, two solutions, F 1 and F 1, are obtained. Multiplying F 1 by F 1 results in a function F that is identical to the solution of the generalized wave equation. Finally, a pair of differential vibration equations are established, whose solution product is identical to the solution of the generalized wave equation. Transforming the Generalized Wave Equation by Means of Matrix The generalized wave equation of free electromagnetic field is F 1 c C F t 1 F = 0, (1) c t where F is the electromagnetic field []. Now we transform Eq.(1) into diagonal matrix. Suppose C = p () and Thus Eq.(1) is rewritten as t = ε. (3) (p pε ε )F = 0. (4)

3 It is required from Eq.(4) that, for arbitrary F, p pε ε = 0. (5) Rewriting Eq.(5), we have p 1 pε 1 pε ε = 0, (6) which is then transformed into ( ) 1 1 p p ε 1 1 ε = 0. (7) Let 1 1 A = 1 (8) 1 whose eigenvalue equation is det(λe A) = 0. (9) Solving Eq.(9), two eigenvalues are obtained. They are λ 1 = 5 (10) and Solving the equation of eigenvector we obtain the eigenvectors as λ = 5. (11) (λe A) x 1 = 0, (1) x α 1 = ( 1 5 ) T (forλ 1 ) (13) 3

4 and α = (1 5 + ) T (forλ ). (14) Then two matrixes 1 1 Y = (15) and 1 5 Y T = (16) are built up in terms of the eigenvectors α 1 in Eq.(13) and α in Eq.(14). And the matrix A in Eq.(8) is transformed into a diagonal matrix Y T AY = = diag(λ 1 Λ ) (17) 5 in virtue of Eq.(15) and Eq.(16). Thus Eq.(4) is transformed into ( ) p p ε F 1 = 0. (18) 5 ε After transformation, the electromagnetic field has changed from the field F to the field F 1. Equation (18) is developed into 5 5(p ε )F 1 10(p + ε )F 1 = 0. (19) Substituting Eq.() and Eq.(3) into Eq.(19) and taking its first bracket, we have F 1 1 c F 1 t = 0, (0) which is the well-known Maxwell wave equation. Taking the second bracket of Eq.(19) results in (p + ε )F 1 = 0, (1) 4

5 which leads to the second form of wave equation. F c F 1 t = 0, () 3 Solutions of Generalized Wave Equation of Free Electromagnetic Field The solution of the free electromagnetic field formulated by Eq.(1) is F = F 0 exp(λkr λωt)exp(iλkr iλωt), (3) where k is the wave number and ω is the wave frequency. In the last section, we transform the the generalized wave equation and obtain the well-known Maxwell wave equation [3] and the the second form of wave equation by matrix transformation. From Maxwell wave equation we have the solution F 1 1 c F 1 t = 0 (4) F 1 = F 10 exp(ikr iωt), (5) where the wave amplitude F 10 is a constant quantity, independent of the wave number k and wave frequency ω. The other form of wave equation is F c F 1 t = 0. (6) By separation of variables of the spatial r and the temporal t [4], F 1 (r, t) = R(r)T (t), (7) 5

6 for the electromagnetic field F 1 in Eq.(6), there is the first form c R(r) = 1 T (t) = λ ω (8) R(r) r T (t) t where ω is supposed to be a constant and the λ a variable. For Eq.(8), there is a solution or rewritten For Eq.(6), there is the second form For Eq.(31), there is a solution or rewritten F 1 = R 0 T 0 exp(±λkr)exp(±iλωt) (9) F 1 = F 0exp(±λkr)exp(±iλωt). (30) c R(r) = 1 T (t) = λ ω. (31) R(r) r T (t) t F 1 = R 0 T 0 exp(±λωt)exp(±iλkr) (3) F 1 = F 0exp(±λωt)exp(±iλkr). (33) Now multiplying Eq.(30) by Eq.(33) results in F = F 1F 1 = F 0 exp(λkr λωt)exp(iλkr iλωt), (34) where k is the wave number and ω is the wave frequency ω. Equation (34) is identical to Eq.(3). By matrix transformation, both wave equations, Eq.(4) and Eq.(6), are derived from the generalized wave equation Eq.(1), thus the general quadratic form becomes a standard quadratic one. The generalized wave equation conforms perfectly with Maxwell wave equation and the second form of wave equation. 6

7 4 Two Vibration Equations and Their Solutions Now let us establish a pair of differential vibration equations for the electromagnetic field. The first vibration equation in terms of the temporal variable t is d F dt + (λω)df dt + λ ω F = 0. (35) In terms of the spatial variable, the second wave equation is F λk F + λ k F = 0. (36) The solution of Eq.(35) is F = F 0 exp( λωt)exp( iλωt). (37) Rewriting Eq.(37) as a function of t, we have T (t) = F 0 exp( λωt)exp( iλωt). (38) In the spheroidal coordinates, Eq.(36) of spherically symmetric becomes For Eq.(39) there is the solution d F dr (λk)df dr + λ k F = 0. (39) F = F 0 exp(λkr)exp(iλkr) (40) or rewritten R(r) = F 0 exp(λkr)exp(iλkr). (41) Evidently, for the free electromagnetic field F, both equations, Eq.(35) and Eq.(36), are vibration equations. Equation (35) is one of the function 7

8 of t and Eq.(36) is of the function of r. Multiplying Eq.(38) by Eq.(41), we arrive at F = R(r)T (t) = F 0 exp(λkr λωt)exp(iλkr iλωt). (4) Equation (4) is identical to Eq.(3). Subtracting Eq.(35) from Eq.(36), we get c F λωc F + λ ω F ( d F dt + λω df dt + λ ω F ) = 0. (43) Equation (43) needs to be rearranged because there are two variables, t and r, included in the equation. Then Eq.(43) is changed to c F F t λω(c F + F t ) + λ ω (F F ) = 0. (44) Considering the continuity equation for passive field, Eq.(44) becomes which is Maxwell wave equation. cf + F t = 0 (45) F 1 c F t = 0, (46) References [1] Gene. H. Golub, Charles F.Ven loan. MATRIX COMPUTATIONS. ISBN rd ed c Hopkins University press. [] X.Z. Zhong. Electromagnetic Field Equation and Field Wave Equation. vixra. org: [84] ,

9 [3] J.D. Jackson. CLASSICAL ELECTRODYNAMICS, 3rd ed. John-Wiley & Sons, Inc., [4] M.D. Weir, R.L. Finney and F.R. Giordano. THOMAS CALCULUS. 10th ed. Publishing as Pearson Education. Inc.,

ROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD

ROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD Progress In Electromagnetics Research Letters, Vol. 11, 103 11, 009 ROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD G.-Q. Zhou Department of Physics Wuhan University

More information

6.6 General Form of the Equation for a Linear Relation

6.6 General Form of the Equation for a Linear Relation 6.6 General Form of the Equation for a Linear Relation FOCUS Relate the graph of a line to its equation in general form. We can write an equation in different forms. y 0 6 5 y 10 = 0 An equation for this

More information

Least squares Solution of Homogeneous Equations

Least squares Solution of Homogeneous Equations Least squares Solution of Homogeneous Equations supportive text for teaching purposes Revision: 1.2, dated: December 15, 2005 Tomáš Svoboda Czech Technical University, Faculty of Electrical Engineering

More information

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and

More information

arxiv:quant-ph/ v1 11 Nov 1995

arxiv:quant-ph/ v1 11 Nov 1995 Quantization of electromagnetic fields in a circular cylindrical cavity arxiv:quant-ph/950v Nov 995 K. Kakazu Department of Physics, University of the Ryukyus, Okinawa 903-0, Japan Y. S. Kim Department

More information

Seminar 6: COUPLED HARMONIC OSCILLATORS

Seminar 6: COUPLED HARMONIC OSCILLATORS Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached

More information

Linearization of Differential Equation Models

Linearization of Differential Equation Models Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking

More information

NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009

NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009 NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and

More information

AA 242B / ME 242B: Mechanical Vibrations (Spring 2016)

AA 242B / ME 242B: Mechanical Vibrations (Spring 2016) AA 242B / ME 242B: Mechanical Vibrations (Spring 206) Solution of Homework #3 Control Tab Figure : Schematic for the control tab. Inadequacy of a static-test A static-test for measuring θ would ideally

More information

A Note on Inverse Iteration

A Note on Inverse Iteration A Note on Inverse Iteration Klaus Neymeyr Universität Rostock, Fachbereich Mathematik, Universitätsplatz 1, 18051 Rostock, Germany; SUMMARY Inverse iteration, if applied to a symmetric positive definite

More information

Modern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur

Modern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur Modern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur Lecture 09 Wave propagation in anisotropic media (Contd.) So, we have seen the various aspects of

More information

Course Code: MTH-S101 Breakup: 3 1 0 4 Course Name: Mathematics-I Course Details: Unit-I: Sequences & Series: Definition, Monotonic sequences, Bounded sequences, Convergent and Divergent Sequences Infinite

More information

Perturbation Theory 1

Perturbation Theory 1 Perturbation Theory 1 1 Expansion of Complete System Let s take a look of an expansion for the function in terms of the complete system : (1) In general, this expansion is possible for any complete set.

More information

Journal of Theoretics

Journal of Theoretics Journal of Theoretics THE MISSED PHYSICS OF SOURCE-FREE MAGNETIC FIELD Author: Bibhas R. De P. O. Box 21141, Castro Valley, CA 94546-9141 BibhasDe@aol.com ABSTRACT Maxwell s equations are shown to have

More information

1 Singular Value Decomposition and Principal Component

1 Singular Value Decomposition and Principal Component Singular Value Decomposition and Principal Component Analysis In these lectures we discuss the SVD and the PCA, two of the most widely used tools in machine learning. Principal Component Analysis (PCA)

More information

Class 15 : Electromagnetic Waves

Class 15 : Electromagnetic Waves Class 15 : Electromagnetic Waves Wave equations Why do electromagnetic waves arise? What are their properties? How do they transport energy from place to place? Recap (1) In a region of space containing

More information

Introduction to Polarization

Introduction to Polarization Phone: Ext 659, E-mail: hcchui@mail.ncku.edu.tw Fall/007 Introduction to Polarization Text Book: A Yariv and P Yeh, Photonics, Oxford (007) 1.6 Polarization States and Representations (Stokes Parameters

More information

Electromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space

Electromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space Electromagnetic Waves 1 1. Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space 1 Retarded Potentials For volume charge & current = 1 4πε

More information

9 Wave solution of Maxwells equations.

9 Wave solution of Maxwells equations. 9 Wave solution of Maxwells equations. Contents 9.1 Wave solution: Plane waves 9.2 Scalar Spherical waves 9.3 Cylindrical waves 9.4 Momentum and energy of the electromagnetic field Keywords: Plane, cylindrical

More information

Eigenvalues and Eigenvectors: An Introduction

Eigenvalues and Eigenvectors: An Introduction Eigenvalues and Eigenvectors: An Introduction The eigenvalue problem is a problem of considerable theoretical interest and wide-ranging application. For example, this problem is crucial in solving systems

More information

Theory of Vibrations in Stewart Platforms

Theory of Vibrations in Stewart Platforms Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a

More information

arxiv: v2 [physics.gen-ph] 20 Mar 2013

arxiv: v2 [physics.gen-ph] 20 Mar 2013 arxiv:129.3449v2 [physics.gen-ph] 2 Mar 213 Potential Theory in Classical Electrodynamics W. Engelhardt 1, retired from: Max-Planck-Institut für Plasmaphysik, D-85741 Garching, Germany Abstract In Maxwell

More information

Coupled Oscillators Monday, 29 October 2012 normal mode Figure 1:

Coupled Oscillators Monday, 29 October 2012 normal mode Figure 1: Coupled Oscillators Monday, 29 October 2012 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l m k θ 2

More information

MathQuest: Differential Equations

MathQuest: Differential Equations MathQuest: Differential Equations Solutions to Linear Systems. Consider the linear system given by dy dt = 4 True or False: Y e t t = is a solution. c False, but I am not very confident Y.. Consider the

More information

Modeling II Linear Stability Analysis and Wave Equa9ons

Modeling II Linear Stability Analysis and Wave Equa9ons Modeling II Linear Stability Analysis and Wave Equa9ons Nondimensional Equa9ons From previous lecture, we have a system of nondimensional PDEs: (21.1) (21.2) (21.3) where here the * sign has been dropped

More information

Copyright 2008 Pearson Education, Inc. Publishing as Pearson Addison-Wesley. Chapter 8 Section 6

Copyright 2008 Pearson Education, Inc. Publishing as Pearson Addison-Wesley. Chapter 8 Section 6 Copyright 008 Pearson Education, Inc. Publishing as Pearson Addison-Wesley Chapter 8 Section 6 8.6 Solving Equations with Radicals 1 3 4 Solve radical equations having square root radicals. Identify equations

More information

Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism

Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger

More information

Quantization of scalar fields

Quantization of scalar fields Quantization of scalar fields March 8, 06 We have introduced several distinct types of fields, with actions that give their field equations. These include scalar fields, S α ϕ α ϕ m ϕ d 4 x and complex

More information

Mathematics with Maple

Mathematics with Maple Mathematics with Maple A Comprehensive E-Book Harald Pleym Preface The main objective of these Maple worksheets, organized for use with all Maple versions from Maple 14, is to show how the computer algebra

More information

Notes on Special Functions

Notes on Special Functions Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.

More information

Part III. Interaction with Single Electrons - Plane Wave Orbits

Part III. Interaction with Single Electrons - Plane Wave Orbits Part III - Orbits 52 / 115 3 Motion of an Electron in an Electromagnetic 53 / 115 Single electron motion in EM plane wave Electron momentum in electromagnetic wave with fields E and B given by Lorentz

More information

Electromagnetic (EM) Waves

Electromagnetic (EM) Waves Electromagnetic (EM) Waves Short review on calculus vector Outline A. Various formulations of the Maxwell equation: 1. In a vacuum 2. In a vacuum without source charge 3. In a medium 4. In a dielectric

More information

Differential of the Exponential Map

Differential of the Exponential Map Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements

More information

H ( E) E ( H) = H B t

H ( E) E ( H) = H B t Chapter 5 Energy and Momentum The equations established so far describe the behavior of electric and magnetic fields. They are a direct consequence of Maxwell s equations and the properties of matter.

More information

Scalar electromagnetic integral equations

Scalar electromagnetic integral equations Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how

More information

arxiv: v1 [physics.class-ph] 8 Apr 2019

arxiv: v1 [physics.class-ph] 8 Apr 2019 Representation Independent Boundary Conditions for a Piecewise-Homogeneous Linear Magneto-dielectric Medium arxiv:1904.04679v1 [physics.class-ph] 8 Apr 019 Michael E. Crenshaw 1 Charles M. Bowden Research

More information

Phonons and lattice dynamics

Phonons and lattice dynamics Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed

More information

HOSTOS COMMUNITY COLLEGE DEPARTMENT OF MATHEMATICS

HOSTOS COMMUNITY COLLEGE DEPARTMENT OF MATHEMATICS HOSTOS COMMUNITY COLLEGE DEPARTMENT OF MATHEMATICS MAT 217 Linear Algebra CREDIT HOURS: 4.0 EQUATED HOURS: 4.0 CLASS HOURS: 4.0 PREREQUISITE: PRE/COREQUISITE: MAT 210 Calculus I MAT 220 Calculus II RECOMMENDED

More information

A Spectral Time-Domain Method for Computational Electrodynamics

A Spectral Time-Domain Method for Computational Electrodynamics A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation

More information

Lecture 4. Alexey Boyarsky. October 6, 2015

Lecture 4. Alexey Boyarsky. October 6, 2015 Lecture 4 Alexey Boyarsky October 6, 2015 1 Conservation laws and symmetries 1.1 Ignorable Coordinates During the motion of a mechanical system, the 2s quantities q i and q i, (i = 1, 2,..., s) which specify

More information

Math 1302, Week 8: Oscillations

Math 1302, Week 8: Oscillations Math 302, Week 8: Oscillations T y eq Y y = y eq + Y mg Figure : Simple harmonic motion. At equilibrium the string is of total length y eq. During the motion we let Y be the extension beyond equilibrium,

More information

In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid.

In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Coupled Oscillators Wednesday, 26 October 2011 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l k θ 2

More information

PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY MOORE PENROSE INVERSE OF BIDIAGONAL MATRICES. I

PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY MOORE PENROSE INVERSE OF BIDIAGONAL MATRICES. I PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Mathematical Sciences 205, 2, p. 20 M a t h e m a t i c s MOORE PENROSE INVERSE OF BIDIAGONAL MATRICES. I Yu. R. HAKOPIAN, S. S. ALEKSANYAN 2, Chair

More information

Lecture 12 Eigenvalue Problem. Review of Eigenvalues Some properties Power method Shift method Inverse power method Deflation QR Method

Lecture 12 Eigenvalue Problem. Review of Eigenvalues Some properties Power method Shift method Inverse power method Deflation QR Method Lecture Eigenvalue Problem Review of Eigenvalues Some properties Power method Shift method Inverse power method Deflation QR Method Eigenvalue Eigenvalue ( A I) If det( A I) (trivial solution) To obtain

More information

Chapter 6: The Rouse Model. The Bead (friction factor) and Spring (Gaussian entropy) Molecular Model:

Chapter 6: The Rouse Model. The Bead (friction factor) and Spring (Gaussian entropy) Molecular Model: G. R. Strobl, Chapter 6 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). R. B. Bird, R. C. Armstrong, O. Hassager, "Dynamics of Polymeric Liquids", Vol. 2, John Wiley and Sons (1977). M. Doi,

More information

Two-Body Problem. Central Potential. 1D Motion

Two-Body Problem. Central Potential. 1D Motion Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of

More information

MEMORANDUM-4. n id (n 2 + n2 S) 0 n n 0. Det[ɛ(ω, k)]=0 gives the Dispersion relation for waves in a cold magnetized plasma: ω 2 pα ω 2 cα ω2, ω 2

MEMORANDUM-4. n id (n 2 + n2 S) 0 n n 0. Det[ɛ(ω, k)]=0 gives the Dispersion relation for waves in a cold magnetized plasma: ω 2 pα ω 2 cα ω2, ω 2 Fundamental dispersion relation MEMORANDUM-4 n S) id n n id n + n S) 0 } n n 0 {{ n P } ɛω,k) E x E y E z = 0 Det[ɛω, k)]=0 gives the Dispersion relation for waves in a cold magnetized plasma: n P ) [

More information

AN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b

AN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b AN ITERATION In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b In this, A is an n n matrix and b R n.systemsof this form arise

More information

Appendix C. Modal Analysis of a Uniform Cantilever with a Tip Mass. C.1 Transverse Vibrations. Boundary-Value Problem

Appendix C. Modal Analysis of a Uniform Cantilever with a Tip Mass. C.1 Transverse Vibrations. Boundary-Value Problem Appendix C Modal Analysis of a Uniform Cantilever with a Tip Mass C.1 Transverse Vibrations The following analytical modal analysis is given for the linear transverse vibrations of an undamped Euler Bernoulli

More information

3.3 Introduction to Infinite Sequences and Series

3.3 Introduction to Infinite Sequences and Series 3.3 Introduction to Infinite Sequences and Series The concept of limits and the related concepts of sequences underscores most, if not all, of the topics which we now call calculus. Here we focus directly

More information

LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES

LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES Figure 3.47 a. Two-mass, linear vibration system with spring connections. b. Free-body diagrams. c. Alternative free-body

More information

Lecture Notes on the Gaussian Distribution

Lecture Notes on the Gaussian Distribution Lecture Notes on the Gaussian Distribution Hairong Qi The Gaussian distribution is also referred to as the normal distribution or the bell curve distribution for its bell-shaped density curve. There s

More information

Systems of Linear ODEs

Systems of Linear ODEs P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here

More information

2 u 1-D: 3-D: x + 2 u

2 u 1-D: 3-D: x + 2 u c 2013 C.S. Casari - Politecnico di Milano - Introduction to Nanoscience 2013-14 Onde 1 1 Waves 1.1 wave propagation 1.1.1 field Field: a physical quantity (measurable, at least in principle) function

More information

Directional Derivatives in the Plane

Directional Derivatives in the Plane Directional Derivatives in the Plane P. Sam Johnson April 10, 2017 P. Sam Johnson (NIT Karnataka) Directional Derivatives in the Plane April 10, 2017 1 / 30 Directional Derivatives in the Plane Let z =

More information

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent

More information

Guide for Ph.D. Area Examination in Applied Mathematics

Guide for Ph.D. Area Examination in Applied Mathematics Guide for Ph.D. Area Examination in Applied Mathematics (for graduate students in Purdue University s School of Mechanical Engineering) (revised Fall 2016) This is a 3 hour, closed book, written examination.

More information

Math 216 Final Exam 24 April, 2017

Math 216 Final Exam 24 April, 2017 Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that

More information

MathQuest: Differential Equations

MathQuest: Differential Equations MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown

More information

Diagonalization by a unitary similarity transformation

Diagonalization by a unitary similarity transformation Physics 116A Winter 2011 Diagonalization by a unitary similarity transformation In these notes, we will always assume that the vector space V is a complex n-dimensional space 1 Introduction A semi-simple

More information

PART 2: INTRODUCTION TO CONTINUUM MECHANICS

PART 2: INTRODUCTION TO CONTINUUM MECHANICS 7 PART : INTRODUCTION TO CONTINUUM MECHANICS In the following sections we develop some applications of tensor calculus in the areas of dynamics, elasticity, fluids and electricity and magnetism. We begin

More information

EIGENVALUE ANALYSIS OF SPHERICAL RESONANT CAVITY USING RADIAL BASIS FUNCTIONS

EIGENVALUE ANALYSIS OF SPHERICAL RESONANT CAVITY USING RADIAL BASIS FUNCTIONS Progress In Electromagnetics Research Letters, Vol. 24, 69 76, 2011 EIGENVALUE ANALYSIS OF SPHERICAL RESONANT CAVITY USING RADIAL BASIS FUNCTIONS S. J. Lai 1, *, B. Z. Wang 1, and Y. Duan 2 1 Institute

More information

SVD and its Application to Generalized Eigenvalue Problems. Thomas Melzer

SVD and its Application to Generalized Eigenvalue Problems. Thomas Melzer SVD and its Application to Generalized Eigenvalue Problems Thomas Melzer June 8, 2004 Contents 0.1 Singular Value Decomposition.................. 2 0.1.1 Range and Nullspace................... 3 0.1.2

More information

Solving Quadratic Equations

Solving Quadratic Equations Solving Quadratic Equations MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: solve quadratic equations by factoring, solve quadratic

More information

The Helmholtz Decomposition and the Coulomb Gauge

The Helmholtz Decomposition and the Coulomb Gauge The Helmholtz Decomposition and the Coulomb Gauge 1 Problem Kirk T. McDonald Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544 (April 17, 2008; updated February 7, 2015 Helmholtz showed

More information

POSITRON SCATTERING BY A COULOMB POTENTIAL. Abstract. The purpose of this short paper is to show how positrons are treated

POSITRON SCATTERING BY A COULOMB POTENTIAL. Abstract. The purpose of this short paper is to show how positrons are treated POSITRON SCATTERING BY A COULOMB POTENTIAL Abstract The purpose of this short paper is to show how positrons are treated in quantum electrodynamics, and to study how positron size affects scattering. The

More information

E = 1 c. where ρ is the charge density. The last equality means we can solve for φ in the usual way using Coulomb s Law:

E = 1 c. where ρ is the charge density. The last equality means we can solve for φ in the usual way using Coulomb s Law: Photons In this section, we want to quantize the vector potential. We will work in the Coulomb gauge, where A = 0. In gaussian units, the fields are defined as follows, E = 1 c A t φ, B = A, where A and

More information

Exercise Set 7.2. Skills

Exercise Set 7.2. Skills Orthogonally diagonalizable matrix Spectral decomposition (or eigenvalue decomposition) Schur decomposition Subdiagonal Upper Hessenburg form Upper Hessenburg decomposition Skills Be able to recognize

More information

REVIEW SESSION. Midterm 2

REVIEW SESSION. Midterm 2 REVIEW SESSION Midterm 2 Summary of Chapter 20 Magnets have north and south poles Like poles repel, unlike attract Unit of magnetic field: tesla Electric currents produce magnetic fields A magnetic field

More information

So far we have derived two electrostatic equations E = 0 (6.2) B = 0 (6.3) which are to be modified due to Faraday s observation,

So far we have derived two electrostatic equations E = 0 (6.2) B = 0 (6.3) which are to be modified due to Faraday s observation, Chapter 6 Maxwell Equations 6.1 Maxwell equations So far we have derived two electrostatic equations and two magnetostatics equations E = ρ ɛ 0 (6.1) E = 0 (6.2) B = 0 (6.3) B = µ 0 J (6.4) which are to

More information

AMC Lecture Notes. Justin Stevens. Cybermath Academy

AMC Lecture Notes. Justin Stevens. Cybermath Academy AMC Lecture Notes Justin Stevens Cybermath Academy Contents 1 Systems of Equations 1 1 Systems of Equations 1.1 Substitution Example 1.1. Solve the system below for x and y: x y = 4, 2x + y = 29. Solution.

More information

THE PHOTON DEMYSTIFIED

THE PHOTON DEMYSTIFIED The photon demystified by Bibhas De THE PHOTON DEMYSTIFIED (wwwchasedesignscom) by Bibhas De (wwwbibhasdecom) Posted 26 October 2009 THE-FOUNDATION-IS-THE-FRONTIER-SERIES Essays on Twenty-first Century

More information

Least Squares Regression

Least Squares Regression Least Squares Regression Chemical Engineering 2450 - Numerical Methods Given N data points x i, y i, i 1 N, and a function that we wish to fit to these data points, fx, we define S as the sum of the squared

More information

THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING

THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING CHINESE JOURNAL OF GEOPHYSICS Vol.47, No.2, 24, pp: 379 386 THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING XUE Guo-Qiang LI Xiu SONG Jian-Ping GUO Wen-Bo

More information

Data Analysis and Manifold Learning Lecture 2: Properties of Symmetric Matrices and Examples

Data Analysis and Manifold Learning Lecture 2: Properties of Symmetric Matrices and Examples Data Analysis and Manifold Learning Lecture 2: Properties of Symmetric Matrices and Examples Radu Horaud INRIA Grenoble Rhone-Alpes, France Radu.Horaud@inrialpes.fr http://perception.inrialpes.fr/ Outline

More information

22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes

22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration

More information

Review of similarity transformation and Singular Value Decomposition

Review of similarity transformation and Singular Value Decomposition Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm

More information

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks. Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the

More information

Ch 5.7: Series Solutions Near a Regular Singular Point, Part II

Ch 5.7: Series Solutions Near a Regular Singular Point, Part II Ch 5.7: Series Solutions Near a Regular Singular Point, Part II! Recall from Section 5.6 (Part I): The point x 0 = 0 is a regular singular point of with and corresponding Euler Equation! We assume solutions

More information

Edexcel New GCE A Level Maths workbook Solving Linear and Quadratic Simultaneous Equations.

Edexcel New GCE A Level Maths workbook Solving Linear and Quadratic Simultaneous Equations. Edexcel New GCE A Level Maths workbook Solving Linear and Quadratic Simultaneous Equations. Edited by: K V Kumaran kumarmaths.weebly.com 1 Solving linear simultaneous equations using the elimination method

More information

Mathematical Constraint on Functions with Continuous Second Partial Derivatives

Mathematical Constraint on Functions with Continuous Second Partial Derivatives 1 Mathematical Constraint on Functions with Continuous Second Partial Derivatives J.D. Franson Physics Department, University of Maryland, Baltimore County, Baltimore, MD 15 Abstract A new integral identity

More information

Electromagnetic energy and momentum

Electromagnetic energy and momentum Electromagnetic energy and momentum Conservation of energy: the Poynting vector In previous chapters of Jackson we have seen that the energy density of the electric eq. 4.89 in Jackson and magnetic eq.

More information

Basic Equations of Elasticity

Basic Equations of Elasticity A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ

More information

MOTT-RUTHERFORD SCATTERING AND BEYOND. Abstract. The electron charge is considered to be distributed or extended in

MOTT-RUTHERFORD SCATTERING AND BEYOND. Abstract. The electron charge is considered to be distributed or extended in MOTT-RUTHERFORD SCATTERING AND BEYOND Abstract The electron charge is considered to be distributed or extended in space. The differential of the electron charge is set equal to a function of the electron

More information

vibrations, light transmission, tuning guitar, design buildings and bridges, washing machine, Partial differential problems, water flow,...

vibrations, light transmission, tuning guitar, design buildings and bridges, washing machine, Partial differential problems, water flow,... 6 Eigenvalues Eigenvalues are a common part of our life: vibrations, light transmission, tuning guitar, design buildings and bridges, washing machine, Partial differential problems, water flow, The simplest

More information

11 Perturbation Theory

11 Perturbation Theory S.K. Saikin Oct. 8, 009 11 Perturbation Theory Content: Variational Principle. Time-Dependent Perturbation Theory. 11.1 Variational Principle Lecture 11 If we need to compute the ground state energy of

More information

Information About Ellipses

Information About Ellipses Information About Ellipses David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view

More information

CLASS OF BI-QUADRATIC (BQ) ELECTROMAGNETIC MEDIA

CLASS OF BI-QUADRATIC (BQ) ELECTROMAGNETIC MEDIA Progress In Electromagnetics Research B, Vol. 7, 281 297, 2008 CLASS OF BI-QUADRATIC (BQ) ELECTROMAGNETIC MEDIA I. V. Lindell Electromagnetics Group Department of Radio Science and Engineering Helsinki

More information

TMA4180 Solutions to recommended exercises in Chapter 3 of N&W

TMA4180 Solutions to recommended exercises in Chapter 3 of N&W TMA480 Solutions to recommended exercises in Chapter 3 of N&W Exercise 3. The steepest descent and Newtons method with the bactracing algorithm is implemented in rosenbroc_newton.m. With initial point

More information

Eigenvalues and Eigenvectors. Review: Invertibility. Eigenvalues and Eigenvectors. The Finite Dimensional Case. January 18, 2018

Eigenvalues and Eigenvectors. Review: Invertibility. Eigenvalues and Eigenvectors. The Finite Dimensional Case. January 18, 2018 January 18, 2018 Contents 1 2 3 4 Review 1 We looked at general determinant functions proved that they are all multiples of a special one, called det f (A) = f (I n ) det A. Review 1 We looked at general

More information

Mathematical Foundations -1- Difference equations. Systems of difference equations. Life cycle model 2. Phase diagram 4. Eigenvalue and eigenvector 5

Mathematical Foundations -1- Difference equations. Systems of difference equations. Life cycle model 2. Phase diagram 4. Eigenvalue and eigenvector 5 Mathematical Foundations -- Difference equations Systems of difference equations Life cycle model Phase diagram 4 Eigenvalue and eigenvector 5 The general two variable model 9 Complex eigenvalues 5 Stable

More information

arxiv: v4 [quant-ph] 9 Jun 2016

arxiv: v4 [quant-ph] 9 Jun 2016 Applying Classical Geometry Intuition to Quantum arxiv:101.030v4 [quant-ph] 9 Jun 016 Spin Dallin S. Durfee and James L. Archibald Department of Physics and Astronomy, Brigham Young University, Provo,

More information

3. Mathematical Properties of MDOF Systems

3. Mathematical Properties of MDOF Systems 3. Mathematical Properties of MDOF Systems 3.1 The Generalized Eigenvalue Problem Recall that the natural frequencies ω and modes a are found from [ - ω 2 M + K ] a = 0 or K a = ω 2 M a Where M and K are

More information

On an Inverse Problem for a Quadratic Eigenvalue Problem

On an Inverse Problem for a Quadratic Eigenvalue Problem International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 1, pp. 13 26 (2017) http://campus.mst.edu/ijde On an Inverse Problem for a Quadratic Eigenvalue Problem Ebru Ergun and Adil

More information

Lecture Presentation. Chapter 15. Chemical Equilibrium. James F. Kirby Quinnipiac University Hamden, CT Pearson Education

Lecture Presentation. Chapter 15. Chemical Equilibrium. James F. Kirby Quinnipiac University Hamden, CT Pearson Education Lecture Presentation Chapter 15 Chemical James F. Kirby Quinnipiac University Hamden, CT The Concept of N 2 O 4 (g) 2 NO 2 (g) Chemical equilibrium occurs when a reaction and its reverse reaction proceed

More information

MECHANICS LAB AM 317 EXP 8 FREE VIBRATION OF COUPLED PENDULUMS

MECHANICS LAB AM 317 EXP 8 FREE VIBRATION OF COUPLED PENDULUMS MECHANICS LAB AM 37 EXP 8 FREE VIBRATIN F CUPLED PENDULUMS I. BJECTIVES I. To observe the normal modes of oscillation of a two degree-of-freedom system. I. To determine the natural frequencies and mode

More information

2.6 The optimum filtering solution is defined by the Wiener-Hopf equation

2.6 The optimum filtering solution is defined by the Wiener-Hopf equation .6 The optimum filtering solution is defined by the Wiener-opf equation w o p for which the minimum mean-square error equals J min σ d p w o () Combine Eqs. and () into a single relation: σ d p p 1 w o

More information

Created by T. Madas VECTOR OPERATORS. Created by T. Madas

Created by T. Madas VECTOR OPERATORS. Created by T. Madas VECTOR OPERATORS GRADIENT gradϕ ϕ Question 1 A surface S is given by the Cartesian equation x 2 2 + y = 25. a) Draw a sketch of S, and describe it geometrically. b) Determine an equation of the tangent

More information

Physics 506 Winter 2004

Physics 506 Winter 2004 Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading

More information

Differentiate Early, Differentiate Often!

Differentiate Early, Differentiate Often! Differentiate Early, Differentiate Often! Robert Dawson Department of Mathematics and Computing Science Saint Mary s University. Halifax, Nova Scotia Canada B3H 3C3 Phone: 902-420-5796 Fax: 902-420-5035

More information