[I,map]=imread('3-22.jpg');
imshow(I,map);
I=double(I);
[Gx,Gy]=gradient(I);       % 計算梯度
G=sqrt(Gx.*Gx+Gy.*Gy);   % 注意是矩陣點乘

J1=G;
figure,imshow(J1,map);    % 第一種圖像增強

J2=I;                   % 第二種圖像增強
K=find(G>=7);
J2(K)=G(K);
figure,imshow(J2,map);

J3=I;                   % 第三種圖像增強
K=find(G>=7);
J3(K)=255;
figure,imshow(J3,map);

J4=I;                   % 第四種圖像增強
K=find(G<=7);
J4(K)=255;
figure,imshow(J4,map);

J5=I;                   % 第五種圖像增強
K=find(G<=7);
J5(K)=0;
Q=find(G>=7);
J5(Q)=255;
figure,imshow(J5,map);