lsm.m

关于系统辨识的一些仿真程序

num=[1 6.5 14 11.5 3];
den=[1 20 120 232 320];
u=mseries(16);% to generate a signal of m series
[z,A,B]=getval(num,den); % to function as a way to get z,A,and B.
x1=length(A);% to get the length of A.
x2=length(B);% to get the length of B.
HL=gmatrix(z,x1,x2);% to get the matrix of HL.
ZL=getzl(z,x1);% to get the matrix ZL.
c1=HL'*HL;
c2=inv(c1);
c3=HL'*ZL;
c=c2*c3;% to calculate parameters.
rz=getrz(c,x1,x2,u) % to get the real values using c mentioned above.
dif=rz-z;% to get the difference between the real values and the ideal values.
i=1:1:16;
plot(i,dif,'red')% to depict the curve of the difference,in red.