Returning several dataset using ode45

조회 수: 1 (최근 30일)
Frederik Bjerregaard
Frederik Bjerregaard 2021년 12월 7일
댓글: Frederik Bjerregaard 2021년 12월 8일
Hi. I have a problem where I want an ode45 function to return 5 dataset using the code
clc; clear; close all;
% Defining parameters and initial conditions
f = 0.35;
r = [0.3:0.3:1.5];
zeta = 0.25;
Theta0 = [0.30, 0];
tspan = [0 10];
% Solver options
opts = odeset('RelTol',1e-6,'AbsTol',[1e-9 1e-9]);
% Solver
[t,theta] = ode45(@(t,theta) odefcn(t,theta,f,r,zeta), tspan, Theta0);
plot(theta(:,1),theta(:,2));
As you can see, r is a parameter in the function and I want to create a plot for each of the values for r. The code is running fine when using a single value for r but breaks when using a vector.

채택된 답변

VBBV
VBBV 2021년 12월 7일
for i = 1:length(r)
[t,theta] = ode45(@(t,theta) odefcn(t,theta,f,r,zeta), tspan,Theta0);
plot(theta(:,1),theta(:,2));hold on;
end
Use a for loop and try it
  댓글 수: 3
Frederik Bjerregaard
Frederik Bjerregaard 2021년 12월 8일
No it still gives me an error. Just to be sure I have incorporated it correctly, was this what you meant?
clc; clear;close all;
% Defining parameters and initial conditions
f = 0.35+0.1;
r = [0.3:0.3:1.5];
zeta = 0.25;
Theta0 = [0.30, 0];
tspan = [0 1000];
% Solver options
opts = odeset('RelTol',1e-6,'AbsTol',[1e-9 1e-9]);
% Solver
for i = 1:length(r)
[t,theta] = ode45(@(t,theta) odefcn(t,theta,f,r,zeta), tspan,Theta0);
plot(theta(:,1),theta(:,2));hold on;
end
I will just include the function file as well for good measure
function dthetadt = odefcn(t,theta,f,r,zeta)
dthetadt = zeros(2,1);
dthetadt(1)=theta(2);
dthetadt(2)=f.*sin(r.*t)-2.*zeta.*theta(2)-sin(theta(1));
end
Frederik Bjerregaard
Frederik Bjerregaard 2021년 12월 8일
I got it working now. Thank you!

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기

태그

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by