Abstract:
Analog electrical and pressure data obtained by an EEG/pressure sensor matrix from the subject are sent to a signal processing module that derives digital EEG data from the electrical signal and pressure data without using sensors attached to a subject&#39;s head. The pressure data are used as a secondary signal to measure the physical orientation of the subject&#39;s head. The physical orientation is used to transform the derived EEG signal to a known coordinate axis (the orientation of the subject&#39;s head) to obtain useful, consistent, and accurate EEG data. The EEG may be used to determine a sleep state of a subject and to awaken the subject while in a particular state.

Description:
CROSS REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims the benefit of Provisional Application No. 61/795,191 filed Oct. 12, 2012. The 61/795,191 application is incorporated by reference herein, in its entirety, for all purposes. 
     
    
     BACKGROUND 
       [0002]    In recent times, sleep has become a research topic of great interest. Sleep is a universal animal behavior that has significant import in modern human development and productivity. As a result, there is great motivation to study sleep, its benefits and the benefits that may be derived from optimizing sleep on an individual basis. 
         [0003]    Electroencephalography (EEG) has been used pervasively in sleep research as a data collection tool due to its advantages in the field. Firstly, EEG creates reliable and reproducible quantitative results in a field that was, previously, largely qualitative. Secondly, it does so non-invasively by measuring electric signals at the scalp. Finally, during sleep, EEG signals have properties and behaviors that are largely invariant among varied populations. From EEG data, researchers have been able to accurately typify sleep stage cycles from EEG signals alone. 
         [0004]    Sleep is known to be an episodic process with various stages, starting from a light state of sleep to a deeper state and back to the light state. When awoken during a light state of sleep, post-waking fatigue and sleep inertia (e.g. the desire to return to a sleeping state while waking) are reduced in physiologically healthy patients. Sleep stages can be typified by brainwave patterns and brainwaves during light sleep are similar to brainwaves while awake. The theta waves (4-7 Hz) that occupy most of Stage 1 and 2 sleep are significantly more analogous to the alpha waves (8-12 Hz) that dominate our awake state than the delta waves (0.5-2 Hz) of deep sleep. Most neural state changes (ex. being asleep to being awake) require some period of time for neurochemical and metabolic elements to direct the change; going from a state of light sleep to an awake state is faster than going from a state of deep sleep to an awake state. Because the transition from light sleep to an awake state is rapid, the duration of sleep inertia and fatigue should be minimized, as light state brainwave patterns are the most analogous to those observed in the awake state. 
       SUMMARY 
       [0005]    Embodiments herein are directed to detecting the sleep state of a human subject and waking the subject when in a light sleep state in order to reduce morning fatigue and improve restfulness. 
         [0006]    In an embodiment, analog electrical and pressure data obtained by an EEG/pressure sensor matrix from the subject are sent to a signal processing module that derives digital EEG data from the electrical signal and pressure data without using sensors attached to a subject&#39;s head. The pressure data are used as a secondary signal to measure the physical orientation of the subject&#39;s head. The physical orientation is used to transform the derived EEG signal to a known coordinate axis (the orientation of the subject&#39;s head) to obtain useful, consistent, and accurate EEG data. 
         [0007]    In an embodiment, a system for acquiring an EEG signal suitable for use in electroencephalography from a subject includes a sensor matrix and processing control unit. 
         [0008]    The sensor matrix includes a plurality of EEG electrodes configured to receive electrical signals from the brain of the subject and a plurality of pressure sensors configured to receive pressure signals indicative of the placement of the head of the subject on the sensor matrix. The processing control unit is configured to perform operations that include obtaining EEG electrode data from the electrical signals received by the sensor matrix, determining a scalar EEG voltage for each EEG electrode, obtaining pressure data from the pressure signals received by the sensor matrix, determining a pressure normal unit vector for each EEG electrode from the pressure data, summing the product of the pressure normal unit vector and the scalar EEG voltage for each EEG electrode to obtain a relative mean EEG vector, and obtaining a scalar absolute EEG signal from a dot product of a unit orientation vector and the relative mean EEG vector. In an embodiment, the unit orientation vector is obtained by determining points of maximum pressure along orthogonal axes using the pressure data, deriving rotation and axial tilt of the subject&#39;s head from attributes of the points of maximum pressure, producing the unit orientation vector from the angles of the rotation and axial tilt. 
         [0009]    In an embodiment, the sensor matrix may be a flexible sheet that is a component of a pillow configured for receiving the head of the subject. 
         [0010]    In another embodiment, the sensor matrix and the processing control unit may be used in conjunction with a memory and an alerting device to awaken a subject. The memory receives the EEG electrode data derived from the electrical signals and receives pressure data derived from the pressure signal. The processing control unit determines a sleep state of the subject from the EEG vector and issues a command to power on the alerting device to awaken the subject while the subject is in or is approaching a light sleep state. In an embodiment, the operations performed by the processing control unit are performed by a computing device selected from the group consisting of a laptop computer, a desktop computer, a portable computing device, a tablet or a smart phone. 
     
    
     
       DESCRIPTION OF THE DRAWINGS 
         [0011]      FIG. 1  is a block diagram illustrating a high level view of a sleep tracking and waking optimization system according to an embodiment. 
           [0012]      FIG. 2A  is a block diagram illustrating a bio-amplifier according to an embodiment. 
           [0013]      FIG. 2B  is an exploded perspective view illustrating components of a pillow assembly incorporating a bio-amplifier according to an embodiment. 
           [0014]      FIGS. 3A ,  3 B and  3 C graphical representations illustrating pressure distributions from three orientations of the head of a subject according to an embodiment. 
           [0015]      FIG. 4A  is a graphical representation illustrating the cross products of the normal vector equation according to an embodiment. 
           [0016]      FIG. 4B  is a graphical representation illustrating a mean normal vector according to an embodiment. 
           [0017]      FIG. 4C  is a contour drawing illustrating unit normal vectors according to an embodiment. 
           [0018]      FIG. 5  is a block diagram illustrating a processing control unit according to an embodiment. 
           [0019]      FIG. 6A  is a flow diagram illustrating a process by which data from an EEG/pressure sensor matrix is processed according to an embodiment. 
           [0020]      FIG. 6B  is a flow diagram illustrating a process by which a sleep state is identified according to an embodiment. 
           [0021]      FIG. 6C  includes graphs illustrating various sleep state signatures. 
           [0022]      FIG. 7  is a system block diagram of a computing device suitable for use with various embodiments. 
       
    
    
     DETAILED DESCRIPTION 
       [0023]      FIG. 1  is a block diagram illustrating a high level view of a sleep tracking and waking optimization system according to an embodiment. 
         [0024]    A sleep tracking and waking optimization system  100  includes a bio-amplifier  102 , a processing control unit  110  and an alerting module  116 . The bio-amplifier  102  includes an EEG/pressure sensor matrix  104 , a raw signal processor  106  and an alerting device  108 . The alerting device may produce an audio signal or it may vibrate to awaken the subject. 
         [0025]    According to an embodiment, the EEG/pressure sensor matrix  104  is a flexible sheet in which both pressure and electrical signal sensors have been embedded. By way of illustration and not by way of limitation, the flexible sheet is embedded on a pillow. As will be explained in more detail below, the EEG/pressure sensor matrix  104  sends analog electrical and pressure data to a raw signal processing module  106  that digitizes the data and stores it in a memory  107 . 
         [0026]      FIG. 2A  is a block diagram illustrating a bio-amplifier according to an embodiment. 
         [0027]    In an embodiment, electrodes  202  are placed on a sheet on which the head of a subject under test rests. Contacts between the fabric near an electrode and the subject&#39;s head acquire electrical signals from the subject&#39;s head. The electrodes  202  feed electrical signals to a raw signal processor  106 . Pressure sensors  204  are also embedded in the flexible sheet and acquire pressure data that is used to orient the subject&#39;s head as more fully described below. 
         [0028]      FIG. 2B  is an exploded perspective view illustrating components of a pillow assembly incorporating a bio-amplifier according to an embodiment. In an embodiment, components of a bio-amplifier  102  are integrated into with a foam portion  118  to form a pillow assembly  120 . The pillow assembly  120  may be used in various embodiments described below. 
         [0029]    Referring again to  FIG. 1 , the digitized electrical and pressure data are communicated from the memory  107  to the processing control unit  110 . By way of illustration and not by way of limitation, the bio-amplifier may communicate the digitized electrical and pressure data to the processing control unit  110  using a wired or wireless link. 
         [0030]    The EEG signal is not acquired from an array of electrodes that are affixed to a subject&#39;s head. Rather, the motion of the subject will constantly change the electrode measurement and the reference state of the EEG signal. 
         [0031]      FIGS. 3A ,  3 B and  3 C illustrate pressure distributions resulting from three orientations of the head of a subject on pillow incorporating a bio-amplifier  102 . In  FIG. 3A , the subject&#39;s head is flat on the pillow. In this orientation, the pressure is distributed evenly. The pressure gradient is described by a parabolic function orthogonal to the axis of the head defined by the middle of the forehead through the chin.  FIGS. 3B and 3C  illustrate distribution functions that are not parabolic. 
         [0032]    Additionally, the EEG signal is a three-dimensional time-varying vector that measures a mean electrical potential of the brain. A two-dimensional electrical map is not sufficient to map to a three dimensional space. 
         [0033]    In an embodiment, pressure is used as an additional variable to create “depth” on the two-dimensional electrode map. The three-dimensional surface defined by pressure is used to solve for the current orientation of the subject&#39;s head, as well as for generating a surface to give electrical signals a three-dimensional orientation. 
         [0034]    In an embodiment, a pressure surface is defined by the following vector function (where x and y are indices and X and Y are real space): 
         [0000]        f ( x,y )=&lt; X,Y,P&gt;={right arrow over (P)}   x,y    
         [0035]    The following equation is used to obtain the unit normal vector for each electrode (which is given index e). 
         [0000]    
       
         
           
             
               
                 
                   
                     Normal 
                      
                     
                         
                     
                      
                     Vector 
                   
                   = 
                     
                    
                   
                     
                       N 
                       → 
                     
                     e 
                   
                 
               
               
                 
                     
                    
                    
                 
               
             
             
               
                 
                   = 
                     
                    
                   
                     
                       
                         ( 
                         
                           
                             
                               P 
                               → 
                             
                             
                               
                                 x 
                                 + 
                                 1 
                               
                               , 
                               y 
                             
                           
                           - 
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               y 
                             
                           
                         
                         ) 
                       
                       × 
                       
                         ( 
                         
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               
                                 y 
                                 + 
                                 1 
                               
                             
                           
                           - 
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               y 
                             
                           
                         
                         ) 
                       
                     
                     + 
                     
                       
                         ( 
                         
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               y 
                             
                           
                           - 
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               
                                 y 
                                 + 
                                 1 
                               
                             
                           
                         
                         ) 
                       
                       × 
                     
                   
                 
               
               
                 
                     
                    
                   □ 
                 
               
             
             
               
                 
                     
                    
                   
                     
                       ( 
                       
                         
                           
                             P 
                             → 
                           
                           
                             
                               x 
                               + 
                               1 
                             
                             , 
                             
                               y 
                               + 
                               1 
                             
                           
                         
                         - 
                         
                           
                             P 
                             → 
                           
                           
                             x 
                             , 
                             
                               y 
                               + 
                               1 
                             
                           
                         
                       
                       ) 
                     
                     + 
                     
                       
                         ( 
                         
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               
                                 y 
                                 + 
                                 1 
                               
                             
                           
                           - 
                           
                             
                               P 
                               → 
                             
                             
                               
                                 x 
                                 + 
                                 1 
                               
                               , 
                               y 
                             
                           
                         
                         ) 
                       
                       × 
                     
                   
                 
               
               
                 
                     
                    
                   □ 
                 
               
             
             
               
                 
                     
                    
                   
                     
                       ( 
                       
                         
                           
                             P 
                             → 
                           
                           
                             x 
                             , 
                             y 
                           
                         
                         - 
                         
                           
                             P 
                             → 
                           
                           
                             
                               x 
                               + 
                               1 
                             
                             , 
                             y 
                           
                         
                       
                       ) 
                     
                     + 
                     
                       
                         ( 
                         
                           
                             
                               P 
                               → 
                             
                             
                               x 
                               , 
                               
                                 y 
                                 + 
                                 1 
                               
                             
                           
                           - 
                           
                             
                               P 
                               → 
                             
                             
                               
                                 x 
                                 + 
                                 1 
                               
                               , 
                               
                                 y 
                                 + 
                                 1 
                               
                             
                           
                         
                         ) 
                       
                       × 
                       
                         
 
                       
                        
                       
                         ( 
                         
                           
                             
                               P 
                               → 
                             
                             
                               
                                 x 
                                 + 
                                 1 
                               
                               , 
                               y 
                             
                           
                           - 
                           
                             
                               P 
                               → 
                             
                             
                               
                                 x 
                                 + 
                                 1 
                               
                               , 
                               
                                 y 
                                 + 
                                 1 
                               
                             
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                     
                    
                   □ 
                 
               
             
           
         
       
     
         [0036]      FIG. 4A  is a graphical representation illustrating the cross products of the normal vector equation according to an embodiment. 
         [0037]    The mean normal vector for the electrode is determined from the following equation: 
         [0000]    
       
         
           
             
               
                 n 
                 → 
               
               e 
             
             = 
             
               
                 
                   N 
                   → 
                 
                 e 
               
               
                  
                 
                   
                     N 
                     → 
                   
                   e 
                 
                  
               
             
           
         
       
     
         [0038]      FIG. 4B  is a graphical representation illustrating the mean normal vector  410  according to an embodiment. 
         [0039]      FIG. 4C  is a contour drawing illustrating unit normal vectors  420 . Each electrode is assigned a scalar voltage as follows: 
         [0000]        V ( e )= V   e    
         [0040]    The relative mean EEG vector is determined as follows: 
         [0000]    
       
         
           
             
               Relative 
                
               
                   
               
                
               Mean 
                
               
                   
               
                
               EEG 
             
             = 
             
               
                 
                   ∑ 
                   e 
                   
                     n 
                     e 
                   
                 
                  
                 
                     
                 
                  
                 
                   
                     V 
                     e 
                   
                    
                   
                     
                       n 
                       → 
                     
                     e 
                   
                 
               
               
                 n 
                 e 
               
             
           
         
       
     
         [0041]    Where n e  is the number of electrodes 
         [0042]    When a subject lies on the EEG/pressure sensor matrix  104 , a two-dimensional contour map of the pressure on the device is generated. This contour may be represented as the following function: 
         [0000]        F ( x,y )= P    
         [0043]    Exemplary contours are illustrated in  FIGS. 3A ,  3 B and  3 C. 
         [0044]    The directional derivative is obtained for the entire vector map in the y direction. All values of y in which the directional derivative passes the 1 st  derivative test for maxima as points y max  are obtained, where j is the unit vector in the y direction: 
         [0000]    
       
         
           
             
               
                 ▽ 
                 
                   j 
                   → 
                 
               
                
               
                 F 
                  
                 
                   ( 
                   
                     x 
                     , 
                     y 
                   
                   ) 
                 
               
             
             = 
             
               
                 ∂ 
                 P 
               
               
                 ∂ 
                 y 
               
             
           
         
       
     
         [0045]    Solve for y max  defined as: 
         [0000]    
       
         
           
             
               
                 ▽ 
                 
                   j 
                   → 
                 
               
                
               
                 F 
                  
                 
                   ( 
                   
                     x 
                     , 
                     
                       y 
                       max 
                     
                   
                   ) 
                 
               
             
             = 
             0 
           
         
       
       
         
           and 
         
       
       
         
           
             
               
                 lim 
                 
                   y 
                   → 
                   
                     y 
                     
                       max 
                       - 
                     
                   
                 
               
                
               
                 
                   ▽ 
                   
                     j 
                     → 
                   
                 
                  
                 
                   F 
                    
                   
                     ( 
                     
                       x 
                       , 
                       
                         y 
                         max 
                       
                     
                     ) 
                   
                 
               
             
             = 
             
               + 
               0 
             
           
         
       
       
         
           and 
         
       
       
         
           
             
               
                 lim 
                 
                   y 
                   → 
                   
                     y 
                     
                       max 
                       + 
                     
                   
                 
               
                
               
                 
                   ▽ 
                   
                     j 
                     → 
                   
                 
                  
                 
                   F 
                    
                   
                     ( 
                     
                       x 
                       , 
                       
                         y 
                         max 
                       
                     
                     ) 
                   
                 
               
             
             = 
             
               - 
               0 
             
           
         
       
     
         [0046]    There exists a function g such that: 
         [0000]        g ( x )= y   max    
         [0047]    The function g(x) is treated as linear so the rotation θ of the head about the z-axis is defined as: 
         [0000]    
       
         
           
             θ 
             = 
             
               
                 tan 
                 
                   - 
                   1 
                 
               
                
               
                 ( 
                 
                   
                     〈 
                     
                       
                         ∂ 
                         g 
                       
                       
                         ∂ 
                         x 
                       
                     
                     〉 
                   
                   
                     for 
                      
                     
                         
                     
                      
                     all 
                      
                     
                         
                     
                      
                     
                       y 
                       max 
                     
                   
                 
                 ) 
               
             
           
         
       
     
         [0048]    where 
         [0000]    
       
         
           
             〈 
             
               
                 ∂ 
                 g 
               
               
                 ∂ 
                 x 
               
             
             〉 
           
         
       
     
         [0000]    represents the average over a set of point. 
         [0049]    All the functions of pressure orthogonal to the linear function g(x) are obtained by reorienting the pressure map and holding the new variables constant. 
         [0000]    
       
         
           
             
               f 
                
               
                 ( 
                 
                   a 
                   , 
                   b 
                 
                 ) 
               
             
             = 
             P 
           
         
       
       
         
           where 
         
       
       
         
           
             a 
             = 
             
               
                 
                   〈 
                   
                     
                       ∂ 
                       g 
                     
                     
                       ∂ 
                       x 
                     
                   
                   〉 
                 
                 
                   for 
                    
                   
                       
                   
                    
                   all 
                    
                   
                       
                   
                    
                   
                     y 
                     max 
                   
                 
               
               * 
               x 
             
           
         
       
       
         
           and 
         
       
       
         
           
             b 
             = 
             
               
                 
                   - 
                   1 
                 
                 
                   
                     〈 
                     
                       
                         ∂ 
                         g 
                       
                       
                         ∂ 
                         x 
                       
                     
                     〉 
                   
                   
                     for 
                      
                     
                         
                     
                      
                     all 
                      
                     
                         
                     
                      
                     
                       y 
                       max 
                     
                   
                 
               
               * 
               y 
             
           
         
       
     
         [0050]    If “a” is held constant to a 0,  a function h of “b” that maps pressure may be obtained: 
         [0000]        f ( a   o   ,b )= h ( b )= P    
         [0051]    When the head is laying flat on the pillow, regardless of its rotation about the z-axis, it spreads pressure evenly across the pillow. This state may be modeled as a parabola orthogonal to the axis of the head. However, if the head is tilted in either direction, the pressure gradient is no longer parabolic.  FIGS. 3B and 3C  illustrate distribution functions that are not parabolic. 
         [0052]    The difference in pressure growth on both sides of the head can be represented as the 2 nd  derivative of both sides of the pressure map. 
         [0053]    The following equation takes the difference of both sides of the pressure curve formed by h(b) and normalizes them to generate a “tilt index”: 
         [0000]    
       
         
           
             
               tilt 
                
               
                   
               
                
               index 
             
             = 
             
               
                 
                   
                     〈 
                     
                       
                         ∂ 
                         h 
                       
                       
                         ∂ 
                         b 
                       
                     
                     〉 
                   
                   + 
                 
                 - 
                 
                   
                     〈 
                     
                       
                         ∂ 
                         h 
                       
                       
                         ∂ 
                         b 
                       
                     
                     〉 
                   
                   - 
                 
               
               
                 
                   
                     〈 
                     
                       
                         ∂ 
                         h 
                       
                       
                         ∂ 
                         b 
                       
                     
                     〉 
                   
                   + 
                 
                 + 
                 
                   
                     〈 
                     
                       
                         ∂ 
                         h 
                       
                       
                         ∂ 
                         b 
                       
                     
                     〉 
                   
                   - 
                 
               
             
           
         
       
     
         [0054]    where + and − indicate which side of the parabola the average is taken from. 
         [0055]    Table 1 is an illustration of a look-up table that maps measured tilt index values to tilt angles of the head of a subject, based on a sample population deemed representative of typical subjects. 
         [0000]    
       
         
               
               
               
             
               
               
               
             
           
               
                   
                 TABLE 1 
               
               
                   
                   
               
               
                   
                 Tilt Index 
                 Φ (tilt angle away from the z-axis) 
               
               
                   
                   
               
             
             
               
                   
               
             
          
           
               
                   
                 −1.5 
                 −45 
               
               
                   
                 −1 
                 −30 
               
               
                   
                 −.5 
                 −15 
               
               
                   
                 0 
                 0 
               
               
                   
                 .5 
                 15 
               
               
                   
                 1 
                 30 
               
               
                   
                 1.5 
                 45 
               
               
                   
                   
               
             
          
         
       
     
         [0056]    In an embodiment, because the function to map tilt index to real head tilt varies significantly from person to person due to head weight, size and shape, a look-up table for a particular subject may be generated by calibrating the pillow for a particular subject. The look up table will then map the tilt index for a given subject to the subject&#39;s actual head. See Table 2 for a sample subject&#39;s pillow calibration. 
         [0000]    
       
         
               
               
               
               
             
               
               
               
               
             
           
               
                   
                 TABLE 2 
               
               
                   
                   
               
               
                   
                 Calibration 
                 Φ 
                 Tilt Index 
               
               
                   
                   
               
             
             
               
                   
               
             
          
           
               
                   
                 Head facing up 
                 0 
                 0 
               
               
                   
                 Head turned 45° to the left 
                 −45 
                 −1.56 
               
               
                   
                 Head turned 45° to the right 
                 45 
                 1.52 
               
               
                   
                   
               
             
          
         
       
     
         [0057]    After the tilt and rotation are determined, a unit orientation vector in spherical coordinates with rotation θ, tilt Φ and radius of 1 is generated. 
         [0058]    This vector is then converted into Cartesian space and then dotted with the relative mean EEG vector to give the absolute scalar EEG signal: 
         [0000]      Absolute Scalar EEG=           r  sin Φ cos θ,  r  sin Φ sin θ,  r  cos Φ         ·Relative Mean EEG
 
         [0059]    This signal is the magnitude projection of the EEG vector onto its correct measurement orientation. This scalar signal is then processed to solve for sleep state. 
         [0060]    Referring again to  FIG. 1 , in an embodiment, the processing control unit  110  includes an EEG vector analysis module  112 , a lookup table  114  and an alerting module  116 . The processing control unit  110  receives the electrical and pressure data from the raw signal processor and uses a lookup table  114  to derive a correctly oriented EEG mean vector. This vector is computed at fixed intervals to produce an EEG mean vector timeseries. 
         [0061]    By way of illustration and not by way of limitation, the processing control unit  110  is incorporated into an alarm clock device, which device also accepts input from the user. In another embodiment, the functions of the processing control unit may be performed by a computing device, such as a computing device illustrated in  FIG. 6 , executing software instructions. By way of illustration and not by way of limitation, the computing device may be a laptop computer, a desktop computer, or a portable computing device such as a tablet or a smart phone. 
         [0062]    The absolute scalar EEG timeseries is provided to the alerting module  116 . The alerting module  116  processes the absolute scalar EEG timeseries to determine a sleep state of the subject under test. The alerting module  116  uses the sleep state to issue an alarm signal to an alerting device  108  in the bio-amplifier  102 . As further illustrated in  FIG. 5  and discussed below, the alerting module may activates an alarm device that is incorporated into a processing control unit or an alarm device that is external to the processing control unit. 
         [0063]    In an embodiment, the sleep state of the subject may be mapped over time to determine a sleep cycle pattern. This pattern may be used to project a future sleep state of the subject. 
         [0064]      FIG. 5  is a block diagram illustrating a processing control unit with added functionality according to an embodiment. 
         [0065]    As described above and illustrated in  FIG. 1 , the processing control unit  110  includes an EEG vector analysis module  112 , a lookup table  114  and an alerting module  116 .  FIG. 5  illustrates a processing control unit  502  that includes additional functional elements that provide user-controllable features. A wireless interface  504  provides communications with the bio-amplifier  102 . The wireless interface may be any short-range protocol including, for example, Bluetooth, Wi-Fi, and Zigbee. The processing control unit  502  may receive signal inputs (filtered and amplified EEG and pressure signals) from the memory  107  of the bio-amplifier  102  and may send alarm signals to the alarm  108  and/or an alarm device  517  incorporated into the processing control unit  502  and/or an alarm device  518  that is remote from the processing control unit  502  using the wireless interface  504 . 
         [0066]    The processing control unit  502  may also include a user interface  506  that receives user selected alarm and timing data. These data are provided to the timer circuit  508  to set a time for display on an LED display device  512 . The data may also be used to establish a preferred wake-up time that is provided to the alerting module  116 . 
         [0067]    In an embodiment, the EEG vector analysis module  112  uses the pressure and EEG data to determine the sleep state of a subject and provides this information to the alerting module  116 . The alerting module assesses the trends of the sleep state of the subject against the preferred wake-up time to determine when to activate the alarm  108  component of the bio-amplifier  102 . Thus, the activation of the alarm  108  may occur before the preferred wake-up time if needed to optimize the wake time of the subject relative to the subject&#39;s sleep state and/or projected sleep state. 
         [0068]      FIG. 6A  is a flow diagram illustrating a process by which data from an EEG/pressure sensor matrix is processed according to an embodiment. 
         [0069]    At Block  602 , EEG electrode signals are acquired from the memory  107  of the bio-amplifier  102  (see,  FIG. 1 ). At Block  604  the signals are band passe filtered. At Block  606  a scalar EEG value is obtained. At Block  612  pressure data are acquired from the memory  107  of the bio-amplifier  102  (see,  FIG. 1 ). Vector operations are performed on the pressure data (Block  614 ) and a pressure normal unit vector is obtained (Block  616 ). 
         [0070]    Operations to determine g (x)=y max  are performed a maximum pressure operation (Block  620 ). The rotation angle of the head  0  is determined (Block  622 ). The tilt index is determined (Block  624 ). The tilt index is used to determine the tilt angle Φ (Block  626 ). The pressure normal unit vector is generated (Block  628 ). 
         [0071]    At Block  632 , the scalar EEG and pressure normal unit vector are multiplied to obtain a relative mean EEG vector. 
         [0072]    At Block  634 , an absolute scalar EEG vector is obtained from the dot product of the relative mean EEG vector and a unit orientation vector. 
         [0073]    At Block  636 , the sleep state of the subject is determined from the absolute scalar EEG vector. 
         [0074]    At Block  638 , the sleep state is provided to an alerting module. 
         [0075]      FIG. 6B  is a flow diagram illustrating a process by which a sleep state is identified according to an embodiment. 
         [0076]    The EEG signal data from the absolute scalar EEG vector are acquired (Block  632 ) and converted to a frequency domain using a Fourier transform (Block  634 ). The amplitude and frequency of the signal is determined (Block  638 ). Each stage of sleep has its own unique signature of frequencies and amplitudes that identify it. The amplitude and frequency signature of the signal that are obtained by Fourier transforming include real and imaginary components which differ based on the phase shift of the incoming signal. 
         [0077]      FIG. 6C  includes graphs illustrating various sleep state signatures. Attributes of the sleep state signatures may be used in at Block  640  to identify a current sleep state of a subject. 
         [0078]    Referring again to  FIG. 6B , a determination is made at Block  640  whether the data are sufficient to identify the sleep stage of the subject. If the sleep stage cannot be identified from the data, that is, if the result of Block  640  is “NO,” the process waits for an updated absolute scalar EEG. In an embodiment, additional data are acquired from the from the memory  107  of the bio-amplifier  102  (see,  FIG. 1 ). In an embodiment, several signal acquisitions may be made and summed to improve the signal to noise ratio. In other embodiment, EEG and pressure data may be recorded for an extended amount of time to mitigate the effects of disturbances to the device and noise. 
         [0079]    If the sleep stage can be identified from the data, that is, if the result of Block  640  is “YES,” at Block  634  the sleep state is provided the alerting module  116  of the processing control unit  110  (see,  FIG. 1 ). The process continues with additional data being acquired from the memory  107  of the bio-amplifier  102  (see,  FIG. 1 ). 
         [0080]    Once an evaluation yields a conclusive result, the algorithm may inform the alerting module of the subject&#39;s current sleep stage. 
         [0081]    In an embodiment, the process illustrated in  FIGS. 6A  occurs several times each second, and  6 B occurs several times each minute, allowing the processing control unit  110  to maintain a strict definition of each stage of sleep. The processing control unit  110  will only alert the alarm component that a change in sleep stage has occurred when it has conclusively made such a determination. Thus, even though the processing control unit  110  may not be able to determine the individual&#39;s stage of sleep upon each signal acquisition, it will still report changes in sleep stage with high accuracy. 
         [0082]    In an embodiment, when the subject is in a deep sleep and is not projected to be in a light sleep state at a time chosen by the subject to be awakened, the alerting device may be signaled to issue an alarm of a type and/or magnitude that will not awaken the subject but is likely to raise the sleep level of the subject from heavy toward light. The non-awakening alarms may be repeated with varying intensity in order to encourage the subject into a lighter level of sleep. At a point when the subject enters a light level of sleep or the time for awakening have been reached, the alerting device may be signaled to issue an alarm of a type and/or magnitude that will awaken the subject. 
         [0083]      FIG. 7  is a system block diagram of a computing device useful to perform functions of a processing control unit. While the computing device  700  is illustrated as a laptop computer, a computing device providing the functional capabilities of the computer device  700  may be implemented as a tablet, a smartphone, a desktop computer or a handheld computer. A typical computing device  700  may include a processor  701  coupled to an electronic display  704 , a speaker  706 , a volatile memory  702  and to a nonvolatile memory, such as a disk drive  922 . When implemented as a laptop computer or desktop computer, the computing device  700  may also include a floppy disc drive, compact disc (CD) or DVD disc drive coupled to the processor  701 . The computing device  700  may include an antenna  710 , a multimedia receiver  712 , a transceiver  718  and/or communications circuitry  920  coupled to the processor  701  for sending and receiving electromagnetic radiation, connecting to a wireless data link, and receiving data from the bio-amplifier  102 . Additionally, the computing device  700  may include network access ports  724  coupled to the processor  701  for establishing data connections with a network (e.g., LAN coupled to a service provider network, etc.). A laptop computer or desktop computer  700  typically also includes a keyboard  714  and a mouse pad  716  for receiving user inputs. 
         [0084]    It will be understood by those skilled in the art that the present invention may be embodied in other specific forms without departing from the scope of the invention disclosed and that the examples and embodiments described herein are in all respects illustrative and not restrictive. Those skilled in the art of the present invention will recognize that other embodiments using the concepts described herein are also possible. Further, any reference to claim elements in the singular, for example, using the articles “a,” “an,” or “the,” is not to be construed as limiting the element to the singular.