大家好,又见面了,我是你们的朋友全栈君。
举例为如下方程 y 1 ′ ( t ) = y 1 ( t − 1 ) y 2 ′ ( t ) = y 1 ( t − 1 ) + y 2 ( t − 0.2 ) y 3 ′ ( t ) = y 2 ( t − 1 ) \begin{aligned} y_1′(t) &= y_1(t-1) \\ y_2′(t) &= y_1(t-1) + y_2(t-0.2) \\ y_3′(t) &= y_2(t-1) \\ \end{aligned} y1′(t)y2′(t)y3′(t)=y1(t−1)=y1(t−1)+y2(t−0.2)=y2(t−1)
%% 测试求解常时滞微分方程
% Author: Zhao-Jichao
% Date: 2021-07-05
clear
clc
%% Define Time-Delay
lags = [1 0.2];
tspan = [0 5];
sol = dde23(@ddefun, lags, @history, tspan);
%% Draw Results
plot(sol.x,sol.y,'-o')
xlabel('Time t');
ylabel('Solution y');
legend('y_1','y_2','y_3','Location','NorthWest');
%% DDE Function
function dydt = ddefun(t,y,Z) % equation being solved
ylag1 = Z(:,1);
ylag2 = Z(:,2);
dydt = [ylag1(1);
ylag1(1)+ylag2(2);
y(2)];
end
%-------------------------------------------
function s = history(t) % history function for t <= 0
s = ones(3,1);
end
%-------------------------------------------
y 1 ′ = y 1 ( t − 0.1 ) y 2 ′ = y 2 ( t − 0.5 ) \begin{aligned} y_1’&=y_1(t-0.1) \\ y_2’&=y_2(t-0.5) \\ \end{aligned} y1′y2′=y1(t−0.1)=y2(t−0.5)
写成了两个并列的形式,方便对比最终的效果。
代码如下
%% 测试求解常时滞微分方程
% Author: Zhao-Jichao
% Date: 2021-07-05
clear
clc
%% Define Time-Delay
lags = [0.1 0.5];
tspan = [0 5];
sol = dde23(@ddefun, lags, @history, tspan);
%% Draw Results
plot(sol.x,sol.y,'-o')
xlabel('Time t');
ylabel('Solution y');
legend('y_1','y_2','y_3','Location','NorthWest');
%% DDE Function
function dydt = ddefun(t,y,Z) % equation being solved
ylag1 = Z(:,1);
ylag2 = Z(:,2);
dydt = [-ylag1(2);
-ylag2(2)];
end
%-------------------------------------------
function s = history(t) % history function for t <= 0
s = ones(2,1);
end
%-------------------------------------------
方程为: y 1 ′ = y 2 ( t − 0.1 ) y 2 ′ = − 2 y 3 ′ = y 4 ( t − 0.5 ) y 4 ′ = − 2 \begin{aligned} y_1′ &= y_2(t-0.1) \\ y_2′ &= -2 \\ y_3′ &= y_4(t-0.5) \\ y_4′ &= -2 \\ \end{aligned} y1′y2′y3′y4′=y2(t−0.1)=−2=y4(t−0.5)=−2
%% 测试求解常时滞微分方程
% Author: Zhao-Jichao
% Date: 2021-07-05
clear
clc
%% Define Time-Delay
lags = [0.1 0.5];
tspan = [0 3];
sol = dde23(@ddefun, lags, @history, tspan);
%% Draw Results
plot(sol.x,sol.y)
xlabel('Time t');
ylabel('Solution y');
legend('y_1','y_2','y_3','y_4','Location','NorthWest');
%% DDE Function
function dydt = ddefun(t,y,Z) % equation being solved
ylag1 = Z(:,1);
ylag2 = Z(:,2);
dydt = [ylag1(2);
-2;
ylag2(4);
-2;];
end
%-------------------------------------------
function s = history(t) % history function for t <= 0
s = ones(4,1);
end
%-------------------------------------------
Ref: 时滞微分方程 — 示例
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/188839.html原文链接:https://javaforall.cn