OMP算法Matlab代码

在MATLAB中，正交匹配追踪（OMP）算法通常用来解决线性最小二乘问题，特别是当问题中的信号是稀疏的。以下是一个简单的示例，展示如何使用MATLAB内置的`omp`函数实现OMP算法：

% 初始化变量
m = length(measurements); % 测量数据长度
n = dimension; % 理论上的信号维度
k = numNonzeros; % 真实的非零系数数量
signal = sparseRandom(n, k); % 创建一个稀疏信号
measurement_matrix = generateMeasurementMatrix(m, n); % 创建随机测量矩阵

% 使用omp函数进行正交匹配追踪
options = struct('NumIterations', max_iterations); % 设置最大迭代次数
[x_hat, support] = omp(measurement_matrix, signal, options);

% "support"就是找到的非零系数索引
% "x_hat"是重建后的信号估计

% 可视化结果
figure;
plot(signal, 'o', 'MarkerSize', 10, 'DisplayName', 'Original Signal');
hold on;
plot(x_hat, 'r', 'DisplayName', 'Recovered Signal (OMP)');
legend('Location', 'SouthEast');