Matlab EVI模型解算

`function result=eivMain(x,y) % 计算x的大小 sizeX=size(x); n=sizeX(:,1); % 生成观测矩阵 a=[ones(n,1) x]; % 构建增广矩阵 n=[a y]; % 计算增广矩阵特征值 [~,a2]=eig(n'*n); diagA2=diag(a2); % 获取最小特征值 minA2=min(diagA2); % 计算求解的中间变量 sizeMinV=size(a); i=eye(sizeMinV(:,2)); i(1,1)=0; temp=a'*a-minA2*i; % 解算待估参数 result=temp\a'*y;`