论坛 › Matlab讨论区 › 求助怎样再加一个大循环？

prince00840 发表于 2009-7-7 14:35

求助怎样再加一个大循环？

clear
for k=1:7
y=;%测量值
A=[1 -1 1 0 0 0 0
    0 1 -1 -1 0 0 0
    0 0 01 -1 0 0
    0 0 0 01 -1 -1];
for i=1:7
    Q(i,i)=0.01*y(i);% 方差协方差矩阵
end
       e= zeros(7,1);
e(k,1) = 1;%e为单位列向量第i个为1，其余为0
r=A*inv(y(k)); %残差
V=A*Q*A';%r的协方差
f=A*e ;   
T=((f'*inv(V)*r).^2)/(f'*inv(V)*f);
Tc=0.248445
    if T(k)>Tc
      disp(k)
   b=inv(f'*inv(V)*f)*(f'*inv(V)*r);         
      y=y-b*e;   %最后一步循环回去不会，怎样把这个新得到y返回到前面 进行新的循环？是不是要增加成二维，请高手帮帮忙
         Ax=-b*e;
          else               
    end
end

dtczhl 发表于 2009-7-7 16:46

y=;%测量值
放到for k=1:7之前不就行了吗？？
不知道楼主是不是这个意思

波波球 发表于 2009-7-7 20:48

把后面一块做成函数,返回值为y可否?

查看完整版本: 求助怎样再加一个大循环？