数码相机标定方法研究

数码相机在计算机视觉中的应用逐渐普及和深入，在应用中数码相机的标定是相当重要的，根据数码相机的特点，提出了一种新的标定方法，针对数码相机的某一状态，先确定其内部参数矩阵，再通过采集的图像术出外部参数矩阵，首先建立了相机成像的几何模型，并将此模型分解成内、外参数矩阵。文中详细介绍了数码相机内部参数矩阵中的各元素 应的物理参数以及求解它们的原理和方法，从而建立内部参数矩阵，推导了求解外部参数矩阵的计算公式，并且介绍了相应的求解方法，再由外部参数矩阵求解出对应的物理参数，通过实验表明这种方法有相当高的精度，且操作方便。     

基于Lyapunov方法的非线性系统自适应观测器设计

针对一类满足L ipsch itz条件的具有未知参数的非线性系统,利用Lyapunov方法对L ipsch itz非线性系统自适应观测器的设计问题进行了研究。基于分析求解代数R iccati方程给出求解问题的不完善性、特征结构配置理论给出设计方法的重特征值限制性、对不同形式的观测器增益矩阵求解方法进行比较,最终选用线性矩阵不等式来改进观测器增益矩阵的选取方法。在观测误差稳定的条件下,得出了基于线性矩阵不等式方程设计状态观测器的增益矩阵,保证系统的状态估计误差收敛到零,并对其进行了仿真研究。结果证明,本文所构造的非线性观测器增益矩阵方法明显优越于其他方法,增强了系统的鲁棒性。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

Lyapunov矩阵方程的一种解法

本文提供了求解Lyapunov矩阵方程的一个新算法。首先,对系数矩阵具有Frobenius形式的方程给出了一种公式解法。然后通过对系数矩阵的约化,将解法推广列一般情形。最后,给出了数值例子,并同有关的方法进行了比较。     

线性方程组的并行算法及在Transputer系统上的实现

在“一种有效的多Transputer系统的并行算法——ABC法”一文的基础上,本文进一步研究将ABC法用于变带宽矩阵线性方程组的求解问题,对线性方程组的系数矩阵采用了逐行一维存储方式,提出了相应的并行Gauss消元法,给出了该算法的效率.分析结果表明,带宽越大方程阶数越高,这种算法的效率就越高。因此本算法适用于高阶的大带宽线性方程组的求解问题. 根据本文的算法,编制了线性方程组的并行求解程序,并分别在一个、二个和四个T414系统上做了若干算例,结果表明本文分析的结论是正确的。     

求解大规模矩阵问题的Krylov子空间方法

求解大规模矩阵问题包括线性方程组和特征值问题等是计算数学和科学工程计算中的重大课题，最近几年，其研究工作取得了许多重大进展。文中给出大型线性方程组和特征值问题Krylov子空间方法若干进展的一个概述，其中包括作者对这些问题的研究成果。涉及的专题包括求解大型线性方程组的共轭梯度法、SYMMLQ算法、MINRES算法、GMRES算法、Lanczos双正交化算法、QMR算法以及这些算法的块格式；求解大对称特征值问题的Lanczos算法和块Lqnczos算法；求解大型非对称特征问题的Lanczos算法、Arnodi算法以及这些算法的推广。讨论求解大规模矩阵问题的加速技术和预处理技术。了一些有待进一步研究的问题。     

基于测量试验数据修正有限元模型质量矩阵

对无阻尼结构系统有限元模型质量矩阵修正问题,以该矩阵修正量的F-范数为目标函数,并以待修正质量矩阵应具有的性质,如满足正交关系,对称性,半正定性和稀疏性作为约束条件,数学上形成带约束的矩阵最佳逼近问题。给出了问题有解的条件,基于循环投影方法,提出了求解矩阵最佳逼近问题的数值方法。数值结果说明了所给方法的有效性。     

阻尼结构振动问题特征值定位的一个定理

本文的主要结果是得出了矩阵B—ωA(其中B为对称正定矩阵,A为Hermite矩阵)的首主子式具有类似于Sturm序列的性质,从而建立了求解阻尼结构振动方程M(?) C(?) Kq=0特征值问题的行列式查找法的理论基础,改正了K.K.Gupta在文献[1]—[5]中的错误提法。     

由三个向量对构造三对角对称矩阵

文章讨论利用给定的三个向量对构造不可约三对角矩阵、Jacobi矩阵和负Jacobi矩阵的反问题.在求解方法中,将已知的-些关系式等价地转化为线性方程组,利用线性方程组有解的条件,得到了所研究问题有惟一解的充要条件,并给出了数值算法和例子.     

用有限元素法计算和多点激振法试验确定某型机翼的动力特性参数

本文采用直接刚度法得到刚度矩阵,用聚缩质量模型得到质量矩阵。应用基于矩阵分解的矩阵迭代法,求解了某型机翼前六阶固有频率、固有振型和广义质量等动力特性参数。对同一机翼采用多点激振法进行了模态试验,测得了前十二阶固有频率、固有振型和几阶广义质量。文中对计算和试验结果进行了比较和分析。     

BEST APPROXIMATION BY NORMAL MATRICES WITH SPECTRAL CONSTRAINTS

考虑在给定谱约束和Ｆｒｏｂｅｎｉｕｓ范数意义下用正规矩阵最佳逼近一个给定复方阵的问题。给出了这个问题可解的充分必要条件，提出了求解这个问题的一个数值算法，并给出了一个数值例子。     

基于细以方法的颗粒增强复合材料弹塑性分析

针对高体积含量的颗粒增强复合材料，提出了一种细观结构模型：将颗粒简化为同质、同尺寸的弹性圆球，两颗粒之间的连接基本简化为一弹塑性短圆柱体，并假设细观应力、应变和塑性区均为轴对称分布。其体和颗粒的变形行为分别简化为类似于弹性地基和弹性半空是，以此为基础，建立了反映一对颗粒法向和切向变形协调关系的两个积分方程。数值求解这两个方程，可得到一对颗粒间基体中细观应力的分布形式，从而建立颗粒对的细观弹塑性本构，采用平均化方法，进一步推导出材料宏观的应力－应变关系。本文利用该模型对一种金属基复合材料的拉伸实验进行了模拟，理论预测与实验结果吻合较。     

求解-维浅水波方程的中心差分格式

使用Haim Nessyahu和Eitan Tadmor给出的求解双曲型守恒律的中心差分交错格式求解浅水波方程,经一维溃坝问题数值实验验证了该格式对于求解浅水波方程的可行性和有效性。     

求解织物热湿耦合方程的控制体-时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方程的算法。首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化连续性的非线性微分方程组为一系列的递归形式的线性代数方程组。文中给出了问题的求解步骤和算例。结论表明该方法的预测结果具有不依赖时间步长的特点。     

求解代数Riccati方程的一种新算法

本文讨论了Hamiltonian矩阵在辛相似变换下的标准形,由此给出代数Riccati方程存在非负定解的一个充分条件,提出了求解代数Riccati方程的一种新的算法。该方法节省运算量,且精度较高,尤其适合于解阶数不太高、系数矩阵为满阵的Riccati方程。最后给出了一个数值例子,并将该方法与其他方法作了比较.     

周期对称三对角矩阵的特征值反问题

本文提出由两个特征值和相应的特征向量构造周期对称三对角矩阵的一类特征值反问题,讨论了这类问题的可解性,给出了这类问题有解的充分必要条件,描述了求解这类问题的数值算法,并且给出了数值例子。     

叶片非线性瞬态响应计算方法与参数选择

在计算叶片的岛撞击响应时，为提高计算精度及降低计算成本，需要为待解的非线性动力方程组制定一个合理的求解方案，这包括选择系数矩阵的算法与非线性动力方程组的解法，以及选择主要计算参数，本文介绍了当使用ＡＤＩＮＡ程序计算平板叶片在冲载荷下的非线性瞬态响应时，对计算方法与计算参数的选择，算例说明了时间步长的重要影响作用。     

求解Jacobi矩阵逆特征值问题的一个新算法

研究以下反问题：给定两组实数｛λｉ｝ｉ＝１＾ｎ，｛μｉ｝ｉ＝１＾ｎ－１，满足如下分隔条件：λｉ〈μｉ〈λｉ＋１，要求构造一个ｎ阶Ｊａｃｏｂｉ矩阵Ｊｎ，使得｛λｉ｝ｉ＝１＾ｎ是Ｊｎ的特征值，｛μｉ｝ｉ＝１＾ｎ－１是Ｊ２，ｎ的特征值，本文利用Ｊａｃｏｂｉ矩阵的性质，导出了求解上述问题的一个算法。     

利用NURB作曲线和曲面的插值

给出了非均匀有理三次Ｂ－样条插值的方法，利用曲线（面）方程的矩阵表达式导出了反求顶点问题的方程组，使方程组的系数矩阵呈三对角型，易于求解。同时，分别推出几种情况下的端点（边界）条件。利用该结果可以使得曲线（面）易于调整，运算速度快，便于处理。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。