run_omp.m

%============================================================
% Paulo R. A. Candido Jr.
% Estimar parametros de localizacao (x,y) via OMP
%============================================================
close all; clc; startup1; tic;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Modelo do sinal
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
UPA = upa(Mx, Mz, ...
    d_x, d_z, elev, lambda, plt_array);
    
[~, ~, Y] = signals_los(UEs, UPA, lambda, L, alpha, ...
    SNR_dB, P_tx, Mx, Mz);

ref = UPA(1,:);  % coordenadas da primeira antena

responsearray_UPA = @(x, y, z) upa_responsearray(Mx, Mz, ...
    elev, d_x, d_z, lambda, x, y, z);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Algoritmo Ortogonal Matching Pursuit
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
K = 3; % máximo de átomos (paths) a selecionar 
yl = Y(:,1);
A = dicionario(x_grid, y_grid, responsearray_UPA);

[est_omp, pseudospec] = omp_2D( ...
    yl, x_grid, y_grid, K, A);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Resultados e analises graficas
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
printOMPResults;