w = 1000*[0.2:.2:3,3.5:1:10];                              %頻率
p = [-1.5728 -1.5747 -1.5772 -1.5790 -1.5816 -1.5852 -1.5877 -1.5922 ...
     -1.5976 -1.6044 -1.6129 -1.6269 -1.6512 -1.6998 -1.8621 1.7252 ...
      1.6124  1.5930  1.5916  1.5708  1.5708  1.5708 ]; %相位
  figure(1);                                            %開一個新的畫圖窗口
  semilogx(w,p,'ro-');                                    %用半對數坐標畫原始數據
  hold on;                                             %保持原來的圖
  semilogx(w,unwrap(p),'b*-');                            %用半對數坐標畫校正后數據
  xlabel('頻率');                                        %添加橫坐標說明
  ylabel('相位');                                        %添加縱坐標說明
  legend('原始相頻特性曲線','矯正后的相頻特性曲線');      %添加圖例