[发明专利]一种灰度值矩阵快速分解重建的正则化方法

摘要

本发明公开了一种灰度值矩阵快速分解重建的正则化方法，包含如下步骤：(1)根据被测场域获取重建所需的边界测量电压b和灵敏度矩阵A；(2)设置初始化参数；(3)更新辅助变量pk+1；(4)更新辅助变量qk+1；(5)更新变量wk+1；(6)更新变量vk+1；(7)更新灰度值gk+1；(8)设置优化参数；(9)更新灰度值uk+1；(10)更新变量κk+1；(11)更新变量mk+1；(12)判断迭代是否符合迭代终止条件uk+1是否收敛，若是则迭代终止，进行下一步操作；若否，设置k＝k+1并跳回第(9)步，继续迭代求解；(13)根据最终求解所得灰度值进行成像。本发明能够有效地降低阶梯伪影、提升背景清晰及提高抗噪性。