最小拍控制器的Matlab仿真

clear all
clc
T0=10;t0=3;ts=1.5;
 n=1;d=[T01];
G=tf（nd‘inputdelay‘t0）;%（控制对象频域传递函数）
HG=c2d（G1.5‘zoh‘）;%广义对象（包括0阶保持器）
Ge=filt（[100-1][1]）;%理想误差传递函数
% Gc=1-Ge;%闭环传递函数
% Dz=（1-Ge）/（Ge*HG）;%最小拍控制器
num=[1 -0.9048];den=[0.09516 0 0 - 0.