Different General Algorithms for Solving Poisson Equation
|
|
- Mariah Gilmore
- 6 years ago
- Views:
Transcription
1 Dffeent Geneal Algotms fo Solvng Posson Equaton Me Yn SUMMARY Te objectve of ts tess s to dscuss te applcaton of dffeent geneal algotms to te soluton of Posson Equaton subject to Dclet bounday condton on a squae doman: u π u Γ 0 sn( π* xsn( π* y wee Ω x, y. I adopt two nds of gds, one wt a coase mes and one wt a fne mes. Accodng to dffeent needs, I dscuss bot te advantages and dsadvantages of tese algotms. Some of te deas mgt be novel and a few mpovements ae offeed. *Undelyng ( J denotes te numbe of unnowns, wee J equals 50 (J and 00 (J, espectvely. INTRODUCTION Fnte Dffeence Metod (FDM s a pmay numecal metod fo solvng Posson Equatons. As electonc dgtal computes ae only capable of andlng fnte data and opeatons, any numecal metod equng te use of computes must fst be dscetzed. Te steps ae as follows: Mes patton s employed and fnte gd nodes ae constucted on a contnuous egon. Dffeentaton opeatos ae dscetzed to tansfe te ognal poblem to fnte sets of lnea equatons. THE ALGORITHM u f (x y Consde Posson Equaton u Γ α(x y wee Ω x, y. Adopt a squae gd wt x and y mes lengts and daw a bunc of paallel lnes to te x * axes. Te ntesectons (, j ae denoted node (, j. y j* Employ second-ode cental dffeence along x and y dectons to substtute fo u xx and u yy, we ave at fve-pont fnte dffeence fomula: u + j u j + u j u j+ u j + u + u Fo smplfcaton, mae j f ( j u and f stand fo mes functons, wee u (x y j u j, f (x y f f (x y. Te above Posson Equaton s ten ewtten as u f. j j j Page
2 Fom Taylo Expanson, we now tat te cutoff eo of s O(. Hee s a setc map fo llustaton: At last, we can set up te appoxmate sets of lnea equatons: Wee L B I I B I O I B I s a (M- OOOOOOO KKKK I B 4 4 B 0 4 s also a (M- matx ( M /. OOOOOOO KKKK 4 L u f matx, I s a (M- dentty matx and As L s a uge spase matx and tee ae at most fve nonzeo elements on evey lne, t can be seen as a band matx wt alf band wdt of N and band wdt of N+. Restcted by compute memoy, t can not be nput by odnay means and MATLAB command "spase" s adopted. Undelyng s te constucton pocess of 50*50 mes: J50; /J; aspase(:(j-^,:(j-^,4*ones(,(j-^,(j-^,(j-^; aspase(j:j^-*j+,:j^-3*j+,-ones(,j^-3*j+,(j-^,(j-^; fo :J^-*J s(-; end fo :J- s(*(j-0; end %sub-dagonal zeo elements ae consdeed wen establsng s a3spase(:j^-*j+,:j^-*j,s,(j-^,(j-^; Aa+a+a'+a3+a3'; Now we tae a loo at ow to solve te sets of lnea equatons wt coeffcent matx A. Undelyng s te man pocess of Steepest Descent (SD Metod: wle (nom>tol Page
3 p A * ; alpa (nom^/(p'*; %fx alpa x x + alpa * ; %fom new teate b - A * x; tete+; ys(tenom; nomnom(; end And te man pocess of Conjugate Gadent (CG Metod: wle (nom>tol q A * p; alpa (nom^/(p'*q; %fx alpa x x + alpa * p; %fom new teate - alpa * q; beta (nom(^/(nom^; %fx beta ; p + beta * p; tete+; yc(tenom; nomnom(; end Wen A s mobd ( K >>, convegence s evolvng vey slowly. To tacle ts poblem, pecondton metod s ntoduced to educe A s condton numbe. Te smplest way s to coose a dagonal matx (Jacoban Pecondton, oweve t mgt not always be effectve. Mspase(:n,:n,0.5*ones(,n,n,n; zm*; pz; wle (nom>tol q A * p; alpa (z'*/(p'*q; %fx alpa x x + alpa * p; %fom new teate - alpa * q; zm*; beta (z'*/(z'*; %fx beta zz; ; p z + beta * p; tete+; yp(tenom; nomnom(; end Tee-dagonal matx M seems to be a nce coce, wc can be solved by Tomas Algotm. adag(dag(a; Page 3
4 adag(dag(a,-,-; Ma+a+a'; %tee-dagonal matx M Ax d b c a b c A O OOOO c n KKKKK a n b n d d d d d M M n n Pefom LU Decomposton of A, u b a l u u b lc.3...n Ax d s tus tansfomed nto Ly d and Ux y. ( u 0 y y d d l y.3...n x x n y u y n n cx u + n.n... It wll be even bette to coose M as te pecondton matx fo SSOR, wc wll educe A s condton numbe to ts squae oot. Especally n te case of w, symmetc GS teaton wll aceve good esults. atl(a; aspase(:n,:n,0.5*ones(,n,n,n; Sa*a; MS*S'; Two-gd metod s a vey popula paallel algotm tese days, wc seves to educe te amount of computaton needed to pocess te data. M u M v M L M u M f M d M Mv M d M Man Pncple: u + L wee v M s te modfe. Messmoot teaton s embedded n-between. (Te oveall sequence of teaton on fne mes, estcton, teaton on te coase mes, ntepolaton and ten fute teaton on te fne mes epesents one cycle of teaton. Geneal teaton fomula s adapted to ts poblem and computaton s fute educed. (J- lnes n te font and bac ae dscussed sepaately and ote lnes ae vewed as a wole. Besdes, sub-dagonal nonzeo elements ae also taen nto consdeaton. educe late computaton Page 4
5 compute (J- lnes n te font ote lnes ae vewed as a wole compute (J- lnes n te bac Pogammng used fo estcton and ntepolaton: Vectos ae tansfomed nto nodes on te gd, wc ae fst aayed ozontally and ten vetcally. Hee s a setc map fo llustaton: Defne estcton opeato I H and ntepolaton opeato I H : I H 6 4 fom fne mes to coase mes Page 5
6 I H 4 fom coase mes to fne mes (dstbuton accodng to wegt coeffcent 4 It s qute stagtfowad to desgn I H : fo ::J- fo j::j- spacng s enlaged u0.065*(a(+,j-+*a(+,j+a(+,j++*a(,j-+4*a(,j+*a(,j++a(-,j-+*a (-,j+a(-,j+; evaluate accodng to wegt coeffcent +; end end Wen dealng wt te desgn of I H, te stuaton becomes complcated and bounday and ntenal ponts need to be consdeed sepaately. Moeove, dffeent elatve postons between fne mes nodes and coase mes nodes wll lead to vaed wegt coeffcents. coase mes nodes fne mes nodes on dffeent ows and dffeent columns fne mes nodes on dffeent ows but same columns fne mes nodes on dffeent columns but same ows fou vetces fne mes nodes on same ows (fst and last column fne mes nodes on dffeent ows (fst and last column fne mes nodes on same columns (fst and last ow fne mes nodes on dffeent columns (fst and last ow mes nodes ae tansfomed nto vectos STRUCTURE PROGRAMMING Data Flow Gap Input Evaluate x Output nom>tol Steepest Descent (SD Metod: Intal value of x0 b Ax α ( (A Eo x x + α Page 6
7 Conjugate Gadent (CG Metod: Eo Intal value of x0 p b Ax α ( (p Ap p α Ap x x + α p + β p β ( ( Pecondton Conjugate Gadent (PCG Metod: Eo Intal value of x0 p b Ax z M α ( (p Ap p + β p z M α Ap x x + α p β (z (z z Page 7
8 Two-gd Metod: Relaxaton on fne mes Modfcaton on coase mes Restcton Eo Modfcaton on fne mes Intepolaton Relaxaton on coase mes Module Stuctue Man M-Fle PCG Algotm Tomas Algotm Two-gd Metod Intepolaton Restcton Sedel PCG Tomas RESULTS AND FURTHER DISCUSSION Caactestcs of te dffeence fomat: u j f j ( x y j D u j α j ( x y j D Suppose u s a functon defned on te doman D j D, we ave: u j 0, ( x y j D ten u j 0, ( x y j D ten max u max u j D D max u max u j D D j j α as a unque soluton, max u j max u j + max u j wee α s te x-sde lengt of D. D D D Page 8
9 Fve-pont fnte dffeence metod conveges wt max u u(x y O( + Numecal esults: (Intal value of x0ones(n, Iteaton Steps Steepest Descent J50 D j CG PCG Revsed PCG j Revsed PCG Two-gd Metod Iteaton Steps Steepest Descent J00 CG PCG Revsed PCG Revsed PCG Two-gd Metod Fo Steepest Descent (SD Metod, Solvng of ϕ ( x (Ax x (b x. Stat fom an abtay A * x b s equvalent to fndng te mnmum pont x and fnd te decton along wc ϕ(x ϕ(x deceases geatest,.e. ϕ( x b Ax. If 0, fnd α tat maes ϕ (x + α least n value. Ts epesents one cycle of teaton. It can be fute poved tat lm x > * * λ λ n * x A b x x A x x A λ + λ n 0.5 A (Au u u Wee λ and λ n ae A s maxmum and mnmum egenvalues, espectvely. Convegence s extemely slow wen λ >> λ n, so ts metod s not of pactcal use. CG s a vaatonal metod. It fst appeaed n te 950 s and ganed populaty wt te development of PCG n te 980 s and s now a pmay metod fo solvng uge spase matces.. In essence, t calls fo ϕ(x mn ϕ(y. Teoetcally, ts metod can fnd te ( y span{p...p } soluton afte at most n teaton steps. Howeve, because of te exstence of oundng eo, ts goal mgt not be aceved so easly. We ave * K * x x A x x A wee K K cond (A. + Wen A s mobd ( K >>, convegence s evolvng vey slowly and pecondton metod s ntoduced to educe A s condton numbe. T Consde Colesy Decomposton M LL and A M N, wee M s symmetc postve defnte and N s as nea zeo as possble. Te smplest way s to coose a dagonal matx (Jacoban Pecondton and my attaced PCG pogamme adopts ts dea. My Revsed PCG cooses tee-dagonal matx M as te Page 9
10 pecondton matx and my Revsed PCG cooses M as te pecondton matx fo SSOR T ( M SS 0.5 wees (D LD. Fo Two-gd Metod: Egenfuncton ϕ x sn(m π* x sn(m π* x m ( And egenvalue λ 0.5*(cos(m π * + cos(m π*. m Undelyng x denotes teaton steps and y denotes eo. J50 Pecson: 0.0 J50 Pecson: 0.00 J00 Pecson: 0.0 J00 Pecson: 0.00 J50 Pecson: 0.0 J50 Pecson: 0.00 Page 0
11 J00 Pecson: 0.0 J00 Pecson: 0.00 Pecson: 0.0 Pecson: 0.00 Multple Gd (MG Metod s nealy always te best and as been wdely used snce 977 fo solvng ellptcal bounday condton poblems. Hee s a setc map fo one cycle of MG teaton. Page
12 Undelyng s a compason table. Dmenson Metod Jacoban Gauss-Sedel SOR MG Two O N 4 O N 4 O N 3 O N Tee O N 5 O N 5 O N 4 O N 3 A LIST OF RELATED PROGRAMS spasepaametes 50*50 mes spasepaametes 00*00 mes steepestdescent Steepest Descent Metod cg CG Metod pcg PCG Metod (Jacoban Pecondton pcgg PCG Metod (tee-dagonal pecondton matx, solved by embedded system pogamme evsedpcg PCG Metod (tee-dagonal pecondton matx, solved by Tomas Algotm Tomas Tomas Algotm evsedpcg PCG Metod (pecondton matx fo SSOR gd Two-gd Metod Page
13 estct Restcton Algotm ntepolate Intepolaton Algotm sedel Messmoot Algotm REFERENCES Yang Huazong, Wang Hu Numecal Metods & C Language Scence Pess 996 Da Jazun, Qu Janxan Numecal Metods fo Dffeental Equatons Souteast Unvesty Pess 00 L Rongua, Feng Guocen Numecal Metods fo Dffeental Equatons People s Educaton Pess 980 L Qngyang, Guan Z, Ba Fengsan Numecal Computaton Tsngua Unvesty Pess 000 Tang Huamn, Hu Janwe Numecal Metods fo Dffeental Equatons Nana Unvesty Pess 990 Wang Moan MATLAB & Scentfc Computaton Publsng House of Electoncs Industy 004 Lu Jnfu, Guan Z Numecal Metods fo Patal Dffeental Equatons Tsngua Unvesty Pess 987 Guan Z, Lu Jnfu Intoducton to Numecal Analyss Hge Educaton Pess 998 Page 3
Slide 1. Quantum Mechanics: the Practice
Slde Quantum Mecancs: te Pactce Slde Remnde: Electons As Waves Wavelengt momentum = Planck? λ p = = 6.6 x 0-34 J s Te wave s an exctaton a vbaton: We need to know te ampltude of te exctaton at evey pont
More informationUNIT10 PLANE OF REGRESSION
UIT0 PLAE OF REGRESSIO Plane of Regesson Stuctue 0. Intoducton Ojectves 0. Yule s otaton 0. Plane of Regesson fo thee Vaales 0.4 Popetes of Resduals 0.5 Vaance of the Resduals 0.6 Summay 0.7 Solutons /
More informationSet of square-integrable function 2 L : function space F
Set of squae-ntegable functon L : functon space F Motvaton: In ou pevous dscussons we have seen that fo fee patcles wave equatons (Helmholt o Schödnge) can be expessed n tems of egenvalue equatons. H E,
More informationCOORDINATE SYSTEMS, COORDINATE TRANSFORMS, AND APPLICATIONS
Dola Bagaoo 0 COORDINTE SYSTEMS COORDINTE TRNSFORMS ND PPLICTIONS I. INTRODUCTION Smmet coce of coodnate sstem. In solvng Pscs poblems one cooses a coodnate sstem tat fts te poblem at and.e. a coodnate
More informationCorrespondence Analysis & Related Methods
Coespondence Analyss & Related Methods Ineta contbutons n weghted PCA PCA s a method of data vsualzaton whch epesents the tue postons of ponts n a map whch comes closest to all the ponts, closest n sense
More informationMultistage Median Ranked Set Sampling for Estimating the Population Median
Jounal of Mathematcs and Statstcs 3 (: 58-64 007 ISSN 549-3644 007 Scence Publcatons Multstage Medan Ranked Set Samplng fo Estmatng the Populaton Medan Abdul Azz Jeman Ame Al-Oma and Kamaulzaman Ibahm
More informationMachine Learning. Spectral Clustering. Lecture 23, April 14, Reading: Eric Xing 1
Machne Leanng -7/5 7/5-78, 78, Spng 8 Spectal Clusteng Ec Xng Lectue 3, pl 4, 8 Readng: Ec Xng Data Clusteng wo dffeent ctea Compactness, e.g., k-means, mxtue models Connectvty, e.g., spectal clusteng
More informationV. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. "Flux": = da i. "Force": = -Â g a ik k = X i. Â J i X i (7.
Themodynamcs and Knetcs of Solds 71 V. Pncples of Ievesble Themodynamcs 5. Onsage s Teatment s = S - S 0 = s( a 1, a 2,...) a n = A g - A n (7.6) Equlbum themodynamcs detemnes the paametes of an equlbum
More informationThermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering
Themodynamcs of solds 4. Statstcal themodynamcs and the 3 d law Kwangheon Pak Kyung Hee Unvesty Depatment of Nuclea Engneeng 4.1. Intoducton to statstcal themodynamcs Classcal themodynamcs Statstcal themodynamcs
More information( ) α is determined to be a solution of the one-dimensional minimization problem: = 2. min = 2
Homewo (Patal Solton) Posted on Mach, 999 MEAM 5 Deental Eqaton Methods n Mechancs. Sole the ollowng mat eqaton A b by () Steepest Descent Method and/o Pecondtoned SD Method Snce the coecent mat A s symmetc,
More informationOptimization Methods: Linear Programming- Revised Simplex Method. Module 3 Lecture Notes 5. Revised Simplex Method, Duality and Sensitivity analysis
Optmzaton Meods: Lnea Pogammng- Revsed Smple Meod Module Lectue Notes Revsed Smple Meod, Dualty and Senstvty analyss Intoducton In e pevous class, e smple meod was dscussed whee e smple tableau at each
More informationPhysics 11b Lecture #2. Electric Field Electric Flux Gauss s Law
Physcs 11b Lectue # Electc Feld Electc Flux Gauss s Law What We Dd Last Tme Electc chage = How object esponds to electc foce Comes n postve and negatve flavos Conseved Electc foce Coulomb s Law F Same
More informationCOMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS
ultscence - XXX. mcocd Intenatonal ultdscplnay Scentfc Confeence Unvesty of skolc Hungay - pl 06 ISBN 978-963-358-3- COPLEENTRY ENERGY ETHOD FOR CURVED COPOSITE BES Ákos József Lengyel István Ecsed ssstant
More informationIntegral Vector Operations and Related Theorems Applications in Mechanics and E&M
Dola Bagayoko (0) Integal Vecto Opeatons and elated Theoems Applcatons n Mechancs and E&M Ι Basc Defnton Please efe to you calculus evewed below. Ι, ΙΙ, andιιι notes and textbooks fo detals on the concepts
More informationScalars and Vectors Scalar
Scalas and ectos Scala A phscal quantt that s completel chaacteed b a eal numbe (o b ts numecal value) s called a scala. In othe wods a scala possesses onl a magntude. Mass denst volume tempeatue tme eneg
More information8 Baire Category Theorem and Uniform Boundedness
8 Bae Categoy Theoem and Unfom Boundedness Pncple 8.1 Bae s Categoy Theoem Valdty of many esults n analyss depends on the completeness popety. Ths popety addesses the nadequacy of the system of atonal
More informationSome Approximate Analytical Steady-State Solutions for Cylindrical Fin
Some Appoxmate Analytcal Steady-State Solutons fo Cylndcal Fn ANITA BRUVERE ANDRIS BUIIS Insttute of Mathematcs and Compute Scence Unvesty of Latva Rana ulv 9 Rga LV459 LATVIA Astact: - In ths pape we
More information19 The Born-Oppenheimer Approximation
9 The Bon-Oppenheme Appoxmaton The full nonelatvstc Hamltonan fo a molecule s gven by (n a.u.) Ĥ = A M A A A, Z A + A + >j j (883) Lets ewte the Hamltonan to emphasze the goal as Ĥ = + A A A, >j j M A
More informationA. Thicknesses and Densities
10 Lab0 The Eath s Shells A. Thcknesses and Denstes Any theoy of the nteo of the Eath must be consstent wth the fact that ts aggegate densty s 5.5 g/cm (ecall we calculated ths densty last tme). In othe
More informationGenerating Functions, Weighted and Non-Weighted Sums for Powers of Second-Order Recurrence Sequences
Geneatng Functons, Weghted and Non-Weghted Sums fo Powes of Second-Ode Recuence Sequences Pantelmon Stăncă Aubun Unvesty Montgomey, Depatment of Mathematcs Montgomey, AL 3614-403, USA e-mal: stanca@studel.aum.edu
More information9/12/2013. Microelectronics Circuit Analysis and Design. Modes of Operation. Cross Section of Integrated Circuit npn Transistor
Mcoelectoncs Ccut Analyss and Desgn Donald A. Neamen Chapte 5 The pola Juncton Tanssto In ths chapte, we wll: Dscuss the physcal stuctue and opeaton of the bpola juncton tanssto. Undestand the dc analyss
More informationON THE FRESNEL SINE INTEGRAL AND THE CONVOLUTION
IJMMS 3:37, 37 333 PII. S16117131151 http://jmms.hndaw.com Hndaw Publshng Cop. ON THE FRESNEL SINE INTEGRAL AND THE CONVOLUTION ADEM KILIÇMAN Receved 19 Novembe and n evsed fom 7 Mach 3 The Fesnel sne
More informationq-bernstein polynomials and Bézier curves
Jounal of Computatonal and Appled Mathematcs 151 (2003) 1-12 q-bensten polynomals and Béze cuves Hall Ouç a, and Geoge M. Phllps b a Depatment of Mathematcs, Dokuz Eylül Unvesty Fen Edebyat Fakültes, Tınaztepe
More informationSOME NEW SELF-DUAL [96, 48, 16] CODES WITH AN AUTOMORPHISM OF ORDER 15. KEYWORDS: automorphisms, construction, self-dual codes
Факултет по математика и информатика, том ХVІ С, 014 SOME NEW SELF-DUAL [96, 48, 16] CODES WITH AN AUTOMORPHISM OF ORDER 15 NIKOLAY I. YANKOV ABSTRACT: A new method fo constuctng bnay self-dual codes wth
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationP 365. r r r )...(1 365
SCIENCE WORLD JOURNAL VOL (NO4) 008 www.scecncewoldounal.og ISSN 597-64 SHORT COMMUNICATION ANALYSING THE APPROXIMATION MODEL TO BIRTHDAY PROBLEM *CHOJI, D.N. & DEME, A.C. Depatment of Mathematcs Unvesty
More informationA Study about One-Dimensional Steady State. Heat Transfer in Cylindrical and. Spherical Coordinates
Appled Mathematcal Scences, Vol. 7, 03, no. 5, 67-633 HIKARI Ltd, www.m-hka.com http://dx.do.og/0.988/ams.03.38448 A Study about One-Dmensonal Steady State Heat ansfe n ylndcal and Sphecal oodnates Lesson
More informationThe Unique Solution of Stochastic Differential Equations With. Independent Coefficients. Dietrich Ryter.
The Unque Soluton of Stochastc Dffeental Equatons Wth Independent Coeffcents Detch Ryte RyteDM@gawnet.ch Mdatweg 3 CH-4500 Solothun Swtzeland Phone +4132 621 13 07 SDE s must be solved n the ant-itô sense
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationMechanics Physics 151
Mechancs Physcs 151 Lectue 18 Hamltonan Equatons of Moton (Chapte 8) What s Ahead We ae statng Hamltonan fomalsm Hamltonan equaton Today and 11/6 Canoncal tansfomaton 1/3, 1/5, 1/10 Close lnk to non-elatvstc
More informationWave Equations. Michael Fowler, University of Virginia
Wave Equatons Mcael Fowle, Unvesty of Vgna Potons and Electons We ave seen tat electons and potons beave n a vey smla fason bot exbt dffacton effects, as n te double slt expement, bot ave patcle lke o
More informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More information24-2: Electric Potential Energy. 24-1: What is physics
D. Iyad SAADEDDIN Chapte 4: Electc Potental Electc potental Enegy and Electc potental Calculatng the E-potental fom E-feld fo dffeent chage dstbutons Calculatng the E-feld fom E-potental Potental of a
More informationAn Approach to Inverse Fuzzy Arithmetic
An Appoach to Invese Fuzzy Athmetc Mchael Hanss Insttute A of Mechancs, Unvesty of Stuttgat Stuttgat, Gemany mhanss@mechaun-stuttgatde Abstact A novel appoach of nvese fuzzy athmetc s ntoduced to successfully
More informationANALYSIS OF AXIAL LOADED PILE IN MULTILAYERED SOIL USING NODAL EXACT FINITE ELEMENT MODEL
Intenatonal Jounal of GEOMATE, Apl, 8 Vol. 4, Issue 44, pp. -7 Geotec., Const. Mat. & Env., DOI: https://do.og/.66/8.44.785 ISS: 86-98 (Pnt), 86-99 (Onlne), Japan AAYSIS OF AXIA OADED PIE I MUTIAYERED
More informationTest 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o?
Test 1 phy 0 1. a) What s the pupose of measuement? b) Wte all fou condtons, whch must be satsfed by a scala poduct. (Use dffeent symbols to dstngush opeatons on ectos fom opeatons on numbes.) c) What
More informationAnalytical and Numerical Solutions for a Rotating Annular Disk of Variable Thickness
Appled Mathematcs 00 43-438 do:0.436/am.00.5057 Publshed Onlne Novembe 00 (http://www.scrp.og/jounal/am) Analytcal and Numecal Solutons fo a Rotatng Annula Ds of Vaable Thcness Abstact Ashaf M. Zenou Daoud
More informationDistinct 8-QAM+ Perfect Arrays Fanxin Zeng 1, a, Zhenyu Zhang 2,1, b, Linjie Qian 1, c
nd Intenatonal Confeence on Electcal Compute Engneeng and Electoncs (ICECEE 15) Dstnct 8-QAM+ Pefect Aays Fanxn Zeng 1 a Zhenyu Zhang 1 b Lnje Qan 1 c 1 Chongqng Key Laboatoy of Emegency Communcaton Chongqng
More information3. A Review of Some Existing AW (BT, CT) Algorithms
3. A Revew of Some Exstng AW (BT, CT) Algothms In ths secton, some typcal ant-wndp algothms wll be descbed. As the soltons fo bmpless and condtoned tansfe ae smla to those fo ant-wndp, the pesented algothms
More informationMachine Learning 4771
Machne Leanng 4771 Instucto: Tony Jebaa Topc 6 Revew: Suppot Vecto Machnes Pmal & Dual Soluton Non-sepaable SVMs Kenels SVM Demo Revew: SVM Suppot vecto machnes ae (n the smplest case) lnea classfes that
More informationChapter 23: Electric Potential
Chapte 23: Electc Potental Electc Potental Enegy It tuns out (won t show ths) that the tostatc foce, qq 1 2 F ˆ = k, s consevatve. 2 Recall, fo any consevatve foce, t s always possble to wte the wok done
More informationAPPLICATIONS OF SEMIGENERALIZED -CLOSED SETS
Intenatonal Jounal of Mathematcal Engneeng Scence ISSN : 22776982 Volume Issue 4 (Apl 202) http://www.mes.com/ https://stes.google.com/ste/mesounal/ APPLICATIONS OF SEMIGENERALIZED CLOSED SETS G.SHANMUGAM,
More informationPattern Analyses (EOF Analysis) Introduction Definition of EOFs Estimation of EOFs Inference Rotated EOFs
Patten Analyses (EOF Analyss) Intoducton Defnton of EOFs Estmaton of EOFs Infeence Rotated EOFs . Patten Analyses Intoducton: What s t about? Patten analyses ae technques used to dentfy pattens of the
More informationRigid Bodies: Equivalent Systems of Forces
Engneeng Statcs, ENGR 2301 Chapte 3 Rgd Bodes: Equvalent Sstems of oces Intoducton Teatment of a bod as a sngle patcle s not alwas possble. In geneal, the se of the bod and the specfc ponts of applcaton
More informationCSJM University Class: B.Sc.-II Sub:Physics Paper-II Title: Electromagnetics Unit-1: Electrostatics Lecture: 1 to 4
CSJM Unvesty Class: B.Sc.-II Sub:Physcs Pape-II Ttle: Electomagnetcs Unt-: Electostatcs Lectue: to 4 Electostatcs: It deals the study of behavo of statc o statonay Chages. Electc Chage: It s popety by
More informationPhysics 202, Lecture 2. Announcements
Physcs 202, Lectue 2 Today s Topcs Announcements Electc Felds Moe on the Electc Foce (Coulomb s Law The Electc Feld Moton of Chaged Patcles n an Electc Feld Announcements Homewok Assgnment #1: WebAssgn
More informationMEI Structured Mathematics. Module Summary Sheets. Numerical Methods (Version B reference to new book)
MEI Matematics in Education and Industy MEI Stuctued Matematics Module Summay Seets (Vesion B efeence to new book) Topic : Appoximations Topic : Te solution of equations Topic : Numeical integation Topic
More informationGroupoid and Topological Quotient Group
lobal Jounal of Pue and Appled Mathematcs SSN 0973-768 Volume 3 Numbe 7 07 pp 373-39 Reseach nda Publcatons http://wwwpublcatoncom oupod and Topolocal Quotent oup Mohammad Qasm Manna Depatment of Mathematcs
More informationPart V: Velocity and Acceleration Analysis of Mechanisms
Pat V: Velocty an Acceleaton Analyss of Mechansms Ths secton wll evew the most common an cuently pactce methos fo completng the knematcs analyss of mechansms; escbng moton though velocty an acceleaton.
More informationThe Forming Theory and the NC Machining for The Rotary Burs with the Spectral Edge Distribution
oden Appled Scence The Fomn Theoy and the NC achnn fo The Rotay us wth the Spectal Ede Dstbuton Huan Lu Depatment of echancal Enneen, Zhejan Unvesty of Scence and Technoloy Hanzhou, c.y. chan, 310023,
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationPhysics 2A Chapter 11 - Universal Gravitation Fall 2017
Physcs A Chapte - Unvesal Gavtaton Fall 07 hese notes ae ve pages. A quck summay: he text boxes n the notes contan the esults that wll compse the toolbox o Chapte. hee ae thee sectons: the law o gavtaton,
More informationIf there are k binding constraints at x then re-label these constraints so that they are the first k constraints.
Mathematcal Foundatons -1- Constaned Optmzaton Constaned Optmzaton Ma{ f ( ) X} whee X {, h ( ), 1,, m} Necessay condtons fo to be a soluton to ths mamzaton poblem Mathematcally, f ag Ma{ f ( ) X}, then
More informationPHY126 Summer Session I, 2008
PHY6 Summe Sesson I, 8 Most of nfomaton s avalable at: http://nngoup.phscs.sunsb.edu/~chak/phy6-8 ncludng the sllabus and lectue sldes. Read sllabus and watch fo mpotant announcements. Homewok assgnment
More informationKhintchine-Type Inequalities and Their Applications in Optimization
Khntchne-Type Inequaltes and The Applcatons n Optmzaton Anthony Man-Cho So Depatment of Systems Engneeng & Engneeng Management The Chnese Unvesty of Hong Kong ISDS-Kolloquum Unvestaet Wen 29 June 2009
More information4 SingularValue Decomposition (SVD)
/6/00 Z:\ jeh\self\boo Kannan\Jan-5-00\4 SVD 4 SngulaValue Decomposton (SVD) Chapte 4 Pat SVD he sngula value decomposton of a matx s the factozaton of nto the poduct of thee matces = UDV whee the columns
More informationPHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle
1 PHYS 705: Classcal Mechancs Devaton of Lagange Equatons fom D Alembet s Pncple 2 D Alembet s Pncple Followng a smla agument fo the vtual dsplacement to be consstent wth constants,.e, (no vtual wok fo
More informationThe Greatest Deviation Correlation Coefficient and its Geometrical Interpretation
By Rudy A. Gdeon The Unvesty of Montana The Geatest Devaton Coelaton Coeffcent and ts Geometcal Intepetaton The Geatest Devaton Coelaton Coeffcent (GDCC) was ntoduced by Gdeon and Hollste (987). The GDCC
More informationEfficiency of the principal component Liu-type estimator in logistic
Effcency of the pncpal component Lu-type estmato n logstc egesson model Jbo Wu and Yasn Asa 2 School of Mathematcs and Fnance, Chongqng Unvesty of Ats and Scences, Chongqng, Chna 2 Depatment of Mathematcs-Compute
More informationChapter Fifiteen. Surfaces Revisited
Chapte Ffteen ufaces Revsted 15.1 Vecto Descpton of ufaces We look now at the vey specal case of functons : D R 3, whee D R s a nce subset of the plane. We suppose s a nce functon. As the pont ( s, t)
More informationDetection and Estimation Theory
ESE 54 Detecton and Etmaton Theoy Joeph A. O Sullvan Samuel C. Sach Pofeo Electonc Sytem and Sgnal Reeach Laboatoy Electcal and Sytem Engneeng Wahngton Unvety 411 Jolley Hall 314-935-4173 (Lnda anwe) jao@wutl.edu
More informationShuai Dong. Isaac Newton. Gottfried Leibniz
Computatonal pyscs Sua Dong Isaac Newton Gottred Lebnz Numercal calculus poston dervatve ntegral v velocty dervatve ntegral a acceleraton Numercal calculus Numercal derentaton Numercal ntegraton Roots
More informationLie Subalgebras and Invariant Solutions to the Equation of Fluid Flows in Toroidal Field. Lang Xia
Le Subalgebas and Invaant Solutons to the Equaton of Flud Flows n Toodal Feld Lang a Emal: langxaog@gmalcom Abstact: Patal dffeental equatons (PDEs), patculaly coupled PDE systems, ae dffcult to solve
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationEnergy in Closed Systems
Enegy n Closed Systems Anamta Palt palt.anamta@gmal.com Abstact The wtng ndcates a beakdown of the classcal laws. We consde consevaton of enegy wth a many body system n elaton to the nvese squae law and
More informationRanks of quotients, remainders and p-adic digits of matrices
axv:1401.6667v2 [math.nt] 31 Jan 2014 Ranks of quotents, emandes and p-adc dgts of matces Mustafa Elshekh Andy Novocn Mak Gesbecht Abstact Fo a pme p and a matx A Z n n, wte A as A = p(a quo p)+ (A em
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 17. a ij x (k) b i. a ij x (k+1) (D + L)x (k+1) = b Ux (k)
STAT 309: MATHEMATICAL COMPUTATIONS I FALL 08 LECTURE 7. sor method remnder: n coordnatewse form, Jacob method s = [ b a x (k) a and Gauss Sedel method s = [ b a = = remnder: n matrx form, Jacob method
More informationEngineering Mechanics. Force resultants, Torques, Scalar Products, Equivalent Force systems
Engneeng echancs oce esultants, Toques, Scala oducts, Equvalent oce sstems Tata cgaw-hll Companes, 008 Resultant of Two oces foce: acton of one bod on anothe; chaacteed b ts pont of applcaton, magntude,
More information4. Linear systems of equations. In matrix form: Given: matrix A and vector b Solve: Ax = b. Sup = least upper bound
4. Lnea systes of eqatons a a a a 3 3 a a a a 3 3 a a a a 3 3 In at fo: a a a3 a a a a3 a a a a3 a Defnton ( vecto no): On a vecto space V, a vecto no s a fncton fo V to e set of non-negatve eal nes at
More informationConjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationMultivariate Ratio Estimator of the Population Total under Stratified Random Sampling
Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng
More informationEffect of a Frequency Perturbation in a Chain of Syntonized Transparent Clocks
Effect of a Fequency Petubaton n a Chan of Syntonzed anspaent Clocs Geoffey M. Gane SAMSUNG Electoncs (Consultant) EEE 80. AVB G 007.03.0 gmgane@comcast.net : Outlne ntoducton ansfe functon fo a chan of
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationMultigrid Methods and Applications in CFD
Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac
More informationN = N t ; t 0. N is the number of claims paid by the
Iulan MICEA, Ph Mhaela COVIG, Ph Canddate epatment of Mathematcs The Buchaest Academy of Economc Studes an CECHIN-CISTA Uncedt Tac Bank, Lugoj SOME APPOXIMATIONS USE IN THE ISK POCESS OF INSUANCE COMPANY
More information1. A body will remain in a state of rest, or of uniform motion in a straight line unless it
Pncples of Dnamcs: Newton's Laws of moton. : Foce Analss 1. A bod wll eman n a state of est, o of unfom moton n a staght lne unless t s acted b etenal foces to change ts state.. The ate of change of momentum
More informationContact, information, consultations
ontact, nfomaton, consultatons hemsty A Bldg; oom 07 phone: 058-347-769 cellula: 664 66 97 E-mal: wojtek_c@pg.gda.pl Offce hous: Fday, 9-0 a.m. A quote of the week (o camel of the week): hee s no expedence
More informationChapter IV Vector and Tensor Analysis IV.2 Vector and Tensor Analysis September 29,
hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe 9, 08 47 hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe 9, 08 48 I. ETOR AND TENOR ANALYI I... Tenso functon th Let A
More informationChapter 3: Theory of Modular Arithmetic 38
Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences
More informationThe intrinsic sense of stochastic differential equations
The ntnsc sense of stochastc dffeental equatons Detch Ryte RyteDM@gawnet.ch Mdatweg 3 CH-4500 Solothun Swtzeland Phone +413 61 13 07 A change of the vaables (establshng a nomal fom of the fowad and bacwad
More informationReading Assignment. Problem Description for Homework #9. Read Chapters 29 and 30.
Reading Assignment Read Chaptes 29 and 30. Poblem Desciption fo Homewok #9 In this homewok, you will solve the inhomogeneous Laplace s equation to calculate the electic scala potential that exists between
More informationThermoelastic Problem of a Long Annular Multilayered Cylinder
Wold Jounal of Mechancs, 3, 3, 6- http://dx.do.og/.436/w.3.35a Publshed Onlne August 3 (http://www.scp.og/ounal/w) Theoelastc Poble of a Long Annula Multlayeed Cylnde Y Hsen Wu *, Kuo-Chang Jane Depatent
More informationClosed-loop adaptive optics using a CMOS image quality metric sensor
Closed-loop adaptve optcs usng a CMOS mage qualty metc senso Chueh Tng, Mchael Gles, Adtya Rayankula, and Pual Futh Klpsch School of Electcal and Compute Engneeng ew Mexco State Unvesty Las Cuces, ew Mexco
More informationImpact of Polarimetric Dimensionality of Forest Parameter Estimation by Means of Polarimetric SAR interferometry
Impact of Polametc Dmensonalty of Foest Paamete Estmaton by Means of Polametc SAR ntefeomety Jun Su Km, Seung-Kuk Lee, Konstantnos Papathanassou, and Iena Hajnsek Geman Aeospace Cente Mcowaves and Rada
More informationTian Zheng Department of Statistics Columbia University
Haplotype Tansmsson Assocaton (HTA) An "Impotance" Measue fo Selectng Genetc Makes Tan Zheng Depatment of Statstcs Columba Unvesty Ths s a jont wok wth Pofesso Shaw-Hwa Lo n the Depatment of Statstcs at
More informationOn Pfaff s solution of the Pfaff problem
Zur Pfaff scen Lösung des Pfaff scen Probles Mat. Ann. 7 (880) 53-530. On Pfaff s soluton of te Pfaff proble By A. MAYER n Lepzg Translated by D. H. Delpenc Te way tat Pfaff adopted for te ntegraton of
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
More informationOn Maneuvering Target Tracking with Online Observed Colored Glint Noise Parameter Estimation
Wold Academy of Scence, Engneeng and Technology 6 7 On Maneuveng Taget Tacng wth Onlne Obseved Coloed Glnt Nose Paamete Estmaton M. A. Masnad-Sha, and S. A. Banan Abstact In ths pape a compehensve algothm
More informationImplementation in the ANSYS Finite Element Code of the Electric Vector Potential T-Ω,Ω Formulation
Implementaton n the ANSYS Fnte Element Code of the Electc Vecto Potental T-Ω,Ω Fomulaton Peto Teston Dpatmento d Ingegnea Elettca ed Elettonca, Unvestà d Cagla Pazza d Am, 0923 Cagla Pegogo Sonato Dpatmento
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More information3.1 Electrostatic Potential Energy and Potential Difference
3. lectostatc Potental negy and Potental Dffeence RMMR fom mechancs: - The potental enegy can be defned fo a system only f consevatve foces act between ts consttuents. - Consevatve foces may depend only
More informationINTRODUCTION. consider the statements : I there exists x X. f x, such that. II there exists y Y. such that g y
INRODUCION hs dssetaton s the eadng of efeences [1], [] and [3]. Faas lemma s one of the theoems of the altenatve. hese theoems chaacteze the optmalt condtons of seveal mnmzaton poblems. It s nown that
More informationOptimal System for Warm Standby Components in the Presence of Standby Switching Failures, Two Types of Failures and General Repair Time
Intenatonal Jounal of ompute Applcatons (5 ) Volume 44 No, Apl Optmal System fo Wam Standby omponents n the esence of Standby Swtchng Falues, Two Types of Falues and Geneal Repa Tme Mohamed Salah EL-Shebeny
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
More informationSmoothed Particle Hydrodynamics
Smooted Patcle Hydodynamcs Applcaton Example Alan Hec Noembe 9, 010 Fluds Nae-Stokes equatons Smooted Patcle Hydodynamcs Smooted Patcle Hydodynamcs Noembe 010 Fluds Lquds, e.g. wate Gasses, e.g. a Plasmas
More informationChapter IV Vector and Tensor Analysis IV.2 Vector and Tensor Analysis September 23,
hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe, 07 47 hapte I ecto and Tenso Analyss I. ecto and Tenso Analyss eptembe, 07 48 I. ETOR AND TENOR ANALYI I... Tenso functon th Let A n n
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationIntroduction. - The Second Lyapunov Method. - The First Lyapunov Method
Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationAccurate Evaluation Schemes for Triangular Domain Integrals
IOSR Jounal of Mechancal and Cvl Engneeng (IOSRJMCE) ISSN : 78-68 Volume, Issue 6 (Sep-Oct 0), PP 38-5 Accuate Evaluaton Schemes fo Tangula Doman Integals Fazana Hussan, M. S. Kam, (Depatment of Mathematc
More information