[zt_zt]

离散化不是问题，c2d(sys,T)就可以了。关键问题是线性化，每走一步就要线性化一次计算量太大。如果你的系统特殊就先吧线性化表示成函数就变成时变系统了。或者采用次优滤波。