matlab练习程序(PCA<SVD>)
clear all;
close all;
clc;
img1=imread('Corner.png');
img2=imread('Corner1.png');
img3=imread('Corner2.png');
img1=imresize(img1,[35 90]); %Matlab的svd不支持太大的数据,我把数据变小了.
img2=imresize(img2,[35 90]);
img3=imresize(img3,[35 90]);
[height width]=size(img2);
a=reshape(double(img1),1,35*90);
b=reshape(double(img2),1,35*90);
c=reshape(double(img3),1,35*90);
ma=mean(a);
mb=mean(b);
mc=mean(c);
a=a-ma;
b=b-mb;
c=c-mc;
x=[a' b' c'];
x=double(x);
%y=x'*x/(height*width);
%[v d]=eig(y);
[u s w]=svd(x); %w相当于排序后的d,用svd就不用求x'*x了
ed=x*w;%v;
img4=zeros(height,width);
img4=reshape(ed(:,1),height,width);
img5=zeros(height,width);
img5=reshape(ed(:,2),height,width);
img6=zeros(height,width);
img6=reshape(ed(:,3),height,width);
imshow(mat2gray(img4));
figure,imshow(mat2gray(img5));
figure,imshow(mat2gray(img6));
参考:
1.http://iiec.cqu.edu.cn/wiki/index.php/SVD%E4%B8%8EPCA%E7%9A%84%E7%93%9C%E8%91%9B
最新文章
- PL/SQL通过免安装客户端连接远端ORACLE数据库
- libev安装与示例程序编译运行
- asp.net判断FileUpload选择的文件是否是图片
- 弹出框四 之toastr.js (完成提示框)
- 情定XMLA,割舍不下的XAML
- bzoj4137 [FJOI2015]火星商店问题
- ural 1114,计数dp
- Android下各个按键对应的key code
- iOS事件机制(二)
- SRM 441(1-250pt, 1-500pt)
- jquery prop()方法 解决全选 不全选 反选 问题 解决执行一次不不能再执行问题
- 【noip 2009】 乌龟棋 记忆化搜索&;动规
- 很污的机器学习:从xhamster网站找到喜欢的片子
- jQuery === 面条式代码?
- 前端的UI设计与交互之布局篇
- Error response from daemon: conflict: unable to remove repository reference 解决方案
- pta编程总结
- MyBatisPlus环境下使用MyBatis的配置类
- window10安装黑苹果
- 初始化Weex项目遇到的问题记录