matlab_code of reconstruction