求教关于多输入最小二乘程序的问题

afei8233

想要对自己的模型进行参数辨识。试着用了下该论坛里的这个程序——
最小二乘法：
fid=fopen('input.txt','r+');
U=fscanf(fid,'%f',a);
fclose(fid);
fid1=fopen('output1.txt','r+');
Y=fscanf(fid1,'%f',a);
fclose(fid1);
U1=U(1:a-1);
U2=U(1:a-2);
Y1=Y(1:a-1);
Y2=Y(1:a-2);
D1=-1*[0;Y1];
D2=-1*[0;0;Y2];
D3=[0;U1];
D4=[0;0;U2];
D=[D1 D2 D3 D4];
Q=inv(D'*D)*D'*Y

好像对多输入但输出的问题不行，举个例子说：y=a1*u1+a2*u2+e，取得3组数据为：y(k)=1,1,3；u1(k)=1,2,1；u2(k)=0,-1,2。求参数a1,a2的最小二乘估计？？怎么写这个程序，请指点以下，谢谢！！

[ 本帖最后由 eight 于 2008-4-14 21:11 编辑 ]

dingd

用1stOpt对这类问题很方便。