Walter D. Pilkey. Dmitry V. Balandin. Nikolai N. Bolotnik. Sergey V. Purtsezov. Introduction. Nomenclature

Size: px
Start display at page:

Download "Walter D. Pilkey. Dmitry V. Balandin. Nikolai N. Bolotnik. Sergey V. Purtsezov. Introduction. Nomenclature"

Transcription

1 Waler D. Pilkey e al Waler D. Pilkey wdp@virgiia.edu Uiversiy of Virgiia Dep. of Mechaical ad Aerospace Eg Charloesville, VA. USA Diry V. Baladi baladi@pk.u.rue.ru izhy ovgorod Sae Uiversiy Dep. of Cop. Maheaics ad Cybereics 6395 izhy ovgorod. Russia ikolai. Boloik boloik@ipe.ru Isiue for Probles i Mechaics of he Russia Acadey of Scieces 956 Moscow. Russia Exreal Disurbace Aalysis for Dyaical Syses wih Ucerai Ipu The philosophy of he exreal disurbace aalysis for dyaical syses wih ucerai ipus is described. This aalysis ivolves solvig opial corol probles i which he ie hisories of he ipus (exeral disurbaces) are regarded as he corol fucios ad a respose easure of he syse serves as he perforace idex. The perforace idex should be axiized or iiized over he disurbaces wihi a prescribed class. Ofe his class is specified by lower ad upper bouds (a corridor) bewee which he values of he disurbace us lie. The axiizaio ad iiizaio probles are referred o as he wors disurbace ad bes disurbace probles, respecively. The soluios of hese probles provide he exree values bewee which he respose easure lies for ay disurbace fro he specified class. The exreal disurbace aalysis is ipora, i paricular, whe desigig sadards for esig devices for he proecio of fragile objecs fro ipac loads. This approach is illusraed for a sigledegree-of-freedo syse ha ca be regarded as a siplified odel of he equipe for sled ess of auoobile resrai syses. Keywords: Ucerai ipu, wors disurbace, bes disurbace, ipac isolaio, sled ess Sergey V. Pursezov pusv@virgiia.edu Uiversiy of Virgiia Dep. of Mechaical ad Aerospace Eg Charloesville, VA. USA Iroducio There are uerous srucural dyaics probles where he loadig is o fully defied. Typically, such probles are reaed as sochasic probles usig he heory of probabiliy (Hass, 5; Iyeger, 97). This paper discusses a aleraive deeriisic approach i which he available iforaio o he loadig is uilized o calculae bouds o possible resposes of he srucure. The lower boud would be he leas axiu respose ha ca occur. The correspodig loadig is referred o as he bes disurbace. The upper boud is he greaes possible peak respose ad he correspodig loadig is he wors disurbace (Koivuiei, 966; O Hara, 973). The respose of he srucure o ay oher loadig fallig wihi he prescribed loadig bouds will lie bewee he exree values of he respose. Oe exaple occurs if a explosive eviroe is characerized i ers of a oal ipulse. This is sufficie iforaio o deerie wha us happe o soe equipe, e.g., iliary equipe, ad he axiu respose of his equipe. Oher exaple occurs i he esig of safey devices for crashig vehicles. For isace, ie-dow syses for wheelchairs are esed o high-speed sleds for which he crash pulse us adhere o sadards ha are defied i ers of a evelope i which he crash pulse us lie (Kag ad Pilkey, 998). A eve ore coo exaple occurs i he sled esig of child seas, where he sled pulse us lie wihi prescribed ie-varyig bouds (Cradall e al., 996). The proposed schee provides upper ad lower bouds for criical resposes. These resposes ca be used o help gauge wheher he prescribed pulse corridor is oo igh or oo broad Preseed a XI DIAME Ieraioal Syposiu o Dyaic Probles of Mechaics, February 8h - March 4h, 5, Ouro Preo. MG. Brazil. Paper acceped: Jue, 5. Techical Ediors: J.R.F. Arruda ad D.A. Rade. o provide effecive sadards. Thus, he procedure here develops a echology o sudy he sesiiviy of dyaic resposes o icopleely prescribed es codiios. We choose o use he sled esig of child seas as a exaple of usig exreal disurbace aalysis o illusrae he priciples of his sesiiviy echology. oeclaure A, B, C, D = coordiaes of he verices of he rapezoid, defied by Eq. (8), (ie, acceleraio), (s, /s ) a = heigh of he rapezoid, defied by Eq. (7), acceleraio, /s c = dapig coefficie of he resrai syse, dapig coefficie, kg/s F = force applied o he sled, force, f = force applied o he duy by he resrai syse, force, g( ) = ipulse respose of he syse, s h = ie bewee saples, ie, s i = saple uber, diesioless j = desig variable uber, diesioless J = peak agiude of he absolue acceleraio of he duy, acceleraio, /s K( ) = fucio, defied by Eq. (), /s k = siffess coefficie of he resrai syse, siffess coefficie, kg/s L = uber of saples, diesioless M = ass of he sled, ass, kg = ass of he duy, ass, kg = uber of desig variables, diesioless p = paraeer defied by Eq. (), diesioless 5 / Vol. XXVIII, o. 4, Ocober-Deceber 6 ABCM

2 Exreal Disurbace Aalysis for Dyaical Syses wih Ucerai Ipu R = sesiiviy of he peak force rasied o he duy relaive o he variaio i he crash pulse, diesioless R d = prescribed value of he sesiiviy of he peak force rasied o he duy relaive o he variaio i he crash pulse, diesioless T = duraio of he crash pulse, ie, s = ie, s u = absolue acceleraio of he duy wih ius sig, acceleraio, /s W = class of he desig variables w, acceleraio, /s w = absolue acceleraio of he sled wih ius sig reaed as a desig variable, acceleraio, /s w ( ) = lower boud of he desig variable w, acceleraio, /s w ( ) = upper boud of he desig variable w, acceleraio, /s x = displacee of he duy relaive o he sled, displacee, z = displacee of he sled relaive o he fixed frae, displacee, Greek Sybols v = velociy chage durig he crash, velociy, /s v = lower boud of he velociy chage durig he crash, velociy, /s v = upper boud of he velociy chage durig he crash, velociy, /s δ ( ) = Dirac dela fucio, diesioless ς = fracio of criical dapig, diesioless λ = siilariy facor, diesioless λ = siilariy facor equal o.4, diesioless µ = reduced ass of he syse, defied by Eq. (), ass, kg τ = variable of iegraio, ie, s τ = paraeers of he rapezoid, defied by Eq. (7), ie, s 4 ψ ( ) = fucio, represeig a rapezoid ad defied by Eq. (7), acceleraio, /s ω = agular udaped aural frequecy of he syse, agular frequecy, rad/s ω D = agular daped aural frequecy of he syse, agular frequecy, rad/s Subscrips i relaive o saple uber j relaive o desig variable uber relaive o he opial soluio of he bes disurbace proble relaive o he lower boud of he desig variable w relaive o he upper boud of he desig variable w Superscrips relaive o he opial soluio of he wors disurbace proble relaive o he upper boud of he desig variable w Basic Coceps of he Bes ad Wors Disurbace Aalyses; Saee of he Opial Corol Probles A sled es is ofe eployed o ivesigae proecive characerisics of safey devices. For exaple, auoobile resrai syses ca be sudied wih sled ess. I paricular, cosider a sled es for a auoobile resrai ha resrics he oio of a occupa o preve a viole coac wih copoes iside he car durig a crash. Typically, he es sled equipe cosiss of a highspeed sled wih a ock-up of he car ierior, a occupa sea, a resrai syse, e.g., sea bel, ad a duy as he occupa. The car ock-up wih ierior ad seabels are rigidly aached o he sled. The sled, duy, ad he sea bels are equipped wih sesors ha easure acceleraios, displacees, srais, forces, ad oher characerisics of he respose of he syse. The sled is acceleraed o a prescribed velociy ad he is subjec o a deceleraio pulse siulaig he crash ipac. The various easured resposes ca be used o for ijury crieria ha ca be copared wih esablished olerace levels. To fix ideas, we cosider a siplified odel of he esig equipe i which he sled ad he duy are regarded as rigid bodies (Fig. ). The sled oves alog a horizoal sraigh lie relaive o a fixed (ierial) referece frae ad he duy oves relaive o he sled alog he sae lie. Le z deoe he displacee (coordiae) of he sled relaive o he fixed frae, x he displacee (coordiae) of he duy relaive o he sled, M he ass of he sled, he ass of he duy, F he force applied o he sled durig he crash siulaio, ad f he force applied o he duy by he resrai syse. The ass of he sea bels is egleced as copared wih he asses of he sled ad he duy. The force f is assued o be a fucio of he paraeers of he relaive oio of he duy, he relaive coordiae x ad he relaive velociy xɺ, i.e., f = f ( x, xɺ ). This fucio akes io accou he elasic ad dissipaive properies of he seabels. The oio of his echaical syse is govered by he siulaeous equaios Mz ɺɺ = F f, ( ɺɺ x ɺɺ z) = f, () which represe ewo s secod law for he sled ad he duy, respecively. The duy is aced upo oly by he resrai syse force f, whereas he sled is aced upo by he ipac force F ad he force f as he sea bel ieracs wih he sled. Solve he syse of Eq. () for xɺɺ ad zɺɺ o obai M µ ɺɺ x = f ( x, xɺ ) F, Mz ɺɺ = F f ( x, xɺ ), µ =. () M M The quaiy µ is called he reduced ass of he syse. I he syse of Eq. (), he equaio goverig he oio of he duy relaive o he sled is idepede of he equaio of oio of he sled relaive o he fixed frae. If oe is ieresed oly i he relaive oio of he duy, ad he ipac force is a kow fucio of ie, i.e., F = F( ), oly he firs relaio of () is used. I he heory of ipac isolaio, he ipac disurbace specified as he force pulse acig o he base (he sled) is referred o as he dyaic disurbace. z M x c k Figure. Sigle-degree-of-freedo odel. J. of he Braz. Soc. of Mech. Sci. & Eg. Copyrigh 6 by ABCM Ocober-Deceber 6, Vol. XXVIII, o. 4 / 53

3 Waler D. Pilkey e al I a uber of cases, i is reasoable o assue ha he ipac disurbace is characerized by he acceleraio of he base, raher ha by he force applied o i. For isace, for he sled ess, he ipac deceleraio of he sled is easured by acceleroeers. The he relaive oio of he body o be proeced (he duy) is govered by he secod relaio of Eq. () where ɺɺ z( ) is reaed as a kow fucio of ie. The ipac disurbace specified as he base acceleraio pulse is referred o as he kieaic disurbace. For boh ypes of disurbaces, he equaio of he relaive oio of he body o be proeced ca be represeed i a uified for as ɺɺ x u( x, xɺ ) = w( ), (3) where u( x, xɺ ) = f ( x, x ɺ )/ µ, w( ) = F( ) / M for he dyaic disurbace ad u( x, xɺ ) = f ( x, xɺ )/, w( ) = ɺɺ z( ) for he kieaic disurbace. I he heory of ipac isolaio, he fucio u( x, xɺ ) is frequely called he characerisic of he ipac isolaor. As a rule, Eq. (3) is subjeced o zero iiial codiios x() =, xɺ () =. (4) These codiios iply ha he objec o be proeced does o ove relaive o he base uil he ipac pulse has bee applied o he base. To copare he disurbaces i ers of he respose of he objec o be proeced, oe should have a uerical easure of his respose. We selec J ( w) = ax u( x( ), xɺ ( )) (5) where x( ) sads for he soluio of he differeial equaio (3) subjec o he iiial codiios of Eq. (4) for a give w = w( ). The quaiy J characerizes he peak agiude of he force rasied o he duy by he resrai syse. (Acually, he peak force is equal o µ J or J for he case of he dyaic or kieaic disurbace, respecively.) The easure J ( w ) is a fucioal of he disurbace, depedig o he eire ie hisory of w. Cosider ow he saee of he wors ad bes disurbace probles for he syse of Eqs. (3) ad (4) wih he respose easure of Eq. (5). Le he disurbaces belog o a prescribed class W, i.e., w W. oe ha W is a se of fucios raher ha a se of values of hese fucios. I ers of he sled es, he class W ca characerize a se of disurbaces ha are represeaive of he es, akig io accou he ieviable iaccuracy i he reproducio of he desired deceleraio pulse. Usually, he class of represeaive disurbaces ivolves a corridor or evelope (Fig. ) i which hese disurbaces us lie, i.e., W = { w( ) : w ( ) w( ) w ( ), [, T ]}, (6) where w ( ) ad w ( ) are prescribed fucios ha defie he lower ad upper bouds of he corridor ad T is a fixed ie. I addiio, he velociy of he sled should chage by a cerai value durig he crash pulse ie T. I he case of he kieaic disurbace, his codiio has he for T v = w( τ ) dτ, (7) where v is a specified posiive quaiy. Usually, he quaiy v is prescribed wih a allowace for a error i he easuree of he velociy of he sled, v v v, (8) where he ierval [ v, v] characerizes he uceraiy i he easuree of he velociy decrease. Sice v is defied as he iegral of w, he iequaliies of Eq. (8) ipose addiioal cosrais o he deceleraio pulse ie hisory. Therefore, he oal class o which he allowed deceleraio pulses us belog is defied by T W = { w( ) : w ( ) w( ) w ( ), v w( τ ) dτ v}. (9) w - () Tie w () Figure. Corridor i which a ipu pulse w() us lie. The wors disurbace proble is ha of he deeriaio of he upper boud of he fucioal J ( w ) for he disurbaces of he class W. Proble (Wors Disurbace Proble). For he syse govered by Eqs. (3) ad (4), fid he disurbace fucio w ( ) belogig o he class W such ha J w = () ( ) ax w W J ( w). The bes disurbace proble is ha of he deeriaio of he lower boud of he fucioal J ( w ) for he disurbaces of he class W. Proble (Bes Disurbace Proble). For he syse govered by Eqs. (3) ad (4), fid he disurbace fucio w( ) belogig o he class W such ha J ( w ) = i J ( w). () w W Fro he aheaical poi of view Probles ad are opial corol probles wih a axiu ype fucioal, sice he perforace idex i hese probles is he axiu of he absolue value of he force rasied o he duy over he ipac respose ie. The wors ad bes disurbace aalyses eable oe o ivesigae he sesiiviy of he resuls of he sled es o he variaio of he ipac pulse wihi he corridor prescribed by he sadards for such ess. If he differece bewee he upper ad lower bouds is large, he es ay have a udesirably low value of he peak force rasied o he duy due o he occurrece of he disurbace close o he bes oe. Such a es ca be decepive whe assessig he qualiy of he resrai syse ad, hece, is illdesiged. I his case, he sadards for he defiig allowable ipac pulses should be revised. Soeies, he sesiiviy of he peak force rasied o he duy relaive o he variaio i he ipac pulse is characerized by he raio 54 / Vol. XXVIII, o. 4, Ocober-Deceber 6 ABCM

4 Exreal Disurbace Aalysis for Dyaical Syses wih Ucerai Ipu R = J w J w () ( ) / ( ). Resrai Syse wih Liear Characerisics; uerical Soluio of he Opial Corol Probles Le he resrai syse have liear elasic ad dapig properies, i.e., f ( x, xɺ ) = cxɺ kx, (3) where c ad k are he dapig ad siffess coefficies, respecively. For a ipac disurbace of he kieaic ype, Eq. (3) becoes ɺɺ x ςω xɺ ω x = w, (4) ( ) where ω is he aural frequecy of he syse ad ς is he dapig raio, / / = ( / ), ς = /(( ) ). (5) ω k c k The crierio J of Eq. (5) (he peak agiude of he occupa s absolue acceleraio) becoes J w ςω xɺ ω x. (6) ( ) = ax ( ) ( ) Eve for he liear syse of Eq. (4), he soluio of he bes ad wors disurbace probles (Probles ad ) requires uerical ehods. To his ed, he coiuous-ie forulaio of hese probles is replaced by a discree-ie approxiaio. The ie axis is discreized wih sep size h. O he ie iervals ( i ) h < ih, i =,,..., he fucio w( ) is assued o be cosa: wi, if ( i ) h < T w( ) =, ( i ) h < ih. (7), if ( i ) h T The cosa paraeers wi, i =,...,, play he role of desig variables whe solvig he opiizaio probles. The uber of hese desig variables is defied as follows: T / h, if { T / h} = =, (8) [ T / h], if { T / h} where he square brackes ad he braces deoe he ieger ad fracioal pars, respecively, of he expressio eclosed. The soluio of Eq. (4) subjec o he iiial codiios of Eq. (4) has he covoluio for where x( ) = g( τ ) w( τ ) dτ, xɺ ( ) = gɺ ( τ ) w( τ ) dτ, (9) ωd exp( ζω)si ωd, if ς <, / g( ) = exp( ςω), if ς =, ωd = ω ς. () ωd exp( ζω)sih ωd, if ς >, The fucio g( ) is referred o as he ipulse respose fucio (fudaeal soluio) for Eq. (4). I saisfies he differeial equaio of (4) wih w( ) = δ ( ), where δ ( ) is he Dirac dela fucio, subjec o he iiial codiios x () = ad xɺ ( ) =. This fucio ca also be defied as he soluio of Eq. (4) wih zero righ-had side ( w( ) ), subjec o he iiial codiios x () = ad x ɺ () =. Subsiue he relaios of Eq. () io Eq. (6) o obai ax ( τ ) ( τ ) τ, ( ) ςω ( ) ω ( ). J = K w d K = gɺ g () The discree-ie approxiaio of he respose easure of Eq. () has he for jh ( j ) h i{ i, } J ( w) = ax p w, i [: L] ij j j = p = K( ih τ ) dτ, j i, w = [ w,..., w ], ij () where L is he uber of discree iervals o he ie axis o which he respose o he ipac pulse is o be calculaed. To obai reliable resuls, i is ecessary o cosider he respose o he ie ierval subsaially exceedig ha of he crash pulse. Sice he crash pulse has a raher shor duraio, he respose easure ca aai is axiu afer he disurbace has ceased o ac. For his reaso, he uber of ie isas a which he respose is easured should exceed he uber of he discreizaio pois i he crash pulse ierval [, T ] ad, hece, L >. The crierio of Eq. () is a fucio of variables w = [ w,..., w ]. The cosrais of Eq. (9) are discreized as follows: w ( ih) w w ( ih), i [: ]; v h w v. (3) i j j= I he discree-ie forulaio, he wors-disurbace proble is reduced o he axiizaio of he fucio of Eq. () ad he bes-disurbace proble o he iiizaio of his fucio, subjec o he cosrais of Eq. (3). These cosrais are liear relaive o he desig variables w,..., w, ad he fucio o be axiized or iiized is a axiu of he absolue values of liear fucios of hese variables. I his case, he cosraied iiizaio (axiizaio) proble ca be reduced o ha of he liear prograig, which subsaially faciliaes he soluio, sice here are rapidly covergig reliable liear prograig algorihs available i os opiizaio sofware. For a ore deailed descripio of he soluio of opiizaio probles siilar o Probles ad o he basis of he liear prograig, see Baladi, Boloik, ad Pilkey (). uerical Exaple Assue he dapig ad siffess coefficies of he resrai syses are: c s k 5 =, =, = 3 kg, (4) J. of he Braz. Soc. of Mech. Sci. & Eg. Copyrigh 6 by ABCM Ocober-Deceber 6, Vol. XXVIII, o. 4 / 55

5 Waler D. Pilkey e al which are reasoably realisic for resrai syses. The ass of he occupa of 3kg correspods o he 6-year-old child duy. Figures 3 o 5 deosrae he resuls of he soluio of Probles ad for he sled deceleraio pulse corridor show i Fig. 3. The corridor correspods o he sadard of HTSA (979) ha specifies es codiios for child resrais. The uceraiy ierval of Eq. (8) for he sled velociy chage is defied as 45.k/h v 48.3 k/h (8 ph v 3 ph), (5) Deceleraio (g) 3 Bes Pulse Wors Pulse Figure 3 shows he ie hisories for he bes ad wors deceleraio pulses. The shaded area i Fig. 3 represes he corridor i which deceleraio pulses us lie. If here were o cosrais such as of Eq. (5) o he velociy chage, boh he bes ad wors disurbace pulses would have swiched bewee he lower ad upper bouds of he corridor. This follows sice he soluio of a liear prograig proble always lies o he boudary of he doai cosraiig he desig variables. Fro he ie hisories show i Fig. 3 i is appare ha boh he bes disurbace ad wors disurbace pulses lie wihi he walls of he corridor o soe ie iervals. Therefore, he velociy chage assues oe of he boudary values prescribed by Eq. (5), specifically, v = 45. k/h for he bes disurbace ad v = 48.3 k/h for he wors disurbace. Figure 4 depics he ie hisory of he force applied o he duy by he sea bel i he case of he wors (solid curve) ad bes (dashed curve) disurbaces. I is appare ha i boh cases he peak agiude of he force is aaied before he ed of he disurbace. This is because he duraio of he pulse is close o he udaped vibraio period of he syse. I accordace wih figs. 3, he duraio of he deceleraio pulse lies bewee.75 ad.9s, while he vibraio period for he paraeers of Eq. (4) is.96s. The ie hisories of he displacee of he duy for he wors ad bes disurbaces are show i Fig. 5. I fac, he curves of Fig. 5 repea hose of Fig. 4 scaled by he siffess coefficie 5 ( / ). This is due o he fac ha he dapig does o have uch effec durig he ie ierval uder cosideraio. The dapig raio ς of Eq. (5) calculaed for he syse wih he paraeers of Eq. (4) equals -.66 s ad, accordigly, he characerisic ie of he dapig equals 5.s, which subsaially exceeds he udaped vibraio period. The sesiiviy raio of Eq. () for he forces rasied o he duy is R =. /7.6 =.34. (6) This reasoably low value would probably be accepable i pracice. Force (k) Displacee () Tie (s) 5 Figure 3. Bes ad wors disurbaces... Bes Wors.5..5 Tie (s) Figure 4. Force applied o he duy by he sea bel Bes Wors.5..5 Tie (s) Figure 5. Displacee of he duy relaive o he sled. Raioal Desig of Sled Tes Sadards I he previous secio, we have illusraed he applicaio of he exreal disurbace aalysis o check he sesiiviy of he respose easure (he force rasied o he occupa) o he variaio of he sled deceleraio pulse wihi a prescribed corridor. This cocep suggess a echique for he raioal desig of such a corridor so ha he sesiiviy raio does o exceed a allowed value. Cosider a exaple. Le he walls of he corridor have a shape of a rapezoid ABCD show i Fig. 6. A geeric rapezoid ca be represeed aalyically by he piecewise liear fucio 56 / Vol. XXVIII, o. 4, Ocober-Deceber 6 ABCM

6 Exreal Disurbace Aalysis for Dyaical Syses wih Ucerai Ipu, (, τ) [, ), a( τ)/ τ, [ τ, τ τ ) [, ), ψ ( ) = a, [ τ τ, τ τ 3) [, ), a( τ τ 4) /( τ 3 τ 4), [ τ τ 3, τ τ 4) [, ),, τ τ 4. (7) The coordiaes of he verices of his rapezoid o he ieacceleraio plae are give by, (, τ ) [, ), a( τ ), [, τ τ λτ ) [, ), τ w ( ) = λa, [ τ λτ, τ λτ 3) [, ), a, τ λτ 4, ( τ λτ 4), [ 3, τ λτ τ λτ 4 ) [, ), ( τ 3 τ 4) (3) A = ( τ,), B = ( τ τ, a), C = ( τ τ, a), D = ( τ τ,). 3 4 (8) I ers of he geoery of he rapezoid, τ 4 ad τ 3 τ are he leghs of he lower ad upper bases, respecively, a is he heigh, / / ( τ a ) ad (( τ τ ) a ) are he leghs of he laeral 4 3 sides AB ad CD, respecively. The paraeer τ characerizes he shif of he rapezoid alog he ie axis. The walls (he upper ad lower bouds) of he corridor for sled ess show i Fig. 3 have he shape of rapezoids wih for he lower boud ad τ = 4s, τ = 9s, τ = 4s, 3 τ = 7s, a=9g 4 τ = s, τ = s, τ = 53s, 3 τ = 9s, a=5g 4 (9) (3) for he upper boud. Wih a high degree of accuracy, he rapezoid of he upper boud is siilar o he rapezoid of he lower boud wih he siilariy facor close o.3. Deceleraio a A B C D τ τ τ τ τ 3 τ τ 4 Tie Figure 6. A rapezoidal boud for he force applied o he duy by he sea bel. We will cofie ourselves o he case where he fucios w ( ) ad w ( ) of Eq. (6) are represeed by siilar rapezoids. Le he lower boud be fixed ad le w ( ) = ψ ( ), where ψ ( ) is defied by Eq. (7). Le he upper boud be described by he fucio where λ ad τ are he paraeers of he upper boud, he siilariy facor ad he ie coordiae of he lef-had boo verex (poi A ) of he respecive rapezoid. For he curve w ( ) of Eq. (3) o lie above he curve w ( ) = ψ ( ) of Eq. (7) i is ecessary ha λ, τ ( λ) τ τ τ. (3) 4 I is appare ha for λ =, he lower ad upper bouds of he corridor ach oe aoher. By varyig he paraeers λ ad τ wihi he doai of Eq. (3) oe ca assure ha he raio R of Eq. () does o exceed a prescribed value. The variaio ca be orgaized i various ways. For exaple, oe ca cosrai he paraeer τ o lie i he iddle of he ierval allowed for his paraeer by Eq. (3), i.e., τ = τ ( λ) τ /. (33) 4 I his case, he variaio is perfored wih respec o oly oe paraeer, λ. As λ icreases, he corridor becoes wider ad, herefore, he bes-disurbace respose easure decreases while he wors-disurbace respose easure icreases. Accordigly, he sesiiviy raio R is a oooically icreasig fucio of λ. The search for he axiu λ ha assures ha he quaiy R does o exceed he prescribed value R is reduced o he soluio of he equaio R( λ ) = R d. This equaio ca be solved by various ehods, for exaple, by he ierval bisecio ehod. For each rial λ, oe should solve he bes disurbace ad wors disurbace probles o calculae R( λ ). Figure 7 preses he curve R( λ ), calculaed for he corridor defied by Eqs. (7), (9), (3), ad (33), copleed by he velociy chage ierval of Eq. (5). The curve becoes ore fla as λ icreases. oe ha his curve begis wih λ = λ.4, raher ha wih λ =, i which case he upper ad lower bouds of he crash pulse corridor coicide. This is because of he cosrai of Eq. (5) o he velociy chage. The lower boud for he velociy chage is 45.k/h, whereas he iegral of he lower boud of he pulse corridor (he area of he rapezoid represeig he respecive wall of he corridor) is 34.9k/h. Therefore, he pulse of he lower boud is uable o decelerae o a coplee sop eve for he leas crash velociy allowed by he ierval of Eq. (5). Moreover, o pulse of he corridor is able o provide a velociy chage of 45.k/h while he area of he upper boud rapezoid is less ha his value. This is he case for λ < λ. For λ = λ, he area of he upper boud rapezoid is equal o 45.k/h. I his case, he oly pulse of he corridor able o decelerae he sled crashig a his velociy o a coplee sop is ha of he upper boud ad, d J. of he Braz. Soc. of Mech. Sci. & Eg. Copyrigh 6 by ABCM Ocober-Deceber 6, Vol. XXVIII, o. 4 / 57

7 Waler D. Pilkey e al accordigly, he wors ad bes disurbaces coicide, i.e., R( λ ) =. R λ Figure 7. Sesiiviy of he peak force rasied o he duy relaive o he variaio i he ipac pulse R versus he siilariy facor λ. Coclusios The exreal disurbace aalysis eables oe o evaluae he sesiiviy of he respose of a dyaical syse o he variaio i he exeral disurbaces wihi a prescribed uceraiy class. To ha ed oe should solve wo opial corol probles, he bes disurbace proble o calculae he lower boud ad he wors disurbace proble o calculae he upper boud of a respose easure. The raio of hese bouds ca serve as he sesiiviy idex. The closer his raio is o uiy, he lower he sesiiviy. The exreal disurbace aalysis is, i paricular, ipora for he validaio of sadards for ipac ess of fragile objecs ad he equipe for he proecio of such objecs. Such sadards cooly prescribe a corridor for ipac pulses o be reproduced o he es faciliy. The es daa are cosidered o be reliable if he spread i he resposes o he ipu pulses wihi he corridor is reasoably sall. This spread is characerized by he sesiiviy raio resulig fro he exreal disurbace aalysis. I addiio, he exreal disurbace aalysis suggess a echique for raioal desig of a corridor for loadig pulses so ha he sesiiviy raio does o exceed a prescribed value. Ackowledges This research was parly suppored by SF (gra BES- 3337), ATO (gra CBP.R.RCLG 988), ad he Russia Foudaio for Basic Research (gras ad 4- -). Refereces Baladi, D.V., Boloik,.., ad Pilkey, W.D.,, Opial Proecio fro Ipac, Shock ad Vibraio, Taylor ad Fracis Publishers, Philadelphia, PA, 437 p. Cradall, J.R., Pilkey, W.D., Kag, W., ad Bass, C.R., 996, Sesiiviy of Occupa Subjec o Prescribed Corridors, Joural of Shock ad Vibraio, Vol. 3, o. 6, pp Hass, M., 5, Applied Fuzzy Ariheic - A Iroducio wih Egieerig Applicaios, Spriger, Berli, 56 p. Iyeger, R.., 97, Wors Ipus ad a Boud o he Highes Peak Saisics of a Class of oliear Syses, Joural of Soud ad Vibraio, Vol. 5, o., pp Kag, W., ad Pilkey, W.D., 998, Crash Siulaios of Wheelchair Occupa Syses i Traspor, Joural of Rehabiliaio Research ad Develope, Vol. 35, o., pp Koivuiei, A.J., 966, Paraeer Opiizaio i Syses Subjec o Wors Bouded Disurbace, IEEE Trasacios o Auoaic Corol, Vol., o. 3, pp aioal Highway Traffic Safey Adiisraio (HTSA), 979 (Ocober, ed.), FMVSS 3 Child Resrai Syses, Vol. 49 CFR O Hara, G.J., 973, Maxiizaio ad Miiizaio of Dyaic Load Facors, Shock ad Vibraio Bullei, Vol. 43, Par, pp / Vol. XXVIII, o. 4, Ocober-Deceber 6 ABCM

Introduction to Mechanical Vibrations

Introduction to Mechanical Vibrations CHAPTER 1 Iroducio o Mechaical Vibraios Vibraio is he oio of a paricle or a body or syse of coeced bodies displaced fro a posiio of equilibriu. Mos vibraios are udesirable i achies ad srucures because

More information

Problems and Solutions for Section 3.2 (3.15 through 3.25)

Problems and Solutions for Section 3.2 (3.15 through 3.25) 3-7 Problems ad Soluios for Secio 3 35 hrough 35 35 Calculae he respose of a overdamped sigle-degree-of-freedom sysem o a arbirary o-periodic exciaio Soluio: From Equaio 3: x = # F! h "! d! For a overdamped

More information

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory Liear Time-Ivaria Sysems (LTI Sysems) Oulie Basic Sysem Properies Memoryless ad sysems wih memory (saic or dyamic) Causal ad o-causal sysems (Causaliy) Liear ad o-liear sysems (Lieariy) Sable ad o-sable

More information

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion BE.43 Tuorial: Liear Operaor Theory ad Eigefucio Expasio (adaped fro Douglas Lauffeburger) 9//4 Moivaig proble I class, we ecouered parial differeial equaios describig rasie syses wih cheical diffusio.

More information

Bernstein Direct Method for Solving. Variational Problems

Bernstein Direct Method for Solving. Variational Problems Ieraioal Maheaical Foru, 5,, o. 48, 35-37 Bersei Direc Mehod for Solvig Variaioal Probles Sadeep Dixi*, Viee K. Sigh**, Ai K. Sigh*, O P. Sigh* *Depare of Applied Maheaics, Isiue of echology, Baaras Hidu

More information

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1) Fourier Series Iroducio I his secio we will sudy periodic sigals i ers o heir requecy is said o be periodic i coe Reid ha a sigal ( ) ( ) ( ) () or every, where is a uber Fro his deiiio i ollows ha ( )

More information

Optimum design of complementary transient experiments for estimating thermal properties

Optimum design of complementary transient experiments for estimating thermal properties Opiu desig of copleeary rasie experies for esiaig heral properies Jaes V. Beck*, Filippo de Moe, Doald E. Aos *Depare of Mechaical Egieerig, Michiga Sae Uiversiy, USA Depare of Idusrial ad Iforaio Egieerig

More information

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation. ecure Phys 375 Eergy Desiy / Eergy Flu / oal Eergy i D Overview ad Moivaio: Fro your sudy of waves i iroducory physics you should be aware ha waves ca raspor eergy fro oe place o aoher cosider he geeraio

More information

On Acoustic Radiation by Rotating Dipole Sources in Frequency Domain

On Acoustic Radiation by Rotating Dipole Sources in Frequency Domain www.ccsee.org/as Moder Applied Sciece Vol. 5, No. 5; Ocober 211 O Acousic Radiaio by Roaig Dipole Sources i Frequecy Doai Zhihog Liu Ceer of Eergy ad Eviroe, QigDao Techological Uiversiy 11 FuShu Road,

More information

CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS

CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS I. Karafyllis ad J. Tsiias Depare of Maheaics, Naioal Techical Uiversiy of Ahes, Zografou Capus 578, Ahes, Greece Eail: jsi@ceral.ua.gr.

More information

E will be denoted by n

E will be denoted by n JASEM ISSN 9-8362 All rigs reserved Full-ex Available Olie a p:// wwwbiolieorgbr/ja J Appl Sci Eviro Mg 25 Vol 9 3) 3-36 Corollabiliy ad Null Corollabiliy of Liear Syses * DAVIES, I; 2 JACKREECE, P Depare

More information

Sliding Mode Control for Robust Stabilization of Uncertain Input-Delay Systems

Sliding Mode Control for Robust Stabilization of Uncertain Input-Delay Systems 98 CSE: he siue of Corol, uoaio ad Syses Egieers, KORE Vol, No, Jue, Slidig Mode Corol for Robus Sabilizaio of Ucerai pu-delay Syses Youg-Hoo Roh ad Ju-Ho Oh bsrac: his paper is cocered wih a delay-depede

More information

Calculus BC 2015 Scoring Guidelines

Calculus BC 2015 Scoring Guidelines AP Calculus BC 5 Scorig Guidelies 5 The College Board. College Board, Advaced Placeme Program, AP, AP Ceral, ad he acor logo are regisered rademarks of he College Board. AP Ceral is he official olie home

More information

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar Ieraioal Joural of Scieific ad Research Publicaios, Volue 2, Issue 7, July 22 ISSN 225-353 EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D S Palikar Depare of Maheaics, Vasarao Naik College, Naded

More information

PLASTIC BUCKLING OF SSSS THIN RECTANGULAR PLATES SUBJECTED TO UNIAXIAL COMPRESSION USING TAYLOR-MACLAURIN SHAPE FUNCTION

PLASTIC BUCKLING OF SSSS THIN RECTANGULAR PLATES SUBJECTED TO UNIAXIAL COMPRESSION USING TAYLOR-MACLAURIN SHAPE FUNCTION I. J. Sruc. & Civil Egg. es. 013 U G Eziefula e al., 013 esearch Paper ISSN 319 6009 www.ijscer.co Vol., No., Noveer 013 013 IJSCE. All ighs eserved PLASTIC BUCKLING OF SSSS THIN ECTANGULA PLATES SUBJECTED

More information

The analysis of the method on the one variable function s limit Ke Wu

The analysis of the method on the one variable function s limit Ke Wu Ieraioal Coferece o Advaces i Mechaical Egieerig ad Idusrial Iformaics (AMEII 5) The aalysis of he mehod o he oe variable fucio s i Ke Wu Deparme of Mahemaics ad Saisics Zaozhuag Uiversiy Zaozhuag 776

More information

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative Malaysia Joural of Maheaical Scieces 7(): 3-7 (3) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Joural hoepage: hp://eispe.up.edu.y/joural Boudary Value Proble for he Higher Order Equaio wih Fracioal Derivaive

More information

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER Maerials Physics ad Mechaics 3 (5) 36-4 Received: March 7 5 METHOD OF THE EQUIVAENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBEM FOR EASTIC DIFFUSION AYER A.V. Zemsov * D.V. Tarlaovsiy Moscow Aviaio Isiue

More information

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis Joural of aerials Sciece ad Egieerig B 5 (7-8 (5 - doi: 765/6-6/57-8 D DAVID PUBLISHING Opimizaio of Roaig achies Vibraios Limis by he Sprig - ass Sysem Aalysis BENDJAIA Belacem sila, Algéria Absrac: The

More information

Introduction to Mobile Robotics Mapping with Known Poses

Introduction to Mobile Robotics Mapping with Known Poses Iroducio o Mobile Roboics Mappig wih Kow Poses Wolfra Burgard Cyrill Sachiss Mare Beewi Kai Arras Why Mappig? Learig aps is oe of he fudaeal probles i obile roboics Maps allow robos o efficiely carry ou

More information

Notes 03 largely plagiarized by %khc

Notes 03 largely plagiarized by %khc 1 1 Discree-Time Covoluio Noes 03 largely plagiarized by %khc Le s begi our discussio of covoluio i discree-ime, sice life is somewha easier i ha domai. We sar wih a sigal x[] ha will be he ipu io our

More information

Electrical Engineering Department Network Lab.

Electrical Engineering Department Network Lab. Par:- Elecrical Egieerig Deparme Nework Lab. Deermiaio of differe parameers of -por eworks ad verificaio of heir ierrelaio ships. Objecive: - To deermie Y, ad ABD parameers of sigle ad cascaded wo Por

More information

Comparison between Fourier and Corrected Fourier Series Methods

Comparison between Fourier and Corrected Fourier Series Methods Malaysia Joural of Mahemaical Scieces 7(): 73-8 (13) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Joural homepage: hp://eispem.upm.edu.my/oural Compariso bewee Fourier ad Correced Fourier Series Mehods 1

More information

Chemistry 1B, Fall 2016 Topics 21-22

Chemistry 1B, Fall 2016 Topics 21-22 Cheisry B, Fall 6 Topics - STRUCTURE ad DYNAMICS Cheisry B Fall 6 Cheisry B so far: STRUCTURE of aos ad olecules Topics - Cheical Kieics Cheisry B ow: DYNAMICS cheical kieics herodyaics (che C, 6B) ad

More information

Big O Notation for Time Complexity of Algorithms

Big O Notation for Time Complexity of Algorithms BRONX COMMUNITY COLLEGE of he Ciy Uiversiy of New York DEPARTMENT OF MATHEMATICS AND COMPUTER SCIENCE CSI 33 Secio E01 Hadou 1 Fall 2014 Sepember 3, 2014 Big O Noaio for Time Complexiy of Algorihms Time

More information

STK4080/9080 Survival and event history analysis

STK4080/9080 Survival and event history analysis STK48/98 Survival ad eve hisory aalysis Marigales i discree ime Cosider a sochasic process The process M is a marigale if Lecure 3: Marigales ad oher sochasic processes i discree ime (recap) where (formally

More information

ES 330 Electronics II Homework 03 (Fall 2017 Due Wednesday, September 20, 2017)

ES 330 Electronics II Homework 03 (Fall 2017 Due Wednesday, September 20, 2017) Pae1 Nae Soluios ES 330 Elecroics II Hoework 03 (Fall 017 ue Wedesday, Sepeber 0, 017 Proble 1 You are ive a NMOS aplifier wih drai load resisor R = 0 k. The volae (R appeari across resisor R = 1.5 vols

More information

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions. Iera. J. Mah. & Mah. Si. Vol. 6 No. 3 (1983) 559-566 559 ASYMPTOTIC RELATIOHIPS BETWEEN TWO HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS TAKA KUSANO laculy of Sciece Hrosh llersy 1982) ABSTRACT. Some asympoic

More information

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4) 7 Differeial equaios Review Solve by he mehod of udeermied coefficies ad by he mehod of variaio of parameers (4) y y = si Soluio; we firs solve he homogeeous equaio (4) y y = 4 The correspodig characerisic

More information

Robust Dynamic Output Feedback Second-Order Sliding Mode Controller for Uncertain Systems

Robust Dynamic Output Feedback Second-Order Sliding Mode Controller for Uncertain Systems Ieraioal Joural of Corol, Auoaio, ad Syses (3 (5:878-884 DOI.7/s555--7-9 ISSN:598-6446 eissn:5-49 hp://www.spriger.co/555 Robus Dyaic Oupu Feedbac Secod-Order Slidig Mode Coroller for cerai Syses Jeag-Li

More information

Lecture 9: Polynomial Approximations

Lecture 9: Polynomial Approximations CS 70: Complexiy Theory /6/009 Lecure 9: Polyomial Approximaios Isrucor: Dieer va Melkebeek Scribe: Phil Rydzewski & Piramaayagam Arumuga Naiar Las ime, we proved ha o cosa deph circui ca evaluae he pariy

More information

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme America Joural of Compuaioal ad Applied Maemaics, (6): 77-8 DOI:.59/.acam.6. Numerical Soluio of Parabolic Volerra Iegro-Differeial Equaios via Bacward-Euler Sceme Ali Filiz Deparme of Maemaics, Ada Mederes

More information

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs America Joural of Compuaioal Mahemaics, 04, 4, 80-88 Published Olie Sepember 04 i SciRes. hp://www.scirp.org/joural/ajcm hp://dx.doi.org/0.436/ajcm.04.4404 Mea Square Coverge Fiie Differece Scheme for

More information

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form Joural of Applied Mahemaics Volume 03, Aricle ID 47585, 7 pages hp://dx.doi.org/0.55/03/47585 Research Aricle A Geeralized Noliear Sum-Differece Iequaliy of Produc Form YogZhou Qi ad Wu-Sheg Wag School

More information

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS Opimal ear Forecasig Alhough we have o meioed hem explicily so far i he course, here are geeral saisical priciples for derivig he bes liear forecas, ad

More information

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

Online Supplement to Reactive Tabu Search in a Team-Learning Problem Olie Suppleme o Reacive abu Search i a eam-learig Problem Yueli She School of Ieraioal Busiess Admiisraio, Shaghai Uiversiy of Fiace ad Ecoomics, Shaghai 00433, People s Republic of Chia, she.yueli@mail.shufe.edu.c

More information

Vibration 2-1 MENG331

Vibration 2-1 MENG331 Vibraio MENG33 Roos of Char. Eq. of DOF m,c,k sysem for λ o he splae λ, ζ ± ζ FIG..5 Dampig raios of commo maerials 3 4 T d T d / si cos B B e d d ζ ˆ ˆ d T N e B e B ζ ζ d T T w w e e e B e B ˆ ˆ ζ ζ

More information

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003 ODEs II, Suppleme o Lecures 6 & 7: The Jorda Normal Form: Solvig Auoomous, Homogeeous Liear Sysems April 2, 23 I his oe, we describe he Jorda ormal form of a marix ad use i o solve a geeral homogeeous

More information

A quadratic convergence method for the management equilibrium model

A quadratic convergence method for the management equilibrium model (IJACSA Ieraioal Joural of Advaced Copuer Sciece ad Applicaios Vol No 9 03 A quadraic covergece ehod for he aagee equilibriu odel Jiayi Zhag Feixia School Liyi Uiversiy Feixia Shadog PRChia Absrac i his

More information

Math 6710, Fall 2016 Final Exam Solutions

Math 6710, Fall 2016 Final Exam Solutions Mah 67, Fall 6 Fial Exam Soluios. Firs, a sude poied ou a suble hig: if P (X i p >, he X + + X (X + + X / ( evaluaes o / wih probabiliy p >. This is roublesome because a radom variable is supposed o be

More information

11. Adaptive Control in the Presence of Bounded Disturbances Consider MIMO systems in the form,

11. Adaptive Control in the Presence of Bounded Disturbances Consider MIMO systems in the form, Lecure 6. Adapive Corol i he Presece of Bouded Disurbaces Cosider MIMO sysems i he form, x Aref xbu x Bref ycmd (.) y Cref x operaig i he presece of a bouded ime-depede disurbace R. All he assumpios ad

More information

Dynamic Model and Analysis of Nonlinear Vibration Characteristic of a Curve-Face Gear Drive Cai, Z. Lin, C. Zhiqin Cai Chao Lin*

Dynamic Model and Analysis of Nonlinear Vibration Characteristic of a Curve-Face Gear Drive Cai, Z. Lin, C. Zhiqin Cai Chao Lin* Srojiški vesik - Joural of Mechaical Egieerig 63(07)3, 6-70 Received for review: 06-07-06 07 Joural of Mechaical Egieerig. All righs reserved. Received revised for: 07-0-05 DOI:0.5545/sv-je.06.3859 Origial

More information

An Exact Solution for the Free Vibration Analysis of Timoshenko Beams

An Exact Solution for the Free Vibration Analysis of Timoshenko Beams www.seipub.org/rap Review of Applied Physics Volue 3, 4 A Exac Soluio for he Free Vibraio Aalysis of Tiosheko Beas Raaza A. Jafari-Talookolaei *, Marya Abedi School of Mechaical Egieerig, Babol Noshirvai

More information

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i) Mah PracTes Be sure o review Lab (ad all labs) There are los of good quesios o i a) Sae he Mea Value Theorem ad draw a graph ha illusraes b) Name a impora heorem where he Mea Value Theorem was used i he

More information

The Importance of Ordering the Number of Lattice Points Inside a Rational Polyhedron Using Generating Functions

The Importance of Ordering the Number of Lattice Points Inside a Rational Polyhedron Using Generating Functions Ieraioal Joural of Copuer Sciece ad Elecroics Egieerig (IJCSEE Volue Issue ( ISSN 48 (Olie he Iporace of Orderig he Nuber of Laice ois Iside a Raioal olyhedro Usig Geeraig Fucios Halil Sopce Absrac I pure

More information

Chapter 3 Moments of a Distribution

Chapter 3 Moments of a Distribution Chaper 3 Moes of a Disribuio Epecaio We develop he epecaio operaor i ers of he Lebesgue iegral. Recall ha he Lebesgue easure λ(a) for soe se A gives he legh/area/volue of he se A. If A = (3; 7), he λ(a)

More information

Paper Introduction. ~ Modelling the Uncertainty in Recovering Articulation from Acoustics ~ Korin Richmond, Simon King, and Paul Taylor.

Paper Introduction. ~ Modelling the Uncertainty in Recovering Articulation from Acoustics ~ Korin Richmond, Simon King, and Paul Taylor. Paper Iroducio ~ Modellig he Uceraiy i Recoverig Ariculaio fro Acousics ~ Kori Richod, Sio Kig, ad Paul Taylor Tooi Toda Noveber 6, 2003 Proble Addressed i This Paper Modellig he acousic-o-ariculaory appig

More information

A Note on Random k-sat for Moderately Growing k

A Note on Random k-sat for Moderately Growing k A Noe o Radom k-sat for Moderaely Growig k Ju Liu LMIB ad School of Mahemaics ad Sysems Sciece, Beihag Uiversiy, Beijig, 100191, P.R. Chia juliu@smss.buaa.edu.c Zogsheg Gao LMIB ad School of Mahemaics

More information

SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 21 Base Excitation Shock: Classical Pulse

SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 21 Base Excitation Shock: Classical Pulse SHOCK AND VIBRAION RESPONSE SPECRA COURSE Ui 1 Base Exciaio Shock: Classical Pulse By om Irvie Email: omirvie@aol.com Iroucio Cosier a srucure subjece o a base exciaio shock pulse. Base exciaio is also

More information

Dynamic h-index: the Hirsch index in function of time

Dynamic h-index: the Hirsch index in function of time Dyamic h-idex: he Hirsch idex i fucio of ime by L. Egghe Uiversiei Hassel (UHassel), Campus Diepebeek, Agoralaa, B-3590 Diepebeek, Belgium ad Uiversiei Awerpe (UA), Campus Drie Eike, Uiversieisplei, B-260

More information

APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY

APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY ZHEN-GUO DENG ad GUO-CHENG WU 2, 3 * School of Mahemaics ad Iformaio Sciece, Guagi Uiversiy, Naig 534, PR Chia 2 Key Laboraory

More information

Transverse Vibrations of Elastic Thin Beam Resting on Variable Elastic Foundations and Subjected to Traveling Distributed Forces.

Transverse Vibrations of Elastic Thin Beam Resting on Variable Elastic Foundations and Subjected to Traveling Distributed Forces. Trasverse Vibraios of Elasic Thi Beam Resig o Variable Elasic Foudaios ad Subjeced o Travelig Disribued Forces. B. Omolofe ad S.N. Oguyebi * Deparme of Mahemaical Scieces, Federal Uiversiy of Techology,

More information

Review Exercises for Chapter 9

Review Exercises for Chapter 9 0_090R.qd //0 : PM Page 88 88 CHAPTER 9 Ifiie Series I Eercises ad, wrie a epressio for he h erm of he sequece..,., 5, 0,,,, 0,... 7,... I Eercises, mach he sequece wih is graph. [The graphs are labeled

More information

NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE

NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE Yugoslav Joural of Operaios Research 8 (2008, Number, 53-6 DOI: 02298/YUJOR080053W NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE Jeff Kuo-Jug WU, Hsui-Li

More information

6.01: Introduction to EECS I Lecture 3 February 15, 2011

6.01: Introduction to EECS I Lecture 3 February 15, 2011 6.01: Iroducio o EECS I Lecure 3 February 15, 2011 6.01: Iroducio o EECS I Sigals ad Sysems Module 1 Summary: Sofware Egieerig Focused o absracio ad modulariy i sofware egieerig. Topics: procedures, daa

More information

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION Malaysia Joural of Mahemaical Scieces 2(2): 55-6 (28) The Soluio of he Oe Species Loka-Volerra Equaio Usig Variaioal Ieraio Mehod B. Baiha, M.S.M. Noorai, I. Hashim School of Mahemaical Scieces, Uiversii

More information

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS PB Sci Bull, Series A, Vol 78, Iss 4, 2016 ISSN 1223-7027 CLOSED FORM EVALATION OF RESTRICTED SMS CONTAINING SQARES OF FIBONOMIAL COEFFICIENTS Emrah Kılıc 1, Helmu Prodiger 2 We give a sysemaic approach

More information

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3 Ieraioal Joural of Saisics ad Aalysis. ISSN 48-9959 Volume 6, Number (6, pp. -8 Research Idia Publicaios hp://www.ripublicaio.com The Populaio Mea ad is Variace i he Presece of Geocide for a Simple Birh-Deah-

More information

Supplementary Figure S1 Characterization of the two-qubit system. a-d, When the nuclear spin is polarized along the

Supplementary Figure S1 Characterization of the two-qubit system. a-d, When the nuclear spin is polarized along the Suppleeary Figure S Characeriaio of he wo-qubi syse. a-d Whe he uclear spi is polaried alog he axis is free precessio sigal abou he axis is odulaed by he relaive agle bewee he codiioal local fields. a-b

More information

Hybrid System Modeling Using Impulsive Differential Equations

Hybrid System Modeling Using Impulsive Differential Equations Hybrid Syse Modelig Usig Ipulsive Differeial Equaios ENDRA JOELIANO Depare of Egieerig Physics Isiu eologi Badug Jl Gaesha 1, Badug 413 INDONESIA Absrac: Hybrid syses have eerged as a rapidly expadig field

More information

State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by

State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by Sae-Space Model I geeral, he dyaic equaio of a luped-paraeer coiuou ye ay be repreeed by x & f x, u, y g x, u, ae equaio oupu equaio where f ad g are oliear vecor-valued fucio Uig a liearized echique,

More information

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are: 3. Iiial value problems: umerical soluio Fiie differeces - Trucaio errors, cosisecy, sabiliy ad covergece Crieria for compuaioal sabiliy Explici ad implici ime schemes Table of ime schemes Hyperbolic ad

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 Seod ad igher Order Liear Differeial Equaios Oober 9, 7 Seod ad igher Order Liear Differeial Equaios Larr areo Mehaial Egieerig 5 Seiar i Egieerig alsis Oober 9, 7 Oulie Reiew las lass ad hoewor ppl aerial

More information

Actuarial Society of India

Actuarial Society of India Acuarial Sociey of Idia EXAMINAIONS Jue 5 C4 (3) Models oal Marks - 5 Idicaive Soluio Q. (i) a) Le U deoe he process described by 3 ad V deoe he process described by 4. he 5 e 5 PU [ ] PV [ ] ( e ).538!

More information

WILD Estimating Abundance for Closed Populations with Mark-Recapture Methods

WILD Estimating Abundance for Closed Populations with Mark-Recapture Methods WILD 50 - Esiaig Abudace for Closed Populaios wih Mark-Recapure Mehods Readig: Chaper 4 of WNC book (especially secios 4. & 4.) Esiaig N is uch ore difficul ha you igh iiially expec A variey of ehods ca

More information

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE Mohia & Samaa, Vol. 1, No. II, December, 016, pp 34-49. ORIGINAL RESEARCH ARTICLE OPEN ACCESS FIED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE 1 Mohia S. *, Samaa T. K. 1 Deparme of Mahemaics, Sudhir Memorial

More information

SUMMATION OF INFINITE SERIES REVISITED

SUMMATION OF INFINITE SERIES REVISITED SUMMATION OF INFINITE SERIES REVISITED I several aricles over he las decade o his web page we have show how o sum cerai iiie series icludig he geomeric series. We wa here o eed his discussio o he geeral

More information

Clock Skew and Signal Representation

Clock Skew and Signal Representation Clock Skew ad Sigal Represeaio Ch. 7 IBM Power 4 Chip 0/7/004 08 frequecy domai Program Iroducio ad moivaio Sequeial circuis, clock imig, Basic ools for frequecy domai aalysis Fourier series sigal represeaio

More information

EGR 544 Communication Theory

EGR 544 Communication Theory EGR 544 Commuicaio heory 7. Represeaio of Digially Modulaed Sigals II Z. Aliyazicioglu Elecrical ad Compuer Egieerig Deparme Cal Poly Pomoa Represeaio of Digial Modulaio wih Memory Liear Digial Modulaio

More information

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws. Limi of a fucio.. Oe-Sided..... Ifiie limis Verical Asympoes... Calculaig Usig he Limi Laws.5 The Squeeze Theorem.6 The Precise Defiiio of a Limi......7 Coiuiy.8 Iermediae Value Theorem..9 Refereces..

More information

1 Notes on Little s Law (l = λw)

1 Notes on Little s Law (l = λw) Copyrigh c 26 by Karl Sigma Noes o Lile s Law (l λw) We cosider here a famous ad very useful law i queueig heory called Lile s Law, also kow as l λw, which assers ha he ime average umber of cusomers i

More information

Discrete-Time Signals and Systems. Introduction to Digital Signal Processing. Independent Variable. What is a Signal? What is a System?

Discrete-Time Signals and Systems. Introduction to Digital Signal Processing. Independent Variable. What is a Signal? What is a System? Discree-Time Sigals ad Sysems Iroducio o Digial Sigal Processig Professor Deepa Kudur Uiversiy of Toroo Referece: Secios. -.4 of Joh G. Proakis ad Dimiris G. Maolakis, Digial Sigal Processig: Priciples,

More information

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

Available online at   J. Math. Comput. Sci. 4 (2014), No. 4, ISSN: Available olie a hp://sci.org J. Mah. Compu. Sci. 4 (2014), No. 4, 716-727 ISSN: 1927-5307 ON ITERATIVE TECHNIQUES FOR NUMERICAL SOLUTIONS OF LINEAR AND NONLINEAR DIFFERENTIAL EQUATIONS S.O. EDEKI *, A.A.

More information

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM Available a hp://pvamu.edu/aam Appl. Appl. Mah. ISSN: 193-9466 Vol. 5, No. Issue (December 1), pp. 575 584 (Previously, Vol. 5, Issue 1, pp. 167 1681) Applicaios ad Applied Mahemaics: A Ieraioal Joural

More information

Extremal graph theory II: K t and K t,t

Extremal graph theory II: K t and K t,t Exremal graph heory II: K ad K, Lecure Graph Theory 06 EPFL Frak de Zeeuw I his lecure, we geeralize he wo mai heorems from he las lecure, from riagles K 3 o complee graphs K, ad from squares K, o complee

More information

[ m] x = 0.25cos 20 t sin 20 t m

[ m] x = 0.25cos 20 t sin 20 t m . x.si ( 5 s [ ] CHAPER OSCILLAIONS x ax (.( ( 5 6. s s ( ( ( xax. 5.7 s s. x.si [] x. cos s Whe, x a x.5. s 5s.6 s x. x( x cos + si a f ( ( [ ] x.5cos +.59si. ( ( cos α β cosαcos β + siαsi β x Acos φ

More information

INVESTMENT PROJECT EFFICIENCY EVALUATION

INVESTMENT PROJECT EFFICIENCY EVALUATION 368 Miljeko Crjac Domiika Crjac INVESTMENT PROJECT EFFICIENCY EVALUATION Miljeko Crjac Professor Faculy of Ecoomics Drsc Domiika Crjac Faculy of Elecrical Egieerig Osijek Summary Fiacial efficiecy of ivesme

More information

Fresnel Dragging Explained

Fresnel Dragging Explained Fresel Draggig Explaied 07/05/008 Decla Traill Decla@espace.e.au The Fresel Draggig Coefficie required o explai he resul of he Fizeau experime ca be easily explaied by usig he priciples of Eergy Field

More information

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY U.P.B. Sci. Bull., Series A, Vol. 78, Iss. 2, 206 ISSN 223-7027 A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY İbrahim Çaak I his paper we obai a Tauberia codiio i erms of he weighed classical

More information

Applications of force vibration. Rotating unbalance Base excitation Vibration measurement devices

Applications of force vibration. Rotating unbalance Base excitation Vibration measurement devices Applicaios of foce viaio Roaig ualace Base exciaio Viaio easuee devices Roaig ualace 1 Roaig ualace: Viaio caused y iegulaiies i he disiuio of he ass i he oaig copoe. Roaig ualace 0 FBD 1 FBD x x 0 e 0

More information

Solutions Manual 4.1. nonlinear. 4.2 The Fourier Series is: and the fundamental frequency is ω 2π

Solutions Manual 4.1. nonlinear. 4.2 The Fourier Series is: and the fundamental frequency is ω 2π Soluios Maual. (a) (b) (c) (d) (e) (f) (g) liear oliear liear liear oliear oliear liear. The Fourier Series is: F () 5si( ) ad he fudameal frequecy is ω f ----- H z.3 Sice V rms V ad f 6Hz, he Fourier

More information

Improvement Over General And Wider Class of Estimators Using Ranked Set Sampling

Improvement Over General And Wider Class of Estimators Using Ranked Set Sampling ITERATIOAL JOURAL OF SIETIFI & TEOLOG RESEAR VOLUME ISSUE 7 AUGUST ISS 77-866 Iprovee Over Geeral Ad ider lass of Esiaors Usi Raked Se Sapli V L Madoara iu Meha Raka Absrac: I his paper Iprovee over eeral

More information

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems Ausralia Joural of Basic ad Applied Scieces, 4(1): 518-57, 1 ISSN 1991-8178 Homoopy Aalysis Mehod for Solvig Fracioal Surm-Liouville Problems 1 A Neamay, R Darzi, A Dabbaghia 1 Deparme of Mahemaics, Uiversiy

More information

Effects of Forces Applied in the Middle Plane on Bending of Medium-Thickness Band

Effects of Forces Applied in the Middle Plane on Bending of Medium-Thickness Band MATEC We of Cofereces 7 7 OI:./ maeccof/77 XXVI R-S-P Semiar 7 Theoreical Foudaio of Civil Egieerig Effecs of Forces Applied i he Middle Plae o Bedig of Medium-Thickess Bad Adre Leoev * Moscow sae uiversi

More information

Recursive Identification of Hammerstein Systems with Polynomial Function Approximation

Recursive Identification of Hammerstein Systems with Polynomial Function Approximation Ieraioal Joural of Maagee ad Fuzzy Syses 7; 3(6): 87-94 hp://www.sciecepublishiggroup.co//ifs doi:.648/.ifs.736. ISSN: 575-4939 (Pri); ISSN: 575-4947 (Olie) Research/echical Noe Recursive Ideificaio of

More information

A Note on Prediction with Misspecified Models

A Note on Prediction with Misspecified Models ITB J. Sci., Vol. 44 A, No. 3,, 7-9 7 A Noe o Predicio wih Misspecified Models Khresha Syuhada Saisics Research Divisio, Faculy of Mahemaics ad Naural Scieces, Isiu Tekologi Badug, Jala Gaesa Badug, Jawa

More information

On stability of first order linear impulsive differential equations

On stability of first order linear impulsive differential equations Ieraioal Joural of aisics ad Applied Mahemaics 218; 3(3): 231-236 IN: 2456-1452 Mahs 218; 3(3): 231-236 218 as & Mahs www.mahsoural.com Received: 18-3-218 Acceped: 22-4-218 IM Esuabaa Deparme of Mahemaics,

More information

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models Oulie Parameer esimaio for discree idde Markov models Juko Murakami () ad Tomas Taylor (2). Vicoria Uiversiy of Welligo 2. Arizoa Sae Uiversiy Descripio of simple idde Markov models Maximum likeliood esimae

More information

12 Getting Started With Fourier Analysis

12 Getting Started With Fourier Analysis Commuicaios Egieerig MSc - Prelimiary Readig Geig Sared Wih Fourier Aalysis Fourier aalysis is cocered wih he represeaio of sigals i erms of he sums of sie, cosie or complex oscillaio waveforms. We ll

More information

Manipulations involving the signal amplitude (dependent variable).

Manipulations involving the signal amplitude (dependent variable). Oulie Maipulaio of discree ime sigals: Maipulaios ivolvig he idepede variable : Shifed i ime Operaios. Foldig, reflecio or ime reversal. Time Scalig. Maipulaios ivolvig he sigal ampliude (depede variable).

More information

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming* The Eighh Ieraioal Symposium o Operaios esearch ad Is Applicaios (ISOA 9) Zhagjiajie Chia Sepember 2 22 29 Copyrigh 29 OSC & APOC pp 33 39 Some Properies of Semi-E-Covex Fucio ad Semi-E-Covex Programmig*

More information

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi Liear lgebra Lecure #9 Noes This week s lecure focuses o wha migh be called he srucural aalysis of liear rasformaios Wha are he irisic properies of a liear rasformaio? re here ay fixed direcios? The discussio

More information

A Complex Neural Network Algorithm for Computing the Largest Real Part Eigenvalue and the corresponding Eigenvector of a Real Matrix

A Complex Neural Network Algorithm for Computing the Largest Real Part Eigenvalue and the corresponding Eigenvector of a Real Matrix 4h Ieraioal Coferece o Sesors, Mecharoics ad Auomaio (ICSMA 06) A Complex Neural Newor Algorihm for Compuig he Larges eal Par Eigevalue ad he correspodig Eigevecor of a eal Marix HANG AN, a, XUESONG LIANG,

More information

Department of Mathematical and Statistical Sciences University of Alberta

Department of Mathematical and Statistical Sciences University of Alberta MATH 4 (R) Wier 008 Iermediae Calculus I Soluios o Problem Se # Due: Friday Jauary 8, 008 Deparme of Mahemaical ad Saisical Scieces Uiversiy of Albera Quesio. [Sec.., #] Fid a formula for he geeral erm

More information

Harmonic excitation (damped)

Harmonic excitation (damped) Harmoic eciaio damped k m cos EOM: m&& c& k cos c && ζ & f cos The respose soluio ca be separaed io par;. Homogeeous soluio h. Paricular soluio p h p & ζ & && ζ & f cos Homogeeous soluio Homogeeous soluio

More information

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu Caper : Time-Domai Represeaios of Liear Time-Ivaria Sysems Ci-Wei Liu Oulie Iroucio Te Covoluio Sum Covoluio Sum Evaluaio Proceure Te Covoluio Iegral Covoluio Iegral Evaluaio Proceure Iercoecios of LTI

More information

N! AND THE GAMMA FUNCTION

N! AND THE GAMMA FUNCTION N! AND THE GAMMA FUNCTION Cosider he produc of he firs posiive iegers- 3 4 5 6 (-) =! Oe calls his produc he facorial ad has ha produc of he firs five iegers equals 5!=0. Direcly relaed o he discree! fucio

More information

INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA

INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA Volume 8 No. 8, 45-54 ISSN: 34-3395 (o-lie versio) url: hp://www.ijpam.eu ijpam.eu INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA A.Arul dass M.Dhaapal

More information

A Two-Level Quantum Analysis of ERP Data for Mock-Interrogation Trials. Michael Schillaci Jennifer Vendemia Robert Buzan Eric Green

A Two-Level Quantum Analysis of ERP Data for Mock-Interrogation Trials. Michael Schillaci Jennifer Vendemia Robert Buzan Eric Green A Two-Level Quaum Aalysis of ERP Daa for Mock-Ierrogaio Trials Michael Schillaci Jeifer Vedemia Rober Buza Eric Gree Oulie Experimeal Paradigm 4 Low Workload; Sigle Sessio; 39 8 High Workload; Muliple

More information

The Eigen Function of Linear Systems

The Eigen Function of Linear Systems 1/25/211 The Eige Fucio of Liear Sysems.doc 1/7 The Eige Fucio of Liear Sysems Recall ha ha we ca express (expad) a ime-limied sigal wih a weighed summaio of basis fucios: v ( ) a ψ ( ) = where v ( ) =

More information