MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL
|
|
- Ralph Booker
- 6 years ago
- Views:
Transcription
1 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1: Specl Issue No 1 MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL M BOTIŞ 1 Astrct: I te ler lyss of structures troug modl superposto fter decouplg te moto equtos, tere must e solved te moto equtos for ec degree of dymc freedom, troug umercl tegrto, y usg te Dumel covoluto tegrl To fd structurl respose to ertquke cto t s lso eeded to determe te mxmum spectrl vlues dsplcemet, speed d ccelerto for y ccelerogrm of te ste; tese vlues re determed y umercl tegrto wt Dumel tegrl I ts pper t s preseted Mtl progrm for tegrtg umercl of Dumel covoluto tegrl wt trpezods metod tt c e used ot for modl lyss d spectrl lyss of structures Key words: Dumel tegrl covoluto, dymc lyss y modl superposto, spectrl lyss ertquke egeerg, dymc respose of te structures for rdom forces 1 Itroducto Dymc lyss of cvl structures s ecessry we ctos re tme-vryg If dymc lyss s ler, modl superposto prcple c e used To use modl superposto prcple we must frst ceve modl decomposto d te mke modl superposto After decomposto modl vrto modes re oted, wc glgt modl deformed spes ssocted wt ec mode of vrto d ertl mss correspodg to ec mode of vrto Troug modl decomposto two mportt results re oted, oe s qulttve d s represeted y te form of ow modes of vrto d oter s qutttvely d dctes level of prtcpto of te ertl msses te movemet of dymc system Te lyss of te vrto mode s very mportt ecuse ec vrto mode s cotruto for te moto of te dymc system After spectrl decomposto, modl superposto c e mde, wc volves composg sgfct vrto modes moto dymc system I cse of system wt degrees of freedom (fg1 dymc equto of moto s: [ M ]{ x} [ C]{ x} [ K]{ x} { F} (1 were: -[M ] ertl mtrx of dymc system; -[C] dmpg mtrx of dymc system; -[K] stffess mtrx of dymc system; 1 Dept of Cvl Egeerg, Trslv Uversty of Brşov
2 46 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1 ẋ -{ } ccelertos vector for te degrees of freedom dymc; -{ẋ} velocty vector for te degrees of freedom dymc; -{x} dsplcemet vector for te degrees of freedom dymc; -{F} force vector for te degrees of freedom dymc Istteous dymc respose o te degree of freedom k s oted y ler superposto of stteous dymc resposes correspodg to te modes of vrto fg: k x ( x ( t ( 1 k, To dscoect equtos of moto te followg ler trsformto s used: k 1 k, ( k, ( 1 x ( x (3 Reltosp etwee turl coordtes dymc system x k (, k 1 d prcpl coordtes or modl (, 1 c e expressed y te modl mtrx s follows: Fg 1 Te system wt degrees of freedom dymc x x x x 1 1 1,1,1 1,1,1 1, 1,,, 1, 1 1, 1, 1, 1 1, 1,,, 1 1 Dsplcemet vector velocty d ccelerto o te dymc degrees of freedom k 1 ecome: Fg Vrtos mode d egevlues for system wt degrees of freedom dymc { } [ ]{ ( }; { } [ ]{ ( }; (4 { } [ ]{ ( } were: -[] modl mtrx -{ ( } dsplcemets vector turl coordtes -{ ( } velocty vector turl coordtes
3 M BOTIŞ: Mtl progrm for te umercl soluto of Dumel covoluto tegrl 47 -{ ( } ccelertos vector turl coordtes Replcg te system equtos of moto (1 rel coordtes wt modl coordtes, te equtos of moto s relesed Due to ortogolty of vectors ssocted vrto modes system of equtos (1 s reduced to decoupled equtos of moto for ec degree of freedom dymc: As te system of dfferetl equtos of moto s decoupled, we c wrte te correspodg equto of moto for ec degree, s follows: F { ( } { ( } { ( } (5 M 1 were: T - F { F( } correspodg equvlet force mode ; T - M M equvlet mss mode ; - correspodg dmpg mode -[ C] [ M ] [ K] Ryleg dmpg mtrx To solve decoupled equtos of moto, Dumel covoluto tegrl c e used Dsplcemet respose oted usg Te covoluto tegrl Dumel tegrl o ec dymc degrees of freedom s: t 1 ( t ( F ( e s ( t d M 1 (6 were, - 1 dmped pulsto If you kow te swer dsplcemets o ec dymc degree of freedom te we c clculte veloctes ccelertos d ertl forces for dymc system: { } [ ]{ ( }; { } [ ]{ ( }; { } [ ]{ ( } Te lyss of te dymcs of system wt degrees of freedom dymcs c e see s solvg te system of equtos s reduced to clculto covoluto tegrl Dumel For ts reso ts pper presets Mtl progrm for umercl tegrto of covoluto tegrl Dumel If te t pply H mpulse of te system, ts s removed from equlrum (7 posto wt velocty x H / m, (fg3 d performs free dmped vrtos descred y te equto : H e m H ( t e m ( t s ( t s ( t for t (8 Fg 3 Free respose dsplcemet due to mpulse fte H If te dymc system rdom dsturg force cts,(fg4 equto
4 48 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1 of moto of te dymc system s: F( m Acto of dsturg force ppled dymc system t rest c e represeted s cotuous elemetry mpulses, o te tme tervl from t t After pplyg elemetry mpulse dsplcemet respose of te dymc system dymc o degree of freedom s: dh ( t d e s ( t d m (9 F( ( t e s ( t d m Dymc dsplcemet respose of te system s ceved y superposto of effects s tegrl mout of moto geerted y elemetry pulse successo o te legt of tme tt force cts [1]: F( ( t e s ( t d m f ( ( t F e s ( t d m Fg 4 Dsplcemet respose oted y superposto of te effects of ftesml mpulses dh If dsturg force s ppled to dymc system wle t s moto due to tl dsplcemet d velocty codtos, te dsplcemet respose s expressed: e [( x x cost s st] f ( ( t F e s ( t d m were, - x tl dsplcemet for t= - x tl velocty for t= It c e see tt f te dymc system wt sgle degree of freedom dymc system to determe dymc respose Dumel covoluto tegrl must e clculted Mtemtcl spects of umercl soluto of te Dumel tegrl To tegrte umerclly covoluto tegrl Dumel, furter s preseted trpezods metod wc s te smplest metod used prctce Trpezods metod s sed o pproxmtg te re delmted y te grp of fucto f :[, ] R wt te re of trpezod If te tervl [, ] s dvded to two sutervls [, x1] d [ x 1, ] te tegrl of fucto f(x c e pproxmted y te relto: dx x1 dx dx [ f ( f ( x1] [ f ( f ( x1] [ f ( f ( x1 f ( ] [ f ( f ( x1 f ( ] 4 To reduce te error of tegrto, te tervl [,] my e dvded sutervls wt dvso ( /, x1
5 M BOTIŞ: Mtl progrm for te umercl soluto of Dumel covoluto tegrl 49 tegrl of fucto f(x c e pproxmted y te relto: dx x1 dx x11 dx [ f ( f ( x1] [ f ( f ( x 1] [ f ( f ( x1 f ( x 1 f ( ] [ f ( f ( x1 f ( x 1 f ( ] Error clcultg troduced y te metod of trpezods depeds o te tegrto step d s te followg expresso o tervl [ x x+]:: 3,, E ( f ( 1 Full error s: 1 3 1,, ET ( E ( f ( 1 I 1 I 1 3 Implemetto of Mtl progrm for clculto Dumel tegrl For umercl clculto of te covoluto tegrl Dumel strts from respose dsplcemet for system wt oe degree of freedom: F( ( t e s( t d m F( t e cos d e st m F( t e s d e cot m t t A( e st B( e cost were, F( A( e cos d m F( B( e s d m To solve te Dumel tegrl trpezods metod wll e used to clculte tegrls A( d B( wt Mtl progrm [] % Solvg covoluto tegrl Dumel %y umercl tegrto metod of %trpezods %Aswer dplcemet, velocty d %ccelerto for dymc sstem oe degree %of freedom % force rmp cler;clf; % Tme step used =; % Mss system wt 1GLD m=175; % Rgdty of te system wt 1GLD k=7; % Frcto of crtcl dmpg u=5/1; % Pulsto system 1GLD wtout %dmpg omeg=sqrt(k/m; % Pulsto system 1GLD wt %dmpg omegd=omegsqrt(1-u^; % Tme tervl tt s computed %respose system wt 1 DOF for %dsplcemet, velocty d ccelerto tmx=14; % Force rmp fucto vrto cure=[ ]; % Numer of pots of fucto curves pt=msze(cure; t=; f=terp1([cure(1,:p]', [cure(,:p]',; y=f1exp(uomeg cos(omegd(1momegd^-1; y=f1exp(uomeg s(omegd(1momegd^-1; % Numer of steps for umercl tegrto ps=tmx/; suplot(,1,1 for =1:ps t=t+; % Clculte te force t tme t fort=terp1([cure(1,:p]', [cure(,:p]', f1=fort5;
6 5 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1 f(=f1; tmp(=t; % Vrto represetto force %fucto tme p=plot(t,f1,'--rs','lewdt',, 'MrkerEdgeColor','', 'MrkerFceColor','', 'MrkerSze',; old o; grd o; % Clculto of vlues y d y %from Dumel covoluto tegrl y(=f11exp(uomeg cos(omegd(1momegd^-1; y(=f11exp(uomeg s(omegd(1momegd^-1; ed % Fl vlues y=[y y(1,:' ]; y=[y y(1,:']; % Metod trpezods t= suplot(,1, for =:ps+1 t=t+; A(1=; B(1=; A(=A(-1+(y(+y(-1/; B(=B(-1+(y(+y(-1/; % Respose moto for rmp sgl %cse 1GLD dymc system d(=(a(exp(-uomeg s(omegd- B(exp(-uomeg cos(omegd(k/5; % Grpcl represetto of te rto of %dymc dsplcemet d sttc %dsplcemet plot(t,d(,'--rs','lewdt',, 'MrkerEdgeColor','r', 'MrkerFceColor','r', 'MrkerSze',; old o; grd o; ed 4 Results d cocluso Te progrm proposed te pper s ee tested te cse of dymcl system wt oe degree of freedom dymc Dymc lyss ws performed for two steps of tegrto = d =1 Te results fter umercl tegrto re preseted fg5 d fg6 Fg 5 Te respose dsplcemet for dymc system wt 1 degree of freedom, tegrto step, Fg 6 Te respose dsplcemet for dymc system wt 1 degree of freedom, tegrto step,1 Te result oted y umercl tegrto ws compred wt te exct soluto, for te tegrto step =, ot error 13,7 % d f tegrto step s = error ws,3% Altoug te progrm preseted s smple, t c e successfully used f te tegrto step s lower t were T/1, were T s ow perod for dymc system Refereces 1 Ifrm, M: Dmc structurlor s gere sesmc Edtur ddctc s pedgogc, 1973 Adelw Kr, Rold B Gueter,: Itroducto to Numercl Metods: A MATLAB Approc CRC Press Tylor &Frcs Group, 1
MTH 146 Class 7 Notes
7.7- Approxmte Itegrto Motvto: MTH 46 Clss 7 Notes I secto 7.5 we lered tht some defte tegrls, lke x e dx, cot e wrtte terms of elemetry fuctos. So, good questo to sk would e: How c oe clculte somethg
More informationIn Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is
Mth Sprg 08 L Approxmtg Dete Itegrls I Itroducto We hve studed severl methods tht llow us to d the exct vlues o dete tegrls However, there re some cses whch t s ot possle to evlute dete tegrl exctly I
More informationThe z-transform. LTI System description. Prof. Siripong Potisuk
The -Trsform Prof. Srpog Potsuk LTI System descrpto Prevous bss fucto: ut smple or DT mpulse The put sequece s represeted s ler combto of shfted DT mpulses. The respose s gve by covoluto sum of the put
More informationNumerical differentiation (Finite difference)
MA umercl Alyss Week 4 To tk d wrte portle progrms we eed strct model of ter computtol evromet. Ftful models do est ut tey revel tt evromet to e too dverse forcg portle progrmmers to lot eve te smplest
More informationChapter Simpson s 1/3 Rule of Integration. ( x)
Cpter 7. Smpso s / Rule o Itegrto Ater redg ts pter, you sould e le to. derve te ormul or Smpso s / rule o tegrto,. use Smpso s / rule t to solve tegrls,. develop te ormul or multple-segmet Smpso s / rule
More informationMathematically, integration is just finding the area under a curve from one point to another. It is b
Numerl Metods or Eg [ENGR 9] [Lyes KADEM 7] CHAPTER VI Numerl Itegrto Tops - Rem sums - Trpezodl rule - Smpso s rule - Rrdso s etrpolto - Guss qudrture rule Mtemtlly, tegrto s just dg te re uder urve rom
More informationChapter 2 Intro to Math Techniques for Quantum Mechanics
Wter 3 Chem 356: Itroductory Qutum Mechcs Chpter Itro to Mth Techques for Qutum Mechcs... Itro to dfferetl equtos... Boudry Codtos... 5 Prtl dfferetl equtos d seprto of vrbles... 5 Itroducto to Sttstcs...
More informationthis is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ]
Atervtves The Itegrl Atervtves Ojectve: Use efte tegrl otto for tervtves. Use sc tegrto rules to f tervtves. Aother mportt questo clculus s gve ervtve f the fucto tht t cme from. Ths s the process kow
More informationChapter 12-b Integral Calculus - Extra
C - Itegrl Clulus Cpter - Itegrl Clulus - Etr Is Newto Toms Smpso BONUS Itroduto to Numerl Itegrto C - Itegrl Clulus Numerl Itegrto Ide s to do tegrl smll prts, lke te wy we preseted tegrto: summto. Numerl
More informationUNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 2007/2008 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 2007 MASA : 3 JAM
UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 7/8 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 7 MASA : 3 JAM KOD KURSUS : KKKQ33/KKKF33 TAJUK : PENGIRAAN BERANGKA ARAHAN :.
More informationCS321. Numerical Analysis
CS Numercl Alyss Lecture 4 Numercl Itegrto Proessor Ju Zg Deprtmet o Computer Scece Uversty o Ketucky Legto, KY 456 6 Octoer 6, 5 Dete Itegrl A dete tegrl s tervl or tegrto. For ed tegrto tervl, te result
More informationOn Several Inequalities Deduced Using a Power Series Approach
It J Cotemp Mth Sceces, Vol 8, 203, o 8, 855-864 HIKARI Ltd, wwwm-hrcom http://dxdoorg/02988/jcms2033896 O Severl Iequltes Deduced Usg Power Seres Approch Lored Curdru Deprtmet of Mthemtcs Poltehc Uversty
More informationCS321. Introduction to Numerical Methods
CS Itroducto to Numercl Metods Lecture Revew Proessor Ju Zg Deprtmet o Computer Scece Uversty o Ketucky Legto, KY 6 6 Mrc 7, Number Coverso A geerl umber sould be coverted teger prt d rctol prt seprtely
More informationIntegration by Parts for D K
Itertol OPEN ACCESS Jourl Of Moder Egeerg Reserc IJMER Itegrto y Prts for D K Itegrl T K Gr, S Ry 2 Deprtmet of Mtemtcs, Rgutpur College, Rgutpur-72333, Purul, West Begl, Id 2 Deprtmet of Mtemtcs, Ss Bv,
More informationi+1 by A and imposes Ax
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 09.9 NUMERICAL FLUID MECHANICS FALL 009 Mody, October 9, 009 QUIZ : SOLUTIONS Notes: ) Multple solutos
More informationunder the curve in the first quadrant.
NOTES 5: INTEGRALS Nme: Dte: Perod: LESSON 5. AREAS AND DISTANCES Are uder the curve Are uder f( ), ove the -s, o the dom., Prctce Prolems:. f ( ). Fd the re uder the fucto, ove the - s, etwee,.. f ( )
More informationModule 2: Introduction to Numerical Analysis
CY00 Itroducto to Computtol Chemtr Autum 00-0 Module : Itroducto to umercl Al Am of the preet module. Itroducto to c umercl l. Developg mple progrm to mplemet the umercl method opc of teret. Iterpolto:
More informationArea and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]
Are d the Defte Itegrl 1 Are uder Curve We wt to fd the re uder f (x) o [, ] y f (x) x The Prtto We eg y prttog the tervl [, ] to smller su-tervls x 0 x 1 x x - x -1 x 1 The Bsc Ide We the crete rectgles
More informationEvaluation of the Singular Stress Field. of the Crack Tip of the Infinite Elastic Plate
ICCM05, 4-7 th July, Auckld, NZ Evluto of the Sgulr Stress Feld of the Crck Tp of the Ifte Elstc lte * J.Y. Lu,, G. L, Z.C.Su 3,. Zhg d Y.Y. Wg Fculty of Ifrstructure Egeerg, Dl Uversty of Techology, Dl
More informationSpectrum Estimation by Several Interpolation Methods
IJCSNS Itertol Jourl o Computer Scece d Network Securty VOL.6 No.A Februry 006 05 Spectrum Estmto by Severl Iterpolto Metods Mbu Isr Oym Ntol College o Tecology Oym-S Tocg 33-0806 JAPAN Summry I ts pper
More informationINTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek
ELM Numerl Alss Dr Murrem Merme INTEROLATION ELM Numerl Alss Some of te otets re dopted from Luree V. Fusett Appled Numerl Alss usg MATLAB. rete Hll I. 999 ELM Numerl Alss Dr Murrem Merme Tod s leture
More informationEFFECT OF CRACK PARAMETERS ON FREE VIBRATIONS OF THE BERNOULLI-EULER BEAM
Jourl of Appled themtcs d Computtol echcs 05 () 67-7 www.mcm.pcz.pl p-in 99-9965 DOI: 0.75/jmcm.05..7 e-in 353-0588 EFFECT OF CRACK PARAETER ON FREE VIBRATION OF THE BERNOULLI-EULER BEA Izbel Zmorsk Dwd
More informationCS473-Algorithms I. Lecture 3. Solving Recurrences. Cevdet Aykanat - Bilkent University Computer Engineering Department
CS473-Algorthms I Lecture 3 Solvg Recurreces Cevdet Aykt - Blket Uversty Computer Egeerg Deprtmet Solvg Recurreces The lyss of merge sort Lecture requred us to solve recurrece. Recurreces re lke solvg
More informationFree Vibration Analysis of Rectangular Plates Using Galerkin-Based Finite Element Method
Itertol Jourl of Mechcl Egeerg ISSN : 77-709 Volume Issue http://.meourl.com/ https://stes.google.com/ste/ourlme/ Free Vbrto Alyss of Rectgulr Pltes Usg Glerk-Bsed Fte Elemet Method Nefft M. Werfll, Abobker
More informationSoo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:
Soo Kg Lm 1.0 Nested Fctorl Desg... 1.1 Two-Fctor Nested Desg... 1.1.1 Alss of Vrce... Exmple 1... 5 1.1. Stggered Nested Desg for Equlzg Degree of Freedom... 7 1.1. Three-Fctor Nested Desg... 8 1.1..1
More information14.2 Line Integrals. determines a partition P of the curve by points Pi ( xi, y
4. Le Itegrls I ths secto we defe tegrl tht s smlr to sgle tegrl except tht sted of tegrtg over tervl [ ] we tegrte over curve. Such tegrls re clled le tegrls lthough curve tegrls would e etter termology.
More informationDESIGNING AND CONSTRUCTION OF ROADS, SUBWAYS, AIRFIELDS, BRIDGES AND TRANSPORT TUNNELS
Issue (), ISSN 5- DESIGNING AND CONSTRUCTION OF ROADS, SUBWAYS, AIRFIELDS, BRIDGES AND TRANSPORT TUNNELS UDC 59.. A. Kozlov THE DEFLECTED MODE OF MULTI COHERENT PRISMATIC CONSTRUCTIE ELEMENTS OF BRIDGE
More informationCURVE FITTING LEAST SQUARES METHOD
Nuercl Alss for Egeers Ger Jord Uverst CURVE FITTING Although, the for of fucto represetg phscl sste s kow, the fucto tself ot be kow. Therefore, t s frequetl desred to ft curve to set of dt pots the ssued
More informationDifferential Method of Thin Layer for Retaining Wall Active Earth Pressure and Its Distribution under Seismic Condition Li-Min XU, Yong SUN
Itertol Coferece o Mechcs d Cvl Egeerg (ICMCE 014) Dfferetl Method of Th Lyer for Retg Wll Actve Erth Pressure d Its Dstrbuto uder Sesmc Codto L-M XU, Yog SUN Key Lbortory of Krst Evromet d Geologcl Hzrd
More informationNumerical Solution of Second order Integro- Differential Equations(Ides) with Different Four Polynomials Bases Functions
Numercl Soluto of Secod order Itegro- Dfferetl Equtos(Ides) wth Dfferet Four Polyomls Bses Fuctos wo O A R M Deprtmet of Mthemtcs Uersty of Ilor Deprtmet of Mthemtcs d Sttstcs he Poly Id Astrct: - I ths
More informationDATA FITTING. Intensive Computation 2013/2014. Annalisa Massini
DATA FITTING Itesve Computto 3/4 Als Mss Dt fttg Dt fttg cocers the problem of fttg dscrete dt to obt termedte estmtes. There re two geerl pproches two curve fttg: Iterpolto Dt s ver precse. The strteg
More informationME 501A Seminar in Engineering Analysis Page 1
Mtr Trsformtos usg Egevectors September 8, Mtr Trsformtos Usg Egevectors Lrry Cretto Mechcl Egeerg A Semr Egeerg Alyss September 8, Outle Revew lst lecture Trsformtos wth mtr of egevectors: = - A ermt
More informationLevel-2 BLAS. Matrix-Vector operations with O(n 2 ) operations (sequentially) BLAS-Notation: S --- single precision G E general matrix M V --- vector
evel-2 BS trx-vector opertos wth 2 opertos sequetlly BS-Notto: S --- sgle precso G E geerl mtrx V --- vector defes SGEV, mtrx-vector product: r y r α x β r y ther evel-2 BS: Solvg trgulr system x wth trgulr
More informationChapter 2 Intro to Math Techniques for Quantum Mechanics
Fll 4 Chem 356: Itroductory Qutum Mechcs Chpter Itro to Mth Techques for Qutum Mechcs... Itro to dfferetl equtos... Boudry Codtos... 5 Prtl dfferetl equtos d seprto of vrbles... 5 Itroducto to Sttstcs...
More informationRoberto s Notes on Integral Calculus Chapter 4: Definite integrals and the FTC Section 2. Riemann sums
Roerto s Notes o Itegrl Clculus Chpter 4: Defte tegrls d the FTC Secto 2 Rem sums Wht you eed to kow lredy: The defto of re for rectgle. Rememer tht our curret prolem s how to compute the re of ple rego
More information19 22 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y Write the given sum or difference as a single integral in
SECTION. THE DEFINITE INTEGRAL. THE DEFINITE INTEGRAL A Clck here for swers. S Clck here for solutos. Use the Mdpot Rule wth the gve vlue of to pproxmte the tegrl. Roud the swer to four decml plces. 9
More information20 23 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y 2 2
SECTION 5. THE DEFINITE INTEGRAL 5. THE DEFINITE INTEGRAL A Clck here for swers. S Clck here for solutos. 7 Use the Mdpot Rule wth the gve vlue of to pproxmte the tegrl. Roud the swer to four decml plces.
More informationScienceDirect. About Verification of Discrete-Continual Finite Element Method of Structural Analysis. Part 2: Three-Dimensional Problems
Avlle ole t wwwscecedrectcom SceceDrect Proced Egeerg 9 (04 4 9 XXIII R-S-P semr heoretcl Foudto of Cvl Egeerg (RSP (FoCE 04 Aout Verfcto of Dscrete-Cotul Fte Elemet Method of Structurl Alyss Prt : hree-dmesol
More informationAnalytical Approach for the Solution of Thermodynamic Identities with Relativistic General Equation of State in a Mixture of Gases
Itertol Jourl of Advced Reserch Physcl Scece (IJARPS) Volume, Issue 5, September 204, PP 6-0 ISSN 2349-7874 (Prt) & ISSN 2349-7882 (Ole) www.rcourls.org Alytcl Approch for the Soluto of Thermodymc Idettes
More informationHopf Bifurcation of Maglev System with Coupled Elastic Guideway
Hopf Bfurcto of Mglev System wth Coupled Elstc udewy No Loghu She, Ho Wg, Dogsheg Zou, Zhzhou Zhg* d Wese Chg Egeerg Reserch Ceter of Mglev Techology, Chgsh Isttute of Techology, Chgsh 7, Hu, PRCh sheloghu@yhoocomc,
More informationITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS
Numercl Alyss for Egeers Germ Jord Uversty ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl soluto of lrge systems of ler lgerc equtos usg drect methods such s Mtr Iverse, Guss
More informationSequences and summations
Lecture 0 Sequeces d summtos Istructor: Kgl Km CSE) E-ml: kkm0@kokuk.c.kr Tel. : 0-0-9 Room : New Mleum Bldg. 0 Lb : New Egeerg Bldg. 0 All sldes re bsed o CS Dscrete Mthemtcs for Computer Scece course
More informationItō Calculus (An Abridged Overview)
Itō Clculus (A Abrdged Overvew) Arturo Ferdez Uversty of Clfor, Berkeley Sttstcs 157: Topcs I Stochstc Processes Semr Aprl 14, 211 1 Itroducto I my prevous set of otes, I troduced the cocept of Stochstc
More informationPreliminary Examinations: Upper V Mathematics Paper 1
relmr Emtos: Upper V Mthemtcs per Jul 03 Emer: G Evs Tme: 3 hrs Modertor: D Grgortos Mrks: 50 INSTRUCTIONS ND INFORMTION Ths questo pper sts of 0 pges, cludg swer Sheet pge 8 d Iformto Sheet pges 9 d 0
More informationStrategies for the AP Calculus Exam
Strteges for the AP Clculus Em Strteges for the AP Clculus Em Strtegy : Kow Your Stuff Ths my seem ovous ut t ees to e metoe. No mout of cochg wll help you o the em f you o t kow the mterl. Here s lst
More informationSt John s College. UPPER V Mathematics: Paper 1 Learning Outcome 1 and 2. Examiner: GE Marks: 150 Moderator: BT / SLS INSTRUCTIONS AND INFORMATION
St Joh s College UPPER V Mthemtcs: Pper Lerg Outcome d ugust 00 Tme: 3 hours Emer: GE Mrks: 50 Modertor: BT / SLS INSTRUCTIONS ND INFORMTION Red the followg structos crefull. Ths questo pper cossts of
More informationNon-Polynomial Spline Method for the Solution of Problems in Calculus of Variations
World Acdemy of Scece Egeerg d Tecology No-Polyoml Sple etod for te Solto of Prolems Clcls of Vrtos. Zre. Hosyr. Sedgt Astrct I ts pper mercl solto sed o opolyoml cc sple fctos s sed for fdg te solto of
More informationGeneral Method for Calculating Chemical Equilibrium Composition
AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt
More informationAPPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS
APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS Pweł Pędzch Jerzy Blcerz Wrsw Uversty of Techology Fculty of Geodesy d Crtogrphy Astrct Usully to pproto of
More informationLeast Squares Method For Solving Integral Equations With Multiple Time Lags
Eg. & Tech. Jourl, Vol.8, No., Lest Squres Method For Solvg Itegrl Equtos Wth Multple Tme Lgs Dr.Suh N. Sheh*, Hyt Adel Al* & Hl Mohmmed Ysee* Receved o:6//9 Accepted o:// Astrct The m purpose of ths work
More informationPubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS
PubH 745: REGRESSION ANALSIS REGRESSION IN MATRIX TERMS A mtr s dspl of umbers or umercl quttes ld out rectgulr rr of rows d colums. The rr, or two-w tble of umbers, could be rectgulr or squre could be
More informationMATRIX AND VECTOR NORMS
Numercl lyss for Egeers Germ Jord Uversty MTRIX ND VECTOR NORMS vector orm s mesure of the mgtude of vector. Smlrly, mtr orm s mesure of the mgtude of mtr. For sgle comoet etty such s ordry umers, the
More informationA Technique for Constructing Odd-order Magic Squares Using Basic Latin Squares
Itertol Jourl of Scetfc d Reserch Publctos, Volume, Issue, My 0 ISSN 0- A Techque for Costructg Odd-order Mgc Squres Usg Bsc Lt Squres Tomb I. Deprtmet of Mthemtcs, Mpur Uversty, Imphl, Mpur (INDIA) tombrom@gml.com
More informationFORMULAE FOR FINITE DIFFERENCE APPROXIMATIONS, QUADRATURE AND LINEAR MULTISTEP METHODS
Jourl o Mtemtcl Sceces: Advces d Alctos Volume Number - Pges -9 FRMULAE FR FINITE DIFFERENCE APPRXIMATINS QUADRATURE AND LINEAR MULTISTEP METHDS RAMESH KUMAR MUTHUMALAI Dertmet o Mtemtcs D G Vsv College
More informationINTRODUCTION ( ) 1. Errors
INTRODUCTION Numercl lyss volves the study, developmet d lyss of lgorthms for obtg umercl solutos to vrous mthemtcl problems. Frequetly umercl lyss s clled the mthemtcs of scetfc computg. Numercl lyss
More informationA METHOD FOR THE RAPID NUMERICAL CALCULATION OF PARTIAL SUMS OF GENERALIZED HARMONICAL SERIES WITH PRESCRIBED ACCURACY
UPB c Bull, eres D, Vol 8, No, 00 A METHOD FOR THE RAPD NUMERAL ALULATON OF PARTAL UM OF GENERALZED HARMONAL ERE WTH PRERBED AURAY BERBENTE e roue o etodă ouă etru clculul rd l suelor rţle le serlor roce
More informationNumerical Differentiation and Integration
Numerl Deretto d Itegrto Overvew Numerl Deretto Newto-Cotes Itegrto Formuls Trpezodl rule Smpso s Rules Guss Qudrture Cheyshev s ormul Numerl Deretto Forwrd te dvded deree Bkwrd te dvded deree Ceter te
More informationChapter 4: Distributions
Chpter 4: Dstrbutos Prerequste: Chpter 4. The Algebr of Expecttos d Vrces I ths secto we wll mke use of the followg symbols: s rdom vrble b s rdom vrble c s costt vector md s costt mtrx, d F m s costt
More information4. Runge-Kutta Formula For Differential Equations. A. Euler Formula B. Runge-Kutta Formula C. An Example for Fourth-Order Runge-Kutta Formula
NCTU Deprme o Elecrcl d Compuer Egeerg Seor Course By Pro. Yo-Pg Ce. Ruge-Ku Formul For Derel Equos A. Euler Formul B. Ruge-Ku Formul C. A Emple or Four-Order Ruge-Ku Formul
More informationPOWERS OF COMPLEX PERSYMMETRIC ANTI-TRIDIAGONAL MATRICES WITH CONSTANT ANTI-DIAGONALS
IRRS 9 y 04 wwwrppresscom/volumes/vol9issue/irrs_9 05pdf OWERS OF COLE ERSERIC I-RIIGOL RICES WIH COS I-IGOLS Wg usu * Q e Wg Hbo & ue College of Scece versty of Shgh for Scece d echology Shgh Ch 00093
More informationEXACT STOCHASTIC ANALYSIS OFCOUPLED BENDING-TORSION BEAMS WITH IN-SPAN SUPPORTS AND MASSES
Meccc de Mterl e delle Strutture Vol. VI 06), o., pp. 5-3 ISS: 035-679X Dprtmeto d Igeger Cvle, Ambetle, Aerospzle, De Mterl DICAM EXACT STOCHASTIC AALYSIS OFCOUPLED BEDIG-TORSIO BEAMS WITH I-SPA SUPPORTS
More informationPGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation
PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad
More informationOutline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III
Itrodcto to Nmercl Alyss Mrc, 9 Nmercl Metods or PDEs Lrry Cretto Meccl Egeerg 5B Semr Egeerg Alyss Mrc, 9 Otle Revew mdterm soltos Revew bsc mterl o mercl clcls Expressos or dervtves, error d error order
More information4. Runge-Kutta Formula For Differential Equations
NCTU Deprme o Elecrcl d Compuer Egeerg 5 Sprg Course by Pro. Yo-Pg Ce. Ruge-Ku Formul For Derel Equos To solve e derel equos umerclly e mos useul ormul s clled Ruge-Ku ormul
More informationLinear Algebra Concepts
Ler Algebr Cocepts Nuo Vscocelos (Ke Kreutz-Delgdo) UCSD Vector spces Defto: vector spce s set H where ddto d sclr multplcto re defed d stsf: ) +( + ) = (+ )+ 5) H 2) + = + H 6) = 3) H, + = 7) ( ) = (
More information10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n
0. Sere I th ecto, we wll troduce ere tht wll be dcug for the ret of th chpter. Wht ere? If we dd ll term of equece, we get whch clled fte ere ( or jut ere) d deoted, for hort, by the ymbol or Doe t mke
More informationComputer Programming
Computer Progrmmg I progrmmg, t s ot eough to be vetve d geous. Oe lso eeds to be dscpled d cotrolled order ot be become etgled oe's ow completes. Hrl D. Mlls, Forwrd to Progrmmg Proverbs b Her F. Ledgrd
More information12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions
HK Km Slghtly moded //9 /8/6 Frstly wrtte t Mrch 5 Itertve Methods er Systems: Guss-Sedel Noler Systems Cse Study: Chemcl Rectos Itertve or ppromte methods or systems o equtos cosst o guessg vlue d the
More informationRendering Equation. Linear equation Spatial homogeneous Both ray tracing and radiosity can be considered special case of this general eq.
Rederg quto Ler equto Sptl homogeeous oth ry trcg d rdosty c be cosdered specl cse of ths geerl eq. Relty ctul photogrph Rdosty Mus Rdosty Rederg quls the dfferece or error mge http://www.grphcs.corell.edu/ole/box/compre.html
More informationAdvanced Algorithmic Problem Solving Le 3 Arithmetic. Fredrik Heintz Dept of Computer and Information Science Linköping University
Advced Algorthmc Prolem Solvg Le Arthmetc Fredrk Hetz Dept of Computer d Iformto Scece Lköpg Uversty Overvew Arthmetc Iteger multplcto Krtsu s lgorthm Multplcto of polyomls Fst Fourer Trsform Systems of
More informationAn Alternative Method to Find the Solution of Zero One Integer Linear Fractional Programming Problem with the Help of -Matrix
Itertol Jourl of Scetfc d Reserch Pulctos, Volue 3, Issue 6, Jue 3 ISSN 5-353 A Altertve Method to Fd the Soluto of Zero Oe Iteger Ler Frctol Progrg Prole wth the Help of -Mtr VSeeregsy *, DrKJeyr ** *
More informationМИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ. Численные методы. Учебно-методическое пособие
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ Нижегородский государственный университет им. Н.И. Лобачевского Численные методы К.А.Баркалов Учебно-методическое пособие Рекомендовано методической
More informationSome Unbiased Classes of Estimators of Finite Population Mean
Itertol Jourl O Mtemtcs Ad ttstcs Iveto (IJMI) E-IN: 3 4767 P-IN: 3-4759 Www.Ijms.Org Volume Issue 09 etember. 04 PP-3-37 ome Ubsed lsses o Estmtors o Fte Poulto Me Prvee Kumr Msr d s Bus. Dertmet o ttstcs,
More informationOn Solution of Min-Max Composition Fuzzy Relational Equation
U-Sl Scece Jourl Vol.4()7 O Soluto of M-Mx Coposto Fuzzy eltol Equto N.M. N* Dte of cceptce /5/7 Abstrct I ths pper, M-Mx coposto fuzzy relto equto re studed. hs study s geerlzto of the works of Ohsto
More information( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.
4.3 - Modal Aalyss Physcal coordates are ot always the easest to work Egevectors provde a coveet trasformato to modal coordates Modal coordates are lear combato of physcal coordates Say we have physcal
More informationFREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM
Joural of Appled Matematcs ad Computatoal Mecacs 04, 3(4), 7-34 FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Ata Cekot, Stasław Kukla Isttute of Matematcs, Czestocowa Uversty of Tecology Częstocowa,
More informationSystems of second order ordinary differential equations
Ffth order dgolly mplct Ruge-Kutt Nystrom geerl method solvg secod Order IVPs Fudzh Isml Astrct A dgolly mplct Ruge-Kutt-Nystróm Geerl (SDIRKNG) method of ffth order wth explct frst stge for the tegrto
More informationQuantization in Dynamic Smarandache Multi-Space
Quatzato Dyamc Smaradache Mult-Space Fu Yuhua Cha Offshore Ol Research Ceter, Beg, 7, Cha (E-mal: fuyh@cooc.com.c ) Abstract: Dscussg the applcatos of Dyamc Smaradache Mult-Space (DSMS) Theory. Supposg
More informationNumerical Analysis Formulae Booklet
Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx
More informationChapter Unary Matrix Operations
Chpter 04.04 Ury trx Opertos After redg ths chpter, you should be ble to:. kow wht ury opertos mes,. fd the trspose of squre mtrx d t s reltoshp to symmetrc mtrces,. fd the trce of mtrx, d 4. fd the ermt
More informationMethod for Hydrodynamic Coupling of Concentric Cylindrical Shells and Beams
ethod for Hydrodymc Couplg of Cocetrc Cyldrcl Shells d Bems tthew D. Syder Westghouse Electrc Compy, LLC Abstrct Dymc lyses of pressurzed wter rector vessels, core support structures, d terls requre ccurte
More informationIntroduction to local (nonparametric) density estimation. methods
Itroducto to local (oparametrc) desty estmato methods A slecture by Yu Lu for ECE 66 Sprg 014 1. Itroducto Ths slecture troduces two local desty estmato methods whch are Parze desty estmato ad k-earest
More informationVibration Analysis Of a Self-Excited Elastic Beam
Vbrto lyss Of Self-Excted Elstc Bem M.. Brró-Mez Deprtmeto de Mterles Uversdd utóom Metropolt zcpotzlco v. S Pblo 80, ol. Reyos-Tmulps.P. 000, Méxco, D.F., MÉXIO bmm@correo.zc.um.mx BSTRT The vbrto behvor
More informationRegression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers
Regresso By Jugl Klt Bsed o Chpter 7 of Chpr d Cle, Numercl Methods for Egeers Regresso Descrbes techques to ft curves (curve fttg) to dscrete dt to obt termedte estmtes. There re two geerl pproches two
More information6.6 Moments and Centers of Mass
th 8 www.tetodre.co 6.6 oets d Ceters of ss Our ojectve here s to fd the pot P o whch th plte of gve shpe lces horzotll. Ths pot s clled the ceter of ss ( or ceter of grvt ) of the plte.. We frst cosder
More informationINTEGRATION TECHNIQUES FOR TWO DIMENSIONAL DOMAINS
IJRET: Itertol Jourl of Reserch Egeerg d Techolog eissn: 9-6 pissn: -78 INTEGRATION TECHNIQUES FOR TWO DIMENSIONAL DOMAINS Logh Peruml Lecturer, Fcult of Egeerg d Techolog, Multmed Uverst, Mlcc, Mls Astrct
More informationChapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II
CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh
More informationNumerical Analysis Topic 4: Least Squares Curve Fitting
Numerl Alss Top 4: Lest Squres Curve Fttg Red Chpter 7 of the tetook Alss_Numerk Motvto Gve set of epermetl dt: 3 5. 5.9 6.3 The reltoshp etwee d m ot e ler. Fd futo f tht est ft the dt 3 Alss_Numerk Motvto
More informationSUM PROPERTIES FOR THE K-LUCAS NUMBERS WITH ARITHMETIC INDEXES
Avlble ole t http://sc.org J. Mth. Comput. Sc. 4 (04) No. 05-7 ISSN: 97-507 SUM PROPERTIES OR THE K-UCAS NUMBERS WITH ARITHMETIC INDEXES BIJENDRA SINGH POOJA BHADOURIA AND OMPRAKASH SIKHWA * School of
More informationCubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem
Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs
More informationThe definite Riemann integral
Roberto s Notes o Itegrl Clculus Chpter 4: Defte tegrls d the FTC Secto 4 The defte Rem tegrl Wht you eed to kow lredy: How to ppromte the re uder curve by usg Rem sums. Wht you c ler here: How to use
More informationAsymptotic Formulas Composite Numbers II
Iteratoal Matematcal Forum, Vol. 8, 203, o. 34, 65-662 HIKARI Ltd, www.m-kar.com ttp://d.do.org/0.2988/mf.203.3854 Asymptotc Formulas Composte Numbers II Rafael Jakmczuk Dvsó Matemátca, Uversdad Nacoal
More informationLecture 3-4 Solutions of System of Linear Equations
Lecture - Solutos of System of Ler Equtos Numerc Ler Alger Revew of vectorsd mtrces System of Ler Equtos Guss Elmto (drect solver) LU Decomposto Guss-Sedel method (tertve solver) VECTORS,,, colum vector
More informationDISCRETE TIME MODELS OF FORWARD CONTRACTS INSURANCE
G Tstsshvl DSCRETE TME MODELS OF FORWARD CONTRACTS NSURANCE (Vol) 008 September DSCRETE TME MODELS OF FORWARD CONTRACTS NSURANCE GSh Tstsshvl e-ml: gurm@mdvoru 69004 Vldvosto Rdo str 7 sttute for Appled
More informationMathematics HL and further mathematics HL formula booklet
Dplom Progrmme Mthemtcs HL d further mthemtcs HL formul boolet For use durg the course d the emtos Frst emtos 04 Publshed Jue 0 Itertol Bcclurete Orgzto 0 5048 Mthemtcs HL d further mthemtcs formul boolet
More informationProceedings of the ASME 2011 International Mechanical Engineering Congress & Exposition IMECE2011 November 11-17, 2011, Denver, Colorado, USA
Proceedgs of the ASME Itertol Mechcl Egeerg ogress & Exosto IMEE ovember -7 Dever olordo USA IMEE-6869 ALORITHM OF LAUHED VEHILES MASS ALULATIO AT THE EARLY STAE OF DESII IMEE-6869 A. AYUTDIOVA Dertmet
More informationMathematics HL and further mathematics HL formula booklet
Dplom Progrmme Mthemtcs HL d further mthemtcs HL formul boolet For use durg the course d the emtos Frst emtos 04 Edted 05 (verso ) Itertol Bcclurete Orgzto 0 5048 Cotets Pror lerg Core 3 Topc : Algebr
More informationAn Analysis of Queueing Model with a Multiple System Governed by a Quasi-Birth Death Process and PH Service and PH Repair
Itertol Jourl of Scece d Reserch (IJSR) ISSN (Ole): 39-764 Idex Copercus Vlue (3): 64 Impct Fctor (4): 56 A Alyss of Queueg Model wth Multple System Govered y Qus-Brth Deth Process d PH Servce d PH Repr
More informationChapter 1. Introduction. Fundamental Concepts. Introduction. Historical background. Historical background. Fundamental Concepts
Chpter udmetl Cocepts Lecture Notes Dr Mohd Afed Uverst Mlys Perls N67 te lemet Alyss Itroducto A or sometmes referred s M, hs ecome powerful tool for umercl soluto of wde rge of egeerg prolems A s computtol
More informationSolutions Manual for Polymer Science and Technology Third Edition
Solutos ul for Polymer Scece d Techology Thrd Edto Joel R. Fred Uer Sddle Rver, NJ Bosto Idols S Frcsco New York Toroto otrel Lodo uch Prs drd Cetow Sydey Tokyo Sgore exco Cty Ths text s ssocted wth Fred/Polymer
More information