MatLAB 中解时滞微分方程 怎么写?
MatLAB中解时滞微分方程sol=dde23(ddefun,tau,history,tspan,options)
ddetun:描述延迟微分方程的句柄
tau=[tau1,tau2,...,taun]
延迟微分方程function dx=ddefun(t,x,z)
我想问的是 这里的Z表示的是什么?是个关于延迟的向量?
例如
function dx=ddefun(t,x,z)
dx=[1-3*x(1)-A-0.2*B^3-tau2(1)
x(3)
4*x(1)-2*x(2)-3*x(3)];
我这里A想表示为x2(t-tau1),B想表示为x1(t-tau2),这个A和B该怎么写?这里Z具体表示的是什么东西?
希望有人能回答,
人气:297 ℃ 时间:2020-06-20 09:33:02
解答
这是matlab中dde23的例子,通过这个例子,应该能看懂dde23个参数的作用.直接复制后边的代码就可以输出图形.%DDEX1 Example 1 for DDE23.% This is a simple example of Wille' and Baker that illustrates the% straig...Z表示x1(t-tau1) x2(t-tau1) x3(t-tau1) ...xn(t-tau1) 这样的n*n一个矩阵? x1(t-tau2) x2(t-tau2) x3(t-tau2) ...xn(t-tau2)................................................................. x1(t-taun) x2(t-taun) x3(t-taun) ...xn(t-taun) 所以ylag1 = Z(:,1)就表示成一列向量[x1(t-tau1)x1(t-tau2) ...x1(t-taun)]'所以ylag1(1)表示的就是x1(t-tau1)我不知道我这么理解对不?Z是初始值,在例子中y是一个3*1的列向量,所以z也是3*1,z的维数看你具体的方程???tau1,tau2....taun 肯定都是初始值。Z是3*1的列向量?那ylag2 = Z(:,2);如何理解?这个不是表示取第二列所有行的元素吗?在ddefun中的z是y相对于t-tau1,t-tau2,t-tau3....t-tauj时的历史数据的集合
推荐
猜你喜欢
- M=N=x=1 rcosA,y=-1 rsinA,
- 1x1-2x2+3x3-4x4******+99x99-100x100+101x101的过程
- 若a,b∈R,且a+b>1,ab>1,则a>1,b>1是真命题或假命题
- 甲、乙两件服装的成本共500元,商店老板为获取利润,决定将甲服装按50%的利润定价,乙服装按40%的利润定价.在实际出售时,应顾客要求,两件服装均按9折出售,这样商店共获利157元,求
- 甲乙两个书柜一共有400本书 ,从甲书柜中拿出50本,从乙书柜中拿出10那么两个书柜中的书相同 ,问甲乙书柜原来有多少本书?
- 英语翻译
- 谁知道:氢气在水中的溶解性?
- 3+2-5*0=?等于多少