时滞系统matlab模板.m文件

clear all;
close all;

 N = 1000;
 ts = 0.01;

 c1=2;        %控制器（SEF）参数 
 x1_1=0;       %被控对象状态变量
 
 Ts=1;                     %延迟时间
 ku=Ts/ts;
 u1_1=zeros（1ku）;        %延迟控制信号


for i =1:N
    t（i） = i*ts;
    
    xd（i） =1;%sin（t（i））;             %参考信号
    
    
    f1（i） = -x1_1;                               %被控对象 
    dx1（i） = f1（i） +  u1_1（1）;
    
    x1（i） = x1_1 + ts*dx1（i）;
    x1_1 = x1（i）;%对象状态更新                                   
 
  
    e1（i） = x1_1 - xd（i）;                 
    u（i）= - c1*e1（i）;
    
%     延迟