%輸入
c1=input（‘請(qǐng)輸入非線性特性c1：‘）;
K1=input（‘請(qǐng)輸入非線性部分的飽和特性K1：‘）;
k=K1/c1;%計(jì)算非線性部分的斜率
A=input（‘請(qǐng)輸入矩陣A：‘）;
B=input（‘請(qǐng)輸入矩陣B：‘）;
C=input（‘請(qǐng)輸入矩陣C：‘）;
T=input（‘請(qǐng)輸入采樣周期T：‘）;

%計(jì)算度Φ、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））;

%輸入初始狀態(tài)和輸入
while?true
????R=input（‘請(qǐng)輸入輸入信號(hào)的幅度R：‘）;
????X=input（‘請(qǐng)輸入初始狀態(tài)X（0）：‘）;
????if?