[发明专利]基于M-稀疏算法的核磁共振回波数据反演方法与装置

权利要求书

1.一种基于M-稀疏算法的核磁共振回波数据反演方法，其特征在于，包括：

基于反演核矩阵和采集的核磁共振回波数据构建目标函数，其中，所述目标函数的解中非零元素的个数小于预设阈值，所述目标函数中不包括正则化项；

使用M-稀疏算法求解所述目标函数，获取所述核磁共振回波数据的反演谱；

所述基于反演核矩阵和采集的核磁共振回波数据构建目标函数，具体包括：

构建如下式所述的目标函数：

其中，所述f≥0表示向量f中的所有元素均大于或等于0，||f||₀表示向量中非零元素的个数，所述N为预设的常数；b为采集的核磁共振回波数据；A是反演核矩阵；

所述使用M-稀疏算法求解所述目标函数，获取所述核磁共振回波数据的反演谱，具体包括：

A、根据公式和公式确定所述Γⁱ⁺¹；

B、判断所述Γⁱ⁺¹与Γⁱ是否相等；

C、若是，则令执行D；

D、计算并判断ξ＜tol是否成立；若是则执行F，若否则执行E；

E、判断当前迭代次数是否大于预设的最大迭代次数；若是则执行F，若否令i＝i+1，返回执行A；

F、若所述fⁱ⁺¹中有负值的项，将所述fⁱ⁺¹中负值的项置0，生成新的fⁱ⁺¹，并将所述新的fⁱ⁺¹确定为所述目标函数的目标解；

其中，所述i为迭代次数，所述gⁱ＝A^T(b-Afⁱ)，所述所述为向量g只保留对应Γⁱ中元素位置的子向量，所述是矩阵A只保留对应所述Γⁱ中元素位置的列的子矩阵，所述所述是向量k中每个元素按从大到小排列后的第M个元素，k为(|k₁|,|k₂|......)，所述tol为预设阈值，f¹＝0，Γ¹＝supp(H_M(A^Tb))，所述supp(p)代表向量p的支撑集；

所述判断所述Γⁱ⁺¹与Γⁱ是否相等之后，所述方法还包括：

G、当所述Γⁱ⁺¹与所述Γⁱ不相等时，判断是否成立；若是，则令执行D；若否，则执行H；

H、根据公式μⁱ＝μⁱ/(d(1-c))更新所述μⁱ，根据公式更新所述直到所述成立为止；

I、令和返回执行所述D；

其中，所述d为常数，且0＜d＜1，所述c为常数。

2.根据权利要求1所述的方法，其特征在于，所述c为0.01。