a=【1 -1.5 0.7】';b=【1 0.5】';%待辨识参数
d=3;na=2;nb=1;N=100;%系统类型
uk=zeros(d+nb,1);%定义输入
yk=zeros(na,1);%定义输出
kxi=sqrt(1)*randn(N,1);%随机噪声
theta=【a(2:na+1);b】;
%% 产生M序列,作为系统的输出
x1=1;x2=1;x3=1;x4=0;S=1;
for k=1:N
M(k)=xor(x3,x4);
IM(k)=xor(M(k),S);
if IM(k)==0
u(k)=-1;
else
u(k)=1;
end
S=not(S);
x4=x3;x3=x2;x2=x1;x1=M(k);
end
for k=1:N
Phi(k,:)=【-yk;uk(d:d+nb)】';
y(k)=Phi(k,:)*theta+kxi(k);
for i=d+nb:-1:2
uk(i)=uk(i-1);
end
uk(1)=u(k);
for i=na:-1:2
yk(i)=yk(i-1);
end
yk(1)=y(k);
end
thetae=inv(Phi'*Phi)*Phi'*y';
参考书:系统辨识理论及MATLAB仿真-刘金琨