
I having problems with MATLAB ode45 solver. For the following code, Matlab is giving me index errors. All I am trying to do is integrate the differential equations in a ode45 function. But I don't see what the problem here is. Can you help me fix this code.
% Initial Conditions
w = [2;-1;1];
T = 30;
I = [150; 400; 400];
EP = [0;0;0;1];
C = eye(3,3);
tspan = 0:20;
% Using ode45 to integrate w dot
options = odeset('RelTol',1e-10,'AbsTol',1e-10);
result = ode45(@K_DDE, tspan, [w';C], options);
% Extracting information from ode solver
t = result.x;
w_ode = result.y(1,:);
C_ode = result.y(2:4,:)
function dwCdt = K_DDE(~,w_C)
w = w_C(1,:);
C = w_C(2:4,:);
% Kinematic Differentia Equations for DCM
dCdt = eye(3,3);
dCdt(1,1) = C(1,2)*w(3) - C(1,3)*w(2);
dCdt(1,2) = C(1,3)*w(1) - C(1,1)*w(3);
dCdt(1,3) = C(1,1)*w(2) - C(1,2)*w(1);
dCdt(2,1) = C(2,2)*w(3) - C(2,3)*w(2);
dCdt(2,2) = C(2,3)*w(1) - C(2,1)*w(2);
dCdt(2,3) = C(2,1)*w(2) - C(2,2)*w(1);
dCdt(3,1) = C(3,2)*w(3) - C(3,3)*w(2);
dCdt(3,2) = C(3,3)*w(1) - C(3,1)*w(3);
dCdt(3,3) = C(3,1)*w(2) - C(3,2)*w(1);
dwCdt = [w'; dCdt];
end

Trending nowThis is a popular solution!
Step by stepSolved in 3 steps with 3 images

- Write a MATLAB function that calculates the approximate value of six(x) using the Maclaurin series approximation from x-0 to xi+1 six(x) = x -+-+. 3 5 The function should accept 3 parameters: value of x, number of significant figures accuracy i.e. n, and the maximum number of iterations. In the function, use ɛ = (0.5 × 10²-" )% in order to continue until the ɛa falls below this criteria. The function should return 3 values: the approximate value of arctan(x) at the end of the program, final ɛg and the number of iterations it took. Submit short snapshots showing the program running.arrow_forwardCan you help me program in MATLAB? I want to plot the magnitude of Modified Rodrigues Parameter (MRP) vs theta. I need 2 plots where theta equals (0,pi) and (pi,2*pi). Also, I need the same plots but for the shadow set of MRP vs theta.arrow_forwardUsing MATLAB : Please explain the process of how you constructed your code and provide screenshots of the output, pictures and grapharrow_forward
- Database System ConceptsComputer ScienceISBN:9780078022159Author:Abraham Silberschatz Professor, Henry F. Korth, S. SudarshanPublisher:McGraw-Hill EducationStarting Out with Python (4th Edition)Computer ScienceISBN:9780134444321Author:Tony GaddisPublisher:PEARSONDigital Fundamentals (11th Edition)Computer ScienceISBN:9780132737968Author:Thomas L. FloydPublisher:PEARSON
- C How to Program (8th Edition)Computer ScienceISBN:9780133976892Author:Paul J. Deitel, Harvey DeitelPublisher:PEARSONDatabase Systems: Design, Implementation, & Manag...Computer ScienceISBN:9781337627900Author:Carlos Coronel, Steven MorrisPublisher:Cengage LearningProgrammable Logic ControllersComputer ScienceISBN:9780073373843Author:Frank D. PetruzellaPublisher:McGraw-Hill Education





