人脸分割 matlab程序

人脸脸部的皮肤具有一定的规律（都落在某个YCbCr表示的一个椭圆区域内），符合此规律的辨识为人脸

程序：（自己写， 但结果颜色很奇怪）

function [] = detection()
clear;
close all;
l = imread('2-2.jpg');
figure(1);
imshow(l);
l=double(l);
k=rgb2ycbcr(l);
imshow(k);
f=zeros(267,400);


for i= 1:267
    for j=1:400
%Y=0.257*l(i,j,1)+0.504*l(i,j,2)+0.098*l(i,j,3);
%Cb=0.148*l(i,j,1)-0.291*l(i,j,2)+0.439*l(i,j,3);
%Cr=0.439*l(i,j,1)-0.368*l(i,j,2)-0.071*l(i,j,3);
Y=k(i,j,1);
Cb=k(i,j,2);
Cr=k(i,j,3);
%Wci(Y)
if(Y<125)
Wcb =23+(Y-16)*(46.97-23)/(125-16);
Wcr =20+(Y-16)*(38.76-20)/(125-16); 
aCb =108+(125-Y)*(118-108)/(125-16);%aCi(Y)(Ci(Y)平均数）
aCr =154-(125-Y)*(154-144)/(125-16);
end
if(Y>188)
Wcb &#