4阶龙格库塔方法离散化Mackey_Glass时间序列

龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法，用于数值求解微分方程。由于此算法精度高，采取措施对误差进行抑制，所以其实现原理也较复杂。4阶龙格库塔方法离散化Mackey_Glass时间序列,

% 4阶龙格库塔方法离散化Mackey_Glass时间序列
function [xt]=Mackey_Glass（T）; 
N=6000/T; 
t=zeros（N1）; 
x=zeros（N1）; 
a=0.2;b=0.1; 
x（1）=1.2; 
t（1）=0; 
h=T;
tau=17; 

for k=1:N-1; 
   t（k+1）=t（k）+h; 
   if t（k）       k1=-b*x（k）; 
       k2=-b*（x（k）+h*k1/2）; 
       k3=-b*（x（k）+k2*h/2）; 
       k4=-b*（x（k）+k3*h）;