预测辨识.m

本人写的预测控制仿真小程序

clc;
clear;
T=6;
g11=poly2tfd(12.8,[16.7 1],0,1);
g12=poly2tfd(6.6,[11 1],0,7);
mod=tfd2mod(T,1,g11,g12);
[A,B,C,D]=mod2ss(mod);
sys=ss(A,B,C,D);
h=tf(sys);
[u,t]=gensig('pulse',7,10,0.1);
x=[u u] ;
t=0:0.1:10;
[y,x1,t1]=lsim(h,x,t);
%[ax,mx,stdx]=autosc(x);
%mx=[0 0] ;
%sx=scal(x,mx,stdx);
n=32;
[xreg,yreg]=wrtreg(x,y,n);
ninput=2;lv=10;
plotopt=2;
%[theta,yres]=mlr(xreg,yreg,ninput,plotopt)
theta=plsr(xreg,yreg,ninput,lv,plotopt)