利用投影算法来计算系统矩阵左乘和右乘
2016-12-21 09:05
357 查看
采用distance driven投影模型
mode 1: 左乘得图像向量
利用系统矩阵进行验证
说明正确!
mode2: 右乘得投影向量
结果不一致,但差异非常小。
这种差异,我实在不知道是怎么造成的。
mode 1: 左乘得图像向量
% 参数设置 para.ht_orig = 0; %切片初始高度,即第一(下层)切片与探测器的间距,单位是切片个数。 para.pv_ind = 14; %投影角序号,如1、3、14等。 para.nvx = 1080; %x轴方向体素个数1080 para.nvy = 2304; %y轴方向体素个数2304 para.nvz = 1; %z轴方向体素个数 % mode = 1; % 右乘 % img_mat = ones(para.nvx*para.nvy*para.nvz,1); mode = 2; % 左乘 img_mat = ones(1,2580480); [sps_sys_mat,myproj_vectori] = voxel_based_ddp2_fun2(img_mat,para,mode);
利用系统矩阵进行验证
>> vector2 = ones(1,2580480); >> rslt2 = vector2*sps_sys_mat; >> isequal(rslt2,myproj_vectori) ans = 1
说明正确!
mode2: 右乘得投影向量
>> vector1 = ones(2488320,1); >> rslt1 = sps_sys_mat*vector1; >> isequal(rslt1,myproj_vectori) ans = 0
结果不一致,但差异非常小。
>> sub=abs(rslt1-myproj_vectori); >> max(sub) ans = 7.1054e-15 >> min(sub) ans = 0
这种差异,我实在不知道是怎么造成的。
相关文章推荐
- 利用投影算法来计算系统矩阵左乘和右乘
- 利用投影算法来计算系统矩阵左乘和右乘
- 利用矩阵计算提升推荐系统的速度
- 利用MPICH2计算矩阵相乘的简单算法
- C++实现:双输入单输出系统:建立一个 感知器网络,实现上述样本的分类。计算出相应的网络权值矩阵w。
- 算法学习之一(二):矩阵计算——LU分解
- Strassen矩阵乘法 + 快速计算乘方的算法 + 矩阵的次幂
- 利用OpenCV计算图像的垂直和水平积分投影
- 利用OpenCV计算图像的垂直和水平积分投影
- 4 利用计算系统的输出
- 算法学习之一:矩阵计算——高斯消元法
- GPU上大规模稀疏矩阵特征值计算高效算法之一——GPU介绍
- 利用高级Java、算法、三角学、分布计算设计自己的智能机器人--(重锤痛击续2)
- GPU上大规模稀疏矩阵特征值计算高效算法之二——稀疏矩阵
- GPU上大规模稀疏矩阵特征值计算高效算法之三——SLEPc测试
- 【算法】 输入n 输出一个n*n的zigzag矩阵 利用c++实现
- 最大子矩阵问题;枚举行的组合,然后利用一维最大子段和的DP算法;
- 数值计算方法与算法:C语言实现利用Gauss消元法解方程组
- HTML5 利用矩阵计算坐标
- 关于Floyd-Warshall算法由前趋矩阵计算出的最短路径反映出了算法的执行过程特性的证明