常见测量矩阵的MATLAB实现
常见测量矩阵的MATLAB实现
下面以文献【吴赟.压缩感知测量矩阵的研究[D]. 西安电子科技大学硕士学位论文,2012】为依据,给出文献中2.2节内容所述的六种测量矩阵MATLAB实现代码,仅为一种参考实现方式,还未验证其正确性。
1、高斯矩阵
以下代码生成的高斯矩阵方差为1,若为改为1/M,只须将除以根号M即可。
- function [ Phi ] = GaussMtx( M,N )
- %GaussMtx Summary of this function goes here
- % Generate Bernoulli matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % Phi -- The Gauss matrix
- %% Generate Gauss matrix
- Phi = randn(M,N);
- %Phi = Phi/sqrt(M);
- end
2、伯努利矩阵
- function [ Phi ] = BernoulliMtx( M,N )
- %BernoulliMtx Summary of this function goes here
- % Generate Bernoulli matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % Phi -- The Bernoulli matrix
- %% (1)Generate Bernoulli matrix(The first kind)
- % 1--P=0.5 -1--P=0.5
- Phi = randi([0,1],M,N);%If your MATLAB version is too low,please use randint instead
- Phi(Phi==0) = -1;
- %Phi = Phi/sqrt(M);
- % %% (2)Generate Bernoulli matrix(The second kind)
- % % 1--P=1/6 -1--P=1/6 0--2/3
- % Phi = randi([-1,4],M,N);%If your MATLAB version is too low,please use randint instead
- % Phi(Phi==2) = 0;%P=1/6
- % Phi(Phi==3) = 0;%P=1/6
- % Phi(Phi==4) = 0;%P=1/6
- % %Phi = Phi*sqrt(3/M);
- end
3、部分哈达玛矩阵
- function [ Phi ] = PartHadamardMtx( M,N )
- %PartHadamardMtx Summary of this function goes here
- % Generate part Hadamard matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % Phi -- The part Hadamard matrix
- %% parameter initialization
- %Because the MATLAB function hadamard handles only the cases where n, n/12,
- %or n/20 is a power of 2
- L_t = max(M,N);%Maybe L_t does not meet requirement of function hadamard
- L_t1 = (12 - mod(L_t,12)) + L_t;
- L_t2 = (20 - mod(L_t,20)) + L_t;
- L_t3 = 2^ceil(log2(L_t));
- L = min([L_t1,L_t2,L_t3]);%Get the minimum L
- %% Generate part Hadamard matrix
- Phi = [];
- Phi_t = hadamard(L);
- RowIndex = randperm(L);
- Phi_t_r = Phi_t(RowIndex(1:M),:);
- ColIndex = randperm(L);
- Phi = Phi_t_r(:,ColIndex(1:N));
- end
4、部分傅里叶矩阵
以下代码生成的是部分傅里叶矩阵,这里要提的一点是,有关“归一化”的概念在网上说法不一,一种说法是向量除以向量各项之和,另一种说法是向量除以向量的2范数(各项平方之和再开方),这里采用后者。
- function [ Phi ] = PartFourierMtx( M,N )
- %PartFourierMtx Summary of this function goes here
- % Generate part Fourier matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % Phi -- The part Fourier matrix
- %% Generate part Fourier matrix
- Phi_t = fft(eye(N,N))/sqrt(N);%Fourier matrix
- RowIndex = randperm(N);
- Phi = Phi_t(RowIndex(1:M),:);%Select M rows randomly
- %normalization
- for ii = 1:N
- Phi(:,ii) = Phi(:,ii)/norm(Phi(:,ii));
- end
- end
5、稀疏随机矩阵
- function [ Phi ] = SparseRandomMtx( M,N,d )
- %SparseRandomMtx Summary of this function goes here
- % Generate SparseRandom matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % d -- The number of '1' in every column,d<M
- % Phi -- The SparseRandom matrix
- %% Generate SparseRandom matrix
- Phi = zeros(M,N);
- for ii = 1:N
- ColIdx = randperm(M);
- Phi(ColIdx(1:d),ii) = 1;
- end
- end
6、托普利兹矩阵和循环矩阵
这里先给出托普利兹矩阵,文中说经过M次循环,这里直接利用MATLAB自带函数生成一个托普利兹方阵再取前M行。
- function [ Phi ] = ToeplitzMtx( M,N )
- %ToeplitzMtx Summary of this function goes here
- % Generate Toeplitz matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % Phi -- The Toeplitz matrix
- %% Generate a random vector
- % %(1)Gauss
- % u = randn(1,2*N-1);
- %(2)Bernoulli
- u = randi([0,1],1,2*N-1);
- u(u==0) = -1;
- %% Generate Toeplitz matrix
- Phi_t = toeplitz(u(N:end),fliplr(u(1:N)));
- Phi = Phi_t(1:M,:);
- end
下面是循环矩阵,仍然利用MATLAB自带函数生成。
- function [ Phi ] = CirculantMtx( M,N )
- %CirculantMtx Summary of this function goes here
- % Generate Circulant matrix
- % M -- RowNumber
- % N -- ColumnNumber
- % Phi -- The Circulant matrix
- %% Generate a random vector
- % %(1)Gauss
- % u = randn(1,N);
- %(2)Bernoulli
- u = randi([0,1],1,N);
- u(u==0) = -1;
- %% Generate Circulant matrix
- Phi_t = toeplitz(circshift(u,[1,1]),fliplr(u(1:N)));
- Phi = Phi_t(1:M,:);
- end
以上都是采用MATLAB自带的toeplitz函数,这需要首先生成一个随机向量,而托普利兹矩阵和循环矩阵的区别也就在于这个随机向量的结构不同,注意代码中有关toeplitz函数的两个输入参数好好体会即可。
常见测量矩阵的MATLAB实现相关推荐
- matlab 测量矩阵,压缩感知中的常见测量矩阵及其MATLAB实现代码
压缩感知中的常见测量矩阵及其MATLAB实现代码 压缩感知中的常见测量矩阵及其MATLAB实现代码 原文链接:http://blog.csdn.net/jbb0523/article/details/ ...
- 压缩感知高斯测量矩阵matlab,压缩感知的常见测量矩阵
题目:压缩感知的常见测量矩阵 下面首先给出十篇参考文献中有关测量矩阵的叙述,然后以一篇硕士论文中对七种常见测量矩阵的描述依据,给出了这七种常见测量矩阵的MATLAB实现代码,以为以后的研究提供一个参考 ...
- 压缩感知的常见测量矩阵
题目:压缩感知的常见测量矩阵 下面首先给出十篇参考文献中有关测量矩阵的叙述,然后以一篇硕士论文中对七种常见测量矩阵的描述依据,给出了这七种常见测量矩阵的MATLAB实现代码,以为以后的研究提供一个参考 ...
- 浅谈压缩感知(三十二):压缩感知的常见测量矩阵
一.参考文献中常见的测量矩阵 [1]喻玲娟,谢晓春.压缩感知介绍[J].电视技术,2008,32(12):16-18. [2]李树涛,魏丹.压缩传感综述[J]. 自动化学报,2009,35(11):1 ...
- 压缩感知高斯测量矩阵matlab,浅谈压缩感知(十七):测量矩阵之有限等距常数RIC的计算...
有限等距常数(RestrictedIsometry Constant, RIC)是与有限等距性质(Restricted IsometryProperty, RIP)紧密结合在一起的一个参数. 一.RI ...
- 压缩感知测量矩阵构造方法研究
前言 这是我本科毕设的题目,现在毕业过了一个多月才想起来把毕设做的整理一下,写这个博客不为别的,只为纪念一下本科最后的时光吧.记得刚拿到这个题目时,满脑子都是黑人问号,压缩感知到底是个什么东东?也正是 ...
- 压缩感知高斯测量矩阵matlab,压缩感知matlab
MATLAB代码如下: 压缩感知仿真实例 压缩感知仿真实例在MATLAB R2001b中的计算结果如下: 原图像 采样率0.7 采样率0.5 采样率0.3 压缩感知仿真实例采用均方误差...... ( ...
- matlab 测量矩阵,急求一个测量矩阵采用分块多项式矩阵时怎样引用的代码!!!
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 各位大佬们,在压缩感知中,当测量矩阵采用分块多项式矩阵时,怎样写代码,下面是采用高斯矩阵时的代码 clc;clear %% 1. 时域测试信号生成 K=7 ...
- 范德蒙德矩阵在MATLAB中怎么表示,Python 之 Python与MATLAB 矩阵操作总结
Python 之 Python与MATLAB 矩阵操作小结 一.线形代数理论基础 线形代数(linear algebra)是数学的一个分支,研究矩阵理论.向量空间.线性变换和有限维线形方程组等内容. ...
最新文章
- nginx已添加systemctl,但不能开机自启动问题
- Using rqt_console and roslaunch
- maven deploy plugin_学习Maven之Maven Surefire Plugin
- mysql 主从同步配置
- sublime支持虚拟环境python2.7
- 未捕获的异常 'NSInternalInconsistencyException'
- jquery ajax xml attribute,获得jQuery ajax和asp.net webmethod xml响应工作
- Spring开发--Bean配置实例讲解
- Seata多微服务互相调用_全局分布式事物使用案例_业务数据库准备---微服务升级_SpringCloud Alibaba工作笔记0059
- 电脑常见故障排除手册(黑屏、死机、重启、蓝屏)
- 网络安全之DOS攻击
- XF660R型号良田高拍仪接口开发,通过图片文件的二进制数据进行图片上传
- 人工智能300年!LSTM之父万字长文:详解现代AI和深度学习发展史
- Docker的volumes踩坑
- deep learning (深度学习)资源整理,非常全面【转】
- Redis -- 主从
- vmware虚拟机上网设置教程(vmware虚拟机设置网络)
- 今晚 8 点,开发者赏金计划正式开启
- 思科网络安全 第十一章测验答案
- Python语法解析器PLY( lex and yacc in Python)