9.2.3 MATLAB代码实现