广义最小二乘算法程序
:@) 下面这个程序是关于广义最小二乘算法的，希望大家共同讨论，多多提出改进意见。%之前先运行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=;
H2=Y;
H3=-1*;
H4=-1*U;
H=;
e=+H*Q;
e1=e(1:a-1);
e2=e(1:a-2);
E1=-1*;
E2=-1*;
E=;
F=inv(E'*E)*E'*e;
%使用估计参数对y、u进行滤波%
Y1=Y(1:a-1);
Y2=Y(1:a-2);
Z1=;
Z2=;
Z=;
Y=Y+Z*F;
U1=U(1:a-1);
U2=U(1:a-2);
Z1=;
Z2=;
Z=;
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*;
D2=-1*;
D3=;
D4=;
D=;
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 编辑 ] 没看懂的说哦
先看到了程序，现在去找书籍，谢谢了
页:
[1]