广义最小二乘算法程序

chengtianmmz

:@) 下面这个程序是关于广义最小二乘算法的，希望大家共同讨论，多多提出改进意见。

%之前先运行input_output和LS程序%
fid=fopen('input.txt','r+');
U=fscanf(fid,'%f',a);
fclose(fid);
fid1=fopen('output1.txt','r+');
Y=fscanf(fid1,'%f',a);
fclose(fid1);
for i=1:a;
%计算系统模型的残差及参数%
y1=Y(3:a);
y2=Y(2:a);
u1=U(2:a);
H1=[y2;0];
H2=Y;
H3=-1*[u1;0];
H4=-1*U;
H=[H1 H2 H3 H4];
e=[y1;0;0]+H*Q;
e1=e(1:a-1);
e2=e(1:a-2);
E1=-1*[0;e1];
E2=-1*[0;0;e2];
E=[E1 E2];
F=inv(E'*E)*E'*e;
%使用估计参数对y、u进行滤波%
Y1=Y(1:a-1);
Y2=Y(1:a-2);
Z1=[0;Y1];
Z2=[0;0;Y2];
Z=[Z1 Z2];
Y=Y+Z*F;
U1=U(1:a-1);
U2=U(1:a-2);
Z1=[0;U1];
Z2=[0;0;U2];
Z=[Z1 Z2];
U=U+Z*F;
Y3=Y(1:a-1);
Y4=Y(1:a-2);
U3=U(1:a-1);
U4=U(1:a-2);
%根据滤波后的输入输出数据再次进行最小二乘参数估计，求Q%
D1=-1*[0;Y3];
D2=-1*[0;0;Y4];
D3=[0;U3];
D4=[0;0;U4];
D=[D1 D2 D3 D4];
QQ=inv(D'*D)*D'*Y;
q=norm(QQ-Q)/norm(Q);
Q=QQ;
if q<10e-14
    i
    q
    F
    Q
    break
end
end
ppp='help me'

[ 本帖最后由 eight 于 2007-6-20 15:24 编辑 ]

ttloseway

没看懂的说哦

ylt9811115

先看到了程序，现在去找书籍，谢谢了