如何给系统加脉冲

大江东去

下面是一个分数阶吸引子，我想在每个方程上乘以一个脉冲函数不知道怎么办，求高手指教：  function [T, Y]=FOChen(parameters, orders, TSim, Y0) % time step: h=0.005; % number of calculated mesh points: n=round(TSim/h); %orders of derivatives, respectively: q1=orders(1); q2=orders(2); q3=orders(3); % constants of Chen's system: b=parameters(1); % binomial coefficients calculation: cp1=1; cp2=1; cp3=1; for j=1:n      c1(j)=(1-(1+q1)/j)*cp1;      c2(j)=(1-(1+q2)/j)*cp2;      c3(j)=(1-(1+q3)/j)*cp3;      cp1=c1(j); cp2=c2(j); cp3=c3(j); end % initial conditions setting:   x(1)=Y0(1); y(1)=Y0(2); z(1)=Y0(3); % calculation of phase portraits /numerical solution/: for i=2:n      x(i)=(35*(y(i-1)-x(i-1))（此处乘以脉冲函数）)*h^q1 - memo(x, c1, i);      y(i)=(-x(i)*z(i-1)-7*x(i)+28*y(i-1))*h^q2 - memo(y, c2, i);      z(i)=(x(i)*y(i)-b*z(i-1))*h^q3 - memo(z, c3, i); end for j=1:n      Y(j,1)=x(j);      Y(j,2)=y(j);      Y(j,3)=z(j); end T=h:h:TSim; [t,y]=FOchen([3], [0.9 0.9 0.9], 100, [2 -1  6]); figure plot3(y(:,1), y(:,3), y(:,2), 'k'); xlabel('z(t)'); ylabel('x(t)'); zlabel('y(t)'); grid; title('chenxiyinzi x,y,z') 下面是脉冲函数：将此函数嵌入到程序的汉字处，谢谢》 function f=f(T) T=0.005:0.005:10; s=0.1; w=2/s; f=square(w*T); plot(T,f);