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

prince00840

clear
for k=1:7
y=[10.5 11.5 2.5 15 9.5 6.5  4.5 ];%测量值
A=[1 -1 1 0 0 0 0
    0 1 -1 -1 0 0 0
    0 0 0  1 -1 0 0
    0 0 0 0  1 -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

y=[10.5 11.5 2.5 15 9.5 6.5  4.5 ];%测量值
放到for k=1:7之前不就行了吗？？
不知道楼主是不是这个意思

波波球

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