Deterministic Modeling and Simulation 1. This problem describes the immune response triggered by encountering foreign antigens....
Fantastic news! We've Found the answer you've been seeking!
Question:
Transcribed Image Text:
Deterministic Modeling and Simulation 1. This problem describes the immune response triggered by encountering foreign antigens. This is called the cytotoxic T cell or CTL response. Suppose that z eliminates infected cells. There are 4 variables: uninfected cells, x, infected cells, y, free virus particles v, and CTL cells z. The equations are: x=λ-dx-Bxv y=Bxv-ay-pyz v=ky-uv ż=c-bz Uninfected cells are produced at a constant rate 2, die at a rate dx and are lost to infection at a rate of Bxv. Infected cells are produced from uninfected cells and free virus at the rate, xv and die at a rate ay and are lost to CTL cells at a rate of pyz. Free virus is produced from infected cells at a rate ky and die at a rate uv. The CTL cells are produced at a rate ofc and die at a rate of b. Parameter values are λ = 1; d = 0.01; a=0.05; p=0.005; p=1; k=50; u = 5; b = 0.05. a). Let c=0. This produces no CTL response, or what happens without an immune system. Simulate the above systems with initial values of x(0)=10 y (0) = 1; v(0)=1; =(0)=0 Determine the steady state value by running the simulation and by running the "trim" function. Check answers with the equations. b). Let c=10. This is a strong immune response. Repeat part (a) and compare the steady-state results via simulation and the trim function. c). Let c=0.05 to simulate a weak immune response. Repeat part (a) and compare the steady- state results via simulation and the trim function. d). For part a, linearize the system using the "linmod" function around the steady-state result for part a. 2. Consider the following differential equation: dt² a. Write the state equations for this system pai-5 b. Add a new state in which x=y-y, where y, is the unit step and represents the desired output of the system. Write the new state equations treating y, as another input. c. Now, pick a vector k = [k,,k,k.] such that if u = kx, the new states will have roots of the characteristic equation at -2, -3, and -4. d. Simulate the system and observe the output y and compare the response for Y₁ = step input to the step response of the system without control. 3. Consider the control system shown below. This system might represent the control of a DC motor. We want to position the shaft to angle determined by, com- The problem is to choose the parameters Kr, the tachometer feedback and Kp, the motor shaft position, so as to minimize the integral squared error between the command input and the controlled output, y. The command input is com = 5 degrees and is shown as a step input of amplitude 5. Subtract1 Subtract amplifier 5 10+1 motor Kr tach pa dot motor int Kp position psi 4. The system of problem1 has the same parameters except that b, p and beta are different, but are between 0 and 1.5 in problem1. The file prob_1_data.mat has the data from a simulation where the first column is time, the second is x, the third is y, the fourth is v and the fifth is z. In the data, the initial condition for z is 0.5 not 0. Write the Matlab code to estimate the unknown parameters, b p and beta. Use a minimum square error criterion between the data and the model output to generate the parameter estimates. Deterministic Modeling and Simulation 1. This problem describes the immune response triggered by encountering foreign antigens. This is called the cytotoxic T cell or CTL response. Suppose that z eliminates infected cells. There are 4 variables: uninfected cells, x, infected cells, y, free virus particles v, and CTL cells z. The equations are: x=λ-dx-Bxv y=Bxv-ay-pyz v=ky-uv ż=c-bz Uninfected cells are produced at a constant rate 2, die at a rate dx and are lost to infection at a rate of Bxv. Infected cells are produced from uninfected cells and free virus at the rate, xv and die at a rate ay and are lost to CTL cells at a rate of pyz. Free virus is produced from infected cells at a rate ky and die at a rate uv. The CTL cells are produced at a rate ofc and die at a rate of b. Parameter values are λ = 1; d = 0.01; a=0.05; p=0.005; p=1; k=50; u = 5; b = 0.05. a). Let c=0. This produces no CTL response, or what happens without an immune system. Simulate the above systems with initial values of x(0)=10 y (0) = 1; v(0)=1; =(0)=0 Determine the steady state value by running the simulation and by running the "trim" function. Check answers with the equations. b). Let c=10. This is a strong immune response. Repeat part (a) and compare the steady-state results via simulation and the trim function. c). Let c=0.05 to simulate a weak immune response. Repeat part (a) and compare the steady- state results via simulation and the trim function. d). For part a, linearize the system using the "linmod" function around the steady-state result for part a. 2. Consider the following differential equation: dt² a. Write the state equations for this system pai-5 b. Add a new state in which x=y-y, where y, is the unit step and represents the desired output of the system. Write the new state equations treating y, as another input. c. Now, pick a vector k = [k,,k,k.] such that if u = kx, the new states will have roots of the characteristic equation at -2, -3, and -4. d. Simulate the system and observe the output y and compare the response for Y₁ = step input to the step response of the system without control. 3. Consider the control system shown below. This system might represent the control of a DC motor. We want to position the shaft to angle determined by, com- The problem is to choose the parameters Kr, the tachometer feedback and Kp, the motor shaft position, so as to minimize the integral squared error between the command input and the controlled output, y. The command input is com = 5 degrees and is shown as a step input of amplitude 5. Subtract1 Subtract amplifier 5 10+1 motor Kr tach pa dot motor int Kp position psi 4. The system of problem1 has the same parameters except that b, p and beta are different, but are between 0 and 1.5 in problem1. The file prob_1_data.mat has the data from a simulation where the first column is time, the second is x, the third is y, the fourth is v and the fifth is z. In the data, the initial condition for z is 0.5 not 0. Write the Matlab code to estimate the unknown parameters, b p and beta. Use a minimum square error criterion between the data and the model output to generate the parameter estimates.
Expert Answer:
Related Book For
Posted Date:
Students also viewed these mathematics questions
-
The missile weighs 40 000 lb. The constant thrust provided by the turbojet engine is T = 15 000 lb. Additional thrust is provided by two rocket boosters B. The propellant in each booster is burned at...
-
From elementary physics, we know that when an object is subjected to a constant acceleration a, the relationship between distance d and time t is given by d = 1/2 at2. Suppose that, during a seek,...
-
From elementary physics, we know that when an object is subjected to a constant acceleration a, the relationship between distance d and time t is given by d = at2. Suppose that, during a seek, the...
-
(a) Write an equation describing a sinusoidal transverse wave traveling on a cord in the positive direction of a y axis with an angular wave number of 60 cm-1, a period of 0.20 s, and an amplitude of...
-
In sports betting, Las Vegas sports books establish winning margins for a team that is favored to win a game. An individual can place a wager on the game and will win if the team bet upon wins after...
-
What should an employee do when he or she discovers that there is an error in a projection? Why do you suggest that action? Would your answer change if the error was not likely to affect other...
-
True or False: A company's beta can be used to accurately forecast its stock price during the coming year.
-
On April 29, 2013, Quality Appliances purchased equipment for $260,000. The estimated service life of the equipment is six years and the estimated residual value is $20,000. Qualitys fiscal year ends...
-
Do a business reporting on the information management of the business focusing of the following aspect: - Incorporating Improved decision-making with information system of a business: - How Systems...
-
After retiring as a physician, Bob Guthrie became an avid downhill skier on the steep slopes of the Utah Rocky Mountains. As an amateur inventor, Bob was always looking for something new. With the...
-
Drinking water treatment plants put a lot of emphasis on monitoring quality of water released to communities for consumption. The presence of indicator organisms such as E. Coli is used a standard...
-
Kim owes Todd a sum of money for money she has borrowed. They cannot agree on the exact amount and have no record of the loans made over a long period of time. Kim doesn't have all the money to pay...
-
Calculate and categorize the Point elasticity for the (inverse)demand curve P=24-2Q at p1=14,p2=12,p3=10
-
7. An investor is considering a capital expenditure for new production facilities. The relevant cash flows are listed below: (12 points) Initial Investment = $280,000 Year 1 2 3 4 5. a. b. Cashflow...
-
If f(x,y)=Sinyx then find dyn enwolf dx
-
A firm's production function is Q=70L06 K03. Its marginal product of labor is thus MP = 42L-0.4 K0.3 and its marginal product of capital is MPK = 21L0.6 K-0.7 a. What returns to scale does this...
-
Create a list of 16 risks including 2 opportunities. In the Category ensure you have at least one risk in each of the following categories: Financial, Technical, Commercial, Execution,...
-
The figure shows a bolted lap joint that uses SAE grade 8 bolts. The members are made of cold-drawn AISI 1040 steel. Find the safe tensile shear load F that can be applied to this connection if the...
-
Without drawing the MOs, state whether the 7r-molecular orbital 6 in 1,3,5,7,9-decapentaene (a 10-carbon conjugated alkene) is symmetric or anti symmetric with respect to the reference plane; is...
-
Using resonance arguments, state which ion or radical within each set is more stable. Explain. CHj HC-C CH2 or HC CH CH CH2
-
Using the Huckel 4n + 2 rule, determine whether each of the following compounds is likely to be aromatic. Explain how you arrived at the -electron count in each case. (a) (b)
-
You are managing the development of a case tracking system project for your large law firm. The requirements phase of the project is almost complete, and preliminary design work has begun. The...
-
At the conclusion of this chapter, the textbook mentions that data and process modeling may eventually become obsolete due to the increasing popularity and usage of object-oriented modeling and...
-
Although data and process models depict the same system with different views, systemi designers must synchronize these different views to make sure that their models are consistent and complete. One...
Study smarter with the SolutionInn App