[h0,h1]=meshgrid([-4:0.2:4]);   %確定濾波器系數矩陣H(n)的范圍J=0.55+h0.^2+h1.^2+2*h0.*h1*cos(pi/8)-sqrt(2)*h0*cos(pi/10)-sqrt(2)*h1*cos(pi/40*9);                                %計算各組系數下的誤差性能J(H)figure(1);surfc(h0,h1,J);                 %繪制誤差性能曲面xlabel('h0');ylabel('h1');zlabel('J(n)');figure(2);contour(h0,h1,J,100);           %繪制曲面的二維等高線xlabel('h0');ylabel('h1');