Question

求大神用matlab的dde23解带时滞的微分方程并画出图来，要完整的程序

求大神用matlab的dde23解带时滞的微分方程并画出图来，要完整的程序可能不一定是dde23，只要解出来就行，毕设的东西，挺急的，路过的大神帮帮忙，谢谢了

Answer 1

请把D，A,k的值补充下，我好写个程序。

追问

D,A就随便定义为2阶矩阵吧，k＝2

比如D=［1，0；0，1］,A=［1，1；1，0］，其实我处理的是几十阶的矩阵，有了程序以后我还要自己改

麻烦了，谢谢大神，程序运行出来，我一定采纳

大神，写出来了吗？

追答

一般问题都是很有针对性，你不把真实的常数写出来，而是随便写个，我也只能随便写了：
function ddemain
clear
clc
lags = 2*ones(4,1);%k=2
sol = dde23(@ddefun,lags,@ddex1hist,[-1,5]);
plot(sol.x,sol.y);
title('dydt=-D*y+A*y(t-k)方程图');
xlabel('time t');
ylabel('solution y');
tint = linspace(0,5,10);
Sint = deval(sol,tint);
hold on
plot(tint,Sint,'o');
function dydt = ddefun(t,y,z)
D=[1 0 0 1]';A=[1 1 1 0]';
ylag=[z(1 ) z( 2) z(3) z(4)]';
dydt =-D.*y+A.*ylag;
end
function S = ddex1hist(t)
S = ones(2,2);
end
end

追问

大神，您给我的程序D,A都设成了单列的矩阵，如果设成方阵怎么做，其它参数怎么改

比如D=［1，0；0，1］,A=［1，1；1，0］

这是我真正的需要处理的方阵，其实有35阶，后面我还有一个62阶和一个115阶的，所以我需要一个模板代码，

大神可不可以帮我改成适用于方阵的，谢谢

追答

可以的，把数据发过来，写个有针对性的程序。

追问

A=rand(35)>0.8;

E=diag(sum(A'));

C= round(rand(35,1));

B = diag(C);

D=B+E;

这是我我的A和D矩阵，35阶的

k就是取常数2

做出来的图类似这种，不过我的是35阶的，所以应该有35条线

这是多智能体研究里面的收敛图像

追答

35阶是35*35个方程，要35*35条曲线，是个大型问题。

追问

对啊，所以当初我先让你做一个2阶的嘛，后面我再上升阶数

可是之前你发给我的程序，A和D都写成了4行1列的了，我不知道该怎么改了

追答

你把矩阵改成行就行了。

追问

怎么改啊，其它参数有没有变动啊，大神把之前的代码改一下吧，注释一下如果改变矩阵哪些地方也需要改

因为我要研究的只有方阵，所以希望大神把之前的代码改为方阵的，并且注释一下改变方阵的话，哪些地方也需要改，谢谢大神了