clear?all;
close?all;
ts=1;

%Plant
sys1=tf（[7.5e005][0.0451.531]‘inputdelay‘0.03）;
dsys1=c2d（sys1ts‘zoh‘）;
[numden]=tfdata（dsys1‘v‘）;
u_1=0.0;u_2=0.0;u_3=0.0;u_4=0.0;u_5=0.0;
y_1=0.0y_2=0.0;
error_1=0.0;error_2=0.0;error_3=0.0;
ei=0;
for?k=1:1:100
time（k）=k*ts;?
rin（k）=sign（sin（0.02*2*pi*k*ts））;??%Tracing?Step?Signal

yo