%--ART2算法迭代格式，乘型ART算法，2005/9/19
P=reshape(p,1,I);
P=P';
F=ones(676,1)*0.2;
for k=0:J
     i_k=mod(k,I)+1;
    for j=1:676
        if W(i_k,:)*F~=0
         F(j)=F(j)*(P(i_k)/(W(i_k,:)*F))^(W(i_k,j)*7);%松弛因子在6—8之間合適
     end
    end
end
%--濾波
for j=1:676
    if F(j)<0
       F(j)=0;
    end
end
%--重建圖像
 F=reshape(F,26,26);
 F=F';
x=-0.5:0.04:0.5;
y=-0.5:0.04:0.5;
[X,Y]=meshgrid(x,y);
figure;
mesh(X,Y,F)