function?With_Delay
clear?all;?close?all;?clc;
Dim=15;q（1:Dim1）=0.9;FinalTime=1.5;tao=0.5;%Dim為方程數(shù)，q為導(dǎo)數(shù)的階數(shù)，F(xiàn)inalTime為t的最大值tao為時(shí)滯
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%求得了作為系數(shù)的a（jn+1）和bi
????for?m=1:Dim
????????a（mN-Ntao+1-n）=（

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????3311??2015-05-02?00:07??With_Delay.m