Code: Difference between revisions

From Control Systems Technology Group
Jump to navigation Jump to search
No edit summary
Line 19: Line 19:


For the complete matlab-script and more information about the output-form look at the [[Experiments]] page.
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
== Sound Control ==

Revision as of 10:44, 30 March 2016

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.

Retrieved from "https://dsdwiki.wtb.tue.nl/index.php?title=Code&oldid=25680"