function FAR_model = get_far_model(N, M, epi)
    FAR_model = zeros(N,M*N);
    %Cn = randperm(M)-1
    for n = 0 : N-1
        Cn = floor(rand()*M);
        for q = 0 : N-1
            for p = 0:M-1
                FAR_model(n+1,q*M+p+1) = exp(1i*2*pi*p/M*Cn+1i*2*pi*q/N*n*(1+Cn*epi));
            end
        end
    end
end