System and method for control of a subject's circadian cycle

• US 20050015122A1
• Filed: 06/03/2004
• Published: 01/20/2005
• Est. Priority Date: 06/03/2003
• Status: Abandoned Application

• Current Assignee: Pulsar Informatics, Inc.
• Original Assignee: Pulsar Informatics, Inc.
• Inventors: Huzmezan, Mihai, Mott, Christopher Grey, Van Wollen, Matthew, Mollicone, Daniel Joseph
• Application Number:
  US10/859,172
• Publication Number:
  US 20050015122A1
• Time in Patent Office: Days
• US Class Current: 607/88
• CPC Class Codes: A61M 2021/0044 by the sight sense
• CPC Class Codes: A61M 21/00 Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
• CPC Class Codes: A61M 2230/06 Heartbeat rate only
• CPC Class Codes: A61M 2230/205 partial oxygen pressure (P-O2)
• CPC Class Codes: A61M 2230/30 Blood pressure A61M2230/04 takes precedence
• CPC Class Codes: A61M 2230/42 Rate
• CPC Class Codes: A61M 2230/50 Temperature
• CPC Class Codes: A61N 5/0618 Psychological treatment
• CPC Class Codes: H05B 47/105 in response to determined parameters
• CPC Class Codes: H05B 47/11 by determining the brightness or colour temperature of ambient light
• CPC Class Codes: Y02B 20/40 Control techniques providing energy savings, e.g. smart controller or presence detection

Abstract

Aspects of the invention provide systems and methods for controllably adjusting the circadian pacemaker cycle of a subject using light (or other stimulus) through application of model-based predictive control techniques. This approach allows the use of closed-loop feedback to compensate for modeling errors, unknown initial conditions and disturbances. It also allows an optimal level of light (or other stimulus) to be generated based on minimization of a cost function. The cost function may incorporate a term associated with tracking errors and a term associated with the amount of light used. The tracking function may be minimized subject to one or more constraints which may include a minimum and maximum amount of light (or other stimulus).

Claims

1. A method for controllably adjusting a circadian cycle of a subject to track a reference trajectory, the method comprising:

providing at least one model representative of a response of a circadian state of the subject to a stimulus control signal; and

generating an optimal stimulus control signal using model predictive control based on the model,

wherein generating the optimal stimulus control signal comprises minimizing a cost function which comprises a cost term which is a function of a tracking error and wherein minimizing the cost function is subject to constraints which comprise minimum and maximum stimulus control signal levels.

2. A method according to claim 1 wherein the cost function comprises a cost term related to an amplitude of the stimulus control signal.

3. A method according to claim 2 wherein generating the optimal stimulus control signal is based, at least in part, on a circadian state feedback signal.

4. A method according to claim 3 wherein the stimulus control signal comprises a light control signal and wherein the method comprises applying the light control signal to determine an intensity of one or more light sources.

5. A method according to claim 4 wherein the at least one model comprises a response predictor subject model and wherein generating the optimal stimulus control signal comprises predicting a free response of the circadian state of the subject starting at a current time and extending out to a control horizon time using the response predictor subject model.

6. A method according to claim 5 wherein predicting the free response comprises assuming that a current value of the circadian state feedback signal represents an initial condition and that the light control signal remains constant between the current time and the control horizon time.

7. A method according to claim 6 wherein the response predictor subject model comprises a mathematical model defined by a plurality of differential equations.

8. A method according to claim 7 wherein the mathematical model is based on a Jewett-Kronauer model.

9. A method according to claim 8 wherein the mathematical model comprises a linearized version of the Jewett-Kronauer model.

10. A method according to claim 5 wherein the at least one model comprises a controller subject model and wherein generating the optimal stimulus control signal comprises determining an optimal series of control moves starting at the current time and extending out to the control horizon time using the controller subject model.

11. A method according to claim 10 wherein determining the optimal series of control moves is based, at least in part, on the free response of the circadian state of the subject, the reference trajectory and the circadian state feedback signal.

12. A method according to claim 11 wherein the controller subject model comprises a mathematical model defined by a plurality of differential equations.

13. A method according to claim 8 wherein the mathematical model is based on a Jewett-Kronauer model.

14. A method according to claim 13 wherein the mathematical model comprises a linearized version of the Jewett-Kronauer model.

15. A method according to claim 10 wherein the circadian state feedback signal is determined by modeling the response of the circadian state of the subject to the light control signal.

16. A method according to claim 15 modeling the response of the subject to the light control signal comprises using a mathematical model based on a Jewett-Kronauer model.

17. A method according to claim 16 wherein the mathematical model comprises a linearized version of the Jewett-Kronauer model.

18. A method according to claim 3 wherein the circadian state feedback signal is determined by modeling the response of the circadian state of the subject to the stimulus control signal.

19. A method according to claim 18 modeling the response of the subject to the stimulus control signal comprises using a mathematical model based on a Jewett-Kronauer model.

20. A method according to claim 19 wherein the mathematical model comprises a linearized version of the Jewett-Kronauer model.

21. A method according to claim 10 wherein the circadian state feedback signal is determined by estimating the circadian state of the subject based, at least in part, on one or more sensed parameters, the sensed parameters relating to the physiology of the subject.

22. A method according to claim 21 wherein the one or more sensed parameters include at least one of:

heart rate, core body temperature, respiration, endocrine function levels, physical activity levels, blood pressure, blood oxygen concentration and skin temperature.

23. A method according to claim 3 wherein the circadian state feedback signal is determined by estimating the circadian state of the subject based, at least in part, on one or more sensed parameters, the sensed parameters relating to the physiology of the subject.

24. A method according to claim 23 wherein the one or more sensed parameters include at least one of:

heart rate, core body temperature, respiration, endocrine function levels, physical activity levels, blood pressure, blood oxygen concentration and skin temperature.

25. A method according to claim 10 wherein predicting the free response is based, at least in part, on a light estimate signal representative of a light intensity experienced by the subject.

26. A method according to claim 25 wherein determining the optimal series of control moves is based, at least in part, on the light estimate signal.

27. A method according to claim 26 comprising estimating the light intensity experienced by the subject using a model to obtain the light estimate signal.

28. A method according to claim 26 comprising sensing the light intensity experienced by the subject using at least one light sensor to obtain the light estimate signal.

29. A method according to claim 13 wherein the stimulus control signal is provided in a domain of a driving input B of the Jewett-Kronauer model and wherein the method comprises converting the stimulus control signal from the domain of the driving input B to the light control signal in a light intensity domain I.

30. A method according to claim 29 comprising converting the minimum and maximum stimulus control signal levels from a light intensity domain I to a domain of a driving input B of the Jewett-Kronauer model.

31. A method according to claim 2 wherein the cost function comprises a first weighting factor associated with the tracking error cost term and a second weighting factor associated with the stimulus control signal cost term.

32. A method according to claim 31 wherein the first and second weighting factors vary over a duration of the circadian cycle of the subject.

33. A method for controllably adjusting a circadian cycle of a subject to track a reference trajectory, the method comprising:

obtaining input information, the input information comprising a signal representative of the reference trajectory, a feedback signal representative of a circadian state of the subject and one or more stimulus level constraints;

predicting a free response of the circadian state of the subject starting at a current time and extending out to a first future time using a response predictor subject model representative of a response of the circadian state of the subject to a stimulus;

determining an optimal series of control moves starting at the current time and extending out to second future time using a controller subject model representative of a response of the circadian state of the subject to a stimulus, the optimal series of control moves based, at least in part, on the free response of the circadian state of the subject, the reference trajectory signal and the current circadian state feedback signal; and

outputting a stimulus control signal to one or more stimulus sources, the stimulus control signal comprising a current one of the series of control moves; and

applying the stimulus control signal to one or more stimulus sources to determine an intensity thereof.

34. A method according to claim 33 wherein the response predictor subject model comprises a mathematical model based on a Jewett-Kronauer model.

35. A method according to claim 33 wherein the controller subject model comprises a mathematical model based on a linearized version of a Jewett-Kronauer model.

36. A system for controllably adjusting a circadian cycle of a subject to track a reference trajectory, the system comprising:

one or more inputs for receiving the reference trajectory and for receiving constraints which comprise minimum and maximum stimulus control signal levels;

a response predictor connected to receive a feedback signal representative of a circadian state of the subject, the response predictor configured to predict a free response of the circadian state of the subject starting at a current time and extending out to a first future time using a response predictor subject model representative of a response of the circadian state of the subject to a stimulus; and

a control sequence generator connected to receive the reference trajectory, the constraints, the feedback signal and the free response of the circadian state of the subject, the control sequence generator configured to determine an optimal series of control moves starting at the current time and extending out to a second future time using a controller subject model representative of a response of the circadian state of the subject to a stimulus, wherein the control sequence generator is connected to output a stimulus control signal to at least one stimulus source, the stimulus control signal comprising a current control move in the optimal series of control moves and the stimulus control signal determining an intensity of the at least one stimulus source.

37. A system according to claim 36 wherein the stimulus source comprises one or more light sources and the stimulus control signal comprises a light control signal representative of a light intensity of the one or more light sources.

38. A system according to claim 37 wherein the response predictor subject model comprises a mathematical model based on a Jewett-Kronauer model.

39. A system according to claim 37 wherein the controller subject model comprises a mathematical model based on a linearized version of a Jewett-Kronauer model.

40. A system according to claim 36 comprising one or more sensors for sensing one or more corresponding physiological parameters of the subject and a circadian state estimator which is connected to receive sensed values of the one or more physiological parameters from the one or more corresponding sensors and to estimate the feedback signal in response thereto.

41. A system according to claim 36 comprising a mathematical model which is connected to receive the stimulus control signal and which is configured to estimate the feedback signal in response thereto.

42. A system according to claim 37 comprising a converter for receiving the constraints which comprise minimum and maximum light control signal levels and converting the constraints from a light intensity domain I to a domain of a Jewett-Kronauer input signal B.

43. A system according to claim 42 comprising a converter for receiving the stimulus control signal in the domain of the Jewett-Kronauer input signal B and converting the stimulus control signal to the light control signal in the light intensity domain I.

44. A system for controllably adjusting a circadian cycle of a subject to track a reference trajectory, the system comprising:

one or more inputs for receiving the reference trajectory and for receiving constraints which comprise minimum and maximum stimulus control signal levels;

means for predicting a free response of a circadian state of the subject starting at a current time and extending out to a first future time based, at least in part, on a feedback signal representative of the circadian state of the subject;

means for determining an optimal series of control moves starting at the current time and extending out to a second future time based, at least in part, on the reference trajectory, the constraints, the feedback signal and the free response of the circadian state of the subject; and

means for applying a stimulus control signal to at least one stimulus source, the stimulus control signal comprising at least a portion of the optimal series of control moves and the stimulus control signal determining an intensity of the at least one stimulus source.

45. A method for altering a phase of a circadian cycle of a subject, the method comprising:

providing a controller subject model representing a response of a circadian state of the subject to a stimulus;

receiving a reference circadian trajectory, the reference circadian trajectory phase shifted from a current circadian cycle of the subject;

determining a series of stimulus control moves predicted by the controller subject model to result in the circadian cycle of the subject changing so as to track the reference circadian trajectory;

applying at least a portion of the series of stimulus control moves to one or more stimulus sources, the one or more stimulus sources providing stimulus which is received by the subject;

wherein determining the series of stimulus control moves comprises applying an optimization process using the controller subject model.

46. A method according to claim 45 comprising using the controller subject model to predict future circadian states of the subject and minimizing a cost function based, at least in part, on differences between the predicted future circadian states of the subject and the reference circadian trajectory.

47. A method according to claim 46 comprising modeling a current state of the circadian cycle of the subject by supplying the stimulus control moves as inputs to a model and supplying the modeled current state as a feedback input to the controller subject model.

48. A method according to claim 45 wherein the one or more stimulus sources comprise one or more light sources.

49. A method according to claim 45 wherein the controller subject model comprises a mathematical model based on a Jewett/Kronauer model.