Software simulation of the EEG

The literature contains many examples of digital procedures for the analytical treatment of electroencephalograms, but there is as yet no standard by which those techniques may be judged or compared. This paper proposes one method of generating an EEG, based on a computer program for Zetterberg's simulation. It is assumed that the statistical properties of an EEG may be represented by stationary processes having rational transfer functions and achieved by a system of software fillers and random number generators.The model represents neither the neurological mechanism response for generating the EEG, nor any particular type of EEG record; transient phenomena such as spikes, sharp waves and alpha bursts also are excluded. The basis of the program is a valid ‘partial’ statistical description of the EEG; that description is then used to produce a digital representation of a signal which if plotted sequentially, might or might not by chance resemble an EEG, that is unimportant. What is important is that the statistical properties of the series remain those of a real EEG; it is in this sense that the output is a simulation of the EEG. There is considerable flexibility in the form of the output, i.e. its alpha, beta and delta content, which may be selected by the user, the same selected parameters always producing the same statistical output. The filtered outputs from the random number sequences may be scaled to provide realistic power distributions in the accepted EEG frequency bands and then summed to create a digital output signal, the ‘stationary EEG’. It is suggested that the simulator might act as a test input to digital analytical techniques for the EEG, a simulator which would enable at least a substantial part of those techniques to be compared and assessed in an objective manner. The equations necessary to implement the model are given. The program has been run on a DEC1090 computer but is suitable for any microcomputer having more than 32 kBytes of memory; the execution time required to generate a 25 s simulated EEG is in the region of 15 s.

Flow through a plateau border of cellular foam

Experimentally measured average velocities through plateau borders of stationary cellular foam, when compared with those calculated with the assumption of rigid Plateau Border walls, show that the assumption of rigid walls severely underestimates the velocities. An analysis of the situation wherein plateau border walls have velocities, as decided by the surface viscosity of the system, is presented here. The plateau border is idealized as a pipe of equilateral triangular cross-section with vertices of the triangle having zero velocity. The pertinent form of Navier-Stoke's equations with inhomogeneous boundary conditions and its solution through a procedure of successive approximations is presented in dimensionless form. The solution reduces to the known solution of slow steady flow through a triangular pipe, when surface viscosity is infinite. Results indicate that the assumption of rigid plateau border walls is valid only when value of the inverse of dimensionless surface viscosity is less than 0.044. Beyond that the assumption severely underestimates the flow and the effect of nonrigidity of the wall must be considered.

Asymptotic behavior of a hierarchical system of learning automata

Learning automata arranged in a two-level hierarchy are considered. The automata operate in a stationary random environment and update their action probabilities according to the linear-reward- -penalty algorithm at each level. Unlike some hierarchical systems previously proposed, no information transfer exists from one level to another, and yet the hierarchy possesses good convergence properties. Using weak-convergence concepts it is shown that for large time and small values of parameters in the algorithm, the evolution of the optimal path probability can be represented by a diffusion whose parameters can be computed explicitly.

Investigation and application of methods for enumerating heterotrophs and Escherichia coli in the air within piggery sheds

Aims: To investigate methods for the recovery of airborne bacteria within pig sheds and to then use the appropriate methods to determine the levels of heterotrophs and Escherichia coli in the air within sheds. Methods and Results: AGI-30 impingers and a six-stage Andersen multi-stage sampler (AMS) were used for the collection of aerosols. Betaine and catalase were added to impinger collection fluid and the agar plates used in the AMS. Suitable media for enumerating E. coli with the Andersen sampler were also evaluated. The addition of betaine and catalase gave no marked increase in the recovery of heterotrophs or E. coli. No marked differences were found in the media used for enumeration of E. coli. The levels of heterotrophs and E. coli in three piggeries, during normal pig activities, were 2Æ2 · 105 and 21 CFU m)3 respectively. Conclusions: The failure of the additives to improve the recovery of either heterotrophs or E. coli suggests that these organisms are not stressed in the piggery environment. The levels of heterotrophs in the air inside the three Queensland piggeries investigated are consistent with those previously reported in other studies. Flushing with ponded effluent had no marked or consistent effect on the heterotroph or E. coli levels. Significance and Impact of the Study: Our work suggests that levels of airborne heterotrophs and E. coli inside pig sheds have no strong link with effluent flushing. It would seem unlikely that any single management activity within a pig shed has a dominant influence on levels of airborne heterotrophs and E. coli

Time dependence of effective slip on textured hydrophobic surfaces

In this paper, we present results on water flow past randomly textured hydrophobic surfaces with relatively large surface features of the order of 50 µm. Direct shear stress measurements are made on these surfaces in a channel configuration. The measurements indicate that the flow rates required to maintain a shear stress value vary substantially with water immersion time. At small times after filling the channel with water, the flow rates are up to 30% higher compared with the reference hydrophilic surface. With time, the flow rate gradually decreases and in a few hours reaches a value that is nearly the same as the hydrophilic case. Calculations of the effective slip lengths indicate that it varies from about 50 µm at small times to nearly zero or “no slip” after a few hours. Large effective slip lengths on such hydrophobic surfaces are known to be caused by trapped air pockets in the crevices of the surface. In order to understand the time dependent effective slip length, direct visualization of trapped air pockets is made in stationary water using the principle of total internal reflection of light at the water-air interface of the air pockets. These visualizations indicate that the number of bright spots corresponding to the air pockets decreases with time. This type of gradual disappearance of the trapped air pockets is possibly the reason for the decrease in effective slip length with time in the flow experiments. From the practical point of usage of such surfaces to reduce pressure drop, say, in microchannels, this time scale of the order of 1 h for the reduction in slip length would be very crucial. It would ultimately decide the time over which the surface can usefully provide pressure drop reductions. ©2009 American Institute of Physics

A comparison of plotless density estimators using Monte Carlo simulation on totally enumerated field data sets

Background: Plotless density estimators are those that are based on distance measures rather than counts per unit area (quadrats or plots) to estimate the density of some usually stationary event, e.g. burrow openings, damage to plant stems, etc. These estimators typically use distance measures between events and from random points to events to derive an estimate of density. The error and bias of these estimators for the various spatial patterns found in nature have been examined using simulated populations only. In this study we investigated eight plotless density estimators to determine which were robust across a wide range of data sets from fully mapped field sites. They covered a wide range of situations including animal damage to rice and corn, nest locations, active rodent burrows and distribution of plants. Monte Carlo simulations were applied to sample the data sets, and in all cases the error of the estimate (measured as relative root mean square error) was reduced with increasing sample size. The method of calculation and ease of use in the field were also used to judge the usefulness of the estimator. Estimators were evaluated in their original published forms, although the variable area transect (VAT) and ordered distance methods have been the subjects of optimization studies. Results: An estimator that was a compound of three basic distance estimators was found to be robust across all spatial patterns for sample sizes of 25 or greater. The same field methodology can be used either with the basic distance formula or the formula used with the Kendall-Moran estimator in which case a reduction in error may be gained for sample sizes less than 25, however, there is no improvement for larger sample sizes. The variable area transect (VAT) method performed moderately well, is easy to use in the field, and its calculations easy to undertake. Conclusion: Plotless density estimators can provide an estimate of density in situations where it would not be practical to layout a plot or quadrat and can in many cases reduce the workload in the field.

Comparison Of If Estimation Methods For Noise Robustness

Non-stationary signal modeling is a well addressed problem in the literature. Many methods have been proposed to model non-stationary signals such as time varying linear prediction and AM-FM modeling, the later being more popular. Estimation techniques to determine the AM-FM components of narrow-band signal, such as Hilbert transform, DESA1, DESA2, auditory processing approach, ZC approach, etc., are prevalent but their robustness to noise is not clearly addressed in the literature. This is critical for most practical applications, such as in communications. We explore the robustness of different AM-FM estimators in the presence of white Gaussian noise. Also, we have proposed three new methods for IF estimation based on non-uniform samples of the signal and multi-resolution analysis. Experimental results show that ZC based methods give better results than the popular methods such as DESA in clean condition as well as noisy condition.

Characteristics of Alfven surface waves along moving cylindrical plasma columns

The dispersive and stability characteristics of Alfven surface waves (ASW) along the boundary of the moving cylindrical plasma column, surrounded by a stationary medium embedded in a parallel magnetic field is studied. The nature of the symmetric and asymmetric modes on the interface parameters is also discussed.

Sleep-related crash characteristics: Implications for applying a fatigue definition to crash reports

Sleep-related (SR) crashes are an endemic problem the world over. However, police officers report difficulties in identifying sleepiness as a crash contributing factor. One approach to improving the sensitivity of SR crash identification is by applying a proxy definition post hoc to crash reports. To identify the prominent characteristics of SR crashes and highlight the influence of proxy definitions, ten years of Queensland (Australia) police reports of crashes occurring in ≥100 km/h speed zones were analysed. In Queensland, two approaches are routinely taken to identifying SR crashes. First, attending police officers identify crash causal factors; one possible option is ‘fatigue/fell asleep’. Second, a proxy definition is applied to all crash reports. Those meeting the definition are considered SR and added to the police-reported SR crashes. Of the 65,204 vehicle operators involved in crashes 3449 were police-reported as SR. Analyses of these data found that male drivers aged 16–24 years within the first two years of unsupervised driving were most likely to have a SR crash. Collision with a stationary object was more likely in SR than in not-SR crashes. Using the proxy definition 9739 (14.9%) crashes were classified as SR. Using the proxy definition removes the findings that SR crashes are more likely to involve males and be of high severity. Additionally, proxy defined SR crashes are no less likely at intersections than not-SR crashes. When interpreting crash data it is important to understand the implications of SR identification because strategies aimed at reducing the road toll are informed by such data. Without the correct interpretation, funding could be misdirected. Improving sleepiness identification should be a priority in terms of both improvement to police and proxy reporting.

Mobile assessment of on-road air pollution and its sources along the East-West Highway in Bhutan

Human exposures in transportation microenvironments are poorly represented by ambient stationary monitoring. A number of on-road studies using vehicle-based mobile monitoring have been conducted to address this. Most previous studies were conducted on urban roads in developed countries where the primary emission source was vehicles. Few studies have examined on-road pollution in developing countries in urban settings. Currently, no study has been conducted for roadways in rural environments where a substantial proportion of the population live. This study aimed to characterize on-road air quality on the East-West Highway (EWH) in Bhutan and identify its principal sources. We conducted six mobile measurements of PM10, particle number (PN) count and CO along the entire 570 km length of the EWH. We divided the EWH into five segments, R1-R5, taking the road length between two district towns as a single road segment. The pollutant concentrations varied widely along the different road segments, with the highest concentrations for R5 compared with other road segments (PM10 = 149 µg/m3, PN = 5.74 × 104 particles/cm-3, CO = 0.19 ppm), which is the final segment of the road to the capital. Apart from vehicle emissions, the dominant sources were road works, unpaved roads and roadside combustion activities. Overall, the highest contributions above the background levels were made by unpaved roads for PM10 (6 times background), and vehicle emissions for PN and CO (5 and 15 times background, respectively). Notwithstanding the differences in instrumentation used and particle size range measured, the current study showed lower PN concentrations compared with similar on-road studies. However, concentrations were still high enough that commuters, road maintenance workers and residents living along the EWH, were potentially exposed to elevated pollutant concentrations from combustion and non-combustion sources. Future studies should focus on assessing the dispersion patterns of roadway pollutants and defining the short- and long-term health impacts of exposure in Bhutan, as well as in other developing countries with similar characteristics.

Resonances, scattering theory, and rigged Hilbert spaces

The problem of decaying states and resonances is examined within the framework of scattering theory in a rigged Hilbert space formalism. The stationary free,''in,'' and ''out'' eigenvectors of formal scattering theory, which have a rigorous setting in rigged Hilbert space, are considered to be analytic functions of the energy eigenvalue. The value of these analytic functions at any point of regularity, real or complex, is an eigenvector with eigenvalue equal to the position of the point. The poles of the eigenvector families give origin to other eigenvectors of the Hamiltonian: the singularities of the ''out'' eigenvector family are the same as those of the continued S matrix, so that resonances are seen as eigenvectors of the Hamiltonian with eigenvalue equal to their location in the complex energy plane. Cauchy theorem then provides for expansions in terms of ''complete'' sets of eigenvectors with complex eigenvalues of the Hamiltonian. Applying such expansions to the survival amplitude of a decaying state, one finds that resonances give discrete contributions with purely exponential time behavior; the background is of course present, but explicitly separated. The resolvent of the Hamiltonian, restricted to the nuclear space appearing in the rigged Hilbert space, can be continued across the absolutely continuous spectrum; the singularities of the continuation are the same as those of the ''out'' eigenvectors. The free, ''in'' and ''out'' eigenvectors with complex eigenvalues and those corresponding to resonances can be approximated by physical vectors in the Hilbert space, as plane waves can. The need for having some further physical information in addition to the specification of the total Hamiltonian is apparent in the proposed framework. The formalism is applied to the Lee–Friedrichs model and to the scattering of a spinless particle by a local central potential. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Liikkumisen ja kokemusten yhteys virtuaaliympäristöissä

The aim of this study was to build a model and analyze how users move in a virtual environment and to explore the experiential dimensions connected with different ways of moving. Due to the lack of previous research on this subject, this was an explorative study. This study also aimed to identify different ways how users move in virtual environments and the background variables connected to them. It was hypothesized that fluent movement in virtual environments is connected to high presence, skill and challenge assessments. Test participants (n = 68) were mostly highly educated young adults. A virtual environment was built using a CAVE -type virtual reality interface. The task was to search for objects that do not belong into a normal house. The participants movement in the virtual house was recorded on a computer. Movement was modelled using a cluster analysis of information entropy based movement measurements, acceleration, amount of stops and time spent being stationary. The experiential dimensions were measured using the EVEQ -questionnaire. We were able to identify four different ways of moving in virtual environments. In respect of background variables, the four groups differed only in the amount of weekly computer usage. However, fluent movement in virtual environments was connected to a high sense of presence. Furthermore, participants who moved fluently in the environment assessed their skills as being high and regarded the use of virtual environment as challenging. The results indicate that different ways of moving affects how people experience virtual environments. Consequently the participants assessment of their skills and level of challenge have an impact on the affective evaluation of the situation at hand. Entropy measures have not been previously applied when studying movement, and in addition the role of movement on the experiential dimensions of virtual environments is an unexplored subject. The movement analysis method introduced here is applicable to other research problems. Finally, this study expands on our knowledge of the special characteristics connected with the experiential dimensions of virtual environments.

A mechanism for the west-north-west movement of monsoon depressions

The monsoon depressions intensify over the Bay of Bengal, move in a west-north-west (WNW) direction and dissipate over the Indian continent. No convincing physical explanation for their observed movement has so far been arrived at, but here, I suggest why the maximum precipitation occurs in the western sector of the depression and propose a feedback mechanism for the WNW movement of the depressions. We assume that a heat source is created over the Bay of Bengal due to organization of cumulus convection by the initial instability. In a linear sense, heating at this latitude (20° N), produces an atmospheric response mainly in the form of a stationary Rossby–gravity wave to the west of the heat source. The low-level vorticity (hence the frictional convergence) and the vertical velocity associated with the steady-state response is such that the maximum moisture convergence (and precipitation) is expected to occur in the WNW sector at a later time. Thus, the heat source moves to the WNW sector at a later time and the feedback continues resulting in the WNW movement of the depressions.

Nonstationary Random Critical Seismic Excitations

A method is presented to find nonstationary random seismic excitations with a constraint on mean square value such that the response variance of a given linear system is maximized. It is also possible to incorporate the dominant input frequency into the analysis. The excitation is taken to be the product of a deterministic enveloping function and a zero mean Gaussian stationary random process. The power spectral density function of this process is determined such that the response variance is maximized. Numerical results are presented for a single-degree system and an earth embankment modeled as shear beam.

Direct measurement of inner prosthesis loading during activities of daily living: Proposal for classification of functional outcome

Background The purpose of this presentation is to outline the relevance of the categorization of the load regime data to assess the functional output and usage of the prosthesis of lower limb amputees. The objectives are • To highlight the need for categorisation of activities of daily living • To present a categorization of load regime applied on residuum, • To present some descriptors of the four types of activity that could be detected, • To provide an example the results for a case. Methods The load applied on the osseointegrated fixation of one transfemoral amputee was recorded using a portable kinetic system for 5 hours. The load applied on the residuum was divided in four types of activities corresponding to inactivity, stationary loading, localized locomotion and directional locomotion as detailed in previously publications. Results The periods of directional locomotion, localized locomotion, and stationary loading occurred 44%, 34%, and 22% of recording time and each accounted for 51%, 38%, and 12% of the duration of the periods of activity, respectively. The absolute maximum force during directional locomotion, localized locomotion, and stationary loading was 19%, 15%, and 8% of the body weight on the anteroposterior axis, 20%, 19%, and 12% on the mediolateral axis, and 121%, 106%, and 99% on the long axis. A total of 2,783 gait cycles were recorded. Discussion Approximately 10% more gait cycles and 50% more of the total impulse than conventional analyses were identified. The proposed categorization and apparatus have the potential to complement conventional instruments, particularly for difficult cases.