Dynamic Ping Optimization for Surveillance in Multistatic Sonar Buoy Networks with Energy Constraints

Size: px
Start display at page:

Download "Dynamic Ping Optimization for Surveillance in Multistatic Sonar Buoy Networks with Energy Constraints"

Transcription

1 Dynamic Ping Optimizatin fr Surveillance in Multistatic Snar Buy Netwrks with Energy Cnstraints Anshu Saksena and I-Jeng Wang Abstract In this paper we study the prblem f dynamic ptimizatin f ping schedule in an active snar buy netwrk deplyed t prvide persistent surveillance f a littral area thrugh multistatic detectin The gal f ping scheduling is t dynamically determine when t ping and which ping surce t engage in rder t achieve the desirable detectin perfrmance Fr applicatins where persistent surveillance is needed fr an extended perid f time it is expected that the energy available at each ping surce is limited relative t the required system lifetime Hence efficient management f pwer cnsumptin fr pinging is imprtant t supprt the required lifetime f the netwrk while maintaining acceptable detectin perfrmance Our apprach t ping ptimizatin is based n the applicatin f apprximate Partially Observable Markv Decisin Prcess (POMDP) techniques such as the rllut algrithm T enable a practical implementatin f the plicy rllut we apply sampling-based techniques based n a simplified mdel that apprximates the detailed multistatic mdel Using high fidelity snar simulatins we evaluate the perfrmance f the prpsed apprach and cmpare it with the greedy technique in terms f detectin perfrmance and system lifetime I INTRODUCTION In this paper we derive a frmulatin f the dynamic ping ptimizatin prblem in multistatic snar buy netwrks using the framewrk f Partially Observable Markv Decisin Prcesses (POMDP) The basic perating cncept f a multistatic snar buy netwrk fr underwater surveillance is t practively ping frm an acustic surce and crrelate the ech returns (specifically time f arrivals and bearings) acrss receivers in the field t detect lcalize and track targets f interest In rder t prvide necessary cverage fr the surveillance area and sufficient diversity in multistatic gemetry fr target lcalizatin multiple receivers and acustic surces are deplyed and distributed thrughut the field The ping ptimizatin prblem cnsidered here is a specific example f sensr management in sensr netwrks which has received significant interest in recent years in generic sensr netwrk settings (see fr example [1] [2] [3]) Our wrk fcuses n the practical applicatin f existing algrithms (in particular the apprximate dynamic prgramming techniques prpsed in [3] [4]) t multistatic snar netwrks where sensing and detectin mdels are mre cmplicated than in settings cnsidered in the literature This wrk was supprted by the ONR Discvery and Inventin Prgram Anshu Saksena is with the Miltn Eisenhwer Research Center Jhns Hpkins University Applied Physics Labratry Laurel MD USA anshusaksena@jhuapledu I-Jeng Wang is with the Miltn Eisenhwer Research Center Jhns Hpkins University Applied Physics Labratry Laurel MD USA i-jengwang@jhuapledu The benefit f ptimizing ping sequences t tracking perfrmance fr multistatic snar systems has been demnstrated thrugh simulatins in [5] where a greedy apprach is presented t ptimize the expected detectin prbability during the search phase when n detectins are cnfirmed and t imprve the track with wrst quality after tracks are established In this paper we fcus n the aspects f ping surce selectin with additinal energy cnstraint during the surveillance phase The surveillance phase is defined as the stage where n specific track with sufficient cnfidence has been established based n the detectin results frm the netwrk f receivers In cntrast t [5] where cmplete absence f cnfirmed detectins is assumed we permit the existence f initial detectins in ur frmulatin The ratinale fr this unique feature is tw-fld: In practical applicatins sparse detectins are expected befre sufficient measurements are available t establish a track (the track initiatin prblem) Furthermre an apprach assuming n detectin has limited applicability since false alarms are inevitable fr mst practical scenaris The separatin f surveillance and tactical (after tracks are established) phases pstulated by existing literature is smewhat artificial In reality bth situatins exist simultaneusly since cntinuing search f new targets is desirable even after tracks are established A frmulatin incrprating detectins during the surveillance phase (such as the ne prpsed here) can likely be extended t address requirements fr bth search and tracking simultaneusly In Sectin II we frmally define all the necessary attributes f a POMDP fr the ping ptimizatin prblem Based n the frmulatin we discuss in Sectin III hw the plicy rllut technique can be applied t sequentially slve fr ping surce selectins by taking a nnmypic view f perfrmance ver a time hrizn int the near future This apprach serves t trade-ff the detectin perfrmance versus energy management requirement We present simulatin results t demnstrate the perfrmance f prpsed algrithms in Sectin IV Thrughut the discussin we assume that there are N s surces and N r receivers deplyed in the field and their psitins are statinary II MATHEMATICAL FORMULATION OF PING OPTIMIZATION PROBLEMS We assume a discrete time mdel and use k = t dente the kth stage with 0th stage denting the starting time We will use time k t dente the beginning f the

2 Sensr detectins Snar Buy Netwrk is the size f the grid The surveillance target state at time k is defined by Fusin Fused Measurements Fig 1 Observable States State Estimatr Psterir distributin f unbservables Cntrller Optimizatin Ping Decisin An illustratin f the POMDP fr ping ptimizatin kth stage Fllwing the standard MDP framewrk we will define the states X(k) the actin a(k) the bservatin Z(k) and its dependency n states and actins the state transitin the immediate reward r(x(k) a(k)) the hrizn H and the terminal cst r T (X(H)) if needed Given the initial state X(0) (r its distributin) the ptimizatin prblem defined by the MDP is given by V H(X(0)) = max E[r (X(0) a(0)) + a(0)a(h 1) + r (X(H 1) a(h 1)) + r T (X(H))] (1) where a(k) 0 N s } is the ping surce selectin at time k with a(k) = 0 crrespnding t the actin f nt pinging frm any surce and the maximizatin is taken ver all feasible sequences f actins Bellman s equatin prvides a characterizatin f the perfrmance f the ptimal decisins VH(X) = max Q(X a) a Q(X a) r(x a) + E X [VH 1(X )] where X is the next state after taking the actin a at state X and Q(X a) is typically referred t as the Q-value f actin a at state X Given that the states are nt fully bservable we cast the prblem as a POMDP using the cncept f the belief state r the cnditinal prbability The basic idea f the POMDP is illustrated in Figure 1 A States and Observatins The states f the ping ptimizatin can be classified int tw types f states: the target states that capture the dynamic behavirs f targets; and the sensr states that characterize the status f the sensr systems 1) Target states and bservatin: Fr the surveillance case where n credible tracks have been established it is difficult t mdel all the uncertainty including the number f targets in the field Instead f directly mdelling the target states we characterize the ptential target presence acrss the field We further discretize the field t arrive at a mre tractable frmulatin Specifically we assume that the field is discretized int a grid (x i y i ) i = 1 N g } where N g where t i (k) = X(k) = [ t 1 (k) t Ng (k) ] T 1 if there is at least 1 target at (x i y i ) at time k 0 therwise Nte that in this mdel we have made tw apprximating assumptins: (1) target lcalizatins are mapped t the grid; and (2) we make n differentiatin in ur mdels between the cases f a single target and multiple targets being at a grid pint Since n specific target tracking infrmatin has been established during the surveillance phase a randmized mdel derived frm a stchastic spatial diffusin prcess determines the randm mvement and arrival f ptential targets We define the state transitin prbability by assuming that the target dynamics are independent f each ther and nt affected by the actin taken: P (X(k + 1) X(k) a(k)) = i P (t i (k + 1) X(k)) where the individual target transitin mdel is given by a Gaussian distributin The target states are nt directly bservable but can be estimated frm the bservatin: where z i (k) = Z(k) = [ z 1 (k) z Ng (k) ] T 1 if a detectin ccurs at (x i y i ) after a(k) 0 therwise The bservatin mdel is specified by the cnditinal prbability P (Z(k) X(k) a(k)) If we assume that the detectin result at the grid pint (x i y i ) with actin a(k) is independent f bth detectin z j (k) and target presence t j (k) at ther grid pints (x j y j ) with j i at time k then P (Z(k) X(k) a(k)) = i P (z i (k) t i (k) a(k)) and P (z i (k) = 1 t i (k) = 1 a(k)) = p d (x i y i a(k)) (2) P (z i (k) = 1 t i (k) = 0 a(k)) = p fa (x i y i a(k)) (3) 2) Sensr states: In the initial frmulatin we assume that all the sensrs are statinary and fcus n mdelling the energy levels at the surces We assume a simple mdel fr the energy states: E(k) = [e 1 (k) e Ns (k)] T e i (k) 0} N where e i (k) is the remaining number f pings the surce i can supprt in the beginning f the time k befre an actin

3 a(k) is taken Assuming that the initial energy state E(0) is given the state transitin is defined by e i (k) 1 if a(k) = i and e i (k) 1 e i (k + 1) = e i (k) therwise In this simple mdel we assume that the energy cnsumptin fr pinging is deterministic Furthermre we assume that the energy state E(k) is directly bservable (that is can be measured withut errrs) B Belief states and estimatin T apply the POMDP technique we define the belief states and derive a recursin that enables us t estimate the belief states sequentially based n the bservatins Z(k) The belief state is defined as the cnditinal prbability f the states given the infrmatin state (all prir bservatins and actins) that is P X(k) Z(0) Z(1) Z(k 1) a(0) a(k 1)} We dente the belief state by P T (k) = [ P T (x 1 y 1 ; k) P T (x 2 y 2 ; k) P T (x Ng y Ng ; k) ] T where P T (x i y i ; k) dentes the cnditinal prbability f target presence at (x i y i ) at time k given all prir bservatins and actins The belief state P T (k) can be estimated and prpagated sequentially by a tw-step recursin Given the current belief state P T (k): 1) Bayes update: Given an actin a(k) and the resulting bservatin Z(k) cmpute the a psteriri belief state P T Z using Bayes rule P T Z(k) (x i y i ; k) = p d P T p d P T +p fa [1 P T ] if z i (k) = 1 (1 p d )P T [1 p d ]P T +[1 p fa ][1 P T ] therwise where p d p fa and P T have been written instead f p d (x i y i a(k)) p fa (x i y i a(k)) and P T (x i y i ; k) respectively fr readability 2) Diffusin: Prpagate the a psteriri belief state accrding t the transitin mdel P T (x i y i ; k+1) = 1 j ( 1 PT Z(k) (x j y j ; k)p ij ) The fllwing diagram illustrates the prpagatin f belief state: P T (k) a(k) Bayes update Z(k) P T Z(k) (k) diffusin P T (k+1) C Feasible actins and plicies An actin that meets all the cnstraints defined in the prblem is called a feasible actin An bvius cnstraint fr the ping ptimizatin prblem is that the energy level at a surce is sufficient t ping That is a i is feasible nly if e i 1 in ur frmulatin Nte that the feasibility f an actin will depend n the current state P T (k) E(k)} P T (k) E(k) } Candidate Actin a Z 1 (k) Z M (k) Fig 2 State Estimatin Mnte Carl Sampling P 1 T (k+1) P M T (k+1) E 1 (k) E M (k) Base Plicy π! a (k+1) An illustratin f the plicy rllut Q 1 Q M We dente the set f feasible actins at time k by A k 0 N s } A (Markv) plicy π is a mapping frm the states P T (k) E(k)} t an actin a(k) A plicy is feasible if π (P T (k) E(k)) A k fr all k = 0 H 1 Tw typical examples f plicy fr the ping ptimizatin prblem are: Fixed sequence plicies: One can define a feasible plicy by repeating the same sequence f ping decisins as lng as the remaining energy level at a surce is sufficient t ping The rund-rbin algrithm falls under this class f plicies Greedy plicies: The greedy plicy is a plicy that selects a decisin amng the feasible actins by maximizing the immediate reward That is π (P T (k) E(k)) = argmax a Ak E [r(p T (k) E(k) a)] III OPTIMIZATION CRITERIA AND POLICY ROLLOUT A Rewards and dynamic prgramming frmulatin Fr detectin perfrmance we define the reward functin using the average detectin prbability acrss the field: r(p T a) = 1 N g PT α (x i y i )p d (x i y i a) (4) N g i=1 where α 1 is a parameter that allws us t cntrl the search strategy T mdel the energy cnstraint we assume a simplified mdel where each ping surce has a pre-specified number f pings thrughut its lifetime Even thugh we d nt explicitly mdel energy related cnsideratins as part f the reward in POMDP it is expected that the energy cnstraint will play a rle in the perfrmance f POMDP if the selected hrizn is lng enugh relative t the available ping resurces With the reward and cnstraint defined abve we frmulate ping ptimizatin as the fllwing stchastic dynamic prgramming prblem: [ H 1 ] max E r (P T (k) a(k)) (5) a(0)a(h 1);a(k) A k k=0 where r is defined in (4) Fr ur experiments we use the frmulatin (5) in a receding hrizn setting The receding

4 hrizn setting can ptentially reduce the impact f mdel inaccuracy n the predictin f future detectin perfrmance in dynamic prgramming B Plicy Rllut Technique T address the cmplexity f the dynamic prgramming frmulated fr the ping ptimizatin prblem we apply the plicy rllut technique [9] The basic idea f the plicy rllut algrithm is t use the perfrmance resulting frm a base plicy π t estimate the Q-value fr making the decisin Given the cmplexity f the underlying mdels it is difficult t directly estimate the perfrmance f a base plicy Hence we apply a sampling-based apprach similar t the nes used in [3] [4] t this estimatin prblem The sampling-based apprach relies n the utcme f multiple Mnte Carl simulatins with the base plicy t estimate the Q-value f each feasible actin Cnsider the traces f states frm M independent Mnte Carl runs by taking an actin a and fllwing a feasible base plicy π afterward: P T (k) E(k)} a Z i (k) P i T (k + 1) E i (k + 1)} a π (k+1) P i T (k + 2) E i (k + 2)} a π (H 1) P i T (H) E i (H)} i = 1 M Nte that the specific actins taken fllwing the base plicy will generally differ amng the traces since they depend n the specific belief and energy state trajectry P i T Ei } btained in each Mnte Carl run Frm each sample path we btain a sample f the Q-value fr the base plicy π as Q i (P T (k) a) = r (P T (k) a) + r ( P i T (k + 1) a π (k + 1) ) + + r ( P i T (H 1) a π (H 1) ) An estimate f the Q-value fr the actin a at the belief state P T (k) dented by Q (P T (k) a) can be btained as the empirical average ver these Mnte Carl runs Fllwing the plicy rllut strategy the actin at time k is chsen based n a (k) = argmax a Ak Q (PT (k) a) In this paper we cnsider the fllwing apprximatin mdel fr simulatins f Z(k) in plicy rllut: z i (k) = 1 with prbability p d P T + p fa (1 P T ) where P T p d and p fa are again used as shrthand fr P T (x i y i ; k) p d (x i y i a(k)) and p fa (x i y i a(k)) respectively fr the sake f readability IV HIGH FIDELITY SIMULATION RESULTS A Snar Mdels In rder t test perfrmance ver a range f envirnmental settings and cnditins numerical mdeling was used t simulate the mtin f a target thrugh a distributed TABLE I TOTAL PROBABILITY OF DETECTION FOR TARGET 1 Mean Std dev Num scenaris Rund-rbin Greedy Rllut H= Rllut H= Rllut H= TABLE II TOTAL PROBABILITY OF DETECTION FOR TARGET 2 Mean Std dev Num scenaris Rund-rbin Greedy Rllut H= Rllut H= Rllut H= sensr netwrk The signal excess fr eches assciated with surce-receiver pairs was determined fr each event in a ping schedule using a physics-based mdel Clutter was established statistically frm a measured characteristic distributin Transmissin lss calculatins are perfrmed using APLNM a range-dependent cupled nrmal mde cde Sund velcity prfiles bathymetry and bttm prperties are extracted frm standard databases Backgrund interference was calculated using a pwer sum f an assumed spectral ambient nise level (NL) and a calculated spectral reverberatin level (RL) The RL was calculated by assuming a cnstant sund speed with an equitime ellipse f reverberant patches (RV B) cnstrained by the signal bandwidth The receiver beam pattern is included directly in establishing RL with a single directivity index applied t the N L Bth N L and RV B are determined in units f intensity and are assumed cnstant ver the bandwidth f the signal Ptential clutter events are established using a unifrm density f events fr each receiver beam (eg ne event per secnd) These events are drawn randmly frm a generalized gamma distributin having mean standard deviatin and skewness parameters fit t represent measured clutter statistics frm at-sea experiments B Evaluatin with Simulated Data Our simulatins are based n a snbuy field cnsisting f 16 receivers arranged in a 4x4 grid and 4 active snar surces lcated within the grid f receivers as shwn in Figure 3 Each surce begins with enugh energy t ping 150 times The first target enters the field at stage 15 and TABLE III SYSTEM LIFETIME Mean Min Max Num scenaris Rund-rbin Greedy Rllut H=

5 the secnd enters at stage 250 We cmpare detectin perfrmance f three surce selectin plicies: the rund-rbin plicy the greedy plicy and plicy rllut using the greedy plicy as the base plicy Plicy rllut ptimized the decisin at each stage based n the average f up t 50 Mnte Carl simulatins f the bservatins At each stage there are 5 pssible actins: ping frm ne f the 4 surces r d nt ping frm any surce We measure surveillance perfrmance by summing the prbability f detecting each target ver each simulatin run A number f different target track scenaris were run using rund-rbin greedy and plicy rllut strategies with their results summarized in Table I fr the first target entering the field and in Table II fr the secnd Plicy rllut was run with hrizn lengths f and 20 stages with the results fr each listed in the table separately Since target 1 is tracked with sufficient energy reserves predictin accuracy ver the hrizn is gverned by the uninfrmed diffusin target mtin mdel the lnger hrizns seem t have verextended the relevance f thse predictins Since target 2 arrived later the end game f limited resurces becmes mre f a factr Lnger hrizns may be beneficial in this case since predictin f energy depletin is accurate even ver lng hrizns While these trends seem reasnable the statistics f these results preclude us frm cncluding that any f these plicies perfrms better r wrse than any ther The system lifetime defined as the time after which all surces have exhausted their energy is shwn in Table III With the exceptin f the first stage when there are n targets believed t be in the area and ur implementatin selects nt t ping the greedy plicy pings at every stage until surce energies are exhausted at stage 601 The rundrbin plicy als pings at every stage and thus als has a lifetime f 600 Hwever lnger lifetimes result frm plicy rllut s anticipatin f the scarcity and subsequent ratining f energy Hwever there was ne scenari that Latitude Time 300 Time 600 Time Time 900 Time 1200 Lngitude Fig 3 Simulated scenari: asterisks indicate acustic surces circles receivers lines shw target tracks and circles indicate the time at which the target is at that pint f the track Targets enter the field frm the Nrthwest turn East and leave Surce 1 Surce 2 Surce 3 Surce Fig 4 Plt f the energy reserves in the number f pings remaining n each surce as a functin f time under the greedy surce selectin plicy Surce 1 Surce 2 Surce 3 Surce Fig 5 Plt f the energy reserves in the number f pings remaining n each surce as a functin f time under the rllut surce selectin plicy with hrizn length f 20 never exhausted its energy fr any f the plicy rllut hrizn lengths This large lifetime dminates the mean f each algrithm but since the number f scenaris increases with increasing hrizn length the effect f the utlying scenari decreases creating an apparent trend f decreasing lifetime Remving this utlier frm each data set prduces the results in Table?? The p-values shwn are fr cmparisn t the lifetime f the greedy algrithm Of curse t be fair an extreme value n the ther end f the distributin shuld als be remved but that wuld nly imprve the statistics ver thse shwn in Table?? Therefre each implementatin f plicy rllut prduces a significantly lnger system lifetime than greedy r rund-rbin a key result fr persistent surveillance withut sacrificing detectin perfrmance significantly We nw briefly present mre detailed data frm ne scenari t illustrate the perfrmance f plicy rllut with hrizn length 20 versus the greedy algrithm Figure 3 shws an example f the simulated tracks that was used t generate the plts presented here Figure 4 shws the energy reserves f each surce as a functin f time fr the greedy plicy and Figure 5 shws this fr plicy rllut giving an idea f hw judiciusly each plicy makes use f each surce

6 V CONCLUSIONS AND FUTURE WORK In this paper we present a mathematical frmulatin f the dynamic ping ptimizatin prblem fr snar buy netwrks with energy cnstraints We develp an implementatin f the sampling-based plicy rllut t ptimize the average detectin perfrmance Even with a much simplified mdel used in plicy rllut the apprximate dynamic prgramming apprach leads t a lnger system lifetime while achieving detectin perfrmance similar t the greedy apprach as illustrated by the experimental results using high-fidelity snar simulatins As an immediate future wrk we will cnsider applying state aggregatin techniques t further imprve the efficiency f plicy rllut We address the energy cnstraint by defining a hrizn lng enugh in stchastic dynamic prgramming One might cnsider explicitly mdelling the energy cnsideratin t enable a receding hrizn frmulatin with shrter hrizn (and hence reducing the cmputatinal csts f plicy rllut) Theretically the mst apprpriate metric fr energy management is the lifetime f the system Hwever directly using the lifetime as the reward will lead t a difficult POMDP A gd alternative is t use a lad-balancing type f metric t capture the fairness f energy cnsumptin acrss the field Given the energy levels E = [e 1 e Ns ] T define the deviatin f energy level at surce i as σi e (E) = e i ē ē where ē is the average amng e i } We can then use the energy deviatins t either define the reward r define a cnstraint N s r(p T E a) = σi e (E ) (6) i=1 σ e i (E ) ρ(e) < 1 fr all i = 1 N s (7) where E is the resulting energy levels after the actin a is taken given current energy levels E and ρ(e) is defined t ensure the existence f a feasible actin while maintaining a desirable level f lad balancing Finally since we cnsider the detectins during the surveillance in ur frmulatin it is straightfrward t incrprate the utcme f tracking algrithms int ur frmulatin with an apprpriate expansin f the state space t include the current target curse and speed etc In essence the state estimatin prcess we discussed in Sectin II-B is a simple Bayesian tracking algrithm with the prbabilistic diffusin as its underlying dynamic mdel Amng the existing tracking appraches the Prbability Hypthesis Density (PHD) filter [10] [11] is perhaps the mst natural ne t cnsider The PHD filter recursins prpagate the intensity f the target states in a way that is very similiar t the prpagatin f the belief state cnsidered in ur frmulatin A key challenge in incrprating tracking utcmes in ping ptimizatin is t devise a perfrmance metric that adequately addresses the peratinal pririty regarding the trade-ff between cntinuing tracking f cnfirmed targets and persistent surveillance f the entire field One pssible apprach is t experiment with the chice f parameter α in the average detectin reward defined in (4) t cntrl the preference f future detectin n lcatins with higher target presence belief where established tracks wuld be VI ACKNOWLEDGEMENTS The authrs gratefully acknwledge the wrk f Jeff Dunne in develping the snar simulatins and the cntributins f Dennis Lucarelli t the develpment f the frmulatin f ur prblem Bth are f the Jhns Hpkins University Applied Physics Labratry REFERENCES [1] V Krishnamurthy Algrithms fr ptimal scheduling f hidden Markv mdel sensrs IEEE Transactins n Signal Prcessing vl 50 n 6 pp [2] C M Kreucher A O Her III K D Kastella and M R Mrelande An infrmatin-based apprach t sensr management in large dynamic netwrks Prceedings f the IEEE vl 95 n 5 May 2007 [3] Y Li L W Krakw E K P Chng and K N Grm Apprximate stchastic dynamic prgramming fr sensr scheduling t track multiple targets Digital Signal Prcessing t appear [4] E K P Chng C Kreucher and A O Her III POMDP apprximatin methds based n heuristics and simulatin Chapter 8 in Fundatins and Applicatins f Sensr Management A O Her D Castann D Cchran and K Kastella Eds Chapter 5 Springer 2008 ISBN: pp [5] D W Krut M A El-Sharkawi W J L Fx and M U Hazen Intelligent ping sequencing fr multistatic snar systems th Internatinal Cnference n Infrmatin Fusin 2006 [6] W A Albersheim A Clsed-Frm Apprximatin t Rbertsn s Detectin Characteristics Prceedings f the IEEE Vl 69 N 7 July 1981 [7] A Saksena L Benmhamed J Dunne D Lucarelli I-J Wang Imprving System-wide Detectin Perfrmance fr Snar Buy Netwrks using In-Netwrk Fusin MILCOM 2007 Unclassified Prceedings Orland FL [8] A C Cn Spatial crelatin f detectins fr impulsive ech ranging snar Jhns Hpkins APL Technical Digest vl 18 n1 pp [9] D P Bertsekas and D A Castañn Rllut algrithms fr stchastic scheduling prblems Jurnal f Heuristics vl 5pp [10] R Mahler Multi-target Bayes filtering via first-rder multi-target mments IEEE Trans Aersp Electrn Syst vl 39 n 4 pp [11] B-N V and W-K Ma The Gaussian mixture prbability hypthesis denisty filter IEEE Transactins n Signal Prcessing vl 54 n11 pp

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) > Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);

More information

You need to be able to define the following terms and answer basic questions about them:

You need to be able to define the following terms and answer basic questions about them: CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f

More information

SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical model for microarray data analysis

SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical model for microarray data analysis SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical mdel fr micrarray data analysis David Rssell Department f Bistatistics M.D. Andersn Cancer Center, Hustn, TX 77030, USA rsselldavid@gmail.cm

More information

the results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must

the results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must M.E. Aggune, M.J. Dambrg, M.A. El-Sharkawi, R.J. Marks II and L.E. Atlas, "Dynamic and static security assessment f pwer systems using artificial neural netwrks", Prceedings f the NSF Wrkshp n Applicatins

More information

Least Squares Optimal Filtering with Multirate Observations

Least Squares Optimal Filtering with Multirate Observations Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical

More information

Chapter 3: Cluster Analysis

Chapter 3: Cluster Analysis Chapter 3: Cluster Analysis } 3.1 Basic Cncepts f Clustering 3.1.1 Cluster Analysis 3.1. Clustering Categries } 3. Partitining Methds 3..1 The principle 3.. K-Means Methd 3..3 K-Medids Methd 3..4 CLARA

More information

, which yields. where z1. and z2

, which yields. where z1. and z2 The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin

More information

Sequential Allocation with Minimal Switching

Sequential Allocation with Minimal Switching In Cmputing Science and Statistics 28 (1996), pp. 567 572 Sequential Allcatin with Minimal Switching Quentin F. Stut 1 Janis Hardwick 1 EECS Dept., University f Michigan Statistics Dept., Purdue University

More information

Methods for Determination of Mean Speckle Size in Simulated Speckle Pattern

Methods for Determination of Mean Speckle Size in Simulated Speckle Pattern 0.478/msr-04-004 MEASUREMENT SCENCE REVEW, Vlume 4, N. 3, 04 Methds fr Determinatin f Mean Speckle Size in Simulated Speckle Pattern. Hamarvá, P. Šmíd, P. Hrváth, M. Hrabvský nstitute f Physics f the Academy

More information

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur Cdewrd Distributin fr Frequency Sensitive Cmpetitive Learning with One Dimensinal Input Data Aristides S. Galanpuls and Stanley C. Ahalt Department f Electrical Engineering The Ohi State University Abstract

More information

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank CAUSAL INFERENCE Technical Track Sessin I Phillippe Leite The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Phillippe Leite fr the purpse f this wrkshp Plicy questins are causal

More information

Optimization Programming Problems For Control And Management Of Bacterial Disease With Two Stage Growth/Spread Among Plants

Optimization Programming Problems For Control And Management Of Bacterial Disease With Two Stage Growth/Spread Among Plants Internatinal Jurnal f Engineering Science Inventin ISSN (Online): 9 67, ISSN (Print): 9 676 www.ijesi.rg Vlume 5 Issue 8 ugust 06 PP.0-07 Optimizatin Prgramming Prblems Fr Cntrl nd Management Of Bacterial

More information

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal

More information

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA Mdelling f Clck Behaviur Dn Percival Applied Physics Labratry University f Washingtn Seattle, Washingtn, USA verheads and paper fr talk available at http://faculty.washingtn.edu/dbp/talks.html 1 Overview

More information

Broadcast Program Generation for Unordered Queries with Data Replication

Broadcast Program Generation for Unordered Queries with Data Replication Bradcast Prgram Generatin fr Unrdered Queries with Data Replicatin Jiun-Lng Huang and Ming-Syan Chen Department f Electrical Engineering Natinal Taiwan University Taipei, Taiwan, ROC E-mail: jlhuang@arbr.ee.ntu.edu.tw,

More information

Determining the Accuracy of Modal Parameter Estimation Methods

Determining the Accuracy of Modal Parameter Estimation Methods Determining the Accuracy f Mdal Parameter Estimatin Methds by Michael Lee Ph.D., P.E. & Mar Richardsn Ph.D. Structural Measurement Systems Milpitas, CA Abstract The mst cmmn type f mdal testing system

More information

MATCHING TECHNIQUES. Technical Track Session VI. Emanuela Galasso. The World Bank

MATCHING TECHNIQUES. Technical Track Session VI. Emanuela Galasso. The World Bank MATCHING TECHNIQUES Technical Track Sessin VI Emanuela Galass The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Emanuela Galass fr the purpse f this wrkshp When can we use

More information

Math Foundations 20 Work Plan

Math Foundations 20 Work Plan Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant

More information

Admin. MDP Search Trees. Optimal Quantities. Reinforcement Learning

Admin. MDP Search Trees. Optimal Quantities. Reinforcement Learning Admin Reinfrcement Learning Cntent adapted frm Berkeley CS188 MDP Search Trees Each MDP state prjects an expectimax-like search tree Optimal Quantities The value (utility) f a state s: V*(s) = expected

More information

A Scalable Recurrent Neural Network Framework for Model-free

A Scalable Recurrent Neural Network Framework for Model-free A Scalable Recurrent Neural Netwrk Framewrk fr Mdel-free POMDPs April 3, 2007 Zhenzhen Liu, Itamar Elhanany Machine Intelligence Lab Department f Electrical and Cmputer Engineering The University f Tennessee

More information

Checking the resolved resonance region in EXFOR database

Checking the resolved resonance region in EXFOR database Checking the reslved resnance regin in EXFOR database Gttfried Bertn Sciété de Calcul Mathématique (SCM) Oscar Cabells OECD/NEA Data Bank JEFF Meetings - Sessin JEFF Experiments Nvember 0-4, 017 Bulgne-Billancurt,

More information

Technical Bulletin. Generation Interconnection Procedures. Revisions to Cluster 4, Phase 1 Study Methodology

Technical Bulletin. Generation Interconnection Procedures. Revisions to Cluster 4, Phase 1 Study Methodology Technical Bulletin Generatin Intercnnectin Prcedures Revisins t Cluster 4, Phase 1 Study Methdlgy Release Date: Octber 20, 2011 (Finalizatin f the Draft Technical Bulletin released n September 19, 2011)

More information

Collocation Map for Overcoming Data Sparseness

Collocation Map for Overcoming Data Sparseness Cllcatin Map fr Overcming Data Sparseness Mnj Kim, Yung S. Han, and Key-Sun Chi Department f Cmputer Science Krea Advanced Institute f Science and Technlgy Taejn, 305-701, Krea mj0712~eve.kaist.ac.kr,

More information

EXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE

EXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE EXPERIMENTAL STUD ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE Tmnbu Gt, Masaaki Ohba, Takashi Kurabuchi 2, Tmyuki End 3, shihik Akamine 4, and Tshihir Nnaka 2

More information

Multiple Source Multiple. using Network Coding

Multiple Source Multiple. using Network Coding Multiple Surce Multiple Destinatin Tplgy Inference using Netwrk Cding Pegah Sattari EECS, UC Irvine Jint wrk with Athina Markpulu, at UCI, Christina Fraguli, at EPFL, Lausanne Outline Netwrk Tmgraphy Gal,

More information

Resampling Methods. Chapter 5. Chapter 5 1 / 52

Resampling Methods. Chapter 5. Chapter 5 1 / 52 Resampling Methds Chapter 5 Chapter 5 1 / 52 1 51 Validatin set apprach 2 52 Crss validatin 3 53 Btstrap Chapter 5 2 / 52 Abut Resampling An imprtant statistical tl Pretending the data as ppulatin and

More information

Hypothesis Tests for One Population Mean

Hypothesis Tests for One Population Mean Hypthesis Tests fr One Ppulatin Mean Chapter 9 Ala Abdelbaki Objective Objective: T estimate the value f ne ppulatin mean Inferential statistics using statistics in rder t estimate parameters We will be

More information

Pattern Recognition 2014 Support Vector Machines

Pattern Recognition 2014 Support Vector Machines Pattern Recgnitin 2014 Supprt Vectr Machines Ad Feelders Universiteit Utrecht Ad Feelders ( Universiteit Utrecht ) Pattern Recgnitin 1 / 55 Overview 1 Separable Case 2 Kernel Functins 3 Allwing Errrs (Sft

More information

ENSC Discrete Time Systems. Project Outline. Semester

ENSC Discrete Time Systems. Project Outline. Semester ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding

More information

THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES

THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES PREFERRED RELIABILITY PAGE 1 OF 5 PRACTICES PRACTICE NO. PT-TE-1409 THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC Practice: Perfrm all thermal envirnmental tests n electrnic spaceflight hardware in a flight-like

More information

CS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007

CS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007 CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is

More information

CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS

CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS 1 Influential bservatins are bservatins whse presence in the data can have a distrting effect n the parameter estimates and pssibly the entire analysis,

More information

Eric Klein and Ning Sa

Eric Klein and Ning Sa Week 12. Statistical Appraches t Netwrks: p1 and p* Wasserman and Faust Chapter 15: Statistical Analysis f Single Relatinal Netwrks There are fur tasks in psitinal analysis: 1) Define Equivalence 2) Measure

More information

ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS

ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS Particle Acceleratrs, 1986, Vl. 19, pp. 99-105 0031-2460/86/1904-0099/$15.00/0 1986 Grdn and Breach, Science Publishers, S.A. Printed in the United States f America ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS

More information

Document for ENES5 meeting

Document for ENES5 meeting HARMONISATION OF EXPOSURE SCENARIO SHORT TITLES Dcument fr ENES5 meeting Paper jintly prepared by ECHA Cefic DUCC ESCOM ES Shrt Titles Grup 13 Nvember 2013 OBJECTIVES FOR ENES5 The bjective f this dcument

More information

and the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:

and the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are: Algrithm fr Estimating R and R - (David Sandwell, SIO, August 4, 2006) Azimith cmpressin invlves the alignment f successive eches t be fcused n a pint target Let s be the slw time alng the satellite track

More information

MATCHING TECHNIQUES Technical Track Session VI Céline Ferré The World Bank

MATCHING TECHNIQUES Technical Track Session VI Céline Ferré The World Bank MATCHING TECHNIQUES Technical Track Sessin VI Céline Ferré The Wrld Bank When can we use matching? What if the assignment t the treatment is nt dne randmly r based n an eligibility index, but n the basis

More information

Reinforcement Learning" CMPSCI 383 Nov 29, 2011!

Reinforcement Learning CMPSCI 383 Nov 29, 2011! Reinfrcement Learning" CMPSCI 383 Nv 29, 2011! 1 Tdayʼs lecture" Review f Chapter 17: Making Cmple Decisins! Sequential decisin prblems! The mtivatin and advantages f reinfrcement learning.! Passive learning!

More information

Enhancing Performance of MLP/RBF Neural Classifiers via an Multivariate Data Distribution Scheme

Enhancing Performance of MLP/RBF Neural Classifiers via an Multivariate Data Distribution Scheme Enhancing Perfrmance f / Neural Classifiers via an Multivariate Data Distributin Scheme Halis Altun, Gökhan Gelen Nigde University, Electrical and Electrnics Engineering Department Nigde, Turkey haltun@nigde.edu.tr

More information

Distributions, spatial statistics and a Bayesian perspective

Distributions, spatial statistics and a Bayesian perspective Distributins, spatial statistics and a Bayesian perspective Dug Nychka Natinal Center fr Atmspheric Research Distributins and densities Cnditinal distributins and Bayes Thm Bivariate nrmal Spatial statistics

More information

AP Statistics Notes Unit Two: The Normal Distributions

AP Statistics Notes Unit Two: The Normal Distributions AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).

More information

NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION

NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION NUROP Chinese Pinyin T Chinese Character Cnversin NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION CHIA LI SHI 1 AND LUA KIM TENG 2 Schl f Cmputing, Natinal University f Singapre 3 Science

More information

Learning to Control an Unstable System with Forward Modeling

Learning to Control an Unstable System with Forward Modeling 324 Jrdan and Jacbs Learning t Cntrl an Unstable System with Frward Mdeling Michael I. Jrdan Brain and Cgnitive Sciences MIT Cambridge, MA 02139 Rbert A. Jacbs Cmputer and Infrmatin Sciences University

More information

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322 ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA by J. C. SPROTT December 4, 1969 PLP N. 3 These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated. They are fr private

More information

On Huntsberger Type Shrinkage Estimator for the Mean of Normal Distribution ABSTRACT INTRODUCTION

On Huntsberger Type Shrinkage Estimator for the Mean of Normal Distribution ABSTRACT INTRODUCTION Malaysian Jurnal f Mathematical Sciences 4(): 7-4 () On Huntsberger Type Shrinkage Estimatr fr the Mean f Nrmal Distributin Department f Mathematical and Physical Sciences, University f Nizwa, Sultanate

More information

Performance Bounds for Detect and Avoid Signal Sensing

Performance Bounds for Detect and Avoid Signal Sensing Perfrmance unds fr Detect and Avid Signal Sensing Sam Reisenfeld Real-ime Infrmatin etwrks, University f echnlgy, Sydney, radway, SW 007, Australia samr@uts.edu.au Abstract Detect and Avid (DAA) is a Cgnitive

More information

Early detection of mining truck failure by modelling its operation with neural networks classification algorithms

Early detection of mining truck failure by modelling its operation with neural networks classification algorithms RU, Rand GOLOSINSKI, T.S. Early detectin f mining truck failure by mdelling its peratin with neural netwrks classificatin algrithms. Applicatin f Cmputers and Operatins Research ill the Minerals Industries,

More information

Module 4: General Formulation of Electric Circuit Theory

Module 4: General Formulation of Electric Circuit Theory Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated

More information

Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart

Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Sandy D. Balkin Dennis K. J. Lin y Pennsylvania State University, University Park, PA 16802 Sandy Balkin is a graduate student

More information

The blessing of dimensionality for kernel methods

The blessing of dimensionality for kernel methods fr kernel methds Building classifiers in high dimensinal space Pierre Dupnt Pierre.Dupnt@ucluvain.be Classifiers define decisin surfaces in sme feature space where the data is either initially represented

More information

Maximum A Posteriori (MAP) CS 109 Lecture 22 May 16th, 2016

Maximum A Posteriori (MAP) CS 109 Lecture 22 May 16th, 2016 Maximum A Psteriri (MAP) CS 109 Lecture 22 May 16th, 2016 Previusly in CS109 Game f Estimatrs Maximum Likelihd Nn spiler: this didn t happen Side Plt argmax argmax f lg Mther f ptimizatins? Reviving an

More information

A New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation

A New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation III-l III. A New Evaluatin Measure J. Jiner and L. Werner Abstract The prblems f evaluatin and the needed criteria f evaluatin measures in the SMART system f infrmatin retrieval are reviewed and discussed.

More information

WRITING THE REPORT. Organizing the report. Title Page. Table of Contents

WRITING THE REPORT. Organizing the report. Title Page. Table of Contents WRITING THE REPORT Organizing the reprt Mst reprts shuld be rganized in the fllwing manner. Smetime there is a valid reasn t include extra chapters in within the bdy f the reprt. 1. Title page 2. Executive

More information

Lab 1 The Scientific Method

Lab 1 The Scientific Method INTRODUCTION The fllwing labratry exercise is designed t give yu, the student, an pprtunity t explre unknwn systems, r universes, and hypthesize pssible rules which may gvern the behavir within them. Scientific

More information

THERMAL TEST LEVELS & DURATIONS

THERMAL TEST LEVELS & DURATIONS PREFERRED RELIABILITY PAGE 1 OF 7 PRACTICES PRACTICE NO. PT-TE-144 Practice: 1 Perfrm thermal dwell test n prtflight hardware ver the temperature range f +75 C/-2 C (applied at the thermal cntrl/munting

More information

Internal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.

Internal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9. Sectin 7 Mdel Assessment This sectin is based n Stck and Watsn s Chapter 9. Internal vs. external validity Internal validity refers t whether the analysis is valid fr the ppulatin and sample being studied.

More information

Intelligent Pharma- Chemical and Oil & Gas Division Page 1 of 7. Global Business Centre Ave SE, Calgary, AB T2G 0K6, AB.

Intelligent Pharma- Chemical and Oil & Gas Division Page 1 of 7. Global Business Centre Ave SE, Calgary, AB T2G 0K6, AB. Intelligent Pharma- Chemical and Oil & Gas Divisin Page 1 f 7 Intelligent Pharma Chemical and Oil & Gas Divisin Glbal Business Centre. 120 8 Ave SE, Calgary, AB T2G 0K6, AB. Canada Dr. Edelsys Cdrniu-Business

More information

Engineering Approach to Modelling Metal THz Structures

Engineering Approach to Modelling Metal THz Structures Terahertz Science and Technlgy, ISSN 1941-7411 Vl.4, N.1, March 11 Invited Paper ngineering Apprach t Mdelling Metal THz Structures Stepan Lucyszyn * and Yun Zhu Department f, Imperial Cllege Lndn, xhibitin

More information

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax .7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical

More information

Aerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison

Aerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison Jurnal f Physics: Cnference Series OPEN ACCESS Aerdynamic Separability in Tip Speed Rati and Separability in Wind Speed- a Cmparisn T cite this article: M L Gala Sants et al 14 J. Phys.: Cnf. Ser. 555

More information

The Sputtering Problem James A Glackin, James V. Matheson

The Sputtering Problem James A Glackin, James V. Matheson The Sputtering Prblem James A Glackin, James V. Mathesn I prpse t cnsider first the varius elements f the subject, next its varius parts r sectins, and finally the whle in its internal structure. In ther

More information

Data Mining: Concepts and Techniques. Classification and Prediction. Chapter February 8, 2007 CSE-4412: Data Mining 1

Data Mining: Concepts and Techniques. Classification and Prediction. Chapter February 8, 2007 CSE-4412: Data Mining 1 Data Mining: Cncepts and Techniques Classificatin and Predictin Chapter 6.4-6 February 8, 2007 CSE-4412: Data Mining 1 Chapter 6 Classificatin and Predictin 1. What is classificatin? What is predictin?

More information

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

Fall 2013 Physics 172 Recitation 3 Momentum and Springs Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.

More information

DECISION BASED COLLABORATIVE OPTIMIZATION

DECISION BASED COLLABORATIVE OPTIMIZATION 8 th ASCE Specialty Cnference n Prbabilistic Mechanics and Structural Reliability PMC2000-217 DECISION BASED COLLABORATIVE OPTIMIZATION X. Gu and J.E. Renaud University f Ntre Dame, Indiana, IN 46556 Renaud.2@nd.edu

More information

APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL

APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL JP2.11 APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL Xingang Fan * and Jeffrey S. Tilley University f Alaska Fairbanks, Fairbanks,

More information

Administrativia. Assignment 1 due thursday 9/23/2004 BEFORE midnight. Midterm exam 10/07/2003 in class. CS 460, Sessions 8-9 1

Administrativia. Assignment 1 due thursday 9/23/2004 BEFORE midnight. Midterm exam 10/07/2003 in class. CS 460, Sessions 8-9 1 Administrativia Assignment 1 due thursday 9/23/2004 BEFORE midnight Midterm eam 10/07/2003 in class CS 460, Sessins 8-9 1 Last time: search strategies Uninfrmed: Use nly infrmatin available in the prblem

More information

Resampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017

Resampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017 Resampling Methds Crss-validatin, Btstrapping Marek Petrik 2/21/2017 Sme f the figures in this presentatin are taken frm An Intrductin t Statistical Learning, with applicatins in R (Springer, 2013) with

More information

Kinetic Model Completeness

Kinetic Model Completeness 5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins

More information

Synchronous Motor V-Curves

Synchronous Motor V-Curves Synchrnus Mtr V-Curves 1 Synchrnus Mtr V-Curves Intrductin Synchrnus mtrs are used in applicatins such as textile mills where cnstant speed peratin is critical. Mst small synchrnus mtrs cntain squirrel

More information

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

Interference is when two (or more) sets of waves meet and combine to produce a new pattern. Interference Interference is when tw (r mre) sets f waves meet and cmbine t prduce a new pattern. This pattern can vary depending n the riginal wave directin, wavelength, amplitude, etc. The tw mst extreme

More information

Lead/Lag Compensator Frequency Domain Properties and Design Methods

Lead/Lag Compensator Frequency Domain Properties and Design Methods Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin

More information

BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS. Christopher Costello, Andrew Solow, Michael Neubert, and Stephen Polasky

BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS. Christopher Costello, Andrew Solow, Michael Neubert, and Stephen Polasky BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS Christpher Cstell, Andrew Slw, Michael Neubert, and Stephen Plasky Intrductin The central questin in the ecnmic analysis f climate change plicy cncerns

More information

A Regression Solution to the Problem of Criterion Score Comparability

A Regression Solution to the Problem of Criterion Score Comparability A Regressin Slutin t the Prblem f Criterin Scre Cmparability William M. Pugh Naval Health Research Center When the criterin measure in a study is the accumulatin f respnses r behavirs fr an individual

More information

^YawataR&D Laboratory, Nippon Steel Corporation, Tobata, Kitakyushu, Japan

^YawataR&D Laboratory, Nippon Steel Corporation, Tobata, Kitakyushu, Japan Detectin f fatigue crack initiatin frm a ntch under a randm lad C. Makabe," S. Nishida^C. Urashima,' H. Kaneshir* "Department f Mechanical Systems Engineering, University f the Ryukyus, Nishihara, kinawa,

More information

BASD HIGH SCHOOL FORMAL LAB REPORT

BASD HIGH SCHOOL FORMAL LAB REPORT BASD HIGH SCHOOL FORMAL LAB REPORT *WARNING: After an explanatin f what t include in each sectin, there is an example f hw the sectin might lk using a sample experiment Keep in mind, the sample lab used

More information

(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f

(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f 1.0 Review f Electrmagnetic Field Thery Selected aspects f electrmagnetic thery are reviewed in this sectin, with emphasis n cncepts which are useful in understanding magnet design. Detailed, rigrus treatments

More information

Comparison of hybrid ensemble-4dvar with EnKF and 4DVar for regional-scale data assimilation

Comparison of hybrid ensemble-4dvar with EnKF and 4DVar for regional-scale data assimilation Cmparisn f hybrid ensemble-4dvar with EnKF and 4DVar fr reginal-scale data assimilatin Jn Pterjy and Fuqing Zhang Department f Meterlgy The Pennsylvania State University Wednesday 18 th December, 2013

More information

COMP 551 Applied Machine Learning Lecture 4: Linear classification

COMP 551 Applied Machine Learning Lecture 4: Linear classification COMP 551 Applied Machine Learning Lecture 4: Linear classificatin Instructr: Jelle Pineau (jpineau@cs.mcgill.ca) Class web page: www.cs.mcgill.ca/~jpineau/cmp551 Unless therwise nted, all material psted

More information

1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp

1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp THE POWER AND LIMIT OF NEURAL NETWORKS T. Y. Lin Department f Mathematics and Cmputer Science San Jse State University San Jse, Califrnia 959-003 tylin@cs.ssu.edu and Bereley Initiative in Sft Cmputing*

More information

An Efficient Load Shedding Scheme from Customer s Perspective

An Efficient Load Shedding Scheme from Customer s Perspective Internatinal Jurnal f Advanced Research in Electrical, Electrnics and Instrumentatin Engineering (An ISO 3297: 2007 Certified Organizatin) Vl. 2, Issue 10, Octber 2013 An Efficient Lad Shedding Scheme

More information

A Novel Electro-thermal Simulation Approach to Power IGBT Modules for Automotive Traction Applications

A Novel Electro-thermal Simulation Approach to Power IGBT Modules for Automotive Traction Applications Special Issue Recent R&D Activities f Pwer Devices fr Hybrid Electric Vehicles 27 Research Reprt A Nvel Electr-thermal Simulatin Apprach t Pwer IGBT Mdules fr Autmtive Tractin Applicatins Takashi Kjima,

More information

Differentiation Applications 1: Related Rates

Differentiation Applications 1: Related Rates Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm

More information

What is Statistical Learning?

What is Statistical Learning? What is Statistical Learning? Sales 5 10 15 20 25 Sales 5 10 15 20 25 Sales 5 10 15 20 25 0 50 100 200 300 TV 0 10 20 30 40 50 Radi 0 20 40 60 80 100 Newspaper Shwn are Sales vs TV, Radi and Newspaper,

More information

Simulation of Line Outage Distribution Factors (L.O.D.F) Calculation for N-Buses System

Simulation of Line Outage Distribution Factors (L.O.D.F) Calculation for N-Buses System Simulatin f Line Outage Distributin Factrs (L.O.D.F) Calculatin fr N-Buses System Rashid H. AL-Rubayi Department f Electrical Engineering, University f Technlgy Afaneen A. Abd Department f Electrical Engineering,

More information

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce

More information

Ecology 302 Lecture III. Exponential Growth (Gotelli, Chapter 1; Ricklefs, Chapter 11, pp )

Ecology 302 Lecture III. Exponential Growth (Gotelli, Chapter 1; Ricklefs, Chapter 11, pp ) Eclgy 302 Lecture III. Expnential Grwth (Gtelli, Chapter 1; Ricklefs, Chapter 11, pp. 222-227) Apcalypse nw. The Santa Ana Watershed Prject Authrity pulls n punches in prtraying its missin in apcalyptic

More information

Pure adaptive search for finite global optimization*

Pure adaptive search for finite global optimization* Mathematical Prgramming 69 (1995) 443-448 Pure adaptive search fr finite glbal ptimizatin* Z.B. Zabinskya.*, G.R. Wd b, M.A. Steel c, W.P. Baritmpa c a Industrial Engineering Prgram, FU-20. University

More information

Performance Comparison of Nonlinear Filters for Indoor WLAN Positioning

Performance Comparison of Nonlinear Filters for Indoor WLAN Positioning Perfrmance Cmparisn f Nnlinear Filters fr Indr WLAN Psitining Hui Wang Andrei Szab Jachim Bamberger Learning Systems Infrmatin and Cmmunicatins Crprate Technlgy Siemens AG Munich Germany Email: {hui.wang.ext.andrei.szab.jachim.bamberger}@siemens.cm

More information

Video Encoder Control

Video Encoder Control Vide Encder Cntrl Thmas Wiegand Digital Image Cmmunicatin 1 / 41 Outline Intrductin Encder Cntrl using Lagrange multipliers Lagrangian ptimizatin Lagrangian bit allcatin Lagrangian Optimizatin in Hybrid

More information

Statistics, Numerical Models and Ensembles

Statistics, Numerical Models and Ensembles Statistics, Numerical Mdels and Ensembles Duglas Nychka, Reinhard Furrer,, Dan Cley Claudia Tebaldi, Linda Mearns, Jerry Meehl and Richard Smith (UNC). Spatial predictin and data assimilatin Precipitatin

More information

Review of Simulation Approaches in Reliability and Availability Modeling

Review of Simulation Approaches in Reliability and Availability Modeling Internatinal Jurnal f Perfrmability Engineering, Vl. 12, N. 4, July 2016, pp. 369-388 Ttem Publisher, Inc., 4625 Stargazer Dr., Plan, Texas 75024, U.S.A Review f Simulatin Appraches in Reliability and

More information

Computational modeling techniques

Computational modeling techniques Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins

More information

COMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d)

COMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d) COMP 551 Applied Machine Learning Lecture 9: Supprt Vectr Machines (cnt d) Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Class web page: www.cs.mcgill.ca/~hvanh2/cmp551 Unless therwise

More information

Building research leadership consortia for Quantum Technology Research Hubs. Call type: Expression of Interest

Building research leadership consortia for Quantum Technology Research Hubs. Call type: Expression of Interest Building research leadership cnsrtia fr Quantum Technlgy Research Hubs Call type: Expressin f Interest Clsing date: 17:00, 07 August 2018 Hw t apply: Expressin f Interest (EI) fr research leaders t attend

More information

Displacement and Deflection Sensitivity of Gas-coupled Laser Acoustic. Detector

Displacement and Deflection Sensitivity of Gas-coupled Laser Acoustic. Detector 1st Internatinal Sympsium n Laser Ultrasnics: Science, echnlgy and Applicatins July 16-18 008, Mntreal, Canada Displacement and Deflectin Sensitivity f Gas-cupled Laser Acustic Detectin James N. CARON

More information

BLAST / HIDDEN MARKOV MODELS

BLAST / HIDDEN MARKOV MODELS CS262 (Winter 2015) Lecture 5 (January 20) Scribe: Kat Gregry BLAST / HIDDEN MARKOV MODELS BLAST CONTINUED HEURISTIC LOCAL ALIGNMENT Use Cmmnly used t search vast bilgical databases (n the rder f terabases/tetrabases)

More information

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu

More information

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018 Michael Faraday lived in the Lndn area frm 1791 t 1867. He was 29 years ld when Hand Oersted, in 1820, accidentally discvered that electric current creates magnetic field. Thrugh empirical bservatin and

More information

Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff

Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeff Reading: Chapter 2 STATS 202: Data mining and analysis September 27, 2017 1 / 20 Supervised vs. unsupervised learning In unsupervised

More information

Ship-Track Models Based on Poisson-Distributed Port-Departure Times

Ship-Track Models Based on Poisson-Distributed Port-Departure Times Naval Research Labratry Washingtn, DC 20375-5320 NRL/FR/7121--06-10,122 Ship-Track Mdels Based n Pissn-Distributed Prt-Departure Times RICHARD HEITMEYER Acustic Signal Prcessing Branch Acustic Divisin

More information