采用最优化方法的张量计算及其应用

来源 :南京师范大学 | 被引量 : 1次 | 上传用户:zhfly6278
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
张量计算的理论与方法广泛应用于科学与工程领域的许多方面,例如:化学计量学、信号处理、数据挖掘和医学工程。本博士论文致力于研究这些实际应用中张量计算问题的内在结构特点,并发展相应的高水平最优化模型与方法。  第二章,我们首先研究医学工程中的一个张量估计问题。在核磁共振成像过程中,确保扩散张量的正定性是非常重要的,因为这个正定性反映的是水分子在复杂的生物组织环境中的扩散(布朗运动)现象。为保持扩散张量的正定性,我们将其转换成一个显式的正半定矩阵约束和一些线性矩阵等式约束。目标函数共有两部分,一是对数线性化的Stejskal-Tanner方程的最小二乘拟合;二是取正半定矩阵的核范数作为正则项,因为我们研究后认为正半定矩阵变量的秩低一点比较好。这样,我们就建立了一个凸二次的半定规划模型。该模型的全局最优解是肯定存在的,且可以通过以下三种算法求得。对原始问题,我们可以用两个成熟的软件:SDPT3和QSDP。然而,在研究了这一实际问题的内在特点之后,我们提出解其等价的对偶问题,算法使用的是增广的拉格朗日交替方向法。数值实验表明,交替方向法的计算速度至少是SDPT3和QSDP的十倍。又由于实际测量的信号有噪音,我们对新模型估计的最优扩散张量的系数与最优目标函数值作了灵敏度分析。对多个纤维交叉的人工实验表明新方法对Rician噪音是稳健的,并且优于现有方法。同时,新方法对纤维方向分布函数的估计与Q-球成像法一样成功。进一步的,对实际的人体大脑数据进行了数值实验,我们发现新算法重建的纤维结构准确表达了真实的神经纤维束的形态。另外,新方法产生的广义扩散张量在任何点处都是正定的,而不考虑正定性约束的最小二乘法是做不到的。  第三章,我们考虑张量分解问题。在信号处理、数据分析和科学计算中,人们常常需要将一个张量分解成一些因子的乘积。求解这类问题的数值方法——如交替最小二乘算法——具有两个关键要素:搜索方向与步长因子。本章就从这两个要素着手。一方面,由于该问题是非线性的,常用的线性化搜索方向不是很有效。所以,我们提出两个高阶的搜索方向。一个称为几何搜索方向,它将最近的两个线性搜索方向作了线性组合。第二个是代数搜索方向,它的构造是通过对三个连续的迭代点进行二次插值而得到的。另一方面,我们考虑沿这些搜索方向进行加强线搜索(精确线搜索)。一个最优的复数步长因子含有两个独立的参数:模与辐角。现有的算法ELSCS是采用交替的方式计算这两个参数的。这样找到的复数步长就有可能不是全局最优的。我们给出了一种同时计算模与辐角这两个参数的新方法,从而可以保证得到的复数步长因子是全局最优的。最后,对带卷积的DS-CDMA混迭数据的盲分离均衡问题进行了数值试验,我们比较了各种搜索方向与步长因子。结果表明新的搜索方向极大的提高了交替最小二乘算法的效率。同时,新的最优步长策略也是很成功的。  第四章,我们对最小二乘意义下高阶张量的最佳秩(一)逼近问题设计了一个快速的优化方法——序列无约束对偶优化算法。该算法的主要思想是把张量最佳秩一逼近的带约束的对偶问题转化成一系列的无约束优化问题。然后,我们设计了一个快速的梯度类方法解这些无约束的优化问题,算法的创新点在于给最速上升方向装备了一个新的初始步长策略,然后采用后退线搜索技术产生新的迭代点。理论分析表明新算法对任意的初始点都可以收敛到对偶问题的一个KKT点。数值实验表明,序列无约束对偶优化算法要优于交替最小二乘算法。
其他文献
据《Scientia Horticulturae》的一篇研究报道(http://dx.doi.org/10.1016/j.scienta.2014.06.007),来自葡萄牙的研究人员研究了生长于地中海的赛美蓉和小粒白麝香酿酒葡萄的
本文对双正则分拆函数的同余性质进行了系统的研究,并证明得到几组双正则分拆函数模5的同余关系式.具体研究工作如下:第一章,给出整数分拆及其同余性质的研究背景和研究意义,简要介绍双正则分拆函数的基本概念,性质和研究现状.第二章,介绍几种双正则分拆函数同余关系的证明方法,如以三次theta函数,Euler五角数定理,Jocabi恒等式等为工具进行证明.第三章,以第二章中介绍到的双正则分拆函数同余关系的证
学位
设K为代数数域,H为K的Hilbert类域.通过Artin映射,Galois群G=Gal(H/K)同构于K的理想类群C(K).域H的被G2固定的子域E=HG2称作K的Hilbert亏格域.由Galois理论  Gal(E/K)(≌)G/G
摘 要:分析单向阀弹簧断裂对高压氨泵安全运行的危害,制定相应的解决方法,解决氨泵长周期运行的问题。  关键词:高压液氨泵;单向阀;对策;导向  一、高压氨泵的简介  宁夏石化公司一化肥尿素装置使用的2台高压液氨泵是德国WORTHINGTON制造的VSE-H型7联立式柱塞泵。高压氨泵运行时,将由合成装置来的25℃、2.5Mpa的液氨升压到16.0Mpa后,送入高压系统进行合成尿素的反应。氨泵尿素装置
组合邻差方法(The method of combinatorial telescoping)是一种证明和发现q-恒等式的有效方法。本文通过对组合邻差方法的深入研究,得到一些新的结果:将组合邻差方法从交错和
保险业是作为金融行业的三大支柱产业之一,保险公司的破产概率自然是研究的重点,在金融保险领域中起着至关重要的作用。破产概率的研究从产生的那一天开始就一直是一个活跃的研
作为金融行业中的一个重要组成部分,保险业在服务社会中的作用日益增多。人寿保险,作为保险公司的主营项目,由于它与社会大众利益的密切联系,引起了普通大众的广泛关注。保险
近年来,最优化方法在物理、化学、金融等领域中得到大量的应用。然而,对于一些实际问题,其目标函数经常是通过计算机模拟出来的,它们的导数信息不可求、不可信或者是非常耗时的。