基于GPU计算平台的Lyapunov方程求解

来源 :哈尔滨工业大学 | 被引量 : 0次 | 上传用户:cangzhe
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Lya p uno v方程在现代控制理论各个分支中都有非常重要的作用。因此对于Lyapuno v方程的求解是很重要的。如今,基于传统CPU或者多核技术已无法快速计算拥有大规模数据的Lyapuno v方程。而图形处理器GPU由于其独特的设计理念,适合并行计算大规模数据,被广泛应用于多个领域。同时简单牛顿迭代算法适合 GPU并行运行,因此基于GPU高性能计算设计牛顿迭代算法求解大规模Ly ap u no v方程是可行的,有意义的。  本文首先基于 GPU高性能计算平台对不同阶数的矩阵相乘和求逆进行计算。实验结果显示随着矩阵阶数增加,并行计算相对于串行计算的加速比也在增大,这充分体现了GPU在并行计算上的强大能力。随后,分别基于CPU串行和GPU并行计算平台,设计了简单牛顿迭代算法和基于 UD分解的牛顿迭代算法,并对不同阶数的Lya p uno v方程进行了数值求解。将计算结果进行对比,数据表明:当矩阵阶数为4096时,两种算法并行相对于串行的加速比分别为34、57,充分体现了GP U高性能计算的高效性。  其次,将简单牛顿迭代算法推广到控制理论中的另一重要方程—Ricatti方程的求解。对比并行计算结果和串行计算结果,数据表明:当矩阵规模达到4096时,并行算法相对于串行算法的计算时间加速比高达280。充分体现了牛顿迭代算法基于GPU对Ricatti方程进行求解的可行性和高效性。  最后,将精确线性搜索和Newton迭代结合起来,构造了一种基于 GPU平台的改进 Newton迭代算法,并应用于不同阶数的Ricatti方程的求解。并行算法与串行算法的计算结果对比,结果表明,当矩阵规模达到4096时,改进牛顿迭代算法相对于串行牛顿迭代算法的计算时间加速比却高达300。同时在使用改进牛顿迭代算法对方程进行求解的过程中,迭代精度相对于串行牛顿迭代算法的迭代精度要高很多,也可以很快的达到收敛条件,表明了改进牛顿迭代算法相对于其它两种算法在精度上有很大的优势。
其他文献
树是无向的、无环、无重边的图,其拉普拉斯多项式为:Λ(G,λ)=det(λIn?L)=∑nk=0(?1)kCkλn?k,Cn?2(T)是T的维纳指数。对于两个n个顶点的树,如果(C0(T1),···, Cn(T1))≤(