5pt-6pt-Li-Hartley求本质矩阵MATLAB代码