How to use my numeric array results as input values for an ODE?

1 次查看（过去 30 天）

显示更早的评论

Thomas Veith 2019-6-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/465507-how-to-use-my-numeric-array-results-as-input-values-for-an-ode

评论： Star Strider 2019-6-4

采纳的回答： Star Strider

在 MATLAB Online 中打开

So I have the following code:

n = 10 ; 
result = zeros(n,5 ); 
for k=1:n ; 
A0=1;P0=29;g=rand;p=rand;B=rand ; 
result(k,:) = [A0,P0,g,p,B ]; 
end;

Which produces a 10x5 numeric array as the result.

And I'd like to use my results from there as the input values to solve this equation ten times:

tspan=[1 7];
[t,x] = ode45(@(t,x) [-g*x(1) + p*x(1); -x(1)*x(2)+ B*x(2)], tspan, [A0 P0]);

Is this possible?

Thanks in advance!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Star Strider 2019-6-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/465507-how-to-use-my-numeric-array-results-as-input-values-for-an-ode#answer_377850

在 MATLAB Online 中打开

Try this:

n = 10 ; 
result = zeros(n,5 ); 
for k=1:n ; 
A0=1;P0=29;g=rand;p=rand;B=rand ; 
result(k,:) = [A0,P0,g,p,B ]; 
end
for k = 1:size(result,1)
    tspan=linspace(0,7,50);
    ODE = @(t,x,g,p,B) [-g*x(1) + p*x(1); -x(1)*x(2)+ B*x(2)];
    [t,x(:,:,k)] = ode45(@(t,x) ODE(t,x,result(k,3),result(k,4),result(k,5)), tspan, [result(k,1) result(k,2)]);
end