Code

Back to main page: PRE2015_3_Groep4

To the input: Input

To the output: Output

to the experiments: Experiments

Determining moment of waking up

From our measurements we got information about respiration during your sleep. Out of this information, we get a global line of of your sleepphase throughout the night. Out of these line, we have to determine what's the best moment to wake someone between a given interval. The best moment to wake someone is just after the REM-sleep. So this is the first or second nREM-sleep stage. (see literature, add reference).

In our algorithm we have developed two methods which can determine at which moment the app should wake you. The first method is looking at the slope of the line. During REM-sleep, the respiration frequency is at its highest. So we can detect the end of the REM-phase with a drop in frequency. Another method is looking at the current frequency relative to the mean frequency of the night. With this method, we can detect a 'good' moment, i.e. when the frequency is below the mean of the past night. But the mean is probably too low, so we have to increase the value with a small factor of the variation.

However, during the REM-sleep, the respiration frequency varies a lot (see literature, add reference), so we need to be sure that a decline, or a value below the mean, is really the end of the REM-sleep and not just a variation during the REM-sleep. This can be done by looking for two or three consecutive decreasing frames. (two/three frames is around 8/12 minutes).

In the coming days, we will test both methods before we choose the best one. Another option is to combine the two methods. So we could call it a good moment of waking up when the line is declining and below the frequency threshold.

For the complete matlab-script and more information about the output-form look at the Experiments page.

Room Brightness Control

function that calculates the brightness of sunlight based on the date/time

 to-report light[day time]
   let psi   51.4408                ;; lat. ehv
   let Gsc   1367                   ;; w/m2 solar constant
   let LC    (15 * 1 - 5.4778) / 15 ;; timezone corrected longitude
   let c     0.21                   ;; fraction of extraterrestrial sunlight that reaches
                                     ; the surface (clear weather
 
   ;; most of the following formulas are from wikipedia, they need to be verified and sourced.
   let Ext   Gsc * (1 + 0.034 * cos(360 / 365 * (day - 2))) ;; day
   let theta 360 * (day - 1) / 365.242
   let d     23.45 * sin(360 * (day + 10) / 365)
   let eqt   0.258 * cos(theta) - 7.416 * sin(theta) - 3.648 * cos(2 * theta) - 9.228 * sin(2 * theta)
   let w     time + eqt / 60 - LC
   let as    asin(cos((w - 12) / 24 * 360) * cos(psi) * cos(d) + sin(psi) * sin(d))
   let Edv   Ext * exp(- c / sin(as)) * sin(90 - as) * (max list 0 as) / as
   report Edv * 93;;lm/w http://physics.ucsd.edu/~tmurphy/papers/lumens-per-watt.pdf
 end

Temperature Control

simple model to calculate the room temperature

 globals [T
          k
          h
          Tr
          Ti
          dT]
 
 to setup
   Clear-all
   set T  T0;; input starting temperature
   set k  -0.03
   set h  0.02
   set Tr 50
   set Ti 12
   set dT 0
   reset-ticks
 end
 
 to Go
     ifelse Heater;; Heater = true when the heater is on
     [set dT h * (Tr - T)] [set dT 0]
     set T T + k * (T - Ti) + dT
     tick
 end