学位论文 > 优秀研究生学位论文题录展示
采用最优化方法的张量计算及其应用
作 者: 陈艳男
导 师: 孙文瑜
学 校: 南京师范大学
专 业: 计算数学
关键词: 扩散张量成像 高阶张量 半定规划 张量分解 最优步长 PARAFAC 最优化方法
分类号: O183.2
类 型: 博士论文
年 份: 2013年
下 载: 85次
引 用: 0次
阅 读: 论文下载
内容摘要
张量计算的理论与方法广泛应用于科学与工程领域的许多方面,例如:化学计量学、信号处理、数据挖掘和医学工程。本博士论文致力于研究这些实际应用中张量计算问题的内在结构特点,并发展相应的高水平最优化模型与方法。第二章,我们首先研究医学工程中的一个张量估计问题。在核磁共振成像过程中,确保扩散张量的正定性是非常重要的,因为这个正定性反映的是水分子在复杂的生物组织环境中的扩散(布朗运动)现象。为保持扩散张量的正定性,我们将其转换成一个显式的正半定矩阵约束和一些线性矩阵等式约束。目标函数共有两部分,一是对数线性化的Stejskal-Tanner方程的最小二乘拟合;二是取正半定矩阵的核范数作为正则项,因为我们研究后认为正半定矩阵变量的秩低一点比较好。这样,我们就建立了一个凸二次的半定规划模型。该模型的全局最优解是肯定存在的,且可以通过以下三种算法求得。对原始问题,我们可以用两个成熟的软件SDPT3和QSDP。然而,在研究了这一实际问题的内在特点之后,我们提出解其等价的对偶问题,算法使用的是增广的拉格朗日交替方向法。数值实验表明,交替方向法的计算速度至少是SDPT3和QSDP的十倍。又由于实际测量的信号有噪音,我们对新模型估计的最优扩散张量的系数与最优目标函数值作了灵敏度分析。对多个纤维交叉的人工实验表明新方法对Rician噪音是稳健的,并且优于现有方法。同时,新方法对纤维方向分布函数的估计与Q-球成像法一样成功。进一步的,对实际的人体大脑数据进行了数值实验,我们发现新算法重建的纤维结构准确表达了真实的神经纤维束的形态。另外,新方法产生的广义扩散张量在任何点处都是正定的,而不考虑正定性约束的最小二乘法是做不到的。第三章,我们考虑张量分解问题。在信号处理、数据分析和科学计算中,人们常常需要将一个张量分解成一些因子的乘积。求解这类问题的数值方法——如交替最小二乘算法——具有两个关键要素:搜索方向与步长因子。本章就从这两个要素着手。一方面,由于该问题是非线性的,常用的线性化搜索方向不是很有效。所以,我们提出两个高阶的搜索方向。一个称为几何搜索方向,它将最近的两个线性搜索方向作了线性组合。第二个是代数搜索方向,它的构造是通过对三个连续的迭代点进行二次插值而得到的。另一方面,我们考虑沿这些搜索方向进行加强线搜索(精确线搜索)。一个最优的复数步长因子含有两个独立的参数:模与辐角。现有的算法ELSCS是采用交替的方式计算这两个参数的。这样找到的复数步长就有可能不是全局最优的。我们给出了一种同时计算模与辐角这两个参数的新方法,从而可以保证得到的复数步长因子是全局最优的。最后,对带卷积的DS-CDMA混迭数据的盲分离均衡问题进行了数值试验,我们比较了各种搜索方向与步长因子。结果表明新的搜索方向极大的提高了交替最小二乘算法的效率。同时,新的最优步长策略也是很成功的。第四章,我们对最小二乘意义下高阶张量的最佳秩一逼近问题设计了一个快速的优化方法——序列无约束对偶优化算法。该算法的主要思想是把张量最佳秩一逼近的带约束的对偶问题转化成一系列的无约束优化问题。然后,我们设计了一个快速的梯度类方法解这些无约束的优化问题,算法的创新点在于给最速上升方向装备了一个新的初始步长策略,然后采用后退线搜索技术产生新的迭代点。理论分析表明新算法对任意的初始点都可以收敛到对偶问题的一个KKT点。数值实验表明,序列无约束对偶优化算法要优于交替最小二乘算法。
|
全文目录
本文创新点 7-8 摘要 8-10 Abstract 10-12 Chapter 1 Fundamentals 12-14 1.1 The concept of tensors 12-13 1.2 The rank of tensors 13-14 Chapter 2 Positive Semidefinite Generalized Diffusion Tensor Imaging 14-42 2.1 Preliminary of diffusion tensor imaging 14-19 2.2 Convex quadratic SDP model 19-23 2.3 State-of-the-art solvers:SDPT3 and QSDP 23-24 2.4 A new solver:ADM 24-27 2.5 Sensitivity analysis 27-30 2.6 Experimental results 30-40 2.6.1 Choosing the regularization parameter μ 30-31 2.6.2 Synthetic experiments 31-35 2.6.3 Human brain study 35-39 2.6.4 Thre solvers comparison 39-40 2.7 Conclusion 40-42 Chapter 3 New ALS Methods for Tensor Decompositions 42-68 3.1 Fundamentals of tensor decompositions 42-45 3.2 Search directions 45-49 3.2.1 Preliminaries and the linear search direction 45-46 3.2.2 Geometric search direction 46-47 3.2.3 Algebraic search direction 47-49 3.3 The optimal step size 49-57 3.4 The new ALS algodthm 57-58 3.5 Simulation results 58-67 3.5.1 The first result:step size 59-63 3.5.2 The second result:directions 63-66 3.5.3 The third result:nearly collinear case 66-67 3.6 Conclusions 67-68 Chapter 4 SUDO for the Tensors' Rank-One Approximation 68-85 4.1 The rank-one approximation to tensors 68-71 4.2 The SUDO Algorithm 71-78 4.2.1 The steepest ascent direction 71-73 4.2.2 Initial step length 73-75 4.2.3 Backtracking line search 75-77 4.2.4 The SUDO algorithm 77-78 4.3 Global convergence 78-81 4.4 Numerical experiments 81-84 4.4.1 The first numerical experiment:Iocal vs. global solution 81-82 4.4.2 The second numerical experiment:convergence rate 82-84 4.5 Conclusions 84-85 Bibliography 85-92 Publications 92 科研项目 92-94 Acknowledgements 94
|
相似论文
- 3.0T高场磁共振扩散张量成像对特发性癫痫的研究,R445.2
- 磁共振扩散张量成像在骨骼肌内肿瘤的应用价值研究,R445.2
- 脑梗死后椎体束华勒变性的功能磁共振研究,R445.2
- 一类非线性切换系统的吸引域估计及控制器设计,TP13
- 半定规划内点算法的搜索方向研究,O221.2
- 半定规划问题的若干算法研究,O221.2
- 解半定规划的两种数值方法,O221.2
- GSM网络优化性能研究,TN929.532
- 热电偶时间常数的测量,TH811
- DTI及DTT技术与NIHSS评分系统在急性脑梗塞患者预后评价中的对比研究,R743.3
- 社会标签推荐张量分解方法研究,TP393.09
- 多波叠前AVA反演中的混沌遗传优化方法研究,P631.44
- 煤矿DTY型挂缆车自动支撑的研究,TD50
- 原发性慢性闭角型青光眼视神经DTI和Bold-fMRI应用,R775
- 海洛因依赖及美沙酮维持治疗者大脑白质微结构的DTI研究,R749.6
- 一个求解非线性半定规划的基于分解的原始对偶内点算法,O221.2
- 城市环境噪声对居住区声环境的影响及优化方法研究,X593
- 扩散张量及纤维束成像在正常人与青光眼患者视通路的临床应用研究,R775
- 磁共振扩散张量成像对中枢神经系统正常老化的定量分析,R741
- 失神经骨骼肌及周围神经的MRI成像研究,R746
- 跨文化翻译中异化策略的使用范围和优化方法,H059
中图分类: > 数理科学和化学 > 数学 > 几何、拓扑 > 向量(矢量)和张量分析 > 张量分析
© 2012 www.xueweilunwen.com
|