www.pudn.com > CTArt.rar > CTArt.m

%--ART算法迭代格式--good 
%function F=CTART(p,W,N,I,J) 
tic 
%p = IMNOISE(p,'gaussian',0,0.000001);%--SNR  30dB 
%p = IMNOISE(p,'gaussian',0,0.0001);  %--SNR  20dB 
P=p(1:I); 
F=ones(J,1).*0.4; 
%F=Z(1:J); 
%F=0.5*F'; 
 
for k=0:10*I 
     i_k=mod(k,I)+1; 
    if sum(W(i_k,:))~=0 
       F=F+0.6*(P(i_k)-W(i_k,:)*F).*W(i_k,:)'/sum(W(i_k,:).^2); 
    end 
    for j=1:J 
             if F(j)<0 
              F(j)=0; 
             end 
         end 
         
       avere(k+1)=sum(abs(Fs-F'))/(fsmax*J); 
       maxe(k+1)=max(abs(Fs-F'))/fsmax; 
       sqrte(k+1)=sqrt((Fs-F')*(Fs-F')'/fssum); 
     
end 
 
%--重建图像 
M=N; 
F=reshape(F,N,M); 
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); 
%--误差 
r=0; 
for c=1:100:50*I 
   r=r+1; 
   avere1(r)=avere(c); 
    maxe1(r)=maxe(c); 
    sqrte1(r)=sqrte(c); 
end 
%     
figure; 
plot(avere1); 
figure; 
plot(maxe1); 
figure; 
plot(sqrte1); 
toc