many vibrating mechanical systems experience undesirable impacts that cause wear or sub - optimal performance .
occasional impacts may be permissible if they result from running components of the system at high speeds for greater efficiency , and some impacts are unavoidable such as those due to random or rare events . in these cases
it is important to have a clear understanding of the dynamical behaviour that impacts may induce .
impacting dynamics is often complicated or chaotic because impacts are highly nonlinear phenomena @xcite .
for instance tubes in heat exchangers vibrate at high fluid velocities and impact against baffles in place to guide the fluid flow .
simple mathematical models of heat exchangers reveal that chaotic dynamics may be created at the onset of recurring impacts @xcite .
rotating cutters spun at high speeds experience repeated contact loss with the material being cut .
the resulting impacts between the cutter and the material may similarly induce chaotic dynamics @xcite .
in contrast , some mechanical systems use impacts to achieve their function . atomic force microscopes measure surface topography and the chemical properties of a sample on the nanoscale by gently hitting the sample with a vibrating cantilever . in this context
it is important to understand the complex impacting dynamics so that the inverse problem of describing the sample can be performed effectively @xcite .
impacts can often be modelled accurately by carefully describing the deformations that components of the system undergo during impacts @xcite .
however , for the purposes of understanding vibro - impacting dynamics , such a modelling approach is too cumbersome and a low degree - of - freedom ode model can be more useful . despite the low - dimensionality ,
such models have been shown to quantitatively match the experimental data of a variety of impacting systems .
examples include a cam - follower system involving occasional detachments between the cam and the follower @xcite , a pendulum experiencing near - instantaneous impacts with a solid wall @xcite , and compliant impacts of a steel block with an elastic beam @xcite .
( 15,3.8 ) ( 0,0 ) , where @xmath0 is a system parameter .
[ fig : grazbifschem ] , title="fig:",height=136 ] ( 5.1,0 ) , where @xmath0 is a system parameter .
[ fig : grazbifschem ] , title="fig:",height=136 ] ( 10.2,0 ) , where @xmath0 is a system parameter .
[ fig : grazbifschem ] , title="fig:",height=136 ] ( .5,3.6)*a * ( 5.6,3.6)*b * ( 10.7,3.6)*c * ( 2,3.6)@xmath1 ( 7.1,3.6)@xmath2 ( 12.2,3.6)@xmath3 ( 4.48,.2)@xmath4 ( 9.58,.2)@xmath4 ( 14.68,.2)@xmath4    using a low degree - of - freedom ode model , the evolution of the system between impacts is tracked in phase space .
periodic behaviour without impacts corresponds to a periodic orbit in phase space that does not reach the switching manifold , @xmath4 , fig .
[ fig : grazbifschem]-a .
physically , @xmath4 corresponds to locations or instances where mechanical components come into contact , or lose contact .
as parameters vary , the system may transition from an impact - free regime to the repeated ( though not necessarily regular ) occurrence of impacts . in phase space ,
the transition occurs when the periodic orbit of the ode model attains an intersection with @xmath4 , fig .
[ fig : grazbifschem]-b .
this is known as a _
grazing bifurcation_.    in this paper we study grazing bifurcations of the three - dimensional piecewise - smooth ode system , @xmath5 where @xmath6 and @xmath7 are smooth functions , @xmath0 is a parameter , and the coordinates @xmath8 are chosen so that @xmath4 is simply the coordinate plane @xmath9 .
we assume that for @xmath1 , there exists an attracting periodic orbit describing non - impacting dynamics located entirely in the region @xmath10 , and that the periodic orbit grazes @xmath4 at the origin when @xmath2 . in the context of vibro - impacting systems
, @xmath10 corresponds to not - in - contact dynamics governed by @xmath6 , and @xmath11 corresponds to in - contact dynamics governed by @xmath7 .
impacts may instead be modelled as instantaneous events with energy loss and velocity reversal , in which case a map is usually defined on the switching manifold to describe the action of an impact @xcite .
theoretical studies of piecewise - smooth and hybrid dynamical systems have led to a useful classification of grazing bifurcations @xcite .
this paper concerns regular grazing bifurcations .
the grazing bifurcation of ( [ eq : f ] ) at @xmath2 is said to be _ regular _ if @xmath12 for all @xmath13 in a neighbourhood of @xmath14 .
this condition arises naturally in mechanical systems with compliant impacts , and implies that @xmath4 is neither attracting nor repelling at any point .    as indicated in fig .
[ fig : grazbifschem]-c , the steady - state dynamics of ( [ eq : f ] ) for @xmath3 is often complicated .
for this reason it is valuable to study the oscillatory dynamics using a return map based on the points on a poincar section .
a normal form for such a map for regular grazing in @xmath15 is the nordmark map @xmath16 where @xmath17 with @xmath18 , and @xmath19 , as determined by the sign of certain coefficients , see  [ sec : nordmark ] .
the coordinates @xmath20 represent points on a poincar section of ( [ eq : f ] ) , with @xmath21 the bifurcation parameter and the grazing bifurcation occurring at @xmath22 . the nordmark map ( [ eq : n ] ) ,
applicable also to models with instantaneous impacts @xcite , includes only the leading order terms of the return map , so is valid only for dynamics close to the grazing bifurcation @xcite .
each iteration of ( [ eq : n ] ) corresponds to one oscillation of ( [ eq : f ] ) near the grazing periodic orbit .
the utility of ( [ eq : n ] ) lies in the fact that the nature of the dynamics can be identified by the location of the corresponding points in the @xmath23-plane .
specifically , ( [ eq : n ] ) is formulated so that if @xmath24 then the oscillation lies entirely in @xmath10 , and if @xmath25 then the oscillation enters @xmath11 .
the non - impacting , attracting periodic orbit shown in fig .
[ fig : grazbifschem]-a corresponds to the fixed point of @xmath26 , @xmath27 with @xmath28 and @xmath29 .
more generally , a periodic orbit of ( [ eq : f ] ) appears as a finite set of points in the @xmath23 plane , as in fig .
[ fig : detbifdiag ] which shows a typical bifurcation diagram of ( [ eq : n ] ) .
this figure shows a period - incrementing cascade , corresponding to several different periodic orbits , and apparently chaotic dynamics , as indicated by a cloud of points . for alternate values of @xmath30 and @xmath31 , the fixed point @xmath32 may bifurcate directly to chaos @xcite . the square - root term in ( [ eq : n ] )
is an artifact of the tangency between the periodic orbit and the switching manifold of ( [ eq : f ] ) at the grazing bifurcation , and is responsible for the distinctive shape of the bifurcation diagram near @xmath22 .
for some vibro - impacting systems it is more appropriate for the return map to be piecewise - linear and either continuous @xcite or discontinuous @xcite .
such maps predict fundamentally different bifurcation structures to those of ( [ eq : n ] ) .
( 12,6 ) ( 0,0 ) ) with @xmath33 , @xmath34 , and @xmath35 .
these parameter values correspond to the vibro - impacting system described in  [ sec : oscillator ] , with @xmath36 .
the fixed point of ( [ eq : n ] ) for @xmath29 , given by ( [ eq : xlyl ] ) , corresponds to an attracting , non - impacting periodic orbit of period @xmath37 .
the map ( [ eq : n ] ) has attracting @xmath38 and @xmath39-cycles for certain values of @xmath40 as shown , as well as an attracting @xmath41-cycle in the approximate range , @xmath42 .
these correspond to periodic orbits of period approximately equal to @xmath43 , for @xmath44 , that experience one impact per period .
there also appears to be a chaotic attractor for the approximate range , @xmath45 .
[ fig : detbifdiag ] , title="fig:",height=226 ] ( 6.2,0)@xmath46 ( 0,3.8)@xmath47    in order to properly explain complicated vibro - impacting dynamics , the effects of randomness and uncertainties needs to be taken into account . mechanical systems are subject to background vibrations and other sources of noise .
experimentally measured parameters involve error , and some physical features are left unmodelled .
for instance , one - degree - of - freedom models do not capture high frequency modes that are often excited by impacts @xcite .    to quantitatively describe stochastic impacting dynamics ,
stochastic averaging methods have proved useful for vibro - impacting systems that experience a wide range of impact velocities @xcite .
if only low - velocity impacts are relevant , then it is useful to study ( [ eq : n ] ) . in his phd thesis @xcite , griffin studied ( [ eq : n ] ) in the presence of additive white noise .
he found that noise blurs bifurcation diagrams and washes out high - period solutions in the same manner as for smooth maps , such as the logistic map @xcite .
recently it was shown that white noise added to ( [ eq : f ] ) translates to additive white noise in ( [ eq : n ] ) , @xcite .
such a noise formulation may be sensible for vibro - impacting systems for which a forcing term or external fluctuations represent a significant source of uncertainty .
however , impact events themselves constitute a substantial source of randomness .
the purpose of this paper is to construct and analyse stochastic versions of ( [ eq : n ] ) for which randomness stems purely from impact events .
we consider three different types of impact noise for ( [ eq : f ] ) using an ornstein - uhlenbeck process with stationary density @xmath48 $ ] , where @xmath49 represents the noise amplitude and @xmath50 is the correlation time .
we first consider uncertainties in @xmath4 , then uncertainties in @xmath7 , and lastly uncertainties in @xmath7 in the white noise ( zero correlation time ) limit .    here
let us indicate the forms of the stochastic maps that we obtain .
coloured noise in @xmath4 leads to random perturbations in both the map and the switching condition as @xmath51 where @xmath52 are gaussian random variables , and @xmath53 is a constant . for coloured noise in the impacting dynamics
we obtain @xmath54 for a particular nonlinear function @xmath55 . for white noise ( @xmath56 ) in the impacting dynamics the map takes the form @xmath57 where @xmath58 and @xmath59 are random variables , and @xmath60 and @xmath61 are nonlinear functions . notice that @xmath62 is stochastic for @xmath24 , whereas @xmath63 and @xmath64 are not .
this is because noise in @xmath4 generates anomalous crossings of @xmath65 .
consequently @xmath62 exhibits stochastic dynamics for @xmath29 , while @xmath63 and @xmath64 do not .
@xmath63 and @xmath64 involve noise terms proportional to @xmath66 , and for this reason exhibit increasing variability for larger values of @xmath40 .    to obtain a more detailed comparison of @xmath62 , @xmath63 and @xmath64 , we first carefully derive ( [ eq : n ] ) in  [ sec : nordmark ] .
in  [ sec : oscillator ] we introduce a prototypical compliant vibro - impacting system to illustrate our results . in  [ sec : randomness ] we add randomness and derive ( [ eq : na2])-([eq : nc2 ] ) .
since ( [ eq : na2])-([eq : nc2 ] ) involve fundamentally different noise terms , they exhibit different sensitivities to the noise amplitude @xmath67 .
therefore we use different values of @xmath67 in each of the different models , in order to make appropriate comparisons . for each @xmath68 ( @xmath69 )
we write @xmath70 and identify the appropriate value @xmath71 , so that @xmath62 , @xmath63 and @xmath64 display roughly the same dynamics for the vibro - impacting system with the parameters of fig .
[ fig : detbifdiag ] and @xmath72 ( chosen for illustration ) and @xmath73 .
this enables us to quantitatively compare @xmath62 , @xmath63 and @xmath64 in  [ sec : dynamics ] .
conclusions are presented in  [ sec : conc ] .
( 15,7.5 ) ( 0,0 ) ) for the three - dimensional piecewise - smooth system ( [ eq : f ] ) .
the solid curve represents an orbit of ( [ eq : f ] ) .
the dashed curves show virtual extensions of this orbit into @xmath11 as governed by @xmath6 . [
fig : grazschem3d ] , title="fig:",height=283 ] ( 14.5,3.25)@xmath74 ( 13.65,5.75)@xmath75 ( 9.7,7.1)@xmath76 ( 1.2,5.8)@xmath77 ( 12.1,1.6)@xmath78 ( 11.3,6.2)@xmath4 ( 1.78,4.92)@xmath79 ( 11.94,3.7)@xmath80 ( 10.7,5.84)@xmath81 ( 8.54,1.12)@xmath82 ( 10.7,2.64)@xmath83 ( 3.02,4.48)@xmath84 ( 11.2,3.5)@xmath85    in this section we derive a return map for the generic deterministic system ( [ eq : f ] ) valid near the grazing bifurcation @xmath2 and provide an explicit coordinate change that transforms the map to ( [ eq : n ] ) . such a derivation is given in @xcite .
we provide an explicit construction of ( [ eq : n ] ) here in order to provide a basis for deriving the stochastic maps in  [ sec : randomness ] . throughout this section
we write @xmath86 for points of ( [ eq : f ] ) in @xmath15 .    as discussed in
[ sec : intro ] , we assume that with @xmath2 the system ( [ eq : f ] ) has a periodic orbit that intersects the origin , but is otherwise contained in @xmath10 .
this implies that @xmath87 for @xmath6 at the origin with @xmath2 ( i.e.  @xmath88 ) . for simplicity
, we assume that we further have @xmath89 for all @xmath90 in a neighbourhood of @xmath14 , which can usually be imposed by an appropriate coordinate change .
this assumption is particularly useful in the case that randomness is present in the switching condition , see  [ sub : s1 ] .
the key to deriving the nordmark map ( [ eq : n ] ) is selecting a convenient poincar section and constructing a _ discontinuity map _ that accounts for the difference between @xmath6 and @xmath7 , that is , the difference between impacting and non - impacting dynamics .
we let @xmath77 represent a generic poincar section of ( [ eq : f ] ) that lies in @xmath10 and intersects the grazing periodic orbit transversally , and let @xmath78 be the coordinate plane @xmath91 , see fig .
[ fig : grazschem3d ] .    given @xmath92
, we define @xmath93 as the next intersection of the forward orbit governed by @xmath6 ( i.e.  ignoring the switching condition at @xmath9 ) with @xmath78 .
if @xmath11 for the point @xmath80 , then @xmath80 does not represent the true intersection of the orbit of ( [ eq : f ] ) with @xmath77 .
nevertheless , we study the return map on @xmath78 ( @xmath94 ) rather than the return map on @xmath77 ( @xmath95 ) , because the two maps are conjugate and the map on @xmath78 has a simpler form than the map on @xmath77 .    to derive the return map on @xmath78
, we must consider three additional points , @xmath81 , @xmath82 , and @xmath83 , see fig .
[ fig : grazschem3d ] .
the points @xmath81 and @xmath82 correspond to the entry and exit locations of the orbit with the impacting region , @xmath11 .
then @xmath80 is obtained by travelling forward from the entry point @xmath81 to @xmath78 using @xmath6 , whilst @xmath83 is obtained by travelling backward from the exit point @xmath82 to @xmath78 using @xmath6 .
the discontinuity map is defined as @xmath96 . if @xmath97 for the point @xmath80 , then @xmath98 .    in order to derive an explicit expression for the discontinuity map
, we consider the three steps , @xmath99 , @xmath100 , and @xmath101 , individually , and expand @xmath6 and @xmath7 about @xmath2 and the origin by writing @xmath102 where @xmath103 , and @xmath104 . in ( [ eq : alphabetagamma ] ) , and throughout this section , @xmath105 is short - hand big - o notation for @xmath106 .
notice that @xmath107 is assumed to be of the same order as @xmath75 , @xmath76 , and @xmath108 , which is appropriate in view of the tangency between the grazing periodic orbit and @xmath4 .
we assume @xmath109 such that the tangency is quadratic for both @xmath6 and @xmath7 , and has the orientation depicted in fig .
[ fig : grazschem3d ] .    from series expansions of the orbits governed by @xmath6 and @xmath7 , see @xcite , we obtain the following formulas for the three steps in the discontinuity map in the case @xmath110 ( writing @xmath111 ) . for @xmath99 , @xmath112 for @xmath100 , @xmath113 and for @xmath101 , @xmath114 by combining ( [ eq : v2w2])-([eq : u4w4 ] ) we get @xmath115 in terms of @xmath116 , @xmath117 where @xmath118 the discontinuity map is then @xmath119    to complete the map ( [ eq : n ] ) , which represents @xmath94 as shown in fig .
[ fig : grazschem3d ] , we must combine @xmath120 with the global return map @xmath121 .
@xmath122 depends on global properties of @xmath6 and is smooth , so we can write @xmath123 for some @xmath124 where each @xmath125
. then the desired return map on @xmath78 is the composition @xmath126 .
finally we apply a coordinate change to convert the map to the normal form ( [ eq : n ] ) that involves only three parameters , @xmath30 , @xmath31 and @xmath127 . under @xmath128 and with higher order terms
omitted , @xmath126 transforms to ( [ eq : n ] ) with @xmath129 since the nordmark map keeps only leading order terms for @xmath90 near @xmath14 , terms that are linear in @xmath47 are omitted since they are dominated by @xmath130 , refer to @xcite for a further discussion . in ( [ eq : coordinatechange ] ) we require
@xmath131 and @xmath132 .
these represent non - degeneracy conditions for the grazing bifurcation .
( 15,5 ) ( 0,0 ) ) .
this system exhibits a _
regular _ grazing bifurcation with a square - root singularity because with @xmath133 and @xmath134 the equations of motion are discontinuous at the grazing point .
[ fig : vibroimpactcompliant ] , title="fig:",height=188 ] ( 6.78,.08)@xmath135 ( 8.44,.08)@xmath136 ( 9.64,.08)@xmath9 ( 3.78,3.38)@xmath137 ( 3.78,1.4)@xmath138 ( 11.58,3.64)@xmath139 ( 11.58,1.18)@xmath140 ( 9.34,.63)@xmath141 ( 6.58,3.88)block ( 11.02,4.12)support ( 10.38,4.64)stop    to motivate and illustrate our results for stochastic nordmark maps , we consider the prototypical vibro - impacting system shown in fig .
[ fig : vibroimpactcompliant ] and studied in @xcite .
this system consists of a harmonically forced one - degree - of - freedom linear oscillator that experiences compliant ( or soft ) impacts with a support , and we use the following non - dimensionalised equations to model the dynamics : @xmath142 here @xmath143 denotes the location of the block , which has the equilibrium position @xmath144 .
a rigid stop prevents the support reaching a position with @xmath10 , and prestresses the support by a distance , @xmath134 .
the constants @xmath137 , @xmath138 , @xmath139 and @xmath140 represent non - dimensionalised spring and damping coefficients for the oscillator and support .
we neglect the mass of the support , ignore energy loss at impacts , and assume that whenever the block is not in contact with the support , the support is located at @xmath9 .
experiments of simple vibro - impacting systems with compliant impacts have shown that low - dimensional models such as ( [ eq : impactoscillator ] ) can quantitatively match the physically observed dynamics near grazing bifurcations @xcite .    here
we treat the forcing amplitude , @xmath145 , as the primary bifurcation parameter , and assume @xmath146 .
the steady - state solution ( behaviour in the limit @xmath147 ) to ( [ eq : impactoscillator ] ) with @xmath10 is @xmath148 when @xmath149 , where @xmath150 the maximum value of @xmath151 over one period is negative , and so @xmath151 is an attracting non - impacting periodic orbit of ( [ eq : impactoscillator ] ) .
the critical value , @xmath152 , is a grazing bifurcation at which @xmath151 has unit amplitude and attains the value @xmath153 at times @xmath154 , for @xmath155 , where @xmath156 and @xmath157 .    to convert ( [ eq : impactoscillator ] ) to the general form ( [ eq : f ] ) , we define @xmath158 here the phase space of ( [ eq : impactoscillator ] ) with ( [ eq : vw ] ) is isomorphic to @xmath159 , rather than @xmath15 , but this does not affect the bifurcation structure near grazing . for ( [ eq : impactoscillator ] ) with ( [ eq : vw ] ) , the coefficients in ( [ eq : alphabetagamma ] ) , which describe the behaviour of the system near the grazing point , are given by @xmath160 and by evaluating ( [ eq : c ] ) with ( [ eq : odecoeffsosc ] ) we obtain @xmath161 in addition , from the general solution to ( [ eq : impactoscillator ] ) we find that the coefficients in the global map ( [ eq : g ] ) are given by @xmath162 [ eq : mapcoeffosc ]
to model noise and uncertainties we use the one - dimensional ornstein - uhlenbeck process @xmath163 where @xmath164 are constants and @xmath165 is a standard brownian motion . given an initial value @xmath166 , at any positive time @xmath167 is a gaussian random variable with mean and variance @xmath168 = \xi_0 \,{\rm e}^{\frac{-t}{\nu } } \ ; , \qquad { \rm var } \left [ \xi(t ) | \xi(0 ) = \xi_0 \right ] = \frac{\ee^2}{2 \nu } \left ( 1 - { \rm e}^{\frac{-2 t}{\nu } } \right ) \;. \label{eq : meanvarxi}\ ] ] in the limit @xmath147 , @xmath169 $ ] , where @xmath170 $ ] denotes the gaussian distribution of mean @xmath46 and variance @xmath171 . the _ correlation time _ of ( [ eq : xi ] ) , defined as @xmath172}{{\rm var } \left [ \xi(0 ) \right ] } \,dt$ ] , with @xmath173 $ ] , is equal to @xmath174 .    in our context , @xmath167 is coloured noise and the parameter @xmath67 governs the size of the noise . unlike white noise , @xmath167 has an inherent time - scale , @xmath174 , and is suitable for modelling various types of uncertainties in mechanical systems , such as background vibrations @xcite . in the white noise limit ,
forcing by @xmath167 becomes a diffusion process @xmath175 .
we first consider the following stochastic perturbation of ( [ eq : f ] ) , @xmath176 where @xmath167 is given by ( [ eq : xi ] ) . in ( [ eq : fa ] ) randomness is present in the switching condition , while evolution between switching events remains deterministic .
we expect ( [ eq : fa ] ) to be applicable to a wide variety of piecewise - smooth systems . for the vibro - impacting system of  [ sec : oscillator ] , @xmath167 may capture uncertainties in the point at which contact between the block and support occurs or is lost . for switched control systems , @xmath167 may correspond to measurement errors that produce variability in evaluations of switching rules @xcite .    here
we consider orbits of ( [ eq : fa ] ) that are close to the grazing periodic orbit of ( [ eq : f ] ) .
orbits of ( [ eq : fa ] ) near grazing only spend short lengths of time in the region @xmath11 while passing near the origin , and for simplicity we suppose that the value of @xmath174 is large compared to such times . in this case
it is reasonable to approximate @xmath167 by a constant while an orbit is near the origin . with this approximation ,
the sum @xmath177 does not switch sign more than twice as the orbit passes near the origin , which substantially simplifies our calculations below .
we let @xmath178 denote the value of @xmath167 during the @xmath179 instance that the orbit of ( [ eq : fa ] ) passes near the origin .
the time between between consecutive traversals near the origin is well - approximated by the period of the grazing periodic orbit , call it @xmath180 . with this approximation , @xmath181 \;.
\label{eq : xiiab}\ ] ] to derive the stochastic version of ( [ eq : n ] ) for ( [ eq : fa ] ) with ( [ eq : xiiab ] ) , we first derive the induced stochastic discontinuity map . here
condition ( [ eq : pipcondition ] ) is useful , as it implies that an orbit governed by @xmath6 attains a local maximum value of @xmath74 at an intersection with @xmath78 . for @xmath182 , we conclude that @xmath183 as the orbit passes near the origin , so that @xmath184 . if instead
@xmath185 , then the discontinuity map @xmath186 is given by ( [ eq : u4w42 ] ) except that @xmath187 appears inside the square root because this quantity represents the distance from the switching condition .
that is , @xmath188 where @xmath189 . by combining ( [ eq : u4w42a ] ) with the global map @xmath122
, applying the coordinate change ( [ eq : coordinatechange ] ) , and dropping higher order terms , we obtain @xmath190 @xmath62 is the stochastic nordmark map corresponding to ( [ eq : fa ] ) . notice that randomness in the switching condition of ( [ eq : fa ] ) has translated to randomness in both the switching condition of ( [ eq : n ] ) and in the image of the map with @xmath25 . in contrast , a piecewise - linear map for which randomness is present purely in the switching condition is studied in @xcite .    in order to fairly compare @xmath62 with other stochastic versions of ( [ eq : n ] ) in  [ sec : dynamics ] , we estimate the effective size of the stochastic contribution for our illustrative parameters values of the prototypical system ( [ eq : impactoscillator ] ) and a representative value of @xmath72 .
this motivates us to express @xmath67 in terms of a scaled parameter @xmath191 , and to obtain comparable stochastic contributions for fixed @xmath192 in the different cases .
the square - root term of @xmath62 is @xmath193 , where @xmath194 .
if @xmath195 is large relative to @xmath196 , then this term is well - approximated by @xmath197 , and so @xmath198 estimates the additive stochastic contribution to @xmath62 . for the impact oscillator with the parameter values of fig .
[ fig : detbifdiag ] ( here @xmath199 and @xmath200 ) and using @xmath201 ( corresponding to the value of @xmath25 in fig .
[ fig : detbifdiag ] for @xmath72 ) , this quantity is approximately @xmath202 . with @xmath203 ( used in  [ sec : dynamics ] ) , the standard derivation of the stochastic contribution is approximately @xmath204 .
therefore , for @xmath191 , where @xmath205 the standard deviation of the stochastic contribution is approximately @xmath206 when @xmath73 .
( 6.8,5.1 ) ( 0,0 ) , ( [ eq : na ] ) , with a numerical solution to ( [ eq : fa ] ) .
the three groups of purple dots were obtained by numerically solving ( [ eq : fa ] ) with ( [ eq : xi ] ) , @xmath203 and @xmath207 , for the vibro - impacting system ( [ eq : impactoscillator ] ) with ( [ eq : vw ] ) and @xmath36 ( as in fig .
[ fig : detbifdiag ] ) and @xmath208 ( which corresponds to @xmath72 ) . more precisely , @xmath209 points on @xmath78 ( labelled @xmath80 in fig .
[ fig : grazschem3d ] ) were obtained by numerically solving ( [ eq : fa ] ) , and these were transformed to @xmath20-coordinates by applying ( [ eq : coordinatechange ] ) and ( [ eq : vw ] ) to produce the purple dots .
the three groups of black dots are @xmath209 iterates of @xmath62 with ( [ eq : xiiab ] ) and parameter values chosen to match the vibro - impacting system ( specifically , @xmath203 , @xmath207 , @xmath210 , @xmath72 , @xmath33 , @xmath34 , @xmath35 , @xmath199 and @xmath200 ) .
the deterministic @xmath38-cycles of ( [ eq : fa ] ) and ( [ eq : na ] ) are shown with triangles . [
fig : comparea_03 ] , title="fig:",height=192 ] ( 3.76,0)@xmath47 ( 0,3.5)@xmath211    to illustrate the accuracy of @xmath62 , fig .  [ fig : comparea_03 ] compares iterates of @xmath62 ( black dots ) with a numerical solution to ( [ eq : fa ] ) for the vibro - impacting system of  [ sec : oscillator ] ( purple dots ) using @xmath207 . for the given parameter values ,
the system has an attracting @xmath38-cycle in the absence of noise .
for this reason , both sets of points are grouped about the @xmath38-cycle .
the two sets of points are slightly separated .
this is because the form of the deterministic nordmark map does not include higher order terms of the true return map , as observed by the separation of the values taken by the deterministic @xmath38-cycles shown in fig .
[ fig : comparea_03 ] .
the size and shape of the spread of the two sets of randomly generated points are similar , as is their location relative to the deterministic values of the map .
this demonstrates that @xmath62 can accurately capture the stochastic dynamics of ( [ eq : fa ] ) .
a more precise characterisation of the accuracy of @xmath62 is beyond the scope of this paper .
next we consider the case where randomness and uncertainty in ( [ eq : f ] ) is associated with @xmath7 . for mechanical systems with impacts , this corresponds to variability in the evolution of the system during an impact .
for simplicity we include noise in only the @xmath75-component of @xmath7 , that is @xmath212 where @xmath167 is given by ( [ eq : xi ] ) . indeed , for the vibro - impacting system of fig .
[ fig : vibroimpactcompliant ] , if noise is incorporated into the force on the block when it is in contact with the support , then the equations of motion may be put in the form ( [ eq : fbc ] ) . with noise added to the @xmath74-component of @xmath7 ( or @xmath6 )
, orbits may cross @xmath153 many times in a short time frame which makes the system substantially more difficult to analyse .
we leave such considerations for future work .    as in  [ sub : s1 ] , we consider near - grazing orbits and assume that the value of @xmath174 is much larger than the time each orbit spends in the region @xmath11 .
in this case it is reasonable to treat @xmath167 as constant while @xmath11 . during the @xmath179 instance that an orbit passes near the origin ,
we approximate @xmath167 by @xmath178 , distributed according to ( [ eq : xiiab ] ) . in this scenario
the three components of the discontinuity map ( [ eq : v2w2])-([eq : u4w4 ] ) are unchanged except that @xmath213 is replaced by @xmath214 in ( [ eq : v3w3 ] ) ( because the @xmath75-component of the system with @xmath11 is given by @xmath215 , see ( [ eq : alphabetagamma ] ) ) . by combining ( [ eq : v2w2])-([eq : u4w4 ] )
we find that for @xmath110 the discontinuity map is given by @xmath216 and therefore the corresponding stochastic nordmark map is @xmath217 notice that with @xmath218 , @xmath63 is identical to ( [ eq :
we can write the stochastic component of @xmath63 as @xmath219 , where @xmath220 . with the parameter values of the impact oscillator ( [ eq : odecoeffsosc ] ) , and @xmath221
, we have @xmath222 .
therefore the noise provides a multiplicative stochastic contribution of approximately @xmath223 , ignoring signs . in order to compare the effect of the noise to the other cases , we write @xmath224 , choosing @xmath225 so that the standard deviation of the stochastic contribution is @xmath206 when @xmath73 . for @xmath201 and @xmath203 ( as in  [ sub : s1 ] ) , the standard deviation of @xmath223 is approximately @xmath226 , so we therefore choose @xmath227    ( 6.8,5.1 ) ( 0,0 ) with a numerical solution to ( [ eq : fbc ] ) .
the three groups of purple dots were obtained by numerically solving ( [ eq : fbc ] ) with ( [ eq : xi ] ) , @xmath203 and @xmath228 , for the vibro - impacting system ( [ eq : impactoscillator ] ) with ( [ eq : vw ] ) using the same parameter values as in fig .
[ fig : comparea_03 ] .
the three groups of black dots are @xmath209 iterates of @xmath63 with ( [ eq : xiiab ] ) and parameter values matching those of the vibro - impacting system ( refer to the caption of fig .
[ fig : comparea_03 ] ) .
the deterministic @xmath38-cycles of ( [ eq : fa ] ) and ( [ eq : na ] ) are shown with triangles .
[ fig : compareb_03 ] , title="fig:",height=192 ] ( 3.76,0)@xmath47 ( 0,3.5)@xmath211    fig .
[ fig : compareb_03 ] compares iterates of @xmath63 to intersections with @xmath78 of a numerical solution to ( [ eq : fbc ] ) using @xmath228 . as expected the two sets of points
are similarly distributed about the deterministic @xmath38-cycle .
lastly we consider ( [ eq : fbc ] ) in the white noise limit , @xmath229 . in this case
( [ eq : fbc ] ) reduces to a diffusion process forced by white noise , specifically @xmath230 is replaced by @xmath175 . by using ( [ eq : alphabetagamma ] ) to expand @xmath7 , ( [ eq : fbc ] ) for @xmath11 may be written as the three - dimensional stochastic differential equation @xmath231 to obtain a stochastic nordmark map corresponding to this scenario , we first derive the stochastic version of the middle component of the discontinuity map ( [ eq : v3w3 ] ) . given an initial point @xmath232 , where @xmath233 and @xmath234 are small , the desired values of @xmath235 and @xmath236 are given by the point @xmath237 of _ first return _ for the stochastic process ( [ eq : frc ] ) to @xmath153 .
first return or first passage problems are an important class of theoretical problems in stochastic calculus with applications traditionally in finance and chemical kinetics @xcite .
we approximate ( [ eq : frc ] ) by keeping only the leading order contributions , i.e. @xmath238 together with @xmath239 . with this approximation
we are able to provide an explicit expression for the joint probability density function of the return location and time .
a formal justification for the omission of the higher order terms is left for future work .
we introduce the change of variables @xmath240 this puts ( [ eq : frc2 ] ) in a standard form studied in @xcite , @xmath241 where @xmath242 specifically , @xmath243 is a diffusion process with constant drift , and @xmath244 .
therefore @xmath245 may be interpreted as _ integrated brownian motion with constant drift_.    we let @xmath246 denote the joint probability density function for the first return of ( [ eq : dpdq ] ) to @xmath247 , at a time @xmath248 , and location @xmath249 . in @xcite , mckean derived an explicit expression for @xmath250 in the case of zero drift by computing the inverse kontorovich - lebedev transform of the corresponding renewal equation . in @xcite , atkinson and clifford extended this result to the case of non - zero drift by applying the radon - nikodyn derivative .
specifically @xmath251 \right ) { \rm erf } \left ( \frac{\sqrt{6 h}}{\sqrt{\varrho r } } \right ) \ ; , \label{eq : f}\ ] ] where @xmath252 is the error function .
the constant @xmath253 governs the shape of @xmath250 .
the limit @xmath254 corresponds to the deterministic case , for which @xmath255 and @xmath256 . with a small value of @xmath257 , @xmath250 is roughly gaussian .
in contrast , the limit @xmath258 corresponds to the case of no drift , as in @xcite , or to the limit @xmath259 . in this limit
the marginal probability density function for @xmath260 is asymptotically proportional to @xmath261 , for large @xmath260 , and so is long - tailed @xcite .    in view of the scaling ( [ eq : pqs ] ) , the stochastic version of ( [ eq : v3w3 ] ) corresponding to ( [ eq : frc ] ) is given by @xmath262 where @xmath58 and @xmath59 have the joint probability density function ( [ eq : f ] ) . by combining ( [ eq : v3w3c ] ) with ( [ eq : v2w2 ] ) and ( [ eq : u4w4 ] ) we obtain @xmath263 which represents the stochastic version of the discontinuity map for points with @xmath110 . then using ( [ eq : u4w4c ] )
we arrive at the following stochastic nordmark map @xmath264 in ( [ eq : nc ] ) we treat each pair @xmath265 as a two - dimensional stochastic random variable with probability density function @xmath266 , where @xmath267 which results from combining ( [ eq : varrho ] ) with @xmath268 , ( [ eq : v2w2 ] ) , and @xmath269 , ( [ eq : coordinatechange ] ) .
we now estimate the size of the stochastic contribution in @xmath64 .
the leading order stochastic component of @xmath64 is @xmath270 , where @xmath271 . with ( [ eq : odecoeffsosc ] ) and @xmath221 , we can write @xmath272 , where the second term represents the multiplicative stochastic contribution , ignoring the sign of @xmath127 , because the deterministic values of @xmath59 and @xmath58 are @xmath273 and @xmath274 respectively .    with a small value of @xmath257 ( @xmath275 is suitable ) , @xmath246 is approximately gaussian because the effective noise amplitude in ( [ eq : dpdq ] ) is small . by ( [ eq : varrho2 ] )
, this approximation is valid when , roughly speaking , @xmath67 is not too large and @xmath195 is not too small .
from ( [ eq : f ] ) we determine the covariance matrix of the gaussian approximation to be @xmath276 and it follows that in this approximation the linear combination @xmath277 has standard deviation @xmath278 .
the standard deviation of the stochastic contribution in @xmath64 is therefore approximately @xmath279 . following the previous cases , we write @xmath280 , and choose @xmath281 so that the standard deviation of the stochastic contribution is @xmath206 when @xmath73 . for @xmath64 ,
this quantity is @xmath282 , using ( [ eq : varrho2 ] ) , the parameter values from fig .
[ fig : detbifdiag ] , and @xmath201
. therefore we let @xmath283 here @xmath284 when @xmath73 , and so for these values the gaussian approximation to ( [ eq : f ] ) is justified .
( 6.8,5.1 ) ( 0,0 ) with a numerical solution to ( [ eq : fbc ] ) .
the three groups of purple dots were obtained by numerically solving ( [ eq : fbc ] ) with ( [ eq : xi ] ) , @xmath229 ( in which case @xmath230 is replaced by @xmath285 ) and @xmath286 , for the vibro - impacting system ( [ eq : impactoscillator ] ) with ( [ eq : vw ] ) using the same parameter values as in fig .
[ fig : comparea_03 ] .
the three groups of black dots are @xmath209 iterates of @xmath64 with ( [ eq : f ] ) and ( [ eq : varrho2 ] ) and parameter values matching those of the vibro - impacting system .
the deterministic @xmath38-cycles of ( [ eq : fa ] ) and ( [ eq : na ] ) are shown with triangles . [
fig : comparec_03 ] , title="fig:",height=192 ] ( 3.76,0)@xmath47 ( 0,3.5)@xmath211    fig .
[ fig : comparec_03 ] compares iterates of @xmath64 to intersections with @xmath78 of a numerical solution to ( [ eq : fbc ] ) using @xmath286 . as with the previous two figures
, this shows that @xmath64 can accurately capture the stochastic dynamics of ( [ eq : fbc ] ) .
in this section we explore the dynamics of the three stochastic nordmark maps , @xmath62 , @xmath63 and @xmath64 , and discuss how the different forms of these maps is evident in their dynamical behaviour .
to briefly summarise , @xmath62 applies to the system with stochastic switching ( [ eq : fa ] ) , whereas @xmath63 and @xmath64 apply to ( [ eq : fbc ] ) . for @xmath62 and @xmath63
it is assumed that the value of @xmath174 ( the correlation time of the noise ) is large relative to the times that orbits spend in @xmath11 . as these times are rarely larger than @xmath287 for the parameter values considered here , we take @xmath203 in @xmath62 and @xmath63 to ensure that the correlation time is large enough .
the values of @xmath178 in @xmath62 and @xmath63 are distributed according to ( [ eq : xiiab ] ) .
@xmath64 corresponds to the limit @xmath56 . in @xmath64 , @xmath58 and @xmath59
are distributed according to ( [ eq : f ] ) and depend on the value of @xmath288 ( [ eq : varrho2 ] ) .    for each @xmath68
we have written @xmath289 where the @xmath71 are given by ( [ eq : tildeee1 ] ) , ( [ eq : tildeee2 ] ) and ( [ eq : tildeee3 ] ) .
these values have been chosen such that for a given value of @xmath192 , the size of the stochastic contribution in @xmath62 , @xmath63 and @xmath64 is roughly the same , at least when @xmath72 . for @xmath46 close to @xmath290 ,
the stochastic contributions are noticeably different for these choices of @xmath71 .    in  [ sub : mudependence ]
we look at stochastic bifurcation diagrams in order to obtain a basic understanding of how the stochastic dynamics differs with the value of @xmath46 . in the subsequent parts of this section we study two - dimensional invariant densities in order to gain a deeper understanding of the dynamics .
( 12,18.4 ) ( 0,12.4 ) ( panel a ) , @xmath63 ( panel b ) , and @xmath64 ( panel c ) , with @xmath33 , @xmath34 , and @xmath35 .
these parameter values correspond to the vibro - impacting system ( [ eq : impactoscillator ] ) , with @xmath36 .
the black dots are iterates of @xmath62 , @xmath63 and @xmath64 with transient points omitted .
the noise amplitudes are given by ( [ eq : tildeee1 ] ) , ( [ eq : tildeee2 ] ) and ( [ eq : tildeee3 ] ) with @xmath73 , and @xmath203 for panels a and b ( panel c corresponds to @xmath229 ) . in each panel the deterministic bifurcation diagram ( fig .
[ fig : detbifdiag ] ) is shown in blue .
[ fig : noisybifdiag ] , title="fig:",height=226 ] ( 0,6.2 ) ( panel a ) , @xmath63 ( panel b ) , and @xmath64 ( panel c ) , with @xmath33 , @xmath34 , and @xmath35 .
these parameter values correspond to the vibro - impacting system ( [ eq : impactoscillator ] ) , with @xmath36 .
the black dots are iterates of @xmath62 , @xmath63 and @xmath64 with transient points omitted .
the noise amplitudes are given by ( [ eq : tildeee1 ] ) , ( [ eq : tildeee2 ] ) and ( [ eq : tildeee3 ] ) with @xmath73 , and @xmath203 for panels a and b ( panel c corresponds to @xmath229 ) . in each panel the deterministic bifurcation diagram ( fig .
[ fig : detbifdiag ] ) is shown in blue .
[ fig : noisybifdiag ] , title="fig:",height=226 ] ( 0,0 ) ( panel a ) , @xmath63 ( panel b ) , and @xmath64 ( panel c ) , with @xmath33 , @xmath34 , and @xmath35 .
these parameter values correspond to the vibro - impacting system ( [ eq : impactoscillator ] ) , with @xmath36 .
the black dots are iterates of @xmath62 , @xmath63 and @xmath64 with transient points omitted .
the noise amplitudes are given by ( [ eq : tildeee1 ] ) , ( [ eq : tildeee2 ] ) and ( [ eq : tildeee3 ] ) with @xmath73 , and @xmath203 for panels a and b ( panel c corresponds to @xmath229 ) . in each panel the deterministic bifurcation diagram ( fig .
[ fig : detbifdiag ] ) is shown in blue .
[ fig : noisybifdiag ] , title="fig:",height=226 ] ( 6.2,12.4)@xmath46 ( 6.2,6.2)@xmath46 ( 6.2,0)@xmath46 ( 0,16.2)@xmath47 ( 0,10)@xmath47 ( 0,3.8)@xmath47 ( 1.7,18.1)*a * ( 1.7,11.9)*b * ( 1.7,5.7)*c *    fig .
[ fig : noisybifdiag ] shows stochastic versions of the bifurcation diagram shown in fig .
[ fig : detbifdiag ] for the three stochastic nordmark maps .
as expected , the noise blurs the bifurcation diagram . in panel
a , which corresponds to the map @xmath62 , for values of @xmath46 very close to zero ( say @xmath291 ) the points are relatively highly spread .
this is because here the deterministic map has an attractor near @xmath292 , so in @xmath62 the sign of @xmath195 is often different to the sign of @xmath293 .
that is , the choice of the half - map of @xmath62 is regularly determined by @xmath178 rather than @xmath195 .
furthermore , as shown in
[ sub : s1 ] , the leading order component of the stochastic contribution to the right half - map of @xmath62 is inversely proportional to @xmath66 , which for very small @xmath46 is large relative to its value for @xmath46 away from zero . in contrast , with @xmath294 say , the underlying attracting @xmath38-cycle is sufficiently far from @xmath292 so that the sign of @xmath195 rarely differs from that of @xmath293 .
the points are randomly distributed    about the @xmath38-cycle due to noise in the right half - map of @xmath62 .
the size of the deviation decreases with increasing @xmath46 , because the strength of attraction of the @xmath38-cycle increases with @xmath46 .    in panels
b and c , which correspond to @xmath63 and @xmath64 respectively , the bifurcation diagrams show no variability for @xmath29 .
this is because for @xmath29 , @xmath32 ( [ eq : xlyl ] ) is a fixed point of @xmath63 and @xmath64 . for @xmath295 , the size of deviations about the underlying attracting @xmath39-cycle increases with @xmath46 .
this is primarily because the coefficient of the @xmath130-term of @xmath63 and @xmath64 is random , and for the @xmath39-cycle this value of @xmath47 increases with @xmath46 . for @xmath294 ,
the size of deviations varies little with @xmath46 because the increased variability caused by a larger value of @xmath47 is balanced by the fact that the strength of attraction of the @xmath38-cycle increases .
panels b and c are similar , suggesting that the value of @xmath174 has little effect on the long - term dynamics , although panel c shows slightly more variability for very small values of @xmath40 .
( 13.8,10.2 ) ( 0,4.7 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
the value of the density is indicated by colour ( dark red  the maximum value of the density ; dark blue  zero ) .
panel b shows the invariant density at three times the noise amplitude as panel a. panel c plots @xmath296  the fraction of instances that points return to @xmath297 in @xmath298 iterations ( [ eq : sigma ] )  against the noise amplitude .
[ fig : invdensitya_03 ] , title="fig:",height=192 ] ( 7,4.7 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
the value of the density is indicated by colour ( dark red  the maximum value of the density ; dark blue
panel b shows the invariant density at three times the noise amplitude as panel a. panel c plots @xmath296  the fraction of instances that points return to @xmath297 in @xmath298 iterations ( [ eq : sigma ] )  against the noise amplitude .
[ fig : invdensitya_03 ] , title="fig:",height=192 ] ( 2.4,0 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
the value of the density is indicated by colour ( dark red  the maximum value of the density ; dark blue
panel b shows the invariant density at three times the noise amplitude as panel a. panel c plots @xmath296  the fraction of instances that points return to @xmath297 in @xmath298 iterations ( [ eq : sigma ] )  against the noise amplitude .
[ fig : invdensitya_03 ] , title="fig:",height=170 ] ( 3.42,10)@xmath73 ( 10.42,10)@xmath299 ( 3.76,4.7)@xmath47 ( .3,8.2)@xmath211 ( 10.76,4.7)@xmath47 ( 7.3,8.2)@xmath211 ( 7.2,0)@xmath192 ( 8,3.44)@xmath300 ( 8.87,1.45)@xmath301 ( 9.33,2)@xmath302 ( .5,9.9)*a * ( 7.5,9.9)*b * ( 2.4,4.2)*c *    ( 13.8,10.2 ) ( 0,4.7 ) with the same parameter values as fig .  [ fig : noisybifdiag ] and @xmath72 .
panel b shows the invariant density with @xmath299 , and panel c is a plot of the fractions @xmath296 ( [ eq : sigma ] ) .
[ fig : invdensityb_03 ] , title="fig:",height=192 ] ( 7,4.7 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
panel b shows the invariant density with @xmath299 , and panel c is a plot of the fractions @xmath296 ( [ eq : sigma ] ) .
[ fig : invdensityb_03 ] , title="fig:",height=192 ] ( 2.4,0 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
panel b shows the invariant density with @xmath299 , and panel c is a plot of the fractions @xmath296 ( [ eq : sigma ] ) .
[ fig : invdensityb_03 ] , title="fig:",height=170 ] ( 3.42,10)@xmath73 ( 10.42,10)@xmath299 ( 3.76,4.7)@xmath47 ( .3,8.2)@xmath211 ( 10.76,4.7)@xmath47 ( 7.3,8.2)@xmath211 ( 7.2,0)@xmath192 ( 8,3.44)@xmath300 ( 6.66,1.64)@xmath301 ( 8.78,2)@xmath302 ( .5,9.9)*a * ( 7.5,9.9)*b * ( 2.4,4.2)*c *    ( 13.8,10.2 ) ( 0,4.7 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
panel b shows the invariant density with @xmath299 , and panel c is a plot of the fractions @xmath296 ( [ eq : sigma ] ) .
[ fig : invdensityc_03 ] , title="fig:",height=192 ] ( 7,4.7 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
panel b shows the invariant density with @xmath299 , and panel c is a plot of the fractions @xmath296 ( [ eq : sigma ] ) .
[ fig : invdensityc_03 ] , title="fig:",height=192 ] ( 2.4,0 ) with the same parameter values as fig .
[ fig : noisybifdiag ] and @xmath72 .
panel b shows the invariant density with @xmath299 , and panel c is a plot of the fractions @xmath296 ( [ eq : sigma ] ) .
[ fig : invdensityc_03 ] , title="fig:",height=170 ] ( 3.42,10)@xmath73 ( 10.42,10)@xmath299 ( 3.76,4.7)@xmath47 ( .3,8.2)@xmath211 ( 10.76,4.7)@xmath47 ( 7.3,8.2)@xmath211 ( 7.2,0)@xmath192 ( 7.3,3.55)@xmath300 ( 8,1.7)@xmath301 ( 8.93,2.03)@xmath302 ( .5,9.9)*a * ( 7.5,9.9)*b * ( 2.4,4.2)*c *    figs .
[ fig : invdensitya_03]-[fig : invdensity_m002 ] show two - dimensional invariant densities of @xmath62 , @xmath63 and @xmath64 . by assuming ergodicity ,
invariant densities were computed on a @xmath303 grid of @xmath47 and @xmath211 values from @xmath304 consecutive iterates of a single orbit with transient points omitted .
let us first explain panel c of figs .
[ fig : invdensitya_03]-[fig : invdensityc_03 ] . given a sample orbit of @xmath62 , @xmath63 or @xmath64 , for each point with @xmath25
we let @xmath305 be the smallest positive integer for which @xmath306 , as in @xcite .
@xmath305 represents the number of iterations required for a return to @xmath307 from the point @xmath308 .
numerically we can compute a large number of values of @xmath305 ( say @xmath309 of them ) . then for each @xmath310
, we let @xmath296 denote the fraction of the @xmath305 that are equal to @xmath298 , i.e. @xmath311 where @xmath312 if @xmath313 , and @xmath314 otherwise .
[ fig : invdensitya_03]-[fig : invdensityc_03 ] correspond to @xmath72 for which there is an underlying attracting @xmath38-cycle . therefore with small noise , @xmath315 , and for all @xmath316 , @xmath317 .
[ fig : invdensitya_03 ] corresponds to the map @xmath62 . in panel
a , the size of the noise is relatively small , so iterates of @xmath62 follow close to the @xmath38-cycle .
the invariant density is well - approximated by a scaled sum of three gaussian densities centred at each point of the @xmath38-cycle .
about the point with @xmath318 , the density is stretched substantially more in @xmath47-direction than in the @xmath211-direction .
this is because points with @xmath318 have likely just undergone an iteration under the right half - map of @xmath62 which is stochastic in the @xmath47-component but not the @xmath211-component .
the stretching around other iterates is then a consequence of iterating under @xmath62 with @xmath319 . with larger values of @xmath192
, it is relatively common for the orbit to return to @xmath307 in a number of iterations other than three , fig .
[ fig : invdensitya_03]-c . for this reason the invariant density displays additional characteristics .
for instance with @xmath299 , the orbit returns to @xmath307 in two iterations almost @xmath320 of the time .
consequently , a substantial part of the invariant density centred roughly about the point of the @xmath38-cycle with @xmath321 , lies in @xmath307 , fig .
[ fig : invdensitya_03]-b . the invariant density in panel b also has a small component with @xmath322 .
this is due to points of the orbit with small values of @xmath307 mapping under the left half - map of @xmath62 due to the noise ( i.e.  returning to @xmath307 in only one iteration ) .
[ fig : invdensityb_03 ] illustrates @xmath63 using the same parameter values . again with small noise
the invariant density is roughly gaussian about each point of the @xmath38-cycle , whereas for relatively large noise iterates often cross into @xmath307 prematurely causing the invariant density to take an irregular shape .
when @xmath299 , points of the orbit that do not return to @xmath307 in three iterations , usually return to @xmath307 in two iterations .
[ fig : invdensityc_03 ] corresponds to @xmath64 and is similar to the previous figure .
this indicates that the correlation time @xmath174 has little effect on these pictures , although the invariant density has a slightly different shape when @xmath299 .
( 13.8,10.4 ) ( 0,5.3)-cycle of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath323 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , using the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath323 .
[ fig : invdensity_001 ] , title="fig:",height=192 ] ( 7,5.3)-cycle of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath323 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , using the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath323 .
[ fig : invdensity_001 ] , title="fig:",height=192 ] ( 0,0)-cycle of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath323 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , using the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath323 .
[ fig : invdensity_001 ] , title="fig:",height=192 ] ( 7,0)-cycle of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath323 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , using the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath323 .
[ fig : invdensity_001 ] , title="fig:",height=192 ] ( 3.76,5.3)@xmath47 ( .3,8.8)@xmath211 ( 10.76,5.3)@xmath47 ( 7.3,8.8)@xmath211 ( 3.76,0)@xmath47 ( .3,3.5)@xmath211 ( 10.76,0)@xmath47 ( 7.3,3.5)@xmath211 ( .6,10.3)*a * ( 7.6,10.3)*b * ( .6,5)*c * ( 7.6,5)*d *    fig .  [ fig : invdensity_001 ] shows invariant densities of @xmath62 , @xmath63 and @xmath64 for parameter values closer to the grazing bifurcation than the previous three figures , specifically @xmath323 . at this value of @xmath46 , there is an underlying attracting @xmath39-cycle , panel a. the noise amplitudes are given by ( [ eq : tildeee1 ] ) , ( [ eq : tildeee2 ] ) and ( [ eq : tildeee3 ] ) with @xmath73 .
recall , these amplitudes were chosen such that the size of the noise response of the three maps is roughly the same for larger values of @xmath46 . here
, however , the size of noise response differs substantially . in panels c and d , which correspond to @xmath63 and @xmath64 respectively ,
the invariant density is approximately a scaled sum of four gaussians about each point of the @xmath39-cycle .
the invariant density in panel d , corresponding to @xmath229 , is noticeably larger than that of panel c.    in panel b , which corresponds to @xmath62 , the noise has a substantial effect because the switching condition of @xmath62 is stochastic , and many points of @xmath62 fall close to @xmath292 .
the invariant density has a small @xmath324-shaped component in @xmath297 corresponding to consecutive points of the orbit mapping under the left half - map of @xmath62 .
the part of the invariant density for @xmath319 and @xmath325 corresponds to images of points under the right half - map of @xmath62 , and is bimodal because the invariant density has roughly two components in @xmath307 .
( 13.8,10.4 ) ( 0,5.3)-cycle and a numerically computed attractor of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath326 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , with the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath326 .
[ fig : invdensity_0145 ] , title="fig:",height=192 ] ( 7,5.3)-cycle and a numerically computed attractor of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath326 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , with the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath326 .
[ fig : invdensity_0145 ] , title="fig:",height=192 ] ( 0,0)-cycle and a numerically computed attractor of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath326 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , with the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath326 .
[ fig : invdensity_0145 ] , title="fig:",height=192 ] ( 7,0)-cycle and a numerically computed attractor of ( [ eq : n ] ) with the same parameter values as fig .
[ fig : detbifdiag ] and @xmath326 .
panels b , c and d show the invariant densities of @xmath62 , @xmath63 and @xmath64 , respectively , with the same parameter values as fig .  [
fig : noisybifdiag ] and @xmath326 .
[ fig : invdensity_0145 ] , title="fig:",height=192 ] ( 3.76,5.3)@xmath47 ( .3,8.5)@xmath211 ( 10.76,5.3)@xmath47 ( 7.3,8.5)@xmath211 ( 3.76,0)@xmath47 ( .3,3.2)@xmath211 ( 10.76,0)@xmath47 ( 7.3,3.2)@xmath211 ( .4,10.3)*a * ( 7.4,10.3)*b * ( .4,5)*c * ( 7.4,5)*d *    ( 13.8,5.1 ) ( 0,0 ) ( [ eq : xlyl ] ) and the attracting @xmath38-cycle of @xmath327 with the same parameter values as fig .
[ fig : detbifdiag ] , except @xmath328 and @xmath329 .
here @xmath330 , @xmath34 and @xmath35 .
panel b shows the invariant density of the corresponding map @xmath62 , with @xmath203 and @xmath331 .
[ fig : invdensity_m002 ] , title="fig:",height=192 ] ( 7,0 ) ( [ eq : xlyl ] ) and the attracting @xmath38-cycle of @xmath327 with the same parameter values as fig .
[ fig : detbifdiag ] , except @xmath328 and @xmath329 .
here @xmath330 , @xmath34 and @xmath35 .
panel b shows the invariant density of the corresponding map @xmath62 , with @xmath203 and @xmath331 .
[ fig : invdensity_m002 ] , title="fig:",height=192 ] ( 3.76,0)@xmath47 ( .3,3.4)@xmath211 ( 10.76,0)@xmath47 ( 7.3,3.4)@xmath211 ( .4,5)*a * ( 7.4,5)*b *    with @xmath326 in fig .
[ fig : detbifdiag ] , there is an attracting @xmath39-cycle and an apparently chaotic attractor .
these are shown in fig .
[ fig : invdensity_0145]-a . as with smooth maps @xcite , in the presence of noise orbits commonly dwell near the attractors for relatively long periods of time , and switch between attractors quickly .
invariant densities of @xmath62 , @xmath63 and @xmath64 are shown in panels b , c and d. in each case the bulk of the density is centred about the two underlying attractors . with white noise ( panel d )
there is no gap in the invariant density around @xmath332 due to randomness in both the @xmath47 and @xmath211-components of @xmath64 .
lastly , fig .
[ fig : invdensity_m002 ] illustrates stochastic dynamics with @xmath29 .
this figure corresponds to @xmath328 ( different to fig .
[ fig : detbifdiag ] ) and @xmath329 at which @xmath327 has an attracting @xmath38-cycle as well as the attracting fixed point @xmath32 ( [ eq : xlyl ] ) .
these are shown in panel a. the dynamics of @xmath63 and @xmath64 for @xmath319 is deterministic , hence @xmath32 is a fixed point of these maps .
given an initial point @xmath333 near the @xmath38-cycle , sample orbits of @xmath63 and @xmath64 eventually reach @xmath32 .
in contrast , @xmath62 has an invariant density concentrated about the two attractors , panel b. the part of the density with @xmath334 and @xmath335 corresponds to points of the orbit repeatedly following the left half - map of @xmath62 ( with @xmath24 and @xmath336 ) .
this paper concerns grazing bifurcations for which the associated dynamics is described by the nordmark map ( [ eq : n ] ) .
the potential influence of randomness and uncertainties on the dynamics of ( [ eq : n ] ) was described in @xcite by studying ( [ eq : n ] ) in the presence of additive white gaussian noise .
such a noise formulation is suitable if the nature of the randomness in the ode system is practically independent to the state of the system , such as if there is a random forcing term .    in this paper
we considered the alternate scenario that randomness is present in the switching condition of the ode system , and in @xmath7  the part of the vector field opposite to the tangential intersection of the grazing periodic orbit .
these cases are especially relevant for vibro - impacting systems for which impact events represent the primary source of uncertainty .
we derived three different stochastic versions of ( [ eq : n ] ) .
these are the maps @xmath62 ( [ eq : na ] ) , which corresponds to a noisy switching condition in the ode system , @xmath63 ( [ eq : nb ] ) , which corresponds to noise in @xmath7 with a large correlation time , and @xmath64 ( [ eq : nc ] ) , which corresponds to white noise in @xmath7 . in each case
the noise is nonlinear and non - additive .
this indicates that some diligence should be taken when formulating stochastic return maps for grazing bifurcations of piecewise - smooth systems .
the stochastic dynamics of @xmath62 , @xmath63 and @xmath64 differs in many ways to that of ( [ eq : n ] ) with additive noise , described in @xcite . for @xmath62 , @xmath63 and @xmath64 ,
dynamics prior to the grazing bifurcation is deterministic , and beyond the grazing bifurcation two - dimensional invariant densities are often skewed dramatically so that they appear almost one - dimensional .    near the grazing bifurcation , @xmath62 exhibits a large noise response relative to @xmath63 and @xmath64 .
this suggests that if experimental data of a physical system shows relatively high variability near a grazing bifurcation , then it is likely to be most appropriate to include randomness in the switching condition of a mathematical model .
invariant densities of @xmath62 near grazing are highly irregular due to the randomness in the switching condition . for @xmath63 and @xmath64 ,
the size of the noise response increases , for most part , with the distance ( in parameter space ) beyond the grazing bifurcation .
the maps @xmath63 and @xmath64 exhibit qualitatively similar invariant densities , which implies that the correlation time has little effect .
indeed the correlation time only influences the short - time dynamics of ( [ eq : fbc ] ) with @xmath11 , and the precise nature of these dynamics has a negligible effect on the invariant densities , which relate to long - time dynamics .
m.  oestreich , n.  hinrichs , k.  popp , and c.j .
budd . analytical and experimental investigation of an impact oscillator . in _ proceedings of the asme 16th biennal conf . on mech .
vibrations and noise .
_ , pages 111 , 1997 .
e.  knobloch and j.b .
effect of noise on discrete dynamical systems with multiple attractors . in m.f .
mcclintock and p.v.e .
moss , editors , _ noise in nonlinear dynamical systems .
theory of noise induced processes in special applications .
_ , volume  2 , pages 6586 . cambridge university press , new york , 1989 .