人脸识别.rar_人脸识别_人脸识别 matlab

clc; clear; phi = zeros(4731,70); fileform = 'D:\湖南大学研究生\研一\课程\图形图像处理作业\人脸识别\训练集GT\*.jpg'; filepathsrc = 'D:\湖南大学研究生\研一\课程\图形图像处理作业\人脸识别\训练集GT\'; file = dir(fileform); for i = 1:length(file) image = imread([filepathsrc,file(i).name]); image = rgb2gray(image); image = imresize(image,[57,83]); phi(:,i) = reshape(image,4731,1); end lable(1:10,1) =1; lable(11:20,1) =2; lable(21:30,1) =3; lable(31:40,1) =4; lable(41:50,1) =5; lable(51:60,1) =6; lable(61:70,1) =7; I = ones(70,70); n = 4731; m = round(0.1*n); q = sqrt(1/m) * randn(m,n); phi = q * phi; fileform = 'D:\湖南大学研究生\研一\课程\图形图像处理作业\人脸识别\测试集GT\*.jpg'; filepathsrc = 'D:\湖南大学研究生\研一\课程\图形图像处理作业\人脸识别\测试集GT\'; file = dir(fileform); count = 0; for i = 1:length(file) x = imread([filepathsrc,file(i).name]); x = rgb2gray(x); x = imresize(x,[57,83]); x = im2double(x); x = reshape(x,4731,1); X = q * x; alph = 0.0001; p = pinv(phi'* phi + alph * I)* phi'* X; r = zeros(1,7); for j = 0:6 r(j+1) = norm(X -phi(:,j*10+1:j*10+10)* p(j*10+1:j*10+10,1))/ norm(p(j*10+1:j*10+10,1)); end min = r(1); num=7; for j = 1:7 if r(j) < min min = r(j); num = j; end end if i > (num-1)*10 && i <= num*10 fprintf(['该图像分类正确\n']); count = count+1; else fprintf(['该图像分类不正确\n']); end end k = count/70