understanding DDE23 function format

Question

alpedhuez 2020-5-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/522327-understanding-dde23-function-format

编辑： alpedhuez 2020-5-4

采纳的回答： Guru Mohanty

在 MATLAB Online 中打开

I look at this example https://www.mathworks.com/help/matlab/math/dde-with-constant-delays.html

DDE23 function has a format of

ddex1de(t,y,Z)

How should one understand Z?

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

请先登录，再进行评论。

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

Answer 1

Guru Mohanty 2020-5-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/522327-understanding-dde23-function-format#answer_430072

在 MATLAB Online 中打开

I understand you are trying to solve system of differential equation using ‘dde23’. To do this you need to do following steps.

Define constant delays.

In this system of equation there are two lags i.e. (t-1) and (t-0.2).

lags = [1 0.2];

2. Define Solution History.

It Defines the first solution from which the solver starts iterations.

        function s = history(t)
            s = ones(3,1);
        end

3. Form the equation.

        function dydt = ddex1de(t,y,Z)
            ylag1 = Z(:,1);
            ylag2 = Z(:,2);
            dydt = [ylag1(1); 
                    ylag1(1)+ylag2(2); 
                    y(2)];
        end

Here In the call to ddex1de,

‘t’ is a scalar indicates the current ‘t’ in the equation

‘y’ is a column vector approximates y(t)

‘Z’ is a column vector approximates y(t – αj) for delay αj= lags(J).

In the following example Solution history is [1;1;1]. So Approximate values of lag ie ‘Z’ will be a matrix of size 3x2.In ‘Z’ Row defines the number of equations and column defines Number of lags.

4. Solve using ‘dde23’.

        tspan = [0 5];
        sol = dde23(@ddefun, lags, @history, tspan);

2 个评论
显示无隐藏无

alpedhuez 2020-5-4

编辑：alpedhuez 2020-5-4

在 MATLAB Online 中打开

Thank you.

1. Would you write down the content of this 3*2 matrix Z?

2. Another question:

function dydt = ddex1de(t,y,Z)
            ylag1 = Z(:,1);
            ylag2 = Z(:,2);
            dydt = [ylag1(1); 
                    ylag1(1)+ylag2(2); 
                    y(2)];
        end

I would think intuitively we define ddex1de by dydt so it should be like

function ddex1de = dydt(t,y,Z)
            ylag1 = Z(:,1);
            ylag2 = Z(:,2);
            dydt = [ylag1(1); 
                    ylag1(1)+ylag2(2); 
                    y(2)];
end

Woud you explain why this is not the case?

Guru Mohanty 2020-5-4

As it is a System of 3 Equations, 'y', 'dydt', 'ylag1' 'ylag2' have dimension 3x1. There are two constant Delay present, So The dimention of 'Z' is 3x2.

请先登录，再进行评论。