A finite difference method for heat equation in the unbounded domain

Size: px
Start display at page:

Download "A finite difference method for heat equation in the unbounded domain"

Transcription

1 Internatona Conerence on Advanced ectronc Scence and Technoogy (AST 6) A nte derence method or heat equaton n the unbounded doman a Quan Zheng and Xn Zhao Coege o Scence North Chna nversty o Technoogy Bejng 44 Chna Abstract. The numerca method o the one-dmensona non-homogeneous heat equaton on an unbounded doman s consdered. Two exact artca boundary condtons are apped on two artca boundares to mt the orgna probem onto a bounded computatona doman. Then the nte derence method s deveoped by usng the method o the reducton o order or the contro equaton and artca boundary condtons. It s proved that the nte derence scheme s stabe and convergent wth the order n space and order 3/ n tme under an energy norm. A non-homogeneous numerca exampe demonstrates the uncondtona stabty and the accuracy o the agorthm. Keywords: non-homogeneous heat equaton; unbounded doman; artca boundary condton; nte derence method; error estmate. Introducton Heat equaton arses rom many eds or exampes the heat transer ud dynamcs astrophyscs nance or other areas o apped mathematcs. For the probem on unbounded domans Han and Huang () n [] and [] derved an exact artca boundary condton obtaned the correspondng numerca souton or heat equatons by FM and FDM n one-dmensona and two-dmensona cases showed the eectveness wth numerca exampes but had no error estmate. Wu and Sun (4) n [3] apped the artca boundary condton n [] to nvestgate the numerca souton o heat equaton on the sem-nnte nterva provng the uncondtona stabty and convergence wth order n space and the order 3/ n tme under an energy norm. In ths paper we consder the oowng Cauchy probem o one-dmensona non-homogeneous heat equaton: ut σ u xx = ( x t ) < x < + < t T u= ( x ) φ ( x < x < + u ( x t ) when x + t T. a Correspondng author : zhengq@ncut.edu.com 6. The authors - Pubshed by Atants Press 387 () () (3)

2 where σ s duson coecent determnng heat conducton rate densty and heat capacty o the matera. The nta term φ ( x) and the source term () xt have compact supports n space. Introducng artca boundares such that { ( )} [ Γ = {( xt ) x = x t T} and Γ r = {( xt ) x = xr t T } su p p xt x xr ] [ T] and su p p { ( x )} [ x x r ]. φ Accordng to the work o Han and Huang n [] and Wu and Sun n [3] we can reduce the nta probem to an nta-boundary probem on a bounded doman: u v σ = t x ( xt) x < x < x r < t T (4) u v = x x < x < x r (5) ux ( ) = φ( x x < x < x r (6) vx ( t) = t ux ( λ) dλ < t T λ t λ vx ( r t) = t ux ( r λ) dλ < t T. λ t λ (7) (8) The nte derence scheme and correspondng error estmaton We dscretze the above probem as oows. Dvde the doman [ x x r ] [ T ] nto an M N mesh wth the spata sze h = ( xr x )/ M n x drecton and the tme step sze τ = T / N n t drecton respectvey where M and N are postve ntegers. Grd ponts ( x t n ) are dened by x = x + h or M and tn = nτ or n N. Moreover denote n be the numerca souton o uxt () at ( x t n ) and ntroduce the notatons: n n n n n n n = ( + ( ( = + δ x = h n n n n n δx = ( ( + + δt = h τ n n n n n M n = ( + ( h ( ). = + = = n Dene the grd uncton: u = ux ( tn M n. The notatons o v are smar. Appyng nte derence method or (4)-(8 we have 388

3 δ tu σδ xv M n v δ xu M n δ tu σδ xv M n v δ xu M n u = φ( x (9) () () () (3) 3 3 where r ( τ + h r h s τ s τ M c are constants ndependent o τ and h. Omttng the sma error terms n (8)-(3 we construct a nte derence scheme: n and c and δ t σδ xv = M n V δ x = M n k V [ ( ) = a a a ] n k k k = k V [ ( ) M = a M a a ] n k k k M = (4) (5) (6) (7) φ( x M. = (8) Theorem. The nte derence scheme (4)-(8) s equvaent to the oowng scheme (9)-(3) wthout ntermedate varabe: ( δt + δt ) σδ x = (9) δ σ k t + a ak ak δx = h k = [ ( ( ) ) ] () δ σ k t + a M an k an km x M + δ = M h M k = [ ( ( ) ) ] () 389

4 = φ( x () where am = ( m+ m m = 3 L. (3) τ We can obtan the stabty and the convergence o the nte derence scheme usng Gronwa nequaty as oows by usng the methods n [4] and [5]. Theorem. Let max ( xt ) 3 c n 3be a postve constant.{ } x x xr t T be the Souton o the derence scheme (4)-(8 uxt () be the souton o probem (4)-(8). Then n τ exp( )( c 3 TL n T / + = L τ τ τ 3 n n T c T σt u exp( ) [ + ( σ + ) L]( τ + h ) n= L [T/ τ] τ τ π where c = max{ c c c3}. 3 Numerca exampe We ntroduce the oowng dentons: L -error: n n h τ = max u M L -order: hτ og ( h τ / h τ L -error: M n n h τ = h [ u ] = L -order: hτ og ( h τ / h τ xampe.consder the oowng nonhomogeneous heat equaton: ( ) xt ut σuxx = σ 4 x e < x < + < t T x ux ( ) = e < x< + uxt ( ) when x + t T. x t whch has the exact souton uxt () = e.the support o φ( x) and () xt s compact on the computatona doman [ x x ] = [ 5 5]. r Tabe. Convergence w.r.t. τ o xampe or σ = and h = τ 3 4. T= T= M N L -error order L -error order N L -error order L -error order e e e e e e e e e e e e e e e e

5 Tabe. Convergence w.r.t. h o xampe or σ = and τ = h 4 3. T= T= M N L -error order L -error order N L -error order L -error order e e e e e e e e e e e e e e e e-4.3 We can see n Tabe that the order o convergence agrees wth the theoretca order 3/ n tme 3 3 durng the derent tme perods snce Theorem gves O( τ + h ) = O( τ ) or h = τ 3 4. Tabe ustrates the theoretca order o convergence snce Theorem gves O( τ + h) = Oh ( ) or τ = h 4 3. In addton the resuts when T= and T= show the uncondtona stabty o the proposed agorthm. 3.5 x Reatve L-error.5.5 r=.5 r= r= t Fgure. Reatve L -error o xampe durng t under derent mesh ratos or σ = and M =. We can aso see n Fgure that the numerca resut s uncondtonay stabe durng the computatona tme perod regardess o the mesh rato στ r =. h 4 Concusons The numerca souton to the non-homogenous heat equaton n the unbounded doman s consdered. Smar to the numerca methods o non-homogenous heat equaton n the haunbounded doman n [] and the nvoved homogenous heat equaton n unbounded doman n [5] we convert the orgna probem to an nta-boundary vaue probem by usng two exact artca boundary condtons. Then the nte derence scheme s obtaned by the method o reducton o order. The stabty and the correspondng error estmate o (9)-(3) are gven n Theorem. A non- 39

6 homogenous exampe demonstrates the uncondtona stabty and the accuracy o the proposed method. Acknowedgement Ths work s supported by Natona Natura Scence Foundaton o Chna. (479). Reerences. H.-D. Han Z.-Y. Huang. A cass o artca boundary condtons or heat equaton n unbounded domans Comput. Math. App. 43 () H.-D. Han Z.-Y. Huang xact and approxmatng boundary condtons or the paraboc probemson unbounded domans Comput. Math. App. 44 () X.-N. Wu Z.-Z. Sun Convergence o derence scheme or heat equaton n unbounded domansusng artca boundary condtons App. Numer. Math. 5 (4) H.-D. Han X.-N. Wu Artca Boundary Method Bejng: Tsnghua nversty Press/SprngerPress. 5. Q. Zheng L. Fan X.-Z. L An artca boundary method or Burgers' equaton n the unboundeddoman CMS: Comput. Mode. ng. Sc. (4)

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory Advanced Scence and Technoogy Letters Vo.83 (ISA 205), pp.60-65 http://dx.do.org/0.4257/ast.205.83.2 Research on Compex etworks Contro Based on Fuzzy Integra Sdng Theory Dongsheng Yang, Bngqng L, 2, He

More information

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory Proceedngs of the 2009 IEEE Internatona Conference on Systems Man and Cybernetcs San Antono TX USA - October 2009 Mutspectra Remote Sensng Image Cassfcaton Agorthm Based on Rough Set Theory Yng Wang Xaoyun

More information

The line method combined with spectral chebyshev for space-time fractional diffusion equation

The line method combined with spectral chebyshev for space-time fractional diffusion equation Apped and Computatona Mathematcs 014; 3(6): 330-336 Pubshed onne December 31, 014 (http://www.scencepubshnggroup.com/j/acm) do: 10.1164/j.acm.0140306.17 ISS: 3-5605 (Prnt); ISS: 3-5613 (Onne) The ne method

More information

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2 Note 2 Lng fong L Contents Ken Gordon Equaton. Probabty nterpretaton......................................2 Soutons to Ken-Gordon Equaton............................... 2 2 Drac Equaton 3 2. Probabty nterpretaton.....................................

More information

Discriminating Fuzzy Preference Relations Based on Heuristic Possibilistic Clustering

Discriminating Fuzzy Preference Relations Based on Heuristic Possibilistic Clustering Mutcrtera Orderng and ankng: Parta Orders, Ambgutes and Apped Issues Jan W. Owsńsk and aner Brüggemann, Edtors Dscrmnatng Fuzzy Preerence eatons Based on Heurstc Possbstc Custerng Dmtr A. Vattchenn Unted

More information

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014 Lecture 16 8/4/14 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 214. Real Vapors and Fugacty Henry s Law accounts or the propertes o extremely dlute soluton. s shown n Fgure

More information

Networked Cooperative Distributed Model Predictive Control Based on State Observer

Networked Cooperative Distributed Model Predictive Control Based on State Observer Apped Mathematcs, 6, 7, 48-64 ubshed Onne June 6 n ScRes. http://www.scrp.org/journa/am http://dx.do.org/.436/am.6.73 Networed Cooperatve Dstrbuted Mode redctve Contro Based on State Observer Ba Su, Yanan

More information

Approximate merging of a pair of BeÂzier curves

Approximate merging of a pair of BeÂzier curves COMPUTER-AIDED DESIGN Computer-Aded Desgn 33 (1) 15±136 www.esever.com/ocate/cad Approxmate mergng of a par of BeÂzer curves Sh-Mn Hu a,b, *, Rou-Feng Tong c, Tao Ju a,b, Ja-Guang Sun a,b a Natona CAD

More information

MARKOV CHAIN AND HIDDEN MARKOV MODEL

MARKOV CHAIN AND HIDDEN MARKOV MODEL MARKOV CHAIN AND HIDDEN MARKOV MODEL JIAN ZHANG JIANZHAN@STAT.PURDUE.EDU Markov chan and hdden Markov mode are probaby the smpest modes whch can be used to mode sequenta data,.e. data sampes whch are not

More information

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks Shengyang Sun, Changyou Chen, Lawrence Carn Suppementary Matera: Learnng Structured Weght Uncertanty n Bayesan Neura Networks Shengyang Sun Changyou Chen Lawrence Carn Tsnghua Unversty Duke Unversty Duke

More information

Journal of Universal Computer Science, vol. 1, no. 7 (1995), submitted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Springer Pub. Co.

Journal of Universal Computer Science, vol. 1, no. 7 (1995), submitted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Springer Pub. Co. Journal of Unversal Computer Scence, vol. 1, no. 7 (1995), 469-483 submtted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Sprnger Pub. Co. Round-o error propagaton n the soluton of the heat equaton by

More information

Nested case-control and case-cohort studies

Nested case-control and case-cohort studies Outne: Nested case-contro and case-cohort studes Ørnuf Borgan Department of Mathematcs Unversty of Oso NORBIS course Unversty of Oso 4-8 December 217 1 Radaton and breast cancer data Nested case contro

More information

Neural network-based athletics performance prediction optimization model applied research

Neural network-based athletics performance prediction optimization model applied research Avaabe onne www.jocpr.com Journa of Chemca and Pharmaceutca Research, 04, 6(6):8-5 Research Artce ISSN : 0975-784 CODEN(USA) : JCPRC5 Neura networ-based athetcs performance predcton optmzaton mode apped

More information

Xin Li Department of Information Systems, College of Business, City University of Hong Kong, Hong Kong, CHINA

Xin Li Department of Information Systems, College of Business, City University of Hong Kong, Hong Kong, CHINA RESEARCH ARTICLE MOELING FIXE OS BETTING FOR FUTURE EVENT PREICTION Weyun Chen eartment of Educatona Informaton Technoogy, Facuty of Educaton, East Chna Norma Unversty, Shangha, CHINA {weyun.chen@qq.com}

More information

Complex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen

Complex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen omplex Varables hapter 8 Integraton n the omplex Plane March, Lecturer: Shh-Yuan hen Except where otherwse noted, content s lcensed under a BY-N-SA. TW Lcense. ontents ontour ntegrals auchy-goursat theorem

More information

n-step cycle inequalities: facets for continuous n-mixing set and strong cuts for multi-module capacitated lot-sizing problem

n-step cycle inequalities: facets for continuous n-mixing set and strong cuts for multi-module capacitated lot-sizing problem n-step cyce nequates: facets for contnuous n-mxng set and strong cuts for mut-modue capactated ot-szng probem Mansh Bansa and Kavash Kanfar Department of Industra and Systems Engneerng, Texas A&M Unversty,

More information

Computational Biology Lecture 8: Substitution matrices Saad Mneimneh

Computational Biology Lecture 8: Substitution matrices Saad Mneimneh Computatonal Bology Lecture 8: Substtuton matrces Saad Mnemneh As we have ntroduced last tme, smple scorng schemes lke + or a match, - or a msmatch and -2 or a gap are not justable bologcally, especally

More information

Optimal Guaranteed Cost Control of Linear Uncertain Systems with Input Constraints

Optimal Guaranteed Cost Control of Linear Uncertain Systems with Input Constraints Internatona Journa Optma of Contro, Guaranteed Automaton, Cost Contro and Systems, of Lnear vo Uncertan 3, no Systems 3, pp 397-4, wth Input September Constrants 5 397 Optma Guaranteed Cost Contro of Lnear

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

More information

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak Thermodynamcs II Department o Chemca ngneerng ro. Km, Jong Hak .5 Fugacty & Fugacty Coecent : ure Speces µ > provdes undamenta crteron or phase equbrum not easy to appy to sove probem Lmtaton o gn (.9

More information

Piecewise Linear Estimate of Attractive Regions for Linear Systems with Saturating Control

Piecewise Linear Estimate of Attractive Regions for Linear Systems with Saturating Control Pecewse Lnear stmate of Attractve Regons for Lnear Systems wth Saturatng Contro Yuzo Ohta Department of Computer and Systems, Kobe Unversty Nada, Kobe 657-8501, Japan ohta@cskobe-uacjp ABSTRACT We show

More information

NON-PARABOLIC INTERFACE MOTION FOR THE 1-D STEFAN PROBLEM Dirichlet Boundary Conditions

NON-PARABOLIC INTERFACE MOTION FOR THE 1-D STEFAN PROBLEM Dirichlet Boundary Conditions Hernandez, E. M., et al.: Non-Parabolc Interace Moton or the -D Stean Problem... THERMAL SCIENCE: Year 07, Vol., No. 6A, pp. 37-336 37 NON-PARABOLIC INTERFACE MOTION FOR THE -D STEFAN PROBLEM Drchlet Boundary

More information

Reliability Sensitivity Algorithm Based on Stratified Importance Sampling Method for Multiple Failure Modes Systems

Reliability Sensitivity Algorithm Based on Stratified Importance Sampling Method for Multiple Failure Modes Systems Chnese Journa o Aeronautcs 3(010) 660-669 Chnese Journa o Aeronautcs www.esever.com/ocate/ca Reabty Senstvty Agorthm Based on Strated Importance Sampng Method or Mutpe aure Modes Systems Zhang eng a, u

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Research Article H Estimates for Discrete-Time Markovian Jump Linear Systems

Research Article H Estimates for Discrete-Time Markovian Jump Linear Systems Mathematca Probems n Engneerng Voume 213 Artce ID 945342 7 pages http://dxdoorg/11155/213/945342 Research Artce H Estmates for Dscrete-Tme Markovan Jump Lnear Systems Marco H Terra 1 Gdson Jesus 2 and

More information

3. Stress-strain relationships of a composite layer

3. Stress-strain relationships of a composite layer OM PO I O U P U N I V I Y O F W N ompostes ourse 8-9 Unversty of wente ng. &ech... tress-stran reatonshps of a composte ayer - Laurent Warnet & emo Aerman.. tress-stran reatonshps of a composte ayer Introducton

More information

A Simple Research of Divisor Graphs

A Simple Research of Divisor Graphs The 29th Workshop on Combnatoral Mathematcs and Computaton Theory A Smple Research o Dvsor Graphs Yu-png Tsao General Educaton Center Chna Unversty o Technology Tape Tawan yp-tsao@cuteedutw Tape Tawan

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

Remark: Positive work is done on an object when the point of application of the force moves in the direction of the force.

Remark: Positive work is done on an object when the point of application of the force moves in the direction of the force. Unt 5 Work and Energy 5. Work and knetc energy 5. Work - energy theore 5.3 Potenta energy 5.4 Tota energy 5.5 Energy dagra o a ass-sprng syste 5.6 A genera study o the potenta energy curve 5. Work and

More information

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by ES 5 (phy 40). a) Wrte the zeroth law o thermodynamcs. b) What s thermal conductvty? c) Identyng all es, draw schematcally a P dagram o the arnot cycle. d) What s the ecency o an engne and what s the coecent

More information

MAXIMUM NORM STABILITY OF DIFFERENCE SCHEMES FOR PARABOLIC EQUATIONS ON OVERSET NONMATCHING SPACE-TIME GRIDS

MAXIMUM NORM STABILITY OF DIFFERENCE SCHEMES FOR PARABOLIC EQUATIONS ON OVERSET NONMATCHING SPACE-TIME GRIDS MATHEMATICS OF COMPUTATION Voume 72 Number 242 Pages 619 656 S 0025-57180201462-X Artce eectroncay pubshed on November 4 2002 MAXIMUM NORM STABILITY OF DIFFERENCE SCHEMES FOR PARABOLIC EQUATIONS ON OVERSET

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

L-Edge Chromatic Number Of A Graph

L-Edge Chromatic Number Of A Graph IJISET - Internatona Journa of Innovatve Scence Engneerng & Technoogy Vo. 3 Issue 3 March 06. ISSN 348 7968 L-Edge Chromatc Number Of A Graph Dr.R.B.Gnana Joth Assocate Professor of Mathematcs V.V.Vannaperuma

More information

6. Stochastic processes (2)

6. Stochastic processes (2) Contents Markov processes Brth-death processes Lect6.ppt S-38.45 - Introducton to Teletraffc Theory Sprng 5 Markov process Consder a contnuous-tme and dscrete-state stochastc process X(t) wth state space

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

6. Stochastic processes (2)

6. Stochastic processes (2) 6. Stochastc processes () Lect6.ppt S-38.45 - Introducton to Teletraffc Theory Sprng 5 6. Stochastc processes () Contents Markov processes Brth-death processes 6. Stochastc processes () Markov process

More information

Delay tomography for large scale networks

Delay tomography for large scale networks Deay tomography for arge scae networks MENG-FU SHIH ALFRED O. HERO III Communcatons and Sgna Processng Laboratory Eectrca Engneerng and Computer Scence Department Unversty of Mchgan, 30 Bea. Ave., Ann

More information

Force = F Piston area = A

Force = F Piston area = A CHAPTER III Ths chapter s an mportant transton between the propertes o pure substances and the most mportant chapter whch s: the rst law o thermodynamcs In ths chapter, we wll ntroduce the notons o heat,

More information

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition) Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes

More information

Endogenous timing in a mixed oligopoly consisting of a single public firm and foreign competitors. Abstract

Endogenous timing in a mixed oligopoly consisting of a single public firm and foreign competitors. Abstract Endogenous tmng n a mxed olgopoly consstng o a sngle publc rm and oregn compettors Yuanzhu Lu Chna Economcs and Management Academy, Central Unversty o Fnance and Economcs Abstract We nvestgate endogenous

More information

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals ECEN 5005 Crystals, Nanocrystals and Devce Applcatons Class 9 Group Theory For Crystals Dee Dagram Radatve Transton Probablty Wgner-Ecart Theorem Selecton Rule Dee Dagram Expermentally determned energy

More information

Predicting Model of Traffic Volume Based on Grey-Markov

Predicting Model of Traffic Volume Based on Grey-Markov Vo. No. Modern Apped Scence Predctng Mode of Traffc Voume Based on Grey-Marov Ynpeng Zhang Zhengzhou Muncpa Engneerng Desgn & Research Insttute Zhengzhou 5005 Chna Abstract Grey-marov forecastng mode of

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

Greyworld White Balancing with Low Computation Cost for On- Board Video Capturing

Greyworld White Balancing with Low Computation Cost for On- Board Video Capturing reyword Whte aancng wth Low Computaton Cost for On- oard Vdeo Capturng Peng Wu Yuxn Zoe) Lu Hewett-Packard Laboratores Hewett-Packard Co. Pao Ato CA 94304 USA Abstract Whte baancng s a process commony

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Chapter 6. Rotations and Tensors

Chapter 6. Rotations and Tensors Vector Spaces n Physcs 8/6/5 Chapter 6. Rotatons and ensors here s a speca knd of near transformaton whch s used to transforms coordnates from one set of axes to another set of axes (wth the same orgn).

More information

Finite Difference Method

Finite Difference Method 7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use

More information

Lowest-Order e + e l + l Processes in Quantum Electrodynamics. Sanha Cheong

Lowest-Order e + e l + l Processes in Quantum Electrodynamics. Sanha Cheong Lowest-Order e + e + Processes n Quantum Eectrodynamcs Sanha Cheong Introducton In ths short paper, we w demonstrate some o the smpest cacuatons n quantum eectrodynamcs (QED), eadng to the owest-order

More information

On a nonlinear compactness lemma in L p (0, T ; B).

On a nonlinear compactness lemma in L p (0, T ; B). On a nonlnear compactness lemma n L p (, T ; B). Emmanuel Matre Laboratore de Matématques et Applcatons Unversté de Haute-Alsace 4, rue des Frères Lumère 6893 Mulouse E.Matre@ua.fr 3t February 22 Abstract

More information

Image Classification Using EM And JE algorithms

Image Classification Using EM And JE algorithms Machne earnng project report Fa, 2 Xaojn Sh, jennfer@soe Image Cassfcaton Usng EM And JE agorthms Xaojn Sh Department of Computer Engneerng, Unversty of Caforna, Santa Cruz, CA, 9564 jennfer@soe.ucsc.edu

More information

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes Numerca Investgaton of Power Tunabty n Two-Secton QD Superumnescent Dodes Matta Rossett Paoo Bardea Ivo Montrosset POLITECNICO DI TORINO DELEN Summary 1. A smpfed mode for QD Super Lumnescent Dodes (SLD)

More information

Analysis of Non-binary Hybrid LDPC Codes

Analysis of Non-binary Hybrid LDPC Codes Anayss of Non-bnary Hybrd LDPC Codes Luce Sassate and Davd Decercq ETIS ENSEA/UCP/CNRS UMR-5 954 Cergy, FRANCE {sassate,decercq}@ensea.fr Abstract Ths paper s egbe for the student paper award. In ths paper,

More information

Associative Memories

Associative Memories Assocatve Memores We consder now modes for unsupervsed earnng probems, caed auto-assocaton probems. Assocaton s the task of mappng patterns to patterns. In an assocatve memory the stmuus of an ncompete

More information

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

CME 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 information

Case Study of Markov Chains Ray-Knight Compactification

Case Study of Markov Chains Ray-Knight Compactification Internatonal Journal of Contemporary Mathematcal Scences Vol. 9, 24, no. 6, 753-76 HIKAI Ltd, www.m-har.com http://dx.do.org/.2988/cms.24.46 Case Study of Marov Chans ay-knght Compactfcaton HaXa Du and

More information

Solution of a nonsymmetric algebraic Riccati equation from a one-dimensional multistate transport model

Solution of a nonsymmetric algebraic Riccati equation from a one-dimensional multistate transport model IMA Journa of Numerca Anayss (2011) 1, 145 1467 do:10.109/manum/drq04 Advance Access pubcaton on May 0, 2011 Souton of a nonsymmetrc agebrac Rccat equaton from a one-dmensona mutstate transport mode TIEXIANG

More information

REAL-TIME IMPACT FORCE IDENTIFICATION OF CFRP LAMINATED PLATES USING SOUND WAVES

REAL-TIME IMPACT FORCE IDENTIFICATION OF CFRP LAMINATED PLATES USING SOUND WAVES 8 TH INTERNATIONAL CONERENCE ON COMPOSITE MATERIALS REAL-TIME IMPACT ORCE IDENTIICATION O CRP LAMINATED PLATES USING SOUND WAVES S. Atobe *, H. Kobayash, N. Hu 3 and H. ukunaga Department of Aerospace

More information

LOW-DENSITY Parity-Check (LDPC) codes have received

LOW-DENSITY Parity-Check (LDPC) codes have received IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 59, NO. 7, JULY 2011 1807 Successve Maxmzaton for Systematc Desgn of Unversay Capacty Approachng Rate-Compatbe Sequences of LDPC Code Ensembes over Bnary-Input

More information

CHAPTER 2. PHOTOABSORPTION AND PHOTOELECTRIC PROCESS IN Si QUANTUM DOTS

CHAPTER 2. PHOTOABSORPTION AND PHOTOELECTRIC PROCESS IN Si QUANTUM DOTS CHAPTER PHOTOABSORPTION AND PHOTOELECTRIC PROCESS IN S QUANTUM DOTS.1 INTRODUCTION At the nanometer sze o a matera comparabe to the characterstc de-broge waveength o the eectron, ts contnuous band structure

More information

3/4 1/2 1/2 3/4 5/4 5/4 3/4 1/2 3/4 3/4 0 1/2 1 3/4 1/2 1 3/4 5/4 1 1/2 3/4 3/4 1/ /2

3/4 1/2 1/2 3/4 5/4 5/4 3/4 1/2 3/4 3/4 0 1/2 1 3/4 1/2 1 3/4 5/4 1 1/2 3/4 3/4 1/ /2 Journal of Computatonal Mathematcs, Vol., No.6, 4, 87{86. THE STRUCTURAL CHARACTERIZATION AND LOCALLY SUPPORTED BASES FOR BIVARIATE SUPER SPLINES ) Zh-qang u (Insttute of Mathematcal Scences, Dalan Unversty

More information

Adaptive and Iterative Least Squares Support Vector Regression Based on Quadratic Renyi Entropy

Adaptive and Iterative Least Squares Support Vector Regression Based on Quadratic Renyi Entropy daptve and Iteratve Least Squares Support Vector Regresson Based on Quadratc Ren Entrop Jngqng Jang, Chu Song, Haan Zhao, Chunguo u,3 and Yanchun Lang Coege of Mathematcs and Computer Scence, Inner Mongoa

More information

Operating conditions of a mine fan under conditions of variable resistance

Operating conditions of a mine fan under conditions of variable resistance Paper No. 11 ISMS 216 Operatng condtons of a mne fan under condtons of varable resstance Zhang Ynghua a, Chen L a, b, Huang Zhan a, *, Gao Yukun a a State Key Laboratory of Hgh-Effcent Mnng and Safety

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

K-Total Product Cordial Labelling of Graphs

K-Total Product Cordial Labelling of Graphs Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Vol. 7, Issue (December ), pp. 78-76 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) K-Total Product Cordal Labellng o Graphs

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty

More information

Journal of Multivariate Analysis

Journal of Multivariate Analysis Journa of Mutvarate Anayss 3 (04) 74 96 Contents sts avaabe at ScenceDrect Journa of Mutvarate Anayss journa homepage: www.esever.com/ocate/jmva Hgh-dmensona sparse MANOVA T. Tony Ca a, Yn Xa b, a Department

More information

Chapter 5 rd Law of Thermodynamics

Chapter 5 rd Law of Thermodynamics Entropy and the nd and 3 rd Chapter 5 rd Law o hermodynamcs homas Engel, hlp Red Objectves Introduce entropy. Derve the condtons or spontanety. Show how S vares wth the macroscopc varables,, and. Chapter

More information

ON THE BEHAVIOR OF THE CONJUGATE-GRADIENT METHOD ON ILL-CONDITIONED PROBLEMS

ON THE BEHAVIOR OF THE CONJUGATE-GRADIENT METHOD ON ILL-CONDITIONED PROBLEMS ON THE BEHAVIOR OF THE CONJUGATE-GRADIENT METHOD ON I-CONDITIONED PROBEM Anders FORGREN Technca Report TRITA-MAT-006-O Department of Mathematcs Roya Insttute of Technoogy January 006 Abstract We study

More information

Optimization of JK Flip Flop Layout with Minimal Average Power of Consumption based on ACOR, Fuzzy-ACOR, GA, and Fuzzy-GA

Optimization of JK Flip Flop Layout with Minimal Average Power of Consumption based on ACOR, Fuzzy-ACOR, GA, and Fuzzy-GA Journa of mathematcs and computer Scence 4 (05) - 5 Optmzaton of JK Fp Fop Layout wth Mnma Average Power of Consumpton based on ACOR, Fuzzy-ACOR, GA, and Fuzzy-GA Farshd Kevanan *,, A Yekta *,, Nasser

More information

Linear Approximation with Regularization and Moving Least Squares

Linear 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 information

Prediction Error of the Multivariate Additive Loss Reserving Method for Dependent Lines of Business

Prediction Error of the Multivariate Additive Loss Reserving Method for Dependent Lines of Business Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness by Mchae Merz and Maro V Wüthrch ABSTRACT Often n non-fe nsurance, cams reserves are the argest poston on the abty

More information

Short-Term Load Forecasting for Electric Power Systems Using the PSO-SVR and FCM Clustering Techniques

Short-Term Load Forecasting for Electric Power Systems Using the PSO-SVR and FCM Clustering Techniques Energes 20, 4, 73-84; do:0.3390/en40073 Artce OPEN ACCESS energes ISSN 996-073 www.mdp.com/journa/energes Short-Term Load Forecastng for Eectrc Power Systems Usng the PSO-SVR and FCM Custerng Technques

More information

SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) M(B) := # ( B Z N)

SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) M(B) := # ( B Z N) SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) S.BOUCKSOM Abstract. The goal of ths note s to present a remarably smple proof, due to Hen, of a result prevously obtaned by Gllet-Soulé,

More information

On the Power Function of the Likelihood Ratio Test for MANOVA

On the Power Function of the Likelihood Ratio Test for MANOVA Journa of Mutvarate Anayss 8, 416 41 (00) do:10.1006/jmva.001.036 On the Power Functon of the Lkehood Rato Test for MANOVA Dua Kumar Bhaumk Unversty of South Aabama and Unversty of Inos at Chcago and Sanat

More information

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

A linear domain decomposition method for partially saturated flow in porous media

A linear domain decomposition method for partially saturated flow in porous media UHasset Computatona Mathematcs Preprnt Seres A near doman decomposton method for partay saturated fow n porous meda Davd Seus, Koondanbha Mtra, Iuu Sorn Pop, Forn Adran Radu and Chrstan Rohde UHasset Computatona

More information

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or

More information

Chapter 7. Potential Energy and Conservation of Energy

Chapter 7. Potential Energy and Conservation of Energy Chapter 7 Potental Energy and Conservaton o Energy 1 Forms o Energy There are many orms o energy, but they can all be put nto two categores Knetc Knetc energy s energy o moton Potental Potental energy

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

More information

Analysis and construction of a family of refinable functions based on generalized Bernstein polynomials

Analysis and construction of a family of refinable functions based on generalized Bernstein polynomials Cheng and Yang Journa of Inequates and Appcatons 6 6:66 DOI.86/s66-6-6- R E S E A R C H Open Access Anayss and constructon of a famy of refnabe functons based on generazed Bernsten poynomas Tng Cheng and

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

More information

Norms, Condition Numbers, Eigenvalues and Eigenvectors

Norms, Condition Numbers, Eigenvalues and Eigenvectors Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b

More information

Econometrics: What's It All About, Alfie?

Econometrics: What's It All About, Alfie? ECON 351* -- Introducton (Page 1) Econometrcs: What's It All About, Ale? Usng sample data on observable varables to learn about economc relatonshps, the unctonal relatonshps among economc varables. Econometrcs

More information

Quantitative Evaluation Method of Each Generation Margin for Power System Planning

Quantitative Evaluation Method of Each Generation Margin for Power System Planning 1 Quanttatve Evauaton Method of Each Generaton Margn for Poer System Pannng SU Su Department of Eectrca Engneerng, Tohoku Unversty, Japan moden25@ececetohokuacjp Abstract Increasng effcency of poer pants

More information

CHAPTER III Neural Networks as Associative Memory

CHAPTER III Neural Networks as Associative Memory CHAPTER III Neural Networs as Assocatve Memory Introducton One of the prmary functons of the bran s assocatve memory. We assocate the faces wth names, letters wth sounds, or we can recognze the people

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

CABLE STRUCTURE WITH LOAD-ADAPTING GEOMETRY

CABLE STRUCTURE WITH LOAD-ADAPTING GEOMETRY Compostes n Constructon Thrd Internatona Conerence Lon, France, Ju 3, CABLE STRUCTURE WITH LOAD-ADAPTING GEOMETRY A.S. Jüch, J.F. Caron and O. Bavere Insttut Naver - Lam, ENPC-LCPC 6 et, avenue Base Pasca

More information

Integral Formula of Minkowski Type and New Characterization of the Wulff Shape

Integral Formula of Minkowski Type and New Characterization of the Wulff Shape Acta athematca Snca, Engsh Seres Apr., 2008, Vo. 24, No. 4, pp. 697 704 Pubshed onne: Apr 5, 2008 DOI: 0.007/s04-007-76-6 Http://www.Actaath.com Acta athematca Snca, Engsh Seres The Edtora Offce of AS

More information

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton

More information

Lecture 3: Shannon s Theorem

Lecture 3: Shannon s Theorem CSE 533: Error-Correctng Codes (Autumn 006 Lecture 3: Shannon s Theorem October 9, 006 Lecturer: Venkatesan Guruswam Scrbe: Wdad Machmouch 1 Communcaton Model The communcaton model we are usng conssts

More information

An Explicit Construction of an Expander Family (Margulis-Gaber-Galil)

An Explicit Construction of an Expander Family (Margulis-Gaber-Galil) An Explct Constructon of an Expander Famly Marguls-Gaber-Gall) Orr Paradse July 4th 08 Prepared for a Theorst's Toolkt, a course taught by Irt Dnur at the Wezmann Insttute of Scence. The purpose of these

More information

A MIN-MAX REGRET ROBUST OPTIMIZATION APPROACH FOR LARGE SCALE FULL FACTORIAL SCENARIO DESIGN OF DATA UNCERTAINTY

A MIN-MAX REGRET ROBUST OPTIMIZATION APPROACH FOR LARGE SCALE FULL FACTORIAL SCENARIO DESIGN OF DATA UNCERTAINTY A MIN-MAX REGRET ROBST OPTIMIZATION APPROACH FOR ARGE SCAE F FACTORIA SCENARIO DESIGN OF DATA NCERTAINTY Travat Assavapokee Department of Industra Engneerng, nversty of Houston, Houston, Texas 7704-4008,

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statstca Mechancs Notes for Lecture 11 I. PRINCIPLES OF QUANTUM STATISTICAL MECHANICS The probem of quantum statstca mechancs s the quantum mechanca treatment of an N-partce system. Suppose the

More information

Labor Share in the Change of Japanese Industrial Structure

Labor Share in the Change of Japanese Industrial Structure Prepared or the 6 th INFORUM Word Conerence hed at European Unersty o eke n North Cyprus, August 3-September 7, 28 abor Share n the Change o Japanese Industra Structure Toshak Hasegaa Proessor o Internatona

More information

Dynamics of a Discrete Predator-Prey System with Beddington-DeAngelis Function Response

Dynamics of a Discrete Predator-Prey System with Beddington-DeAngelis Function Response Apped Mathematcs 3 389-394 http://dxdoorg/46/am46 Pshed Onne Apr (http://wwwscrporg/jorna/am) Dynamcs of a Dscrete Predator-Prey System wth Beddngton-DeAnges Fncton Response Qn Fang Xaopng L * Mey Cao

More information

Numerical integration in more dimensions part 2. Remo Minero

Numerical integration in more dimensions part 2. Remo Minero Numerca ntegraton n more dmensons part Remo Mnero Outne The roe of a mappng functon n mutdmensona ntegraton Gauss approach n more dmensons and quadrature rues Crtca anass of acceptabt of a gven quadrature

More information