Info

이 질문은 마감되었습니다. 편집하거나 답변을 올리려면 질문을 다시 여십시오.

can anyone help me with "ode45" ?

조회 수: 1 (최근 30일)
ADNAN KIRAL
ADNAN KIRAL 2018년 11월 25일
마감: MATLAB Answer Bot 2021년 8월 20일
Hi
i do have a function like;
function [xdot,NOD]=random_func(t,x)
.
.
.
end
and a "ode45" like
[t,x]=ode45(@random_func,tspan,x0,cn3);
I would like to take "NOD" out given in my function;
How can I get it in a ".mat" ? is that possible? Or I need to use the Runge Kutta method manually?
Please a bit urgent, thanks in advance. If you want a clear understanding, I can upload a draft code.
  댓글 수: 2
James Tursa
James Tursa 2018년 11월 26일
편집: James Tursa 2018년 11월 26일
Note that you cannot do anything random inside a derivative function that affects the derivative calculation. This violates an assumption that ode45 makes, will severly confuse it, and will cause it to either fail outright or give you a wrong answer.
ADNAN KIRAL
ADNAN KIRAL 2018년 11월 27일
i see, will consider that. thanks for your reply.

답변 (1개)

Star Strider
Star Strider 2018년 11월 25일
Probably the easiest way would be to integrate your function first, the using the solved values for ‘t’ and ‘x’, run it again, recover the ‘NOD’ output, and save it.
Example —
function [xdot,NOD]=random_func(t,x)
xdot(1,:) = t-x(1); % Create Equation
xdot(2,:) = t+x(2); % Create Equation
NOD = rand; % Create ‘NOD’
end
tspan = [0 1];
x0 = [1;1];
[T,X] = ode45(@random_func, tspan, x0);
figure
plot(T, X)
grid
NODv = zeros(numel(T),1); % Preallocate
for k1 = 1:numel(T)
[~,NOD(k1,:)] = random_func(T(k1),X(k1,:)); % Calculate ‘NOD’, Save As Vector
end
fprintf(1, '%.4f,%.4f\n', [T NOD]') % Save Data
I simply printed out the ‘NOD’ vector here (alternatively saving it as a ‘.csv’ file with the same statement although with an appropriate file identification number created by fopen), since I did not want to create the file and then have to delete it. Use the appropriate save syntax to create your ‘.mat’ file.
This illustrates the idea. You will likely need to experiment with it and tweak it to make it work with your functions.
  댓글 수: 2
ADNAN KIRAL
ADNAN KIRAL 2018년 11월 26일
thanks for reply.
Star Strider
Star Strider 2018년 11월 26일
My pleasure.

제품


릴리스

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by