
Elements Of Electromagnetics
7th Edition
ISBN: 9780190698614
Author: Sadiku, Matthew N. O.
Publisher: Oxford University Press
expand_more
expand_more
format_list_bulleted
Question
Can you use MATLAB to help me solve part (d)
![We will consider a linear system and a nonlinear system under uncertainty, each expressed in the form
of a set of stochastic differential equation (SDE) as follows:
dx = (Ax+ Bu)dt + Gdw,
dx = f(x, u, t)dt+ Gdw,
(1)
(2)
where is the state, u is the control, and dw is a differential increment of standard Brownian motion, i.e.,
E[dw] 0 and E [dw(t)dw(t)] = dt. I.
Problem Set 9 Linear Stochastic Process
In this problem, we consider the linear SDE, Eq. (1), with a very simple model where x = R², u = [0,0]T
(no control), and dw R². The matrices A, B, and G are given as follows:
01
A=02x2, B=02x2, G
-69
(3)
where σp E R represents the degree of the uncertainty, and let us take σ₁ = 2 and 02 3. Assume that the
initial state is deterministic and (t = 0) = [0,0]. Take the following steps to simulate the given SDE for
t = [0, 1]:
(a): Consider the increments of w between each time interval t€ [tk,tk+1). Derive the analytical expression
of Awk using w~N(02, 12), where Awk = w(tk+1) - w(tk).
(b): Generate and plot the time history of Awk, Vk, with Atk = tk+1-tk = 10-3, Vk, for sample number
M = 20; include 3-0 bounds in the plot and discuss the consistency with the Monte Carlo result.
Hint: an increment of standard Brownian motion is white Gaussian noise.
(c): Derive the approximate continuous-time EoM from Eq. (1) by assuming that the same noise is applied
to the system over the interval t = [tk, tk+1) while satisfying the increment Awk derived in (a).
(d): Perform Monte Carlo simulation (again M = 20) by propagating the linear SDE with the approx-
imated Brownian motion, and show the time history of each element of a over time; include 3-0
bounds (i.e., ±30) in the plot and discuss the consistency.](https://content.bartleby.com/qna-images/question/ad0d55fe-d83b-4711-86a1-cee8ecea510f/dec12682-e91a-4c9a-a50d-9d070c9829ac/iw1z5m_thumbnail.png)
Transcribed Image Text:We will consider a linear system and a nonlinear system under uncertainty, each expressed in the form
of a set of stochastic differential equation (SDE) as follows:
dx = (Ax+ Bu)dt + Gdw,
dx = f(x, u, t)dt+ Gdw,
(1)
(2)
where is the state, u is the control, and dw is a differential increment of standard Brownian motion, i.e.,
E[dw] 0 and E [dw(t)dw(t)] = dt. I.
Problem Set 9 Linear Stochastic Process
In this problem, we consider the linear SDE, Eq. (1), with a very simple model where x = R², u = [0,0]T
(no control), and dw R². The matrices A, B, and G are given as follows:
01
A=02x2, B=02x2, G
-69
(3)
where σp E R represents the degree of the uncertainty, and let us take σ₁ = 2 and 02 3. Assume that the
initial state is deterministic and (t = 0) = [0,0]. Take the following steps to simulate the given SDE for
t = [0, 1]:
(a): Consider the increments of w between each time interval t€ [tk,tk+1). Derive the analytical expression
of Awk using w~N(02, 12), where Awk = w(tk+1) - w(tk).
(b): Generate and plot the time history of Awk, Vk, with Atk = tk+1-tk = 10-3, Vk, for sample number
M = 20; include 3-0 bounds in the plot and discuss the consistency with the Monte Carlo result.
Hint: an increment of standard Brownian motion is white Gaussian noise.
(c): Derive the approximate continuous-time EoM from Eq. (1) by assuming that the same noise is applied
to the system over the interval t = [tk, tk+1) while satisfying the increment Awk derived in (a).
(d): Perform Monte Carlo simulation (again M = 20) by propagating the linear SDE with the approx-
imated Brownian motion, and show the time history of each element of a over time; include 3-0
bounds (i.e., ±30) in the plot and discuss the consistency.
Expert Solution

This question has been solved!
Explore an expertly crafted, step-by-step solution for a thorough understanding of key concepts.
Step by stepSolved in 2 steps with 5 images

Knowledge Booster
Similar questions
- Can you use MATLAB?arrow_forwardWe will consider a linear system and a nonlinear system under uncertainty, each expressed in the form of a set of stochastic differential equation (SDE) as follows: = da (Ax+ Bu)dt + Gdw, dx = f(x,u,t)dt+Gdw, (1) (2) where x is the state, u is the control, and dw is a differential increment of standard Brownian motion, i.e., E[dw] = 0 and E [dw(t)dw(t)] = dt-1. In this problem, we consider the linear SDE, Eq. (1), with a very simple model where x = R², u = [0,0] (no control), and dw R². The matrices A, B, and G are given as follows: A=02x2, B=02x2; G = [000] (3) where σp Є R represents the degree of the uncertainty, and let us take σ₁ = 2 and σ2 = 3. Assume that the initial state is deterministic and e(t = 0) = [0,0]. Take the following steps to simulate the given SDE for 1 € [0, 1]: (a): Consider the increments of w between each time interval [+1). Derive the analytical expression of Aw using w~N(02, 12), where Aw, w(tk+1) — w(tk). (c): Derive the approximate continuous-time EoM from…arrow_forwardCan you help me by providing the MATLAB code?arrow_forward
- 3. The relationship between arterial blood flow and blood pressure in a single artery satisfies the following first-order differential equation: dP(t) + dt RC mmHg (cm³/s) P(t) = where Qin is the volumetric blood flow, R is the peripheral resistance, and C is arterial compliance (all constant). Qin-60 cm³/s and the initial arterial pressure is 6 mmHg. Also, assume R = 4 and C= 0.4- Oin cm³ mmHg (a) Find the transient solution Ptran(t) for the arterial pressure. The unit for P(t) is mmHg. (b) Determine the steady-state solution Pss(t) for the arterial pressure. (c) Determine the total solution P(t) assuming that the initial arterial pressure is 0.arrow_forwardThe relative displacement u(t) of a single-storey shear building subjected to an earthquake ground motion is represented by the following second-order linear ordinary differential equation: d?u dt2 du + c + ku = a, (t) m dt where m, c, and k are the mass (kg), damping constant (Ns/m), and stiffness of the structure (N/m), respectively. Meanwhile, a,(t) is the function of earthquake ground acceleration. Suppose the building with a mass of 2000 kg and supported by columns of combined stiffness of 32 x 103 N/m is subjected to earthquake with ground acceleration given by the following function: ag(t) = 36000 cos 2t Find the equation of the displacement, u(t), given the damping is not installed to the building. Then, find how much the building is displaced from t =30 seconds to t =60 seconds of earthquake.arrow_forwardCan you help me by providing the solutions in MATLAB code?arrow_forward
- The stress profile shown below is applied to six different biological materials: Log Time (s] The mechanical behavior of each of the materials can be modeled as a Voigt body. In response to o,= 20 Pa applied to each of the six materials, the following responses are obtained: 2 of Maferial 6 Material 5 0.12 0.10 Material 4 0.08 Material 3 0.06 0.04 Material 2 0.02 Material 1 (a) Which of the materials has the highest Young's Modulus (E)? Why? Log Time (s) (b) Using strain value of 0.06, estimate the coefficient of viscosity (n) for Material 6. Stress (kPa) Strainarrow_forwardThe governing equation of motion for a base motion system is given by (assume the units are Newton) mä(t)+ci(t)+kx(t) =cY@, cos(@,t) +kY sin(@t) %3D Given that m = 180 kg, c = 30 kg/s, Y = 0.02 m, and о, — 3.5 rad/s %3| 1. Use Excel or Matlab to find the largest value of the stiffness, k that makes the transmissibility ratio less than 0.85 2. Using the value of the stiffness obtained in part (1), determine the transmitted force to the base motion system using Matlab or Excel. 3. Display and discuss your results.arrow_forwardManufacturing process with two variables x1,x2 described by the empirical model: y=bo +b1 x1 + b2 x2 + b12 x1 x2 + b3 (x1)^2 +b4 (x2)^2 please refere to the image attachedarrow_forward
- The pressure rise, Ap across a centrifugal pump from a given manufacture can be expected to depend on the angular velocity of the impeller w, the Diameter D, of the impeller, the volume flow-rate Q, and the density of the fluid, p. By using the method of repeating variables show that Др ρω- D wD³ A model pump having an impeller diameter of 0.200 m is tested in the laboratory using water. The pressure rise when tested at an angular velocity of 407 rad/sec is shown on the graph. What would be the pressure rise for a geometrically similar pump with an impeller diameter of 0.30 m used to pump water operating at an angular velocity of 807 rad/sec and at a flow rate of 0.070 m³/s? Sol: 40 30 (kPa) Apm 10 0 € 0.02 Q Model data (₁ = 40 rad/s Dm = 20 cm 0.04 0.06 Qm (cubic meters per sec) 0.08arrow_forwardi need help with all parts in this questionarrow_forwardlail - Ahmed Amro Hussein Ali A x n Course: EN7919-Thermodynamic x Homework Problerns(First Law)_S x noodle/pluginfile.php/168549/mod_resource/content/1/Homework%20Problems%28First%20Law%29_SOLUTIONS.pdf UTIONS.pdf 4 / 4 100% Problem-4: When a system is taken from a state-a to a state-b, the figure along path a-c-b, 84 kJ of heat flow into the system, and the system does 32 kJ of work. (i) How much will the heat that flows into the system along the path a-d-b be, if the work done is 10.5 kJ? (Answer: 62.5 kJ) (ii) When the system is returned from b to a along the curved path, the work done on the system is 21 kJ. Does the system absorb or liberate heat, and how much? (Answer:-73 k]) (iii) If Ua = 0 andU, = 42kJ , find the heat absorbed in the processes ad and db. (Answer: 52.5 kJ, 10 kJ)arrow_forward
arrow_back_ios
SEE MORE QUESTIONS
arrow_forward_ios
Recommended textbooks for you
- Elements Of ElectromagneticsMechanical EngineeringISBN:9780190698614Author:Sadiku, Matthew N. O.Publisher:Oxford University PressMechanics of Materials (10th Edition)Mechanical EngineeringISBN:9780134319650Author:Russell C. HibbelerPublisher:PEARSONThermodynamics: An Engineering ApproachMechanical EngineeringISBN:9781259822674Author:Yunus A. Cengel Dr., Michael A. BolesPublisher:McGraw-Hill Education
- Control Systems EngineeringMechanical EngineeringISBN:9781118170519Author:Norman S. NisePublisher:WILEYMechanics of Materials (MindTap Course List)Mechanical EngineeringISBN:9781337093347Author:Barry J. Goodno, James M. GerePublisher:Cengage LearningEngineering Mechanics: StaticsMechanical EngineeringISBN:9781118807330Author:James L. Meriam, L. G. Kraige, J. N. BoltonPublisher:WILEY

Elements Of Electromagnetics
Mechanical Engineering
ISBN:9780190698614
Author:Sadiku, Matthew N. O.
Publisher:Oxford University Press

Mechanics of Materials (10th Edition)
Mechanical Engineering
ISBN:9780134319650
Author:Russell C. Hibbeler
Publisher:PEARSON

Thermodynamics: An Engineering Approach
Mechanical Engineering
ISBN:9781259822674
Author:Yunus A. Cengel Dr., Michael A. Boles
Publisher:McGraw-Hill Education

Control Systems Engineering
Mechanical Engineering
ISBN:9781118170519
Author:Norman S. Nise
Publisher:WILEY

Mechanics of Materials (MindTap Course List)
Mechanical Engineering
ISBN:9781337093347
Author:Barry J. Goodno, James M. Gere
Publisher:Cengage Learning

Engineering Mechanics: Statics
Mechanical Engineering
ISBN:9781118807330
Author:James L. Meriam, L. G. Kraige, J. N. Bolton
Publisher:WILEY