Dynamics of Structures

Similar documents
ECE 102 Engineering Computation

Introduction to Matrix Algebra

FREE Download Study Package from website: &

Statistics for Financial Engineering Session 1: Linear Algebra Review March 18 th, 2006

Lesson 4 Linear Algebra

Handout #2. Introduction to Matrix: Matrix operations & Geometric meaning

Addendum. Addendum. Vector Review. Department of Computer Science and Engineering 1-1

Chapter System of Equations

Ch. 12 Linear Bayesian Estimators

Section 2.2. Matrix Multiplication

Taylor Polynomials. The Tangent Line. (a, f (a)) and has the same slope as the curve y = f (x) at that point. It is the best

[ 20 ] 1. Inequality exists only between two real numbers (not complex numbers). 2. If a be any real number then one and only one of there hold.

Autar Kaw Benjamin Rigsby. Transforming Numerical Methods Education for STEM Undergraduates

Basic Maths. Fiorella Sgallari University of Bologna, Italy Faculty of Engineering Department of Mathematics - CIRAM

Lecture 2: Matrix Algebra

Mathematical Notations and Symbols xi. Contents: 1. Symbols. 2. Functions. 3. Set Notations. 4. Vectors and Matrices. 5. Constants and Numbers

The total number of permutations of S is n!. We denote the set of all permutations of S by

Section 7.3, Systems of Linear Algebraic Equations; Linear Independence, Eigenvalues, Eigenvectors (the variable vector of the system) and

Elementary Linear Algebra

Notes 17 Sturm-Liouville Theory

FOURIER SERIES PART I: DEFINITIONS AND EXAMPLES. To a 2π-periodic function f(x) we will associate a trigonometric series. a n cos(nx) + b n sin(nx),

CS 331 Design and Analysis of Algorithms. -- Divide and Conquer. Dr. Daisy Tang

DETERMINANT. = 0. The expression a 1. is called a determinant of the second order, and is denoted by : y + c 1

Inner Product Spaces (Chapter 5)

Eigenfunction Expansion. For a given function on the internal a x b the eigenfunction expansion of f(x):

1 4 6 is symmetric 3 SPECIAL MATRICES 3.1 SYMMETRIC MATRICES. Defn: A matrix A is symmetric if and only if A = A, i.e., a ij =a ji i, j. Example 3.1.

Surds, Indices, and Logarithms Radical

1. (25 points) Use the limit definition of the definite integral and the sum formulas to compute. [1 x + x2

Riemann Integral Oct 31, such that

Unit 1. Extending the Number System. 2 Jordan School District

BRILLIANT PUBLIC SCHOOL, SITAMARHI (Affiliated up to +2 level to C.B.S.E., New Delhi)

Chapter 2. LOGARITHMS

z line a) Draw the single phase equivalent circuit. b) Calculate I BC.

Fast Fourier Transform 1) Legendre s Interpolation 2) Vandermonde Matrix 3) Roots of Unity 4) Polynomial Evaluation

Frequency-domain Characteristics of Discrete-time LTI Systems

Infinite Series Sequences: terms nth term Listing Terms of a Sequence 2 n recursively defined n+1 Pattern Recognition for Sequences Ex:

Week 13 Notes: 1) Riemann Sum. Aim: Compute Area Under a Graph. Suppose we want to find out the area of a graph, like the one on the right:

Chapter 2 Infinite Series Page 1 of 9

Linear Programming. Preliminaries

SM2H. Unit 2 Polynomials, Exponents, Radicals & Complex Numbers Notes. 3.1 Number Theory

Discrete Mathematics I Tutorial 12

Chapter Real Numbers

MATH 104: INTRODUCTORY ANALYSIS SPRING 2009/10 PROBLEM SET 8 SOLUTIONS. and x i = a + i. i + n(n + 1)(2n + 1) + 2a. (b a)3 6n 2

M3P14 EXAMPLE SHEET 1 SOLUTIONS

MTH 146 Class 16 Notes

EXPONENTS AND LOGARITHMS

( ) 2 3 ( ) I. Order of operations II. Scientific Notation. Simplify. Write answers in scientific notation. III.

Accuplacer Elementary Algebra Study Guide

lecture 16: Introduction to Least Squares Approximation

Reference : Croft & Davison, Chapter 12, Blocks 1,2. A matrix ti is a rectangular array or block of numbers usually enclosed in brackets.

Numerical Methods (CENG 2002) CHAPTER -III LINEAR ALGEBRAIC EQUATIONS. In this chapter, we will deal with the case of determining the values of x 1

For all Engineering Entrance Examinations held across India. Mathematics

Project 3: Using Identities to Rewrite Expressions

Vectors. Vectors in Plane ( 2

Abel Resummation, Regularization, Renormalization & Infinite Series

Exponents and Radical

Schrödinger Equation Via Laplace-Beltrami Operator

SOLUTION OF SYSTEM OF LINEAR EQUATIONS. Lecture 4: (a) Jacobi's method. method (general). (b) Gauss Seidel method.

MATH 104: INTRODUCTORY ANALYSIS SPRING 2008/09 PROBLEM SET 10 SOLUTIONS. f m. and. f m = 0. and x i = a + i. a + i. a + n 2. n(n + 1) = a(b a) +

Lincoln Land Community College Placement and Testing Office

0 otherwise. sin( nx)sin( kx) 0 otherwise. cos( nx) sin( kx) dx 0 for all integers n, k.

=> PARALLEL INTERCONNECTION. Basic Properties LTI Systems. The Commutative Property. Convolution. The Commutative Property. The Distributive Property

INTEGRATION TECHNIQUES (TRIG, LOG, EXP FUNCTIONS)

Review of Sections

MATRIX ALGEBRA, Systems Linear Equations

ALGEBRA II CHAPTER 7 NOTES. Name

1.3 Continuous Functions and Riemann Sums

A GENERAL METHOD FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS: THE FROBENIUS (OR SERIES) METHOD

ALGEBRA. Set of Equations. have no solution 1 b1. Dependent system has infinitely many solutions

2 Dirac delta function, modeling of impulse processes. 3 Sine integral function. Exponential integral function

Definite Integral. The Left and Right Sums

Matrix. Definition 1... a1 ... (i) where a. are real numbers. for i 1, 2,, m and j = 1, 2,, n (iii) A is called a square matrix if m n.

Section 11.5 Notes Page Partial Fraction Decomposition. . You will get: +. Therefore we come to the following: x x

( a n ) converges or diverges.

Content: Essential Calculus, Early Transcendentals, James Stewart, 2007 Chapter 1: Functions and Limits., in a set B.

GRAPHING LINEAR EQUATIONS. Linear Equations. x l ( 3,1 ) _x-axis. Origin ( 0, 0 ) Slope = change in y change in x. Equation for l 1.

Solutions to Problem Set 7

Physics of Semiconductor Devices Vol.10

Matrix Algebra Notes

Chapter 7 Infinite Series

Linear Algebra. Lecture 1 September 19, 2011

Numbers (Part I) -- Solutions

Student Success Center Elementary Algebra Study Guide for the ACCUPLACER (CPT)

Certain sufficient conditions on N, p n, q n k summability of orthogonal series

CH 19 SOLVING FORMULAS

Thomas J. Osler Mathematics Department Rowan University Glassboro NJ Introduction

Algebra II, Chapter 7. Homework 12/5/2016. Harding Charter Prep Dr. Michael T. Lewchuk. Section 7.1 nth roots and Rational Exponents

Lecture 3: A brief background to multivariate statistics

a f(x)dx is divergent.

AIEEE CBSE ENG A function f from the set of natural numbers to integers defined by

Using Quantum Mechanics in Simple Systems Chapter 15

Chapter 5. The Riemann Integral. 5.1 The Riemann integral Partitions and lower and upper integrals. Note: 1.5 lectures

PROGRESSIONS AND SERIES

Inverse Matrix. A meaning that matrix B is an inverse of matrix A.

We will begin by supplying the proof to (a).

Numerical Methods. Lecture 5. Numerical integration. dr hab. inż. Katarzyna Zakrzewska, prof. AGH. Numerical Methods lecture 5 1

M.A. (ECONOMICS) PART-I PAPER - III BASIC QUANTITATIVE METHODS

CH 20 SOLVING FORMULAS

Solutions to RSPL/1. log 3. When x = 1, t = 0 and when x = 3, t = log 3 = sin(log 3) 4. Given planes are 2x + y + 2z 8 = 0, i.e.

EXERCISE a a a 5. + a 15 NEETIIT.COM

Transcription:

UNION Dymis of Strutures Prt Zbigiew Wójii Je Grosel Projet o-fied by Europe Uio withi Europe Soil Fud

UNION Mtries Defiitio of mtri mtri is set of umbers or lgebri epressios rrged i retgulr form with m s d olums. mtri is deoted by uderlied letter or bold letter. The bbrevitio dim deotes the dimesio of mtri. y elemet of mtri is desribed by two subsripts ij, where the first oe, i, deotes the umber of the d the seod, j - the umber of the olum. The elemet ij is pled o the itersetio of i-th d j-th olum. == m m m m m m m m, dim = m i m j ij mj i m i-th : i j-th olum: j elemet: ij every elemet of mtri : [] ij = i, j ij Types of mtries dim = m retgulr mtri dim B= mtri dim C= m olum mtri, vetor Projet o-fied by Europe Uio withi Europe Soil Fud

UNION dim D= squre mtri dim E= m d [E] ij = zero mtri,, m dim F= d f ij = for ij digol mtri, f f f f f f f digf, f,..., f dim I = d i ij = for ij d i ii = iity mtri I = I = = I = e.g. I, I Equlity of mtries =B [] ij =[B] ij Slr multiplitio B = [B] ij = [] ij ; - slr Mtri dditio Mtri dditio is possible oly if the dimesios of mtries re equl. C=+B [C] ij =[] ij +[B] ij dim = dim B = dim C Property of mtri dditio + B = B + + B + C = ( + B) + C = + (B + C) + = = Projet o-fied by Europe Uio withi Europe Soil Fud

UNION Mtri multiplitio B = C multiplitio is possible if dim C= mp where dim = m, dim B= p. The result of multiplitio is equl to ij = i b j = i b j i b j i b j... i b j The ft tht the multiplitio B is possible does ot me tht the multiplitio B is possible. Emple: = ; B= ; C = B=?, 9 8 7??????? B = C = 9 8 8 7, but B does ot eist! 7 Projet o-fied by Europe Uio withi Europe Soil Fud

UNION Hee, the followig phrses re used to desribe multiplitio B: mtri is right multiplied by mtri B or mtri is postmultiplied by mtri B; mtri B is left multiplied by mtri or mtri B is premultiplied by mtri. Properties of mtri multiplitio I geerl the mtri multiplitio is ot ommuttive. B C = ( B) C = (B C ) ( B) = () B = (B) ( + B) C = C + B C = B C + C C ( + B) = C + C B = C B + C = = I = I = From the results B = oe ot derive tht =, B = or B =. Mtri multiplitio is ommuttive oly for digol mtries with the sme dimesios: {} {b} = {b} {} = {b} dim{} =dim{b} Emple: =, {b} = dig (, ), {b}=? {b} =? 9 {b}=. It is olum multiplitio, i.e. the i-th olum of mtri is multiplied by the i-th elemet of digol {b} {b} = [ j b j ] {b} =. It is multiplitio, i.e. the j-th of mtri is multiplied by the j-th elemet of digol {b} {b} = [b i i ]. Projet o-fied by Europe Uio withi Europe Soil Fud

UNION Mtri trspositio = B T dim = m, dim B= m d [] ij = [B] ji Emple: T B= ; B T = = B T B = T = ( T ) T (+B+C) T = T +B T +C T ( B C) T = C T B T T If = T the is symmetri mtri, [] ij = [] ji If = T the is sew-symmetri mtri, [] ij = -[] ji, [] ii = Every digol mtri is lso symmetri mtri {} T = {}. Tre of mtri dim =, tr =... tr T = tr tr (+B+C) = tr + tr B + tr C tr ( ) = (tr ) tr ( B) = tr (B ) tr ( B C) = tr (B C ) = tr (C B ) tr ( C B ) Projet o-fied by Europe Uio withi Europe Soil Fud

UNION Emple: tr 7 Determit of mtri ermit of mtri is deoted by. The ermit is slr d is defied solely for the squre mtri dim =. The esiest wy to defie the ermit is reursio. For mtri with dimesios, the ermit is equl i vlue to the oly elemet of this mtri. For mtri with dimesios, the ermit is equl to: For mtri with dimesios, the rule of Srrus be used. Two first olums should be writte o the right side of the mtri (or two first s -- uder the mtri). The produt of the elemets rossed with blue lie re dded, d the oes rossed with red lie re subtrted: For mtri with dimesios dim = is defied by ermits of mtries with dimesios (-)(-). The mior M ij is mtri obtied from mtri by ersig the i-th d the j-th olum. The oftor ij of the elemet ij is ermit of the mior M ij multiplied by (- ) i+j. ij = (-) i+j M ij The ermit my be lulted by pplitio of the Lple epsio log y i-th Projet o-fied by Europe Uio withi Europe Soil Fud

Projet o-fied by Europe Uio withi Europe Soil Fud UNION = i i M i = i i or log y j-th = j j M j = j j Emple? Usig the Lple epsio log the rd olum d the rule of Srrus, the ermit be esily lulted: 7 8 ) ( ) ( ddig olum/ multiplied by y umber to other olum/ does ot hge the vlue of the ermit. This property mes it possible to lulte the ermit esily. I the emple below: () the d multiplied by - is dded to the th. fter this, there is oly oe ozero elemet i the rd olum () the Lple epsio is pplied log the rd olum () the st multiplied by - is dded to the d d the st multiplied by is dded to the rd. () the Lple epsio is pplied log the st olum d the ermit be lulted esily for mtri.

UNION () 8 ( ) () () 8 ( 8) 7 () Some properties of ermit: T = ( ) = ( ) ( B C) = B C {} =... is lled sigulr mtri or degeerte mtri if =. is lled o-sigulr mtri (odegeerte mtri or ivertible mtri) if. Iversio of mtri squre mtri B with the sme dimesios s mtri is lled iverse to it whe: B = B = I The symbol - is used for mtri iverse to. Usig the properties of the ermit: =, = ( ) = ( ) - = - = I ( - ) = - = I = Projet o-fied by Europe Uio withi Europe Soil Fud

Projet o-fied by Europe Uio withi Europe Soil Fud UNION - = / - eist oly if ( is o-sigulr mtri). Mtri - by derived usig formul: - = dj / where dj is djugte (or djoit) mtri i.e. trsposed mtri of oftors ij. dj = T = Emple: =, = 9++ = dj = T = 9 9 T - = dj / =,, 7,, y mtri ould be trsformed ito iity mtri I by performig the followig trsformtios (multiplitio of the /olum by slr, ddig s/olums). By pplyig the sme trsformtios to the iity mtri oe obti mtri iverse to.

Projet o-fied by Europe Uio withi Europe Soil Fud UNION =, I = rd d d st rd st (step) =, I (step) = / elemet d d d st st (step ) =, I (step ) = / d rd d rd st st rd rd (step ) = =I, I (step ) =,, 7,, = - Verifitio of: - =?, - =?,, 7,,,, 7,, Vetors Mtries with oly oe olum re ofte lled vetors. Vetors re deoted by overliig = ol (,,... ) = = [,,... ] T Two vetors re orthogol if y y T T

UNION T T Two vetors re orthogol with weight if y y. The ormliztio of vetor is result of dividig elemets of this vetor by umber. The most ommo method of ormliztio is dividig vetor ) by its mimum bsolute vlue or ) by the Eulide orm. Emple: 7 ) m i mimum orm 7 T ) - Eulide orm ( 7) 9 Normlized vetor orm = / ) ol,, - orm, orm 7 7 7 ) ol, orm, 9 9 9 i orm System of lier equtios The system of lier equtios with uows (,,, ) ould be writte dow with mtri ottio......... y y y y Projet o-fied by Europe Uio withi Europe Soil Fud

UNION system of lier equtios is lled homogeous system if y =, d ohomogeous system if y. The eistee of solutio depeds o the type of the system (homogeous or ot) d o the properties of mtri (if it is sigulr mtri or ot).. ohomogeous system ( y ).. solutio eists = - y.. = solutio does ot eist. homogeous system ( y = ).. solutio eists d is trivil =.. = there is ifiite umber of solutios = dj v, v - y rbitrry vetor, for v =ol (,,...,,... ) the solutio is y o-zero olum of mtri dj ; geerlly the solutio is lier ombitio of olums of mtri dj. Eigeproblem Numbers i re eigevlues of squre mtri if ozero vetor equtio w = i w. Vetors w = w i i wi re lled eigevetors. w i eists whih fulfils the w w = ( I) w = ; w ( I) = j i ij i j ( I) =... hrteristi equtio Projet o-fied by Europe Uio withi Europe Soil Fud

UNION The solutios of this hrteristi equtio,,..., re eigevlues. Emple: Clulte eigevlues of mtri =, hrteristi equtio:, 9 solutios, 7 The eigevetor is y o-zero olum of the mtri djoit to mtri ( i I) = ; B =( I) =, dj B = ormliztio, w, = 7; B =( I) =, dj B = ormliztio, w y produt of w is lso eigevetor, where - y umber. i w, w, W = w w =, mtri W, the i-th olum of whih is equl to w i, is lled eigemtri. Usig system Mthemti for mtri lultios I the system Mthemti, mtries re defied s lists of lists e.g. lists of s. The system gives output i the sme form. I order to see the mtri i the form of tble, oe eeds to use the ommd TrditiolForm or //MtriForm. Projet o-fied by Europe Uio withi Europe Soil Fud

UNION I[]:=,,,,,,,, Out[]=,,,,,,,, I[]:= Out[]=,,,,,,,, I[]:= MtriForm Out[]//MtriForm= I[]:= TrditiolForm Out[]//TrditiolForm= I[]:= Out[]= v,,,, I[]:= MtriFormv Out[]//MtriForm= I order to trspose the mtri, the ommd Trspose should be used. The trspositio of vetor with oe olum yields mtri with oe oly. I[7]:= Out[7]= Trspose,,,,,,,, I[8]:= TrditiolFormTrspose Out[8]//TrditiolForm= I[9]:= Trsposev MtriForm Out[9]//MtriForm= For mtries, the multiplitio symbol. (dot) is used. Projet o-fied by Europe Uio withi Europe Soil Fud

UNION I[]:=. MtriForm Out[]//MtriForm= 8 8 8 I[]:=.v MtriForm Out[]//MtriForm= 8 I[]:= Trsposev. MtriForm Out[]//MtriForm= 8 The ommd Det is used for lultig the ermit, Tr is used for lultig the tre Iverse is used for lultig the iversio of the mtri. I[]:= Det Out[]= I[]:= TrditiolFormIverse Out[]//TrditiolForm= 9 I[]:= Tr Out[]= 7 Projet o-fied by Europe Uio withi Europe Soil Fud

UNION I[7]:= TrditiolFormIverse. Out[7]//TrditiolForm= I[8]:= TrditiolFormIityMtri Out[8]//TrditiolForm= I[9]:= DigolMtri,, MtriForm Out[9]//MtriForm= The ommd IityMtri be used i order to delre iity mtri, d the ommd DigolMtri to delre digol mtri. The ommd Eigevlues yields the eigevlues of the mtri, rrged from the gretest to the smllest. The ommd Eigevetors yields the eigevetors of the mtri. These eigevetors re rrged i s! The first orrespods to the gretest eigevlue. The ommd Eigesystem ombies two previous ommds. I[]:=,,, ; MtriForm Out[]//MtriForm= I[7]:= Out[7]= Eigevlues 7, I[8]:= Eigevetors MtriForm Out[8]//MtriForm= I[9]:= Eigesystem MtriForm Out[9]//MtriForm= 7,, Projet o-fied by Europe Uio withi Europe Soil Fud

Projet o-fied by Europe Uio withi Europe Soil Fud UNION Problems The followig mtries re give: =, B =, C =, D =, E =, F = 7, G = 7, U =, V =, X =, L = dig (, ), M = dig (,, ), N = dig (,,, ).. Evlute: () +, (b), () B, (d) B C T, (e) E F, (f) D V, (g) X T G, (h) D M, (i) N G, (j) V T V, () V V T.. Evlute: () tr G, (b) C, () D, (d) G, (e) (L C L).. Evlute: () dj B, (b) dj N, () dj (U U T ).. Evlute eigevlues d eigevetors of mtries: (), (b) D, () G.. Solve the system of equtios: y z y z y. Fid the vlue of for whih otrivil solutio v of these equtios eists: (K B ) v, Fid vetors v d ormlize them () K =, B = m m, fid the solutio for m d or ssume =, m=, (b) K =, B = m m m, =, m=,

Projet o-fied by Europe Uio withi Europe Soil Fud UNION swers to problems. () + =, (bd) = B = B C T = 8, (e) E F = 8 8 9, (f) D V =, (g) X T G = 7 9 (h) D M =, (i) N G = 8 8, (j) V T V =[], () V V T = 9,. () tr G =, (b) C= 8 () D = - (d) G = (e) (L C L) = =,. () dj B, (b) dj N = dig (8,,, ) () dj (U U T ) = dj =. () : = =,77, = + =,88,. W =,,. D = ; D()= D I; D() = = 7. -++8 - = =,9 =,7 = 7,9

UNION 8. eigemtri W =,8,,,88,8 ;, 9. eg. dj D( ) = dj,,, =,,8 9,, 7, 9 8, T,,,8 9, olums' ormliztio, =, 7, 9 8, 8,7,, 8,7,8,,8,,8,. G = 7, The hrteristi equtio + + = =,89 ; =,,987 i; =, +,987 i; =8,7. Eigemtri W =,,9,7,7,88 i,,i,,77i,7,88 i,,i,,77i,,8,,. y z y z y System mtri = Projet o-fied by Europe Uio withi Europe Soil Fud

Projet o-fied by Europe Uio withi Europe Soil Fud UNION dj =,,,, 8 ormliztio olums' z y - y umber. () K =, B = m m, (K B ) = m + m = = ( ) / m=,8 /m = ( + ) / m =,8 /m W =,8,8 =,8,8 ; v - y olum from W. (b) K =, B = m m m, =, m=, -+ - = =,98 =, =,7 v - y olum from,8,,,8,,8