ECE Spring Prof. David R. Jackson ECE Dept. Notes 20

Similar documents
LAPLACE TRANSFORM AND TRANSFER FUNCTION

ECE Spring Prof. David R. Jackson ECE Dept.

THE SINE INTEGRAL. x dt t

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Phys1112: DC and RC circuits

Echocardiography Project and Finite Fourier Series

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

f t te e = possesses a Laplace transform. Exercises for Module-III (Transform Calculus)

EECE 301 Signals & Systems Prof. Mark Fowler

Lab 10: RC, RL, and RLC Circuits

ECE Spring Prof. David R. Jackson ECE Dept. Notes 39

Differential Equations

CHAPTER 6: FIRST-ORDER CIRCUITS

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Math Final Exam Solutions

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

8. Basic RL and RC Circuits

EE243 Advanced Electromagnetic Theory Lec # 13: Waveguides and sources

The following report makes use of the process from Chapter 2 in Dr. Cumming s thesis.

10. State Space Methods

6.2 Transforms of Derivatives and Integrals.

Section 7.4 Modeling Changing Amplitude and Midline

F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

The Arcsine Distribution

Some Basic Information about M-S-D Systems

Chapter 2. First Order Scalar Equations

EECE 301 Signals & Systems Prof. Mark Fowler

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

Section 3.8, Mechanical and Electrical Vibrations

EE100 Lab 3 Experiment Guide: RC Circuits

Solutions to Assignment 1

5.1 - Logarithms and Their Properties

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

2. Nonlinear Conservation Law Equations

Basilio Bona ROBOTICA 03CFIOR 1

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Solutions from Chapter 9.1 and 9.2

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

Math 333 Problem Set #2 Solution 14 February 2003

EXERCISES FOR SECTION 1.5

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

= ( ) ) or a system of differential equations with continuous parametrization (T = R

3, so θ = arccos

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

1 Evaluating Chromatograms

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

Question 1: Question 2: Topology Exercise Sheet 3

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

Matlab and Python programming: how to get started

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Chapter 1 Fundamental Concepts

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

4.1 - Logarithms and Their Properties

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

Signal and System (Chapter 3. Continuous-Time Systems)

Let us start with a two dimensional case. We consider a vector ( x,

ECE 2100 Circuit Analysis

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Final Spring 2007

Q1) [20 points] answer for the following questions (ON THIS SHEET):

Traveling Waves. Chapter Introduction

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

RC, RL and RLC circuits

ME 391 Mechanical Engineering Analysis

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Chapter 3 Boundary Value Problem

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Chapter 8 The Complete Response of RL and RC Circuits

Západočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

6.003: Signal Processing

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

ADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations

EE 301 Lab 2 Convolution

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

5.2. The Natural Logarithm. Solution

ON THE DEGREES OF RATIONAL KNOTS

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

IE1206 Embedded Electronics

Second Order Linear Differential Equations

15. Vector Valued Functions

SOLUTIONS TO ECE 3084

ECE Spring Prof. David R. Jackson ECE Dept. Notes 37

Roller-Coaster Coordinate System

The average rate of change between two points on a function is d t

Theory of! Partial Differential Equations!

Ordinary dierential equations

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

Y 0.4Y 0.45Y Y to a proper ARMA specification.

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Transcription:

ECE 6345 Spring 5 Prof. David R. Jacson ECE Dep. Noes

Overview In his se of noes we apply he SDI mehod o invesigae he fields produced by a pach curren. We calculae he field due o a recangular pach on op of a subsrae. We examine he pole and branch poin singulariies in he complex plane. We examine he pah of inegraion in he complex plane.

y Pach Fields W x J sx π x xˆ cos L L Find E x (x,y,) z J sx ε r L h x 3

From Noes 9 we have ( π ) Pach Fields (con.) + + jx ( x + y y ) E G z J e d d (,, ) (, ) x xx x y sx x y x y where ~ G xx x V i ( z) y V TE i ( z) For he pach curren, we have: L cos x π W J sx ( x, y ) LW sinc y π L x 4

Pach Fields (con.) From he TEN we have V T i T () () Zin () T Yin () Y jy co h ( ) T T z T denoes or TE z / / z h T Z V / T T Z ( ) + - x y Amp ( z or TE z ) z Y Y Y Y TE TE ωε z ωε z z ωµ z ωµ 5

Pach Fields (con.) Define he denominaor erm as T T D ( ) Y () in so ha D ( ) Y jy co h ( ) z D ( ) Y jy co h ( ) TE TE TE z 6

Pach Fields (con.) We hen have ~ G xx + x y TE D D The final form of he elecric field a he inerface is hen E ( xy,,) + y x ( ) + π D D jx ( x + yy) J e d d x TE (, ) sx x y x y 7

Polar Coordinaes Use he following change of variables: y d d x y d dφ ( is also ofen called ρ ) φ x + π + ( ) jx ( x + y y ) jx ( x + y y ) e d d e d dφ ( ) x y π / + 4 cos cos ( ) ( x) ( ) yddφ x y cosφ sinφ y x Advanage: The poles and branch poins are locaed a a fixed posiion in he complex plane. z / / z 8

Hence Polar Coordinaes (con.) + y x jx ( x + y y ) E ( x, y, ) J (, TE ) e d d ( π ) + D D x sx x y x y π / ( ) ( ) Ex xy,, J, cos sin sx φ φ φ TE π + D ( ) ( ) D This is in he following general form: cos ( x) cos( ) yddφ x y π / (, φ) E F d dφ x 9

Poles Poles occur when eiher of he following condiions are saisfied: D ( ) ( ) D TE ( ) ( TE ) p p z : D ( ) Y jy co h z

+ V ( ) V ( ) Y + I ( ) I ( ) Y Hence Y Y Poles (con.) This coincides wih he well-nown Transverse Resonance Equaion (TRE) for deermining he characerisic equaion of a guided mode. (e.g, SW mode) V I + ( ) V ( ) + ( ) I( ) Y Y (Kirchhoff s laws) so ha ( ) jy co h Y z SW: β z β z / / Z z Z I( z ) + V( z) -

Comparison: Poles (con.) Poles in plane ( ) Y jy co h z TRE (surface-wave mode) ( ) Y jy co h z Y Y ωε z ωε z Y Y ωε z ωε z z / p z / p β z β z / / (A similar comparison holds for he TE case)

Poles (con.) Hence, we have he conclusion ha p TE p β β TE Tha is, he poles are locaed a he wavenumbers of he guided modes (he surface-wave modes). Noe: In mos pracical cases, here is only a surface-wave mode. 3

Poles (con.) The complex plane hus has poles on he real axis a he wavenumbers of he surface waves. Im + x y β β Re 4

Pah of Inegraion The pah avoids he poles by going above hem. Lossy case Im β C Re Lossy case Im This pah can be used for numerical compuaion. h R L R β C Re 5

Pah of Inegraion (con.) The pah avoids he poles by going above hem. Lossless case Im h R L R C Re h. 5 R L. R (ypical choices) Pracical noe: If h R is oo small, we are oo close o he pole. If h R is oo large, here is oo much round-off error due o exponenial growh in he sin and cos funcions. 6

Branch Poins To explain why we have branch poins, consider he funcion: wih ( ) co ( ) D Y jy h z ωε ωε j co z z ( h) z z z ( ) ( ) / / Noe: There are no branch cus for z (he funcion D is an even funcion of z ). If z z D ( ) changes (We need branch cus for z.) 7

Branch Poins (con.) z ( ) / ( ) ( ) / / + ( ) ( ) / / j + ( ) / ( ) ( ) j / Noe: The represenaion of he square roo of as j is arbirary here. Im ( ) Re 8

Branch Poins (con.) ( ) / ( ) ( ) j z / jφ jφ ( ) j e e / / Branch cus are necessary o preven he angles from changing by π : Noe: The shape of he branch cus is arbirary, bu verical cus are shown here. Im φ φ Re Noe: The branch cus should no cross he real axis when here is loss in he air (he inegrand mus be coninuous). 9

Branch Poins (con.) We obain he correc signs for z if we choose he following branches: ( ) ( ) π /< Arg < 3 π / ( ) 3 π /< Arg < π / z j( + ) Im z jφ jφ ( ) j e e z ( ) / / + j( + ) z The wave is hen eiher decaying or ougoing in he air region when we are on he real axis. Re

Branch Poins (con.) The wavenumber z is hen uniquely defined everywhere in he complex plane: z jφ jφ ( ) j e e / / Example φ.9π φ 5 π /4 Im φ π /4 φ.π π /< φ < 3 π / 3 π / < φ < π / φ Arg ( ) ( ) ( ) φ Arg Re φ.π φ 3 π /4 φ π /4 φ.π

Riemann Surface The Riemann surface is a pair of complex planes, conneced by ramps (where he branch cus used o be). The angles (and hence he funcion) change coninuously over he surface. All possible values of he funcion are found on he surface.

Riemann Surface (con.) Riemann surface for z / Top shee π < φ < π Boom shee y r φ z j re φ x π < φ < 3π Top view MATLAB : π < φ π Noe: A horizonal branch cu has been arbirarily chosen. A ramp now exiss where he branch cu used o be. 3

3D view Riemann Surface (con.) Riemann surface for z / x Top y B D Boom D B y B D D Side view B Top view x 4

Riemann Surface The Riemann surface can be consruced for he wavenumber funcion. / ( ) ( ) ( ) / / ( ) j z Example: Im We go couner-clocwise around he branch poin a saring on he op shee on he real axis, and end up bac where we sared. Top shee o 3 o 45 Re φ π /4 φ π /6 Boom shee φ π /4+ π φ π /6 5

Sommerfeld Branch Cus Sommerfeld branch cus are a convenien choice for heoreical purposes (discussed more in ECE 634): Im ( ) on branch cu z Im z ( ) / Im ( ) < z (everywhere in complex plane) Re Assume : z j( + ) 6

Sommerfeld Branch Cus (con.) Im z ( ) / Im ( ) < z (everywhere in complex plane) Re z is real here, for a lossless air. < < Pracical noe: If we give he air a small amoun of loss, we can simply chec o mae sure ha Im( z ) <. Noe: The branch poins move off of he axes for a lossy air. 7

Sommerfeld Branch Cus (con.) The Riemann surface wih Sommerfeld branch cus: Im Im ( ) < z Im ( ) > z (op shee) (boom shee) Re Noe: Surface wave poles mus lie on he op shee, and leay-wave poles mus lie on he boom shee. 8

Pah of Inegraion where Im π / ( φ ) Ex I dφ ( φ) (, ) φ I F d Top shee Im ( ) < z h R L R β C Re 9