非线性特征值问题出现在量子点数值模拟、计算电磁场、粘弹性结构的振动分析、时滞系统的稳定性分析等许多领域,其研究具有重要的理论意义和应用价值。本文研究求解非线性特征值问题的数值方法,并将围道积分方法应用于判定n维时滞系统的稳定性,主要创新工作如下:提出了求解非线性特征值问题的逐次m次近似方法,给出了该方法的收敛性结果。为解决逐次m次近似方法产生的大规模多项式特征值问题,结合求解多项式特征值问题的部分正交投影方法,发展了求解非线性特征值问题的一个有效的逐次m次近似方法。为降低基于奇异值分解的Newton方法的计算量,采用逆迭代技术,每步近似计算最小奇异值以及相应的左右奇异向量,发展了改进的Newton方法,并证明了改进的Newton方法具有局部二次收敛性。研究了确定某一区域内非线性特征值问题的特征值代数重数的方法。将幅角原理推广到有理数情形,给出了求非线性特征值问题在某一区域内所有特征值及相应代数重数的数值方法;为了计算多项式特征值问题在一个半开平面的所有特征值,基于分式线性映射,建立了一个多项式特征值问题在半开平面内的所有特征值与求另一个多项式特征值问题在开圆盘中所有特征值的关系,提出了计算多项式特征值问题在半开平面内所有特征值的围道积分方法。对n维线性多时滞系统,利用围道积分方法,提出了计算该系统最右端特征值的数值方法,发展了判定n维线性多时滞系统在正平衡点处局部稳定性的新方法;对n维非线性多时滞Lotka-Volterra系统,给出了判定Lotka-Volterra系统在平衡点处局部稳定性的数值方法。数值结果说明了本文所给求解非线性特征值问题的数值方法都是有效的。