%輸入
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?