FIRST-ORDER FILTERS GENERALIZED TO THE FRACTIONAL DOMAIN
|
|
- Simon Austin
- 6 years ago
- Views:
Transcription
1 Journal of ircuits, Systems, an omputers Vol. 7, No. (8) c Worl Scientific Publishing ompany FIST-ODE FILTES GENEALIZED TO THE FATIONAL DOMAIN A. G. ADWAN,,A.M.SOLIMAN, ana.s.elwakil, Department of Engineering Mathematics, Faculty of Engineering, airo University, Egypt Department of Electronics an ommunications, Faculty of Engineering, airo University, Egypt Department of Electrical an omputer Engineering, University of Sharjah, P. O. Box 77, Emirates ahmegom@yahoo.com asoliman@ieee.org elwakil@ieee.org evise 8 July 7 Traitional continuous-time filters are of integer orer. However, using fractional calculus, filters may also be represente by the more general fractional-orer ifferential equations in which case integer-orer filters are only a tight subset of fractional-orer filters. In this work, we show that low-pass, high-pass, ban-pass, an all-pass filters can be realize with circuits incorporating a single fractance evice. We erive expressions for the pole frequencies, the quality factor, the right-phase frequencies, an the halfpower frequencies. Examples of fractional passive filters supporte by numerical an PSpice simulations are given. Keywors: Analog filters; fractional calculus; fractional-orer circuits; filter esign.. Introuction Filter esign is one of the very few areas of electrical engineering for which a complete esign theory exists. Whether passive or active, filters necessarily incorporate inuctors an capacitors, the total number of which ictates the filter orer. However, an inuctor or capacitor is not but a special case of the more general so-calle fractance evice; which is an electrical element whose impeance in the complex frequency omain is given by Z(s) =as α Z(j)=a α e j(πα/). 5 For the special case of α = this element represents an inuctor while for α = itrepresents a capacitor. In the range <α<, this element may generally be consiere to represent a fractional-orer inuctor while for the range <α<, it may be consiere to represent a fractional-orer capacitor. At α =, it represents the well-known frequency-epenent negative resistor (FDN). Although a physical 55
2 56 A. G. awan, A. M. Soliman & A. S. Elwakil fractance evice oes not yet exist in the form of a single commercial evice, it may be emulate via higher-orer passive or L trees, as escribe in efs. 3 for simulation purposes. a However, very recently the authors of ef. 6 have escribe an emonstrate a single apparatus which preforms as a physical fractional-orer capacitor. It is not easy to reconstruct the apparatus escribe in ef. 6, but it is an inication that a fractance evice might soon be commercially available. Hence, it is important to generalize the filter esign theory to the fractional-orer omain. This work is a contribution in this irection. Fractional calculus is the fiel of mathematics which is concerne with the investigation an application of erivatives an integrals of arbitrary (real or complex) orer. 7 The iemann Liouville efinition of a fractional erivative of orer α is given by m t f(τ) D α Γ(m α) t f(t) := m α+ m τ m <α<m, (t τ) () m f(t) α = m. tm A more physical interpretation of a fractional erivative is given by the Grünwal Letnikov approximation b D α f(t) ( t) α m j= Γ(j α) f((m j) t), () Γ( α)γ(j +) where t is the integration step. Applying the Laplace transform is wiely use to escribe electronic circuits in the complex frequency s-omain. Hence, applying the Laplace transform to Eq. (), assuming zero initial conitions, yiels 7 L { α t f(t)} = sα F (s), (3) where α t f(t) = f (t)/t with zero initial conitions. In this work, we explore the characteristics of a filter which inclues a single fractance evice. We erive expressions for the filter center frequency an quality factor an also for the halfpower an right-phase frequencies. PSpice simulations of three passive fractionalorer filters are shown using higher-orer emulation trees. 3 Finally, impeance an frequency scaling are iscusse in the case of fractional-orer filters. It is worth noting that a fractional-orer Wien brige oscillator, which inclues two equalvalue equal-orer fractional capacitors, was stuie in ef.. a The use of a high integer-orer transfer function to emulate a fractional-orer transfer function whose orer is less than was also explaine in hap. 3 of ef. 5. This approximation is base on a Boe-plot approximation but oes not imply equivalence in the state space. b Numerical simulations in this paper are carrie out using a backwar ifference metho base on the Grünwal Letnikov approximation.
3 First-Orer Filters Generalize to the Fractional Domain 57. Single Fractional Element Filters The general transfer function of a filter with one fractional element is T (s) = bsβ + s α + a. (4) From the stability point of view, this system is stable if an only if a>anα< while it will oscillate if an only if a>anα = ; otherwise it is unstable. The location of the poles is important to etermine the filter center frequency o an its quality factor Q. From Eq. (4) it is seen that the poles in the s-plane are locate at s p = a /α e ±j(m+)π/α,m I +. The possible range of the angle in physical s-plane is θ s <πan hence there are no poles in the physical s-plane for α<. For <α<, there are only two poles locate at s, = a /α e ±j(π/α). omparing with a classical secon-orer system whose poles are locate s, = ( o /Q) ± j o (/4Q )= o e ±jδ where δ =cos ( /Q), it can be seen that o an Q are given, respectively, by o = a /α, Q = cos(π/α). (5) From the above equation, Q is negative for α<, an positive for α. At α = (classical first-orer filter), it is seen that Q =.5 as expecte. It is important now to efine the following critical frequencies () m is the frequency at which the magnitue response has a maximum or a minimum an is obtaine by solving the equation (/) T (j) =m =. () h is the half-power frequency at which the power rops to half the passban power, i.e., T (j h )=(/ ) T (j passban ). (3) rp is the right-phase frequency at which the phase T (j rp )=± π/... Fractional-orer low-pass filter onsier the fractional low-pass filter (FLPF) whose transfer function is T FLPF (s) = s α + a. (6) The magnitue an phase of this transfer function are T FLPF (j) = α +a α cos(απ/) + a, (7) T FLPF (j)= tan α sin(απ/) α cos(απ/) + a. The important critical frequencies for this FLPF are foun as m = o ( cos(απ/)) /α, rp = o /( cos(απ/)) /α,an h = o ( +cos (απ/) cos(απ/)) /α, where o is as given by Eq. (5). From these expressions it is seen that both m an rp exist only if α>, in agreement with the quality-factor expression (5). Table summarizes the magnitue an phase values at some important
4 58 A. G. awan, A. M. Soliman & A. S. Elwakil Table. Magnitue an phase values at important frequencies for the FLPF. = T FLPF (j) T FLPF (j) a o απ a cos(απ/4) 4 απ ( α)π m a sin(απ/) h tan sin(απ/) a cos(.5απ)+ p +cos (απ/) rp a cot(απ/) π frequencies. Figure (a) plots the values of m, h,an rp (normalize with respect to o ) for ifferent values of the fractional-orer α. Figure (b) is a plot of transfer function magnitue (normalize with respect to T (j) ) at ifferent critical frequencies versus α. Finally, Fig. (c) shows numerically simulate filter magnitue an phase responses for two ifferent cases α =.4 anα =.6, respectively. Note the peaking in the magnitue response at α =.6 which is classically only possible to observe in a secon-orer (α =)low-passfilter... Fractional-orer high-pass filter onsier the fractional high-pass filter (FHPF) whose transfer function is T FHPF (s) = bsα s α + a. (8) The magnitue an phase of this transfer function are T FHPF (j) = b α α +a α cos(απ/) + a, T FHPF (j)= απ + T FLPF(j). The important critical frequencies are foun as m =( a/cos(απ/)) /α, rp = ( a cos(απ/)) /α,an h = o [cos(απ/) + +cos (απ/)] /α. Figure plots the filter magnitue an phase responses for the two cases α =.4 anα =.6. (9).3. Fractional-orer ban-pass filter onsier the fractional system T FBPF (s) = bsβ s α + a, ()
5 First-Orer Filters Generalize to the Fractional Domain rp o h o m o α (a).5.5 T(j ) h T(j) T(j ) o T(j) T(j ) m T(j) T(j ) rp T(j) α (b).5 T(j).5 α =.4 Phase(T(j)) 5 5 α =.4 h 4 h 4 T(j).5.5 α =.6 Phase(T(j)) 5 5 α =.6 h 4 h 4 (c) Fig.. Numerical simulation for the FLPF representing (a) normalize critical frequencies versus α, (b) normalize magnitue response versus α, an (c) filter magnitue an phase responses at α =.4 anα =.6 assuming a = =4. whose magnitue an phase functions are T FBPF (j) = b β α +a α cos(απ/) + a, T FBPF (j)= π T FLPF(j). ()
6 6 A. G. awan, A. M. Soliman & A. S. Elwakil.5 5 T(j).5 α =.4 Phase(T(j)) 5 α =.4 4 h 4 h T(j).5.5 α =.6 Phase(T(j)) 5 5 α =.6 4 h 4 h Fig.. Magnitue an phase response of the FHPF when a =4anb =forα =.4 an α =.6, respectively. Table. Magnitue an phase values at important frequencies for the FBPF. = T FBPF (j) T FBPF (j) o ba β/α a cos(απ/4) π βπ απ 4 b (β α) (β α)π Table summarizes the magnitue an phase values at important frequencies for this system. It is seen from Table that for β<α lim T (j) =an hence the filter can be a ban-pass filter (BPF). For β = α lim T (j) = b which makes the filter a HPF. The maxima frequency m is equal to o (X) /α where X is given by X = cos(απ/)[(β α)+ α +4β(α β)tan (απ/)] (α β). ()
7 First-Orer Filters Generalize to the Fractional Domain 6 The case α =β yiels m = o. It is easy to see that there is always a maximum point in the magnitue response if α>β. Figure 3(a) shows the value of X versus α for ifferent ratios of β an α. Figure 3(b) shows the magnitue response for the filter when a = b =. Note from Fig. 3(b) that the center frequency o is not β=α X 5 4 β=.7α 3 β=.5α β=.α (a) α T(j) α = β =.3 T(j) α = β =.5.. o 4 o 4 T(j) α = β =.8 T(j) α = β =.. o 4 o 4 (b) Fig. 3. Numerical simulations for the FBPF (a) values of X versus α an (b) magnitue response for ifferent values of α an β (a = b =).
8 6 A. G. awan, A. M. Soliman & A. S. Elwakil necessarily equal to m (where the maxima occurs); which is significantly ifferent from what is known in integer-orer filters..4. Fractional-orer all-pass filter onsier the fractional all-pass filter system T APF (s) = b(sα a) s α + a. (3) The magnitue an phase of this system are, respectively, α a T FAPF (j) = b α cos(απ/) + a α +a α cos(απ/) + a, T FAPF (j)= b (α cos(απ/) a)+j α sin(απ/) ( α cos(απ/) + a)+j α sin(απ/). Table 3 summarizes the magnitue an phase values at important frequencies for this filter. It is seen here that m = rp = o an that at this frequency a minima occurs if α<anamaximaoccursifα> while the magnitue remains flat when α = (classical integer-orer all-pass filter). The half-power frequency is given by h = o [ cos(απ/) + 4cos (απ/) ] /α. Figure 4 shows the magnitue an phase responses for the two ifferent cases α =.4 anα =.6. (4) 3. PSpice Simulations Passive filters are chosen for simulations. In all simulations, we fix the fractance evice to a fractional capacitor of orer α =.4,, or.6. Figure 5(a) is the structure propose in ef. to simulate a fractional capacitor of orer.5 (Y F = F s.5 ; F = /) while the circuit in Fig. 5(b), propose in ef. 3, is use to simulate a fractional capacitor of arbitrary orer α (Y F = F s α ). To realize α =.4, for example, we nee n = 3 branches with n+ / n.5686 an n+ / n To simulate a capacitor of orer α =.6, a floating GI circuit 3 is use. The input impeance of a GI is Z i = Z Z Z 3 /Z 4 Z 5 ;takingz 3 = Z 4 =, Z = Z =/s F,anZ 4 =/s.4 F results in Z i =/s.6 F. Figures 5(c) an 5() show a passive FLPF an PSpice simulation of its magnitue response in three cases α =.4,, an.6, respectively. Figures 5(e) an 5(f) show a passive FHPF an its PSpice magnitue response while Figs. 5(g) an 5(h) Table 3. Magnitue an phase values at important frequencies for the FAPF. = T FAPF (j) T FAPF (j) b π απ π o b tan 4 b
9 First-Orer Filters Generalize to the Fractional Domain 63 T(j).5.5 α =.4 Phase(T(j)) 5 5 α =.4 4 o o 4 3 T(j).5.5 α =.6 Phase(T(j)) 5 5 α = o 4 o Fig. 4. FAPF magnitue an phase responses when α =.4 anα =.6 (a =4,b =). show a passive FAPF an its PSpice simulation results for the same three values of α. Note the peaking in the response for α> which is expecte to increase as the filter approaches α =, i.e., secon-orer filter. We have constructe the fractional low-pass filter of Fig. 5(c) in the laboratory using one of the fractional capacitors onate by the authors of ef. 6. This fractional capacitor has α.5 up to approximately 3 khz. esults are shown in Fig. 6 compare with a normal capacitor. 4. Scaling Impeance an frequency scaling can be use to ajust the filter component values or operating frequency. A fractional-orer filter is similar to an integer-orer one in terms of impeance scaling. However, for frequency scaling an assuming all critical frequencies are to be scale by a factor λ in which the new frequencies equal λ times the ol ones, then the components may be scale accoring to the following equations new = λ α or new = λ α ol. (5)
10 64 A. G. awan, A. M. Soliman & A. S. Elwakil F F p n n (a) (b) + + V i α V o (c) () α + + V i _ V o _ (e) (f) Fig. 5. (a) ealization of a fractional capacitor of orer.5, (b) realization of a fractional capacitor of orer α<, 3 (c) FLPF circuit with transfer function T (s) =/(s α + a), = a = / = 4, () PSpice simulation of the FLPF ( =6.74 kω, = F =37µF), (e) FHPF circuit with transfer function T (s) =s α /(s α + a), =, a =/ = 4, (f) PSpice simulation of the FHPF ( =6.74 kω, = F =37µF), (g) FAPF circuit with transfer function T (s) = (/)(s α a)/(s α + a), a =/ = 4, an (h) PSpice simulation of the FAPF ( =6.74 kω, = F =37µF).
11 First-Orer Filters Generalize to the Fractional Domain 65 + V i _ V o + _ α (g) (h) Fig. 5. (ontinue) NOrmal apacitance Prob apacitance Low-Pass filter Output (mv) 5 5 Frequency (Hz) Fig. 6. Experimental results of a fractional-orer low-pass filter with α.5 comparewitha normal capacitor. 5. onclusion In this work, we have generalize classical first-orer filter networks to be of fractional-orer. We have shown simulation results for filters of orer <α<
12 66 A. G. awan, A. M. Soliman & A. S. Elwakil an <α<. It is clear that more flexibility in shaping the filter response can be obtaine via a fractional-orer filter. It is also clear that the ban-pass filter, classically known to be realizable only through a secon-orer system, can actually be realize by a fractional filter of orer <α<. eferences. M. Nakagawa an K. Sorimachi, Basic characteristics of a fractance evice, IEIE Trans. Funam. Electron. ommun. omput. Sci. E75 (99) K. Saito an M. Sugi, Simulation of power-law relaxations by analog circuits: Fractal istribution of relaxation times an non-integer exponents, IEIE Trans. Funam. Electron. ommun. omput. Sci. E76 (993) M. Sugi, Y. Hirano, Y. F. Miura an K. Saito, Simulation of fractal immittance by analog circuits: An approach to the optimize circuits, IEIE Trans. Funam. Electron. ommun. omput. Sci. E8 (999) A. Abbisso,. aponetto, L. Fortuna an D. Porto, Non-integer-orer integration by using neural networks, Proc. Int. Symp. ircuits an System ISAS, Vol. 3 (), pp P. Arena,. aponetto, L. Fortuna an D. Porto, Nonlinear Noninteger Orer ircuits an Systems (Worl Scientific, ). 6. K. Biswas, S. Sen an P. Dutta, ealization of a constant phase element an its performance stuy in a ifferentiator circuit, IEEE Trans. ircuits Syst.-II 53 (6) K. B. Olham an J. Spanier, Fractional alculus (Acaemic Press, New York, 974). 8. S. G. Samko, A. A. Kilbas an O. I. Marichev, Fractional Integrals an Derivatives: Theory an Application (Goron & Breach, 987). 9. K. S. Miller an B. oss, An Introuction to the Fractional alculus an Fractional Differential Equations (John Wiley & Sons, 993).. T. T. Hartley an. F. Lorenzo, Initialization, conceptualization, an application in the generalize fractional calculus, National Aeronautics an Space Aministration (NASA/TP ) (998).. W. Ahma,. El-Khazali an A. S. Elwakil, Fractional-orer Wien-brige oscillator, Electron. Lett. 37 ().. D. Matignon, Stability results in fractional ifferential equations with applications to control processing, Proc. Multiconf. omputational Engineering in Systems an Application IMIS IEEE-SM, Vol. (996), pp A. Buak, Passive an Active Network Analysis an Synthesis (Wavelan Press Inc., Illinois, 99).
Fractional Order Oscillator with Independent Control of Phase and Frequency
Fractional Order Oscillator with Independent Control of Phase and Frequency Lobna A. Said Mathematics Dept., Faculty of IET, German University in Cairo (GUC), Egypt l.a.said@ieee.org Ahmed H. Madian Electronics
More informationChapter 6. Electromagnetic Oscillations and Alternating Current
hapter 6 Electromagnetic Oscillations an Alternating urrent hapter 6: Electromagnetic Oscillations an Alternating urrent (hapter 31, 3 in textbook) 6.1. Oscillations 6.. The Electrical Mechanical Analogy
More informationModeling time-varying storage components in PSpice
Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationGeneralized-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal
Commun. Theor. Phys. (Beijing, China) 44 (25) pp. 72 78 c International Acaemic Publishers Vol. 44, No. 1, July 15, 25 Generalize-Type Synchronization of Hyperchaotic Oscillators Using a Vector Signal
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationDesign and Application of Fault Current Limiter in Iran Power System Utility
Australian Journal of Basic an Applie Sciences, 7(): 76-8, 13 ISSN 1991-8178 Design an Application of Fault Current Limiter in Iran Power System Utility M. Najafi, M. Hoseynpoor Department of Electrical
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationOrdinary Differential Equations: Homework 1
Orinary Differential Equations: Homework 1 M. Gameiro, J.-P. Lessar, J.D. Mireles James, K. Mischaikow January 12, 2017 2 Chapter 1 Motivation 1.1 Exercises Exercise 1.1.1. (Frictionless spring) Consier
More informationUncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control
Int. J. of Computers, Communications & Control, ISSN 1841-9836, E-ISSN 1841-9844 Vol. VI (011, No. 3 (September, pp. 418-47 Uncertain Fractional Orer Chaotic Systems Tracking Design via Aaptive Hybri Fuzzy
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationInvestigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model
ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More informationA Quantitative Analysis of Coupling for a WPT System Including Dielectric/Magnetic Materials
Progress In Electromagnetics Research Letters, Vol. 72, 127 134, 2018 A Quantitative Analysis of Coupling for a WPT System Incluing Dielectric/Magnetic Materials Yangjun Zhang *, Tatsuya Yoshiawa, an Taahiro
More informationAn Interdisciplinary Topic
: An Interdisciplinary Topic Ahmed Elwakil T S P 2 10 00 1 6 A H M E D E L W A K I L 1 Outline Fractional Calculus History and Definitions Fractional Trigonometric Identities s-domain, Stability and Impulse
More informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationAdjustable Fractional-Delay Filters Utilizing the Farrow Structure and Multirate Techniques
Ajustable Fractional-Delay Filters Utilizing the Farrow Structure an Multirate Techniques Håkan Johansson (hakanj@isy.liu.se)* an Ewa Hermanowicz (hewa@eti.pg.ga.pl)** *Division of Electronics Systems,
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationFractional Geometric Calculus: Toward A Unified Mathematical Language for Physics and Engineering
Fractional Geometric Calculus: Towar A Unifie Mathematical Language for Physics an Engineering Xiong Wang Center of Chaos an Complex Network, Department of Electronic Engineering, City University of Hong
More informationCompact MOS-RC Voltage-Mode Fractional- Order Oscillator Design
dspace.vutbr.cz Compact MOS- Voltage-Mode Fractional- Order Oscillator Design KARTCI, A.; HERENCSÁR, N.; KOTON, J.; PSYCHALINOS, C. Proceedings of the 07 3 European Conference on Circuit Theory and Design
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More information6.003 Homework #7 Solutions
6.003 Homework #7 Solutions Problems. Secon-orer systems The impulse response of a secon-orer CT system has the form h(t) = e σt cos(ω t + φ)u(t) where the parameters σ, ω, an φ are relate to the parameters
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationA LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM
Teor Imov r. ta Matem. Statist. Theor. Probability an Math. Statist. Vip. 81, 1 No. 81, 1, Pages 147 158 S 94-911)816- Article electronically publishe on January, 11 UDC 519.1 A LIMIT THEOREM FOR RANDOM
More informationChapter Primer on Differentiation
Capter 0.01 Primer on Differentiation After reaing tis capter, you soul be able to: 1. unerstan te basics of ifferentiation,. relate te slopes of te secant line an tangent line to te erivative of a function,.
More informationx = c of N if the limit of f (x) = L and the right-handed limit lim f ( x)
Limit We say the limit of f () as approaches c equals L an write, lim L. One-Sie Limits (Left an Right-Hane Limits) Suppose a function f is efine near but not necessarily at We say that f has a left-hane
More informationDesigning Information Devices and Systems II Spring 2019 A. Sahai, J. Roychowdhury, K. Pister Midterm 1: Practice
EES 16B Designing Information Devices an Systems II Spring 019 A. Sahai, J. Roychowhury, K. Pister Miterm 1: Practice 1. Speaker System Your job is to construct a speaker system that operates in the range
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationNonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain
Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationUNIT 4:Capacitors and Dielectric
UNIT 4:apacitors an Dielectric SF7 4. apacitor A capacitor is a evice that is capable of storing electric charges or electric potential energy. It is consist of two conucting plates separate by a small
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationSimilar Operators and a Functional Calculus for the First-Order Linear Differential Operator
Avances in Applie Mathematics, 9 47 999 Article ID aama.998.067, available online at http: www.iealibrary.com on Similar Operators an a Functional Calculus for the First-Orer Linear Differential Operator
More informationTime Headway Requirements for String Stability of Homogeneous Linear Unidirectionally Connected Systems
Joint 48th IEEE Conference on Decision an Control an 8th Chinese Control Conference Shanghai, PR China, December 6-8, 009 WeBIn53 Time Heaway Requirements for String Stability of Homogeneous Linear Uniirectionally
More informationLecture 6: Control of Three-Phase Inverters
Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationCouncil for Innovative Research
ISSN: 347-3487 A solution of Fractional Laplace's equation y Moifie separation of variales ABSTRACT Amir Pishkoo,, Maslina Darus, Fatemeh Tamizi 3 Physics an Accelerators Research School (NSTRI) P.O. Box
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationConsider for simplicity a 3rd-order IIR filter with a transfer function. where
Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationChapter 2 Derivatives
Chapter Derivatives Section. An Intuitive Introuction to Derivatives Consier a function: Slope function: Derivative, f ' For each, the slope of f is the height of f ' Where f has a horizontal tangent line,
More informationFrIP6.1. subject to: V 0 &V 1 The trajectory to the target (x(t)) is generated using the HPFbased, gradient dynamical system:
Proceeings of the 45th IEEE onference on Decision & ontrol Manchester Gran Hyatt Hotel San Diego, A, USA, December 3-5, 006 Agile, Steay Response of Inertial, onstraine Holonomic Robots Using Nonlinear,
More informationBasic IIR Digital Filter Structures
Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient
More informationMath 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v
Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through
More information. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp
. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an ID Control Theory Chapter 11: Feeback an ID Control Theory I. Introuction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works
More informationPhysics 4B. Chapter 31: Questions: 2, 8, 12 Exercises & Problems: 2, 23, 24, 32, 41, 44, 48, 60, 72, 83. n n f
Physics 4B Solutions to hapter 1 HW hapter 1: Questions:, 8, 1 Exercises & Probles:,, 4,, 41, 44, 48, 60, 7, 8 Question 1- (a) less; (b) greater Question 1-8 (a) 1 an 4; (b) an Question 1-1 (a) lea; (b)
More informationThe Sokhotski-Plemelj Formula
hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationDerivative of a Constant Multiple of a Function Theorem: If f is a differentiable function and if c is a constant, then
Bob Brown Math 51 Calculus 1 Chapter 3, Section Complete 1 Review of the Limit Definition of the Derivative Write the it efinition of the erivative function: f () Derivative of a Constant Multiple of a
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationExam 2 Answers Math , Fall log x dx = x log x x + C. log u du = 1 3
Exam Answers Math -, Fall 7. Show, using any metho you like, that log x = x log x x + C. Answer: (x log x x+c) = x x + log x + = log x. Thus log x = x log x x+c.. Compute these. Remember to put boxes aroun
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationAnalytical accuracy of the one dimensional heat transfer in geometry with logarithmic various surfaces
Cent. Eur. J. Eng. 4(4) 014 341-351 DOI: 10.478/s13531-013-0176-8 Central European Journal of Engineering Analytical accuracy of the one imensional heat transfer in geometry with logarithmic various surfaces
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationTIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS
TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical
More information3.2 Differentiability
Section 3 Differentiability 09 3 Differentiability What you will learn about How f (a) Might Fail to Eist Differentiability Implies Local Linearity Numerical Derivatives on a Calculator Differentiability
More informationModule 2. DC Circuit. Version 2 EE IIT, Kharagpur
Moule 2 DC Circuit Lesson 9 Analysis of c resistive network in presence of one non-linear element Objectives To unerstan the volt (V ) ampere ( A ) characteristics of linear an nonlinear elements. Concept
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationx f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.
Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?
More informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More information