gaborcons.m

gabor algorithm in matlab

clear all;


M=input('enter M');     % M - No of translations 
N=input('enter N');     % N - No of Modulations
x=input('elements of signal are');

L=M*N;                  % L-Length of the Transform 
x=x.';

%----------------------------Gaussian Function---------------------------
%according to wexler and raz
for k=1:L
    h(k)=sqrt(sqrt(2)/N)*exp(-pi*((k-L/2)/N)^2);
end;


%----------------------------Fourier Matrix-----------------------------
w=exp((2*pi*j)/N);
for i=0:N-1
    for j=0:N-1
        W(i+1,j+1)=w^(i*j);          %verified(no normalization)
    end;
end;

E=[];
for a=0:M-1;
E=blkdiag(E,W);
end

E=(1/sqrt(N))*E;
%-----------------------------caculatin H-------------------------------

for count=1:L
    
    i=count;
    j=rem(count,N);
    if j==0
        j=N;
    end;
    for p=1:M
        if i>L
            i=i-L;
        end;
        if j>L
            j=j-L;
        end;
        
        H(i,j)=h(count);
        i=i+N;
        j=j+N;
    end;
    
           
end;

%---------------------------finally gabor-------------------------------
inE=inv(E);
inH=inv(H);
a=inv(H*E)*x;

%-----------------------------inverse gabor------------------------------
signal=H*E*a;