完成公式和MATLAB程序