function f=spectremat(A,R);

%           f=spectremat(A,R);
%
%Calcule, étant donné un tableau de signaux lignes 
%             A=[:,-R:R/M:R-R/M] 
%le tableau où chaque ligne a été remplacée par son spectre approché aux
%points omega=k pi/T, k=-M:1:M-1. 


[M,MM]=size(A);
MMM=fix(MM/2);
w=-ones(M,MM);
for k=1:2:MM
  w(:,k)=-w(:,k);
end
AA=A';
UU=fft(AA);
U=UU';
U=w.*U;
U1=U(:,1:MMM);
U2=U(:,MMM+1:MM);
U=[U2,U1];
U=(2*R/MM)*U;
f=U;