Consider the following modified logistic equation with a threshold dV dt (1) where t represents time,...
Fantastic news! We've Found the answer you've been seeking!
Question:
Transcribed Image Text:
Consider the following modified logistic equation with a threshold dV dt (1) where t represents time, V is the dependent variable (e.g., voltage) and the function F(V) is given by = F(V) == -r V V V (₁ - 1) (1 - 1) 1 T K F(V) (2) The parameters r, K,T and I are constants, r, K,T > 0 and T < K. The parameter r represents the rate constant (inverse of the time constant for the unbiased case I app 0). The units depend on the problem at hand. If the model is intended for the dynamics of a neuron, then the units of t are milliseconds (ms) and the units of V are millivolts (mV). This equation has many applications. In particular, it can be used to describe the transi- tion of a system from a resting state (Vrest, lower equilibrium Veq) to and activated state (Vact, higher equilibrium Veg), when they exist, through a threshold (Vth, intermediate equilibrium Veg). Vact is also referred to as the saturation value V sat. In the context of neuronal systems, eq. (1) can be used as a (very simple) model of ac- tion potential generation (but not termination) where the occurrence of an action potential + I. = is determined by the transition from Vrest to Vact. In this case, the parameter I is interpreted as the applied current (Iapp) to the neuron (an external DC current that is controlled by the experimenter/modeler). The description of the termination of an action potential is not de- scribed by this model, and requires additional mechanisms and additional model complexity (e.g., a system of two or more differential equations). The existence and values of the equilibria (Vrest, Vth and Vsat) depend on the model parameters (r, T, K and I). For I = 0, eq (1) has three equilibria (Veg): Vrest = 0, Vth = T and Vsat = K. Because of this, T and K are referred as the threshold and carrying capacity, respectively, in the literature. For other values of I, the activity attributes (Vrest, Vth and Vsat are not longer explicitly given by a single parameter (T, K) and the attributes may even cease to exist for certain ranges of values of I (all other parameters fixed). 1. Write a code to solve numerically the ODE (1)-(2) (or adapt the template code provided in the course website). The simulation output for each set of parameter values must be (a) A graph of the solution V(t). (b) The equilibrium value(s) Veq = limt→∞ V(t) (c) A graph of F as a function of V. • Each simulation should be run long enough (large enough value of t) so that V(t) reaches values close enough to Veq, but not too long so the changes in V(t) are clearly shown. . Plot the two graphs as two panels in the same graph. . Use correct labels for the axes. • Use large enough fonts (suggested: "fontsize" - 24) 2. Consider the following parameter values: r = 1, T = 0.25, K = 1. Perform simulations as described above for V(0) = 0.01 and three values of I: I = 0, I = 0.05, I = 0.1. 3. (Undergraduate level) Investigate the behavior of the model in terms of the model pa- rameters (e.g., plot graphs of the dependence of the attributes of activity as a function of the model parameters). Consider the following modified logistic equation with a threshold dV dt (1) where t represents time, V is the dependent variable (e.g., voltage) and the function F(V) is given by = F(V) == -r V V V (₁ - 1) (1 - 1) 1 T K F(V) (2) The parameters r, K,T and I are constants, r, K,T > 0 and T < K. The parameter r represents the rate constant (inverse of the time constant for the unbiased case I app 0). The units depend on the problem at hand. If the model is intended for the dynamics of a neuron, then the units of t are milliseconds (ms) and the units of V are millivolts (mV). This equation has many applications. In particular, it can be used to describe the transi- tion of a system from a resting state (Vrest, lower equilibrium Veq) to and activated state (Vact, higher equilibrium Veg), when they exist, through a threshold (Vth, intermediate equilibrium Veg). Vact is also referred to as the saturation value V sat. In the context of neuronal systems, eq. (1) can be used as a (very simple) model of ac- tion potential generation (but not termination) where the occurrence of an action potential + I. = is determined by the transition from Vrest to Vact. In this case, the parameter I is interpreted as the applied current (Iapp) to the neuron (an external DC current that is controlled by the experimenter/modeler). The description of the termination of an action potential is not de- scribed by this model, and requires additional mechanisms and additional model complexity (e.g., a system of two or more differential equations). The existence and values of the equilibria (Vrest, Vth and Vsat) depend on the model parameters (r, T, K and I). For I = 0, eq (1) has three equilibria (Veg): Vrest = 0, Vth = T and Vsat = K. Because of this, T and K are referred as the threshold and carrying capacity, respectively, in the literature. For other values of I, the activity attributes (Vrest, Vth and Vsat are not longer explicitly given by a single parameter (T, K) and the attributes may even cease to exist for certain ranges of values of I (all other parameters fixed). 1. Write a code to solve numerically the ODE (1)-(2) (or adapt the template code provided in the course website). The simulation output for each set of parameter values must be (a) A graph of the solution V(t). (b) The equilibrium value(s) Veq = limt→∞ V(t) (c) A graph of F as a function of V. • Each simulation should be run long enough (large enough value of t) so that V(t) reaches values close enough to Veq, but not too long so the changes in V(t) are clearly shown. . Plot the two graphs as two panels in the same graph. . Use correct labels for the axes. • Use large enough fonts (suggested: "fontsize" - 24) 2. Consider the following parameter values: r = 1, T = 0.25, K = 1. Perform simulations as described above for V(0) = 0.01 and three values of I: I = 0, I = 0.05, I = 0.1. 3. (Undergraduate level) Investigate the behavior of the model in terms of the model pa- rameters (e.g., plot graphs of the dependence of the attributes of activity as a function of the model parameters).
Expert Answer:
Related Book For
Posted Date:
Students also viewed these accounting questions
-
Describe inventory turnover and how it can be used to manage inventory.
-
Describe how a SWOT analysis can be used to conduct a first-pass assessment of whether an idea is likely to become a viable business opportunity.
-
Describe several methods that can be used to relax constrained resources.
-
Using the following information, determine the activity rate for the quality inspections activity. Budgeted Activity Cost Pool Activity Setups Estimated Activity Base 1,000 $60,000 Purchase orders...
-
What changes must occur in a production environment to effectively implement JIT? Why are these changes necessary?
-
What is the symbol of an ion with 16 protons, 18 neutrons, and 18 electrons? What is the symbol for an ion that has 16 protons, 16 neutrons, and 18 electrons?
-
If your course involves exercises with a CFD code, study the manual to determine whether the discretization uses staggered or colocated grids. Does the manual say anything about the exact mass...
-
Soto Corporations balance sheet indicates that the company has $300,000 invested in operating assets. During 2011, Soto earned operating income of $45,000 on $600,000 of sales. Required a. Compute...
-
Consider the following two mutually exclusive projects: Year Cash Flow (A) Cash Flow (B) 0 -$ 420,000 -$ 37,500 1 46,500 19,900 2 59,500 13,800 3 4 76,500 16,100 535,000 12,900 The required return on...
-
Lowes Companies, Inc., is a home improvement company offering a range of products for maintenance, repair, remodeling, and decorating. During the recovery phase since the financial crisis of 2008,...
-
The decomposition of hydrogen iodide 2HI(g) H(g) + I2(g) has a rate constant of 9.51x10- L/mols at 500K with an activation energy Ea of 176kJ/mol. What is the rate constant at 600K? J 1)-(5)- 8.314...
-
What is the maximum pressure a cement block of 100 lb can generate if it is 2*3*8 in?
-
A sentimental consulting company is in charge of putting people who could have an affinity in contact so that they can meet. For this they have a program in C++ (thousands of lines of code) that is...
-
You will draft a client-server system for writing/storing notes. Each note has an ID (integer), a title (short text), and a body (long text). A client can has three kinds of requests to the server:...
-
We can link root cause analysis to a culture of safety. In today's world of quality, organizations that embrace a "no blame" policy are better equipped to open communication when things go wrong. ...
-
In the interval [-7, 7] the function d,, is defined: Sn(x) = n for 0 for = x <1/(2n) x>1/(2n) (a) Find the Fourier series of dn(x) as a Fourier cosine and sine expansion. (b) Simplify your result (a)...
-
Premium Amortization On the first day of the fiscal year, a company issues a $3,500,000, 8%, 6-year bond that pays semiannual interest of $140,000 ($3,500,000 8% x 2), receiving cash of $4,038,534....
-
Software Solution is family-owned business that has been in operation for more than 15 year. The board of directors is comprised of mainly family members, plus a few professionals such as an...
-
The current I flowing through a semiconductor diode is given by I = Is(eqvd/kt 1) where Is = 10-12 A is the saturation current, q = 1.6 Ã 10-19 C is the elementary charge value, k = 1.38...
-
The ladder of a fire truck can be elevated (increase of angle ), rotated about the z axis (increase of angle ), and extended (increase of r). Initially the ladder rests on the truck ( > = 0, = 0,...
-
According to special relativity, a rod of length L moving at velocity v will shorten by an amount 8, given by: = L(1 - 1 v2/c2 where c is the speed of light (about 300 x 106 m/s). Consider a rod of...
-
Yale Corporation issued to Zap Corporation \(\$ 60,000,8 \%\) (cash interest payable semiannually on June 30 and December 31) 10 -year bonds dated and sold on January 1, 2020. Assume that the company...
-
Yale Corporation issued to Zap Corporation \(\$ 60,000,8 \%\) (cash interest payable semiannually on July 1 and January 1) 10 -year bonds dated and sold on January 1,2020 . If the bonds were sold at...
-
Lacey Corp. issued a three-year, \(\$ 5,000\) note with an \(8 \%\) stated rate to Hayley Co. on January 1, 2020, and received cash of \(\$ 5,000\). The note requires semiannual interest payments on...
Study smarter with the SolutionInn App