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

转载于:https://www.cnblogs.com/yangxiao99/p/4751245.html

matlab练习程序(PCASVD)相关推荐

  1. matlab多元回归程序,多元回归程序MATLAB程序

    <多元回归程序MATLAB程序>由会员分享,可在线阅读,更多相关<多元回归程序MATLAB程序(45页珍藏版)>请在人人文库网上搜索. 1.程序MATLAB多元回归程序matl ...

  2. 无源定位之时差估计的精确时差估计算法(ETDE)及MATLAB实现程序

    精确时差估计算法(ETDE)及MATLAB实现程序 算法原理 算法总结 性能分析 实验结果 算法原理 假设两接收站分别接收的带噪信号为 {x(kT)=s(kT)+ε1(kT)y(kT)=s(kT−D) ...

  3. matlab潮流程序,IEEE33节点matlab潮流程序.doc

    IEEE33节点matlab潮流程序 Bus [1,0,0 ; 2,100, 60; 3,90,40; 4,120,80; 5,60,30; 6,60,20; 7, 200, 100 ; 8,200, ...

  4. matlab有意思程序,matlab有意思的小程序

    10个C++趣味小程序,很有意思的.VIP专享文档 VIP专享文档是百度文库认... 现在很多人使用微信的时间已经非常长了,他们注册的微信号往上可能已经是5年前的事情了,正是由于不少使用者在这个过程当 ...

  5. matlab程序里nbus=x的意思,MATLAB潮流程序(IEEE14 直角坐标 牛拉法)

    MATLAB潮流程序(IEEE14 直角坐标牛拉法)clear baseMVA=100; %功率基值 %%读Data1中数据 load Data1.txt Bus=Data1(:,1); %节点号 V ...

  6. 怎样在matlab q-q图上读出斜率,Q分解法潮流计算matlab小程序

    <Q分解法潮流计算matlab小程序>由会员分享,可在线阅读,更多相关<Q分解法潮流计算matlab小程序(18页珍藏版)>请在人人文库网上搜索. 1.Q分解法潮流计算matl ...

  7. matlab期权风险评估算法,使用 MATLAB 应用程序根据期权价格估算风险中性密度 (risk-neutral density, RND)...

    我们为什么选择 MATLAB 我们决定创建 MATLAB 应用程序估算 RND,因为我们希望研究成果得到广泛应用,而不是仅限于在量化经济学领域或是熟悉 MATLAB 编程的和金融研究人员使用. 我们发 ...

  8. matlab密度特征,关于用MATLAB应用程序估算风险中性密度的性能分析和介绍

    由于场外期权合约的买卖在交易双方间私下进行而非通过公开市场,因而可能很难确定合约的价格有利于买方还是卖方.为对这些合约进行定价,金融分析师往往依据看涨期权或看跌期权价格估算出风险中性密度(RND)值. ...

  9. 10 竞争神经网络与SOM神经网络matlab参考程序

    %% I. 清空环境变量 clear all clc %% II. 训练集/测试集产生 %% % 1. 导入数据 load water_data.mat %% % 2. 数据归一化 attribute ...

  10. 多自由度有阻尼matlab,多自由度阻尼系统固有振型的MATLAB求解程序

    多自由度阻尼系统固有振型的MATLAB求解程序function [v,w,zeta]=vbr_sf(m,d,k) %vbr_sf  vbr_sf(m,d,k) %        [v,w,zeta]= ...

最新文章

  1. 使.NET使用数据库连接池 .
  2. 测视力距离5米还是3米_视力表只能看到0.3,相当于近视度数200度?
  3. PHP中define定义的常量如何在smarty模板调用
  4. SQL Server 聚合函数 (方差和标准差)
  5. 服务器内存怎么看型号,怎么查看服务器内存的型号
  6. 铁道部正在研究高铁月票 重申上座率120%准确可靠
  7. 阿米洛键盘失灵_改变静电容键盘手感单一限制,阿米洛静电容机械轴V2上手
  8. 依图在实时音视频中语音处理的挑战丨RTC Dev Meetup
  9. JavaScript学习攻略
  10. android开发笔记之锁屏界面未读短信未接来电提醒(android 4.4)
  11. 宇称不守恒说了啥?杨振宁和李政道的发现有多大意义
  12. POI之图片导出到PPT简单实例
  13. Linux下xmms
  14. UDS安全访问之27服务SecurityAccess
  15. U盘目录穿越获取车机SHELL(内含模拟环境可上手实践)
  16. OKEx调整上币规则,波多野结衣出席“AVH”发布会 | 区块链日报
  17. Dell XPS13 9350 安装win7,遇到的坑
  18. 有限差分法在地震波二维正演模拟中的应用
  19. 解决新能源充电难!小桔充电牵手南网电动,共建充电场站
  20. Linux scp命令 | 还在为两个 ip 服务器之间如何传输项目 发愁吗

热门文章

  1. 程序员才能看懂,看到第18张终于忍不住笑喷了。
  2. Linux中Shell的命令替换用法笔记
  3. 2017已过半,这半年,你累吗?
  4. 多态的概念、对象上下转型、多态的应用、异常(异常概念、异常分类、java异常处理机制、try...catch...finally、throw和throws、自定义异常)
  5. Linux系统文件编程(1)
  6. sql多层嵌套别名无效_SQL之复杂查询
  7. 怎样开发一个 Node.js 命令行工具包
  8. 你应该知道的 Node 基础知识
  9. Realtime Ray Tracing RenderMan Point Cloud
  10. SQL Server CONVERT() 日期转换为新数据类型的 通用函数