Code: Difference between revisions
Line 8: | Line 8: | ||
== Determining moment of waking up == | == Determining moment of waking up == | ||
From our measurements we | 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 | 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 | 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 | 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 | 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. | ||
Revision as of 17:19, 19 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.
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