【人脸识别】基于PCA实现ORL人脸识别附matlab代码和报告

1 简介

人脸识别技术先进,应用广泛。借助PCA算法,利用MATLAB GUI可以简单操作,通过对待识别图像的预处理即可提高识别率。本文首先对相关概念进行了阐述,对工作原理进行了介绍,具体对基于PCA算法人脸识别的MATLAB实现进行了解析。

2 部分代码

function [neednum,average_face, immin, newVT] = newVT(imdata)%此过程为标准PCA流程%计算平均脸并显示average_face=mean(imdata,2); %按行求平均mean(a,2)  按列mean(a)Average_face=reshape(average_face,112,92);%将[112*92,1]的脸灰度数据转成[112,92]figure;subplot(1,1,1);imshow(Average_face,[]);%imshow(I,[]) 显示灰度图像 I,根据 I 中的像素值范围对显示进行转换。title(strcat('40*5张训练样本的平均脸'));%步骤一:图像预处理,去中心化,即减去平均值train_num=5;immin=zeros(112*92,40*train_num);for i=1:40*train_num      immin(:,i) = imdata(:,i)-average_face;end%步骤二:计算协方差矩阵W=immin'*immin; %n*d x d*n= n*n 较小%步骤三:对协方差矩阵W进行特征值分解,计算特征向量与特征值(向量)[V,D]=eig(W);%步骤四:对特征值进行排序,选取贡献率大于前85%的特征脸[D_sort,index] = sort(diag(D),'descend');%降序排列SumAllFaceEigenValue=sum(D_sort);NowFaceEigenValue=0;    %选取累计贡献大于85%的前neednum个特征脸for i=1:size(D_sort,1)NowFaceEigenValue=NowFaceEigenValue+D_sort(i);neednum=i;if(NowFaceEigenValue>SumAllFaceEigenValue*0.85)%累计贡献率达到85%以上即可    break;endendV_sort = V(:,index);VT=immin*V_sort; %dxn*nxk=d*kfor i=1:40*train_num    VT(:,i)=VT(:,i)/norm(VT(:,i));%归一化处理end% need_num=neednum;for i=1:neednumnewVT(i,:)=VT(:,i);%取前neednum个特征值endend

3 仿真结果

4 参考文献

[1]刘向东. 基于PCA算法人脸识别的MATLAB实现[J]. 电脑知识与技术:学术版, 2016(4X):2.

博主简介:擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真,相关matlab代码问题可私信交流。

部分理论引用网络文献,若有侵权联系博主删除。

  • 0
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是基于PCA和BP神经网络的MATLAB人脸识别代码实现,其中包括数据预处理、特征提取、模型训练和测试等步骤。 1. 数据预处理 首先需要准备训练数据和测试数据。数据集可以使用公开的人脸数据库,如Yale人脸数据库、ORL人脸数据库等。这里以Yale人脸数据库为例,该数据库包含15个人的165张灰度图像,每个人有11张不同表情的图像。代码如下: ```matlab clear all; clc; % 读取数据 dataDir = 'yalefaces'; imgList = dir(fullfile(dataDir,'*.*')); imgNum = length(imgList); imgSize = [243, 320]; % 图像大小 imgData = zeros(imgSize(1)*imgSize(2), imgNum); for i = 1:imgNum img = imread(fullfile(dataDir, imgList(i).name)); img = imresize(img, imgSize); imgData(:,i) = img(:); end % 数据归一化 imgData = double(imgData); imgData = imgData - mean(imgData, 2); % 减去均值 imgData = imgData ./ std(imgData, 0, 2); % 归一化 ``` 2. 特征提取 接下来,使用PCA方法对数据进行降维,提取出最重要的特征。代码如下: ```matlab % PCA降维 [U,S,V] = svd(imgData, 'econ'); eigVals = diag(S).^2; energy = cumsum(eigVals) / sum(eigVals); thres = find(energy >= 0.99, 1); U = U(:,1:thres); feaData = U.' * imgData; ``` 3. 模型训练 使用BP神经网络对特征进行分类。首先,将数据集分为训练集和测试集,代码如下: ```matlab % 数据集分割 trainNum = 10; % 每个人的训练样本数 testNum = 11 - trainNum; % 每个人的测试样本数 trainData = zeros(size(feaData,1), trainNum*15); trainLabel = zeros(15, trainNum*15); testData = zeros(size(feaData,1), testNum*15); testLabel = zeros(15, testNum*15); for i = 1:15 idx = (i-1)*11+1:i*11; trainData(:,(i-1)*trainNum+1:i*trainNum) = feaData(:,idx(1:trainNum)); trainLabel(i,(i-1)*trainNum+1:i*trainNum) = 1; testData(:,(i-1)*testNum+1:i*testNum) = feaData(:,idx(trainNum+1:end)); testLabel(i,(i-1)*testNum+1:i*testNum) = 1; end ``` 然后,搭建BP神经网络模型并进行训练。代码如下: ```matlab % BP神经网络训练 net = feedforwardnet([20,10]); net.trainFcn = 'trainlm'; net.trainParam.show = 50; net.trainParam.epochs = 1000; net.trainParam.goal = 1e-5; net.trainParam.lr = 0.01; [net, tr] = train(net, trainData, trainLabel); ``` 4. 模型测试 最后,使用测试数据对模型进行测试,并计算识别准确率。代码如下: ```matlab % BP神经网络测试 testOutput = net(testData); [~, testPred] = max(testOutput); [~, testTarget] = max(testLabel); accuracy = sum(testPred == testTarget) / length(testTarget); fprintf('Accuracy: %.2f%%\n', accuracy*100); ``` 完整代码如下: ```matlab clear all; clc; % 读取数据 dataDir = 'yalefaces'; imgList = dir(fullfile(dataDir,'*.*')); imgNum = length(imgList); imgSize = [243, 320]; % 图像大小 imgData = zeros(imgSize(1)*imgSize(2), imgNum); for i = 1:imgNum img = imread(fullfile(dataDir, imgList(i).name)); img = imresize(img, imgSize); imgData(:,i) = img(:); end % 数据归一化 imgData = double(imgData); imgData = imgData - mean(imgData, 2); % 减去均值 imgData = imgData ./ std(imgData, 0, 2); % 归一化 % PCA降维 [U,S,V] = svd(imgData, 'econ'); eigVals = diag(S).^2; energy = cumsum(eigVals) / sum(eigVals); thres = find(energy >= 0.99, 1); U = U(:,1:thres); feaData = U.' * imgData; % 数据集分割 trainNum = 10; % 每个人的训练样本数 testNum = 11 - trainNum; % 每个人的测试样本数 trainData = zeros(size(feaData,1), trainNum*15); trainLabel = zeros(15, trainNum*15); testData = zeros(size(feaData,1), testNum*15); testLabel = zeros(15, testNum*15); for i = 1:15 idx = (i-1)*11+1:i*11; trainData(:,(i-1)*trainNum+1:i*trainNum) = feaData(:,idx(1:trainNum)); trainLabel(i,(i-1)*trainNum+1:i*trainNum) = 1; testData(:,(i-1)*testNum+1:i*testNum) = feaData(:,idx(trainNum+1:end)); testLabel(i,(i-1)*testNum+1:i*testNum) = 1; end % BP神经网络训练 net = feedforwardnet([20,10]); net.trainFcn = 'trainlm'; net.trainParam.show = 50; net.trainParam.epochs = 1000; net.trainParam.goal = 1e-5; net.trainParam.lr = 0.01; [net, tr] = train(net, trainData, trainLabel); % BP神经网络测试 testOutput = net(testData); [~, testPred] = max(testOutput); [~, testTarget] = max(testLabel); accuracy = sum(testPred == testTarget) / length(testTarget); fprintf('Accuracy: %.2f%%\n', accuracy*100); ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

matlab科研助手

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值