离散相似法仿真控制系统

%输入
c1=input（‘请输入非线性特性c1：‘）;
K1=input（‘请输入非线性部分的饱和特性K1：‘）;
k=K1/c1;%计算非线性部分的斜率
A=input（‘请输入矩阵A：‘）;
B=input（‘请输入矩阵B：‘）;
C=input（‘请输入矩阵C：‘）;
T=input（‘请输入采样周期T：‘）;

%计算度Φ、G
syms t;
%EAT=simplify（expm（（T-t）*A））;
GT1=double（int（expm（（T-t）*A）*Bt0T））;
%GT=GT1*B;
Fai=simplify（expm（T*A））;

%输入初始状态和输入
while true
    R=input（‘请输入输入信号的幅度R：‘）;
    X=input（‘请输入初始状态X（0）：‘）;
    if