理论与试验建模第三章.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 2、最小二乘法求解 a,b * function [num,den]=lsqident(u,y,m,n) %u，y输入输出数列，m，n输入输出数列的阶次 M=length(u); % 获得观测序列u,y的长度 Y=y(n+1:M); for i=1: length(Y) A(i,:)=[-y(n+i-1:-1:i),u(m+i:-1:i)]; End % 建立A,Y阵,A*X=Y,X未知 X= A\Y % A 左除 Y num=X(n+1:length(X)) % 返回差分方程系数 bi den=[1 X(1:n)]