function MPC2()
num1=5.72;den1=[60 1];g1=poly2tfd(num1,den1,0,14);
num2=1.52;den2=[25 1];g2=poly2tfd(num2,den2,0,15);
%將MPC傳遞函數模型轉化為MPC狀態空間模型
mod=tfd2mod(7,1,g1,g2);
%將MPC狀態空間模型轉化為通用狀態空間模型
[A,B,C,D]=mod2ss(mod);
%將通用狀態空間模型轉換LTI對象的狀態空間模型
sys=ss(A,B,C,D);
h=tf(sys);
%獲得脈沖信號x
[u,t]=gensig('pulse',8,10,0.1);
x=[u u];
%求解LTI對象的單位脈沖響應y
t=0:0.1:10;
[y,x1,t1]=lsim(h,x,t);
%輸入脈沖信號x的歸一化處理
[ax,mx,stdx]=autosc(x);
mx=[0 0];
sx=scal(x,mx,stdx);
%生成用于線性回歸計算的輸入/輸出數據矩陣
n=35;
[xreg,yreg]=wrtreg(sx,y,n);
%基于多變量最小二乘法的脈沖響應模型辨識
ninput=2;
plotopt=2;
[theta,yres]=mlr(xreg,yreg,ninput,plotopt)