资源简介
脉冲时滞分数阶多智能体动力系统的matlab仿真,差分方法,有例子,详尽注释,也可自行更改系数矩阵。
代码片段和文件信息
function With_Delay
clear all; close all; clc;
Dim=15;q(1:Dim1)=0.9;FinalTime=1.5;tao=0.5;%Dim为方程数,q为导数的阶数,FinalTime为t的最大值tao为时滞
h=0.0001;t=-tao+h:h:FinalTime;N=length(t)-1;u=-tao+h:h:FinalTime+h;Ntao=tao/h;
x=zeros(DimN);
for i=1:4%X(t)的初值[-tao+h0] 前12行,即i=1,2,3,4
for j=1:3
for ti=1:Ntao+1
x((i-1)*3+jti)=cos((i+j)*(ti-Ntao-1)*h);
end
end
end
x(131:Ntao+1)=1;
x(141:Ntao+1)=2;
x(151:Ntao+1)=3;
% for i=1:3%X(t)的初值[-tao+h0] 后最3行,即i=0
% for ti=1:Ntao+1
% x(i+12ti)=sin(i*(ti-Ntao-1)*h);
% end
% end
% x(:Ntao+1)=1;
%%
a=zeros(DimN-Ntao);b=zeros(DimN+1-Ntao); a0=zeros(DimN+1-Ntao);
T1=(h.^q./q)./gamma(q);
T2=(h.^q)./gamma(q+2);
for n=1:N-Ntao%求得了作为系数的a(jn+1)和bi
for m=1:Dim
a(mN-Ntao+1-n)=(
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 3311 2015-05-02 00:07 With_Delay.m
评论
共有 条评论