Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs.
|
|
- Aileen Smith
- 6 years ago
- Views:
Transcription
1 Handou # 6 (MEEN 67) Numercal Inegraon o Fnd Tme Response of SDOF mechancal sysem Sae Space Mehod The EOM for a lnear sysem s M X DX K X F() () X X X X V wh nal condons, a 0 0 ; 0 Defne he followng varables, and wre EOM () as wo frs-order Eqs. o o o Y X; Y X () M Y DY KY F() & Y Y (3) whch can be wren n marx form as Y 0 Y 0 Y M K M D Y M F (4) Or, Y AY b (5) Y X 0 0 wh Y ; ; Y X A M K M D b M F Ths s known as he sae-space formulaon. Eq. (5) s o be negraed numercally wh nal condon vecor Y X V T o o o MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03
2 If he appled load s NOT a funcon of me, hen an equlbrum sae s defned afer a very long me as Y 0 YE A bo (6) Compuaonal sofware such as Mahcad, Mapple, Mahemaca, Malab, ec has bul-n funcons or commands o perform he numercal negraon of equaons se n he form Y AYb, even when sysem s nonlnear,.e. A=A(Y). A few words abou numercal negraon mehods Typcal numercal negraon mehods nclude a) Euler (smple) mehod b) Fourh and Ffh-Order Runge-Kua negraors, c) Rosenbrock Mehod see references on page, d) Adams Predcor Correcor Mehods e) Average Acceleraon and Wlson-θ (Implc) Mehods In mos mehods, he selecon of an adequae me sep s crucal for numercally sable and accurae resuls. (a)-(b) are favored by he young naes no numercal compung and because of her ready avalably n modern compuaonal sofware. (c) (d) are more modern (mplc) mehods wh auomaed nermedae reszng of he me sep whle performng he negraon. Mehods (e) have long been favored by srucural mechanc analyss when negrang Mulple DOF (lnear) sysems All mehods suffer from defcences when nonlneares are apparen hus forcng exremely small me seps and he ensung MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03
3 cos wh los of numercal compung (me). (Memory) Sorage appears no o be an ssue anymore. Sae-space mehod for MDOF sysems. Recall he EOMS for a lnear sysem are where MU+DU+KU () =F() (7) U,U, and U are he vecors of generalzed dsplacemen, F s he vecor of velocy and acceleraon, respecvely; and () generalzed (exernal forces) acng on he sysem. M,D,K represen he marces of nera, vscous dampng and sffness coeffcens, respecvely. Defne he followng varables, Y=U ; Y=U and wre EOM (7) as a se of n-frs-order Eqs. () (8) MY +DY +KY =F & Y =Y (9) whch can be wren n marx form as Y 0 I Y 0 = Y -M K -M DY -M F (0) Or, Y AY b () U 0 I 0 wh Y ; A ; b U -M K -M D M F () The marces are square wh n-rows = n columns, whle he vecors are n- rows. MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 3
4 and nal condons Y( 0) Uo U o. A s a n x n marx. I s he nxn deny marx, and 0 s a nxn marx full of zeroes. Condons for a good numercal negraor In general a numercal negraon scheme should a) reproduce EOM as me sep 0 b) provde, as wh physcal model, bounded soluons for any sze of me sep,.e. mehod should be sable c) reproduce he physcal response wh fdely and accuracy. The numercal negraon reles n represenng me dervaves of a funcon wh an algebrac approxmaon, for example T dx x x lm d 0 x x x x x ~ Eq. above s exac only f 0 (3) Numercal negraon mehods are usually dvded no wo caegores, mplc and explc. Consder he ODE x fx, (4) In an explc numercal scheme, he ODE s represened n erms of known values a a pror me sep,.e. x x f, (5) x, MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 4
5 whle n an mplc numercal scheme x x f (6) x, Explc numercal schemes are condonally sable. Tha s, hey provde bounded numercal soluons for (very) small me seps. For example, T where Tn and n n K M are he naural perod and naural frequency of he sysem, respecvely. The resrcon on he me sep s oo severe when analyzng sff sysems,.e. hose wh large naural frequences. n cr (7) Implc numercal schemes are uncondonally sable,.e. do no mpose a resrcon on he sze of he me sep. (However, accuracy may be compromsed f oo large me seps are used). ANALOGY beween numercal schemes and a fler A few words of wsdom released n class MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 5
6 The average acceleraon mehod for numercal negraon of SDOF equaon: M X DX K X F() () Consder a change of hnkng frame by defnng Arhmec ~ Connuum (funcon) X X, F F,, X X F F (7) Wre Eq. () a wo mes, M X DX K X F (8a) a M X DX K X F (8b) a Subrac (b) from (a) o oban: where M X DX KX F (9) X X X FF F, X X X, ec Noe ha known quanes a = are X, X, X. (0) Ths numercal mehod s exremely popular among he srucural dynamcs communy. Is exenson o MDOF sysems wll be shown laer. The oher favore mehod, he Wlson-θ scheme, wll also be gven n laer lecures (MDOF sysems). MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 6
7 Now, assume he acceleraon s consan whn he me nerval for 0 () X a se as an average value a X X. The velocy and dsplacemen follow from negraon of Eq. (6) whn he me nerval,.e. X X a (a) X X X a (b) Acceleraon velocy X X X a X X X X a + + dsplacemen X X X X a X + MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 7
8 A he end of he me nerval, he velocy and dsplacemen equal X X X a (3a) X X X X a (3b) And he dfferences n velocy and dsplacemen re X X X X X X X X a X X X X X X a X X X X X 4 X X X 4 (4a) (5b) 4 (6a) from (5b), X X X X and no (4a) X X X 4 X X X X MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 8
9 X X X (6b) Noe ha n Eqs (6), X, X, depend on he known values obaned a he pror me sep,.e. X, X and he unknown X. Thus, replace X, X no he dfference equaon (9), M X DX KX F 4 M X X X D X X KX F Rearrangng erms leads o where K X F (7) * * * 4 K K D M (8a) * 4 F FM X D M X (8b) are known as pseudo dynamc sffness and dynamc force, respecvely The recpe for he numercal negraon o fnd he sysem me response s MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 9
10 A me, known varables are X, X (curren sae) X M F DX K X () fnd from EOM: () form pseudo sffness and forcng funcons, K *, F Eqs. (8), (3) Calculae * X K F, and X X X ; (4) X XX, X X X a + (5) Increase me o + and reurn o sep () from The average acceleraon mehod s an mplc mehod,.e. numercally sable and conssen. The dsadvanage s ha requres memory 3 o sore X, X, X, X, F. 3 A non-ssue n he s cenury MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 0
11 Average acceleraon mehods for numercal negraon of a nonlnear sysem Consder he sysem wh EOM where gx, X 3 o o 3 M X g X, X F( ) (30) s a nonlnear funcon, for example g X, X g k X k X F sgn( X ) As wh he lnear sysem, evaluae Eq. (30) a wo mes (closely spaced): a M X g X, X F (3a) a M X g X, X F (3b) Subrac (b) from (a) o oban: where (3) M X g g F X X X F F F,,,, g g X X g g X X (33) A Taylor seres expanson of he nonlnear funcon a gves g g g g X X O X, X X X X, X 0 X (34) defne local lnearzed sffness and dampng coeffcens as MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03
12 Hence, K g X ; D g X, X X X, X g g K X D X (35) and he dfference Eq. (3) becomes lnear M X D X K X F (36) Eq. (35) s formally dencal o he one devsed for a lnear sysem. Thus, he numercal reamen s smlar, excep ha a each me sep, lnearzed sffness and dampng coeffcens need be calculaed. The recpe s hus dencal; however wh he apparen nonlneary, he mehod does no guaranee sably for (oo) large me seps. The recpe for he numercal negraon o fnd he sysem me response s A me, known varables are X, X (curren sae) () fnd from EOM: X M F g X, X (a) fnd local (lnearzed) sffness and dampng coeffcens, (K, D ) from eq. (35) (b) form pseudo sffness and forcng funcons, K *, F from MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03
13 * 4 K K D M ; * 4 F F M X D M X (3) Calculae * X K F, and X X X ; (4) Se X XX, X X X a + (5) Increase me o + and reurn o sep () Ofen regarded as an ar, numercal compung s n acualy a well esablshed branch of mahemacs. References The followng are a mus! Press, W.H., Flannery, B.P., Teukolsky, S.A., and Veerlng, W.T., 986 ( s edon), Numercal Recpes, The Ar of Scenfc Compung, Cambrdge Unversy Press, UK. Bahe, J-K, 98 ( s ed), 007 or laes, Fne Elemen Procedures, Prence Hall. Pche, R., and P. Nevalanen, 999, Varable Sep Rosenbrock Algorhm for Transen Response of Damped Srucures, Proc. IMechE, Vol. 3, par C, Paper C Rosenbrock, H.H., Some General Implc Processes for he Numercal Soluon of Dfferenal Equaons, MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 3
14 Appendx A 4 Numercal Inegraon Usng Modfed Euler s Mehod I s ofen dffcul o solve (exacly) a nonlnear dfferenal equaon. Numercal negraon s hen employed o oban resuls (predcons of moon). The Modfed Euler Mehod s one ype of numercal negraon scheme. Solve for q() governed by M qcq KqQ( q, q, ) (A.) wh nal condons q0and q 0a 0. In Eqn. (A.), Qq (, q, ) may conan erms ha are nonlnear n Le V dq d (A.) and wre Eq. (A.) as a sysem of TWO frs order dfferenal equaons,.e. C K Q( q, V, ) V V q, M M M q V (A.3) Defne K C n and (A.4) M K M 4 Ths Appendx s ncluded because mos young engneerng sudens have learned only abou Euler s mehod. Hence, he appendx complemens her educaon.. However, I encourage you o abandon he usage of hs poor mehod. MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 4
15 as he naural frequency and vscous dampng rao, respecvely. Wh he noed defnons Eqs. (A.3) become dv Q( q, V, ) V F ( q, V, ) nv n q d M dq q F ( q, V, ) V d (A.5) Noe ha F and F are he slopes of he (V vs ) and (q vs ) curves, respecvely. q() V() q+ V q V+ + + Le he numercal approxmaons (arhmec values) be V V,and q q (A.6) ; and s a suably small me sep for numercal negraon. where 0,,... F q+ In Euler s numercal scheme, approxmae he me dervaves (or slopes) as: q V F V+ V V f, f F q, V,, q q f, f F q, V, for =,.. (A.7) MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 5
16 wh nal condons q0and q 0a 0 0 Eq. (7) offer a recursve relaon o calculae he numercal (arhmec) values of he varables V and q a successve mes ). The regular Euler mehod s frs order wh a runcaon error of order ( ). A modfed Euler mehod (second order accurae) wh error order ( 3 ) follows: (a) Compue prelmnary esmaes of V, q as V V F q, V,, q q F q, V, (A8.a) (b) Use hese prelmnary esmaes o oban mproved slopes as f F q, V,, f F q, V,, (A8.b) (c) Defne average slopes as per f f f f f f (A8.c) (c) and oban new esmaons usng he averaged slopes,.e. MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 6
17 V V f q q f, (A8.d) (d) Repea seps (a) hrough (d) a each me sep, ; 0,,..., and sarng wh he nal condons q and q a The sze of he me sep s very mporan o oban accurae and numercally sable resuls. If s oo large, hen numercal predcons wll be n error and very lkely show unsable (oscllang) resuls. If s oo small, hen he numercal mehod wll be slow and cosly snce he number of operaons ncreases accordngly. In pracce, a me sep of he order =T n /60, where T n s he naural perod of moon gven as (/ω n ). Euler s mehod s mos mes NOT a good choce o perform he numercal negraon of lnear or nonlnear ODES. Alas, s wdely used by rooke engneers and engneerng sudens because s easy o mplemen. Ofen enough, however, predcons can be wrong and msleadng. I call Euler s mehod a BRUTE FORCE approach, MEEN 67 HD 6 Numercal Inegraon for Tme Response: SDOF sysem L. San Andrés 03 7
Handout # 13 (MEEN 617) Numerical Integration to Find Time Response of MDOF mechanical system. The EOMS for a linear mechanical system are
Handou # 3 (MEEN 67) Numercal Inegraon o Fnd Tme Response of MDOF mechancal sysem The EOMS for a lnear mechancal sysem are MU+DU+KU =F () () () where U,U, and U are he vecors of generalzed dsplacemen,
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at
More informationFTCS Solution to the Heat Equation
FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence
More informationChapter 2 Linear dynamic analysis of a structural system
Chaper Lnear dynamc analyss of a srucural sysem. Dynamc equlbrum he dynamc equlbrum analyss of a srucure s he mos general case ha can be suded as akes no accoun all he forces acng on. When he exernal loads
More informationIncluding the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.
Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample
More information[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5
TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres
More information3. OVERVIEW OF NUMERICAL METHODS
3 OVERVIEW OF NUMERICAL METHODS 3 Inroducory remarks Ths chaper summarzes hose numercal echnques whose knowledge s ndspensable for he undersandng of he dfferen dscree elemen mehods: he Newon-Raphson-mehod,
More informationCH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal
More information( ) () we define the interaction representation by the unitary transformation () = ()
Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger
More informationP R = P 0. The system is shown on the next figure:
TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems
Swss Federal Insue of Page 1 The Fne Elemen Mehod for he Analyss of Non-Lnear and Dynamc Sysems Prof. Dr. Mchael Havbro Faber Dr. Nebojsa Mojslovc Swss Federal Insue of ETH Zurch, Swzerland Mehod of Fne
More informationJohn Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany
Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy
More informationCONSISTENT EARTHQUAKE ACCELERATION AND DISPLACEMENT RECORDS
APPENDX J CONSSTENT EARTHQUAKE ACCEERATON AND DSPACEMENT RECORDS Earhqake Acceleraons can be Measred. However, Srcres are Sbjeced o Earhqake Dsplacemens J. NTRODUCTON { XE "Acceleraon Records" }A he presen
More information2/20/2013. EE 101 Midterm 2 Review
//3 EE Mderm eew //3 Volage-mplfer Model The npu ressance s he equalen ressance see when lookng no he npu ermnals of he amplfer. o s he oupu ressance. I causes he oupu olage o decrease as he load ressance
More informationNotes on the stability of dynamic systems and the use of Eigen Values.
Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon
More informationHEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD
Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,
More informationJ i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.
umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal
More informationOn One Analytic Method of. Constructing Program Controls
Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna
More informationVolatility Interpolation
Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local
More informationApproximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy
Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae
More informationSolution in semi infinite diffusion couples (error function analysis)
Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of
More informationLet s treat the problem of the response of a system to an applied external force. Again,
Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem
More informationGENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim
Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran
More informationLecture 18: The Laplace Transform (See Sections and 14.7 in Boas)
Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on
More informationResponse of MDOF systems
Response of MDOF syses Degree of freedo DOF: he nu nuber of ndependen coordnaes requred o deerne copleely he posons of all pars of a syse a any nsan of e. wo DOF syses hree DOF syses he noral ode analyss
More informationOrdinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s
Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class
More informationCubic Bezier Homotopy Function for Solving Exponential Equations
Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.
More information. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.
Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons
More informationOn computing differential transform of nonlinear non-autonomous functions and its applications
On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,
More informationRobust and Accurate Cancer Classification with Gene Expression Profiling
Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem
More informationChapter Lagrangian Interpolation
Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and
More informationFall 2010 Graduate Course on Dynamic Learning
Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure
More information. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.
Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as
More informationPhysical Simulation Using FEM, Modal Analysis and the Dynamic Equilibrium Equation
Physcal Smulaon Usng FEM, Modal Analyss and he Dynamc Equlbrum Equaon Paríca C. T. Gonçalves, Raquel R. Pnho, João Manuel R. S. Tavares Opcs and Expermenal Mechancs Laboraory - LOME, Mechancal Engneerng
More informationOutline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model
Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon
More informationImplementation of Quantized State Systems in MATLAB/Simulink
SNE T ECHNICAL N OTE Implemenaon of Quanzed Sae Sysems n MATLAB/Smulnk Parck Grabher, Mahas Rößler 2*, Bernhard Henzl 3 Ins. of Analyss and Scenfc Compung, Venna Unversy of Technology, Wedner Haupsraße
More informationA NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION
S19 A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION by Xaojun YANG a,b, Yugu YANG a*, Carlo CATTANI c, and Mngzheng ZHU b a Sae Key Laboraory for Geomechancs and Deep Underground Engneerng, Chna Unversy
More informationPHYS 705: Classical Mechanics. Canonical Transformation
PHYS 705: Classcal Mechancs Canoncal Transformaon Canoncal Varables and Hamlonan Formalsm As we have seen, n he Hamlonan Formulaon of Mechancs,, are ndeenden varables n hase sace on eual foong The Hamlon
More informationFI 3103 Quantum Physics
/9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon
More informationOnline Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading
Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng
More informationPendulum Dynamics. = Ft tangential direction (2) radial direction (1)
Pendulum Dynams Consder a smple pendulum wh a massless arm of lengh L and a pon mass, m, a he end of he arm. Assumng ha he fron n he sysem s proporonal o he negave of he angenal veloy, Newon s seond law
More informationComb Filters. Comb Filters
The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of
More informationMEEN Handout 4a ELEMENTS OF ANALYTICAL MECHANICS
MEEN 67 - Handou 4a ELEMENTS OF ANALYTICAL MECHANICS Newon's laws (Euler's fundamenal prncples of moon) are formulaed for a sngle parcle and easly exended o sysems of parcles and rgd bodes. In descrbng
More informationIn the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!
ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal
More informationNew M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)
Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor
More informationLecture 6: Learning for Control (Generalised Linear Regression)
Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson
More informationVEHICLE DYNAMIC MODELING & SIMULATION: COMPARING A FINITE- ELEMENT SOLUTION TO A MULTI-BODY DYNAMIC SOLUTION
21 NDIA GROUND VEHICLE SYSTEMS ENGINEERING AND TECHNOLOGY SYMPOSIUM MODELING & SIMULATION, TESTING AND VALIDATION (MSTV) MINI-SYMPOSIUM AUGUST 17-19 DEARBORN, MICHIGAN VEHICLE DYNAMIC MODELING & SIMULATION:
More informationUNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION
INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he
More informationM. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria
IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund
More informationChapter 6: AC Circuits
Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.
More informationHow about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?
lmcd Lnear ransformaon of a vecor he deas presened here are que general hey go beyond he radonal mar-vecor ype seen n lnear algebra Furhermore, hey do no deal wh bass and are equally vald for any se of
More informationChapters 2 Kinematics. Position, Distance, Displacement
Chapers Knemacs Poson, Dsance, Dsplacemen Mechancs: Knemacs and Dynamcs. Knemacs deals wh moon, bu s no concerned wh he cause o moon. Dynamcs deals wh he relaonshp beween orce and moon. The word dsplacemen
More informationWiH Wei He
Sysem Idenfcaon of onlnear Sae-Space Space Baery odels WH We He wehe@calce.umd.edu Advsor: Dr. Chaochao Chen Deparmen of echancal Engneerng Unversy of aryland, College Par 1 Unversy of aryland Bacground
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours
NATONAL UNVERSTY OF SNGAPORE PC5 ADVANCED STATSTCAL MECHANCS (Semeser : AY 1-13) Tme Allowed: Hours NSTRUCTONS TO CANDDATES 1. Ths examnaon paper conans 5 quesons and comprses 4 prned pages.. Answer all
More informationDynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005
Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc
More informationEcon107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)
Econ7 Appled Economercs Topc 5: Specfcaon: Choosng Independen Varables (Sudenmund, Chaper 6 Specfcaon errors ha we wll deal wh: wrong ndependen varable; wrong funconal form. Ths lecure deals wh wrong ndependen
More informationBifurcation & Chaos in Structural Dynamics: A Comparison of the Hilber-Hughes-Taylor-α and the Wang-Atluri (WA) Algorithms
Bfurcaon & Chaos n Srucural Dynamcs: A Comparson of he Hlber-Hughes-Taylor-α and he Wang-Alur (WA) Algorhms Absrac Xuechuan Wang*, Wecheng Pe**, and Saya N. Alur* Texas Tech Unversy, Lubbock, TX, 7945
More information(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function
MACROECONOMIC THEORY T J KEHOE ECON 87 SPRING 5 PROBLEM SET # Conder an overlappng generaon economy le ha n queon 5 on problem e n whch conumer lve for perod The uly funcon of he conumer born n perod,
More informationTHE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS
THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he
More informationExistence and Uniqueness Results for Random Impulsive Integro-Differential Equation
Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal
More information10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :
. A. IUITS Synopss : GOWTH OF UNT IN IUIT : d. When swch S s closed a =; = d. A me, curren = e 3. The consan / has dmensons of me and s called he nducve me consan ( τ ) of he crcu. 4. = τ; =.63, n one
More informationV.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS
R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon
More informationPart II CONTINUOUS TIME STOCHASTIC PROCESSES
Par II CONTINUOUS TIME STOCHASTIC PROCESSES 4 Chaper 4 For an advanced analyss of he properes of he Wener process, see: Revus D and Yor M: Connuous marngales and Brownan Moon Karazas I and Shreve S E:
More informationRobustness Experiments with Two Variance Components
Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference
More informationMohammad H. Al-Towaiq a & Hasan K. Al-Bzoor a a Department of Mathematics and Statistics, Jordan University of
Ths arcle was downloaded by: [Jordan Unv. of Scence & Tech] On: 05 Aprl 05, A: 0:4 Publsher: Taylor & Francs Informa Ld Regsered n England and ales Regsered umber: 07954 Regsered offce: Mormer House, 37-4
More informationComparison of Differences between Power Means 1
In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,
More informationRelative controllability of nonlinear systems with delays in control
Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.
More informationCHAPTER 10: LINEAR DISCRIMINATION
CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g
More informationAdvanced time-series analysis (University of Lund, Economic History Department)
Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng
More informationLecture 9: Dynamic Properties
Shor Course on Molecular Dynamcs Smulaon Lecure 9: Dynamc Properes Professor A. Marn Purdue Unversy Hgh Level Course Oulne 1. MD Bascs. Poenal Energy Funcons 3. Inegraon Algorhms 4. Temperaure Conrol 5.
More informationOn the numerical treatment ofthenonlinear partial differentialequation of fractional order
IOSR Journal of Mahemacs (IOSR-JM) e-iss: 2278-5728, p-iss: 239-765X. Volume 2, Issue 6 Ver. I (ov. - Dec.26), PP 28-37 www.osrjournals.org On he numercal reamen ofhenonlnear paral dfferenalequaon of fraconal
More informationTransient Response in Electric Circuits
Transen esponse n Elecrc rcus The elemen equaon for he branch of he fgure when he source s gven by a generc funcon of me, s v () r d r ds = r Mrs d d r (')d' () V The crcu s descrbed by he opology equaons
More information( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model
BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng
More informationDEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL
DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA
More informationLi An-Ping. Beijing , P.R.China
A New Type of Cpher: DICING_csb L An-Png Bejng 100085, P.R.Chna apl0001@sna.com Absrac: In hs paper, we wll propose a new ype of cpher named DICING_csb, whch s derved from our prevous sream cpher DICING.
More informationNovel Technique for Dynamic Analysis of Shear-Frames Based on Energy Balance Equations
Novel Technque for Dynamc Analyss of Shear-Frames Based on Energy Balance Equaons Mohammad Jall Sadr Abad, Mussa Mahmoud,* and Earl Dowell 3. Ph.D. Suden, Deparmen of Cvl Engneerng, Shahd Rajaee Teacher
More informationA Deterministic Algorithm for Summarizing Asynchronous Streams over a Sliding Window
A Deermnsc Algorhm for Summarzng Asynchronous Sreams over a Sldng ndow Cosas Busch Rensselaer Polyechnc Insue Srkana Trhapura Iowa Sae Unversy Oulne of Talk Inroducon Algorhm Analyss Tme C Daa sream: 3
More informationGear System Time-varying Reliability Analysis Based on Elastomer Dynamics
A publcaon of CHEMICAL ENGINEERING TRANSACTIONS VOL. 33, 013 Gues Edors: Enrco Zo, Pero Barald Copyrgh 013, AIDIC Servz S.r.l., ISBN 978-88-95608-4-; ISSN 1974-9791 The Ialan Assocaon of Chemcal Engneerng
More informationF-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction
ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or
More informationDensity Matrix Description of NMR BCMB/CHEM 8190
Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon Alernae approach o second order specra: ask abou x magnezaon nsead of energes and ranson probables. If we say wh one bass se, properes vary
More informationDepartment of Economics University of Toronto
Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of
More informationGraduate Macroeconomics 2 Problem set 5. - Solutions
Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K
More informationEP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES
EP Queung heory and eleraffc sysems 3rd lecure Marov chans Brh-deah rocess - Posson rocess Vora Fodor KTH EES Oulne for oday Marov rocesses Connuous-me Marov-chans Grah and marx reresenaon Transen and
More informationControl Systems. Mathematical Modeling of Control Systems.
Conrol Syem Mahemacal Modelng of Conrol Syem chbum@eoulech.ac.kr Oulne Mahemacal model and model ype. Tranfer funcon model Syem pole and zero Chbum Lee -Seoulech Conrol Syem Mahemacal Model Model are key
More informationHidden Markov Models Following a lecture by Andrew W. Moore Carnegie Mellon University
Hdden Markov Models Followng a lecure by Andrew W. Moore Carnege Mellon Unversy www.cs.cmu.edu/~awm/uorals A Markov Sysem Has N saes, called s, s 2.. s N s 2 There are dscree meseps, 0,, s s 3 N 3 0 Hdden
More informationAvailable online at ScienceDirect. Procedia Engineering 153 (2016 )
Avalable onlne a www.scencedrec.com ScenceDrec Proceda Engneerng 5 (6 ) 747 75 XXV Polsh Russan Slova Semnar Theorecal Foundaon of Cvl Engneerng Vbraon of embedded foundaon a mul-layered base ang no accoun
More informationEpistemic Game Theory: Online Appendix
Epsemc Game Theory: Onlne Appendx Edde Dekel Lucano Pomao Marcano Snscalch July 18, 2014 Prelmnares Fx a fne ype srucure T I, S, T, β I and a probably µ S T. Le T µ I, S, T µ, βµ I be a ype srucure ha
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ
More informationPerformance Analysis for a Network having Standby Redundant Unit with Waiting in Repair
TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen
More informationIMPACT MODELING BY MANIFOLD APPROACH IN EXPLICIT TRANSIENT DYNAMICS
COMPDYN III ECCOMS Themac Conference on Compuaonal Mehods n Srucural Dynamcs and Earhquake Engneerng M. Papadrakaks, M. Fragadaks, V. Plevrs (eds.) Corfu, Greece, 5 8 May IMPCT MODELING Y MNIFOLD PPROCH
More informationBayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance
INF 43 3.. Repeon Anne Solberg (anne@f.uo.no Bayes rule for a classfcaon problem Suppose we have J, =,...J classes. s he class label for a pxel, and x s he observed feaure vecor. We can use Bayes rule
More informationScattering at an Interface: Oblique Incidence
Course Insrucor Dr. Raymond C. Rumpf Offce: A 337 Phone: (915) 747 6958 E Mal: rcrumpf@uep.edu EE 4347 Appled Elecromagnecs Topc 3g Scaerng a an Inerface: Oblque Incdence Scaerng These Oblque noes may
More informationDensity Matrix Description of NMR BCMB/CHEM 8190
Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon If we say wh one bass se, properes vary only because of changes n he coeffcens weghng each bass se funcon x = h< Ix > - hs s how we calculae
More informationCS286.2 Lecture 14: Quantum de Finetti Theorems II
CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4
CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped
More informationFirst-order piecewise-linear dynamic circuits
Frs-order pecewse-lnear dynamc crcus. Fndng he soluon We wll sudy rs-order dynamc crcus composed o a nonlnear resse one-por, ermnaed eher by a lnear capacor or a lnear nducor (see Fg.. Nonlnear resse one-por
More informationCalculus Chapter 1 Introduction to Calculus
Inroducon o Calculus Cal 1-3 Calculus Chaper 1 Inroducon o Calculus CHAPER 1 CALCULUS INRODUCION O hs chaper, whch replaces Chaper 4 n Physcs 2, s nended for sudens who have no had calculus, or as a calculus
More information