How to run function within for loop

Question

Mirlan Karimov 2019-3-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/451605-how-to-run-function-within-for-loop

评论： Mirlan Karimov 2019-3-22

The code given below solves the delayed type differential equation. I want to run this code for >10,000 different combinations of

and

. That is,

and

;

function [sol]=system_simulation
po=1.4; do=8; tau=0.5; xth=0.2; vth=0.2;ksi=0.25; 
 
option=odeset('RelTol', 1e-4, 'AbsTol', 1e-4);
sol = dde23(@syst,[tau],@hystory,[0 100],option);
figure(1), plot(sol.x,sol.y(1,:));
% --------------------------------------------------------------------------
    function s = hystory(t)
        % Constant history function for inverted pendulum.
        s = [0.4;0.1];
    end
% --------------------------------------------------------------------------
    function xprime = syst(t,x,Z)
        xlag=Z(:,1);
        xprime=zeros(2,1);
        xprime(1)=x(2);
        xprime(2)=x(1)-p*xlag(1)-d*xlag(2);
                 
    end
end

2 个评论
显示无隐藏无

Geoff Hayes 2019-3-21

Mirlan - you want to iterate over po and do but nowhere do you reference them. How should they be used? (It is also doubtful you want to plot the results for all 10000 iterations so please discuss what you want to save/store from each iteration.

Mirlan Karimov 2019-3-22

No plotting, I just want to store the data of sol.y(1,:) at each time iteration. Oh sorry it is p and d. There are some additional steps that turn po and do to p and d respectively depending on conditions but it has nothing to do now so I removed them.

请先登录，再进行评论。

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

Answer 1

Jan 2019-3-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/451605-how-to-run-function-within-for-loop#answer_366838

编辑：Jan 2019-3-22

在 MATLAB Online 中打开

I prefer providing parameters by anonymous functions instead of nested function. See Answers: Anonymous function for params . Then:

pv = 0.5:0.1:14;
dv = 0.5:0.1:14;
Result = cell(numel(pv), numel(dv));
ip = 0;
for p = pv
    ip = ip + 1;
    id = 0;
    for d = dv
        sol = dde23(@(t,x,Z) syst(t,x,Z, p, d), tau, ...
                    @hystory, [0, 100], option);
        id  = id + 1;
        Result{ip, id} = sol.y(1,:);
    end
end