1.一种主从因素综合分析的转子堆叠装配优化方法，其特征在于，包括步骤：基于转子堆叠误差传播规律和不平衡量叠加机理建立转子堆叠精度预测模型的优化方法，基于所述转子堆叠精度预测模型的优化方法得到初步优化装配相位；基于所述初步优化装配相位，分别对每相邻两级转子之间的装配进行装配过盈量分析，得到离散过盈量函数；对所述离散过盈量函数做离散傅里叶变换，得到离散频谱能量函数；根据所述离散频谱能量函数构造优化目标函数；根据所述优化目标函数以及所述初步优化装配相位得到最终的优化装配相位；所述初步优化装配相位为：式(1)中，N为转子的总级数，m为相邻两级转子之间装配相位的个数，θ n 为第n级转子与第n+1级转子之间的初步优化装配相位，θ表示一组转子的装配相位参数， 为每相邻两种装配相位之间的相位差；基于所述初步优化装配相位，分别对每相邻两级转子之间的装配进行装配过盈量分析，得到离散过盈量函数具体为：测量每级转子的前安装边止口径向跳动数据和后安装边止口径向跳动数据；根据所述前安装边止口径向跳动数据、后安装边止口径向跳动数据以及所述初步优化装配相位建立转子止口表面形貌变化曲线；令相邻两级转子安装边的所述转子止口表面形貌变化曲线的最小二乘圆圆心重合，对于第i级转子装配，在 范围内，根据配合止口的表面形貌变化曲线积分计算过盈量或间隙量，得到离散过盈量函数；所述离散过盈量函数为：式(2)中，δ i (nθ 1 )表示第i级转子与第i+1级转子之间的装配在 范围内的离散过盈量；ρ 1 (θ)，ρ 2 (θ)分别代表外柱面和内柱面对应的转子止口表面形貌变化曲线在以相配合两级转子表面形貌曲线最小二乘圆圆心为原点建立的极坐标系下的函数关系，其形式与第i级转子与第i+1级转子之间的装配相位θ i 相关， 为常数，n＝0，1，...，m-1为离散过盈量函数的自变量；对所述离散过盈量函数做离散傅里叶变换，得到离散频谱能量函数具体为：对所述离散过盈量函数δ i (nθ 1 )在 范围内做离散傅里叶变换，得到离散频谱能量函数Δ i (kω 1 )，所述Δ i (kω 1 )为公式(3)所示：式(3)中，k＝0，1，...，n-1为离散频谱能量函数的自变量，j表示虚数单位， 为常数，i＝1，2，...，N-1；根据所述离散频谱能量函数构造优化目标函数具体为：取能量临界值E，当Δ i (1)>E时，认为接触应力分布不均主要由频率为ω＝ω 1 的部分引起，此装配相位不可取，当Δ i (1)<E时，认为接触应力分布不均主要由频率ω＞ω 1 的部分引起，以高频能量的最大值为优化目标，构造优化目标函数f i (θ i )；所述优化目标函数f i (θ i )如下：式(4)中，f i (θ i )为分段函数，当|Δ i (1)|为离散频谱能量函数Δ i (kω 1 )在k＝1时的模，|Δ i (1)|＞E时取-|Δ i (1)|，当|Δ i (1)|≤E时取Δ i (kω 1 )在k≠1时模的最大值，θ i 表示第i级转子与第i+1级转子之间的装配相位，由于离散过盈量函数的表达式中ρ 1 (θ)，ρ 2 (θ)的形式依赖于θ i 的取值，因此改变θ i 能够得到不同的优化目标函数值。
2.根据权利要求1所述的主从因素综合分析的转子堆叠装配优化方法，其特征在于，对第i级转子装配，取 范围内的装配相位分别计算优化目标函数f i (θ)值，当f i (θ)值不全为负时，取f i (θ)值最小的装配相位θ作为最终优化装配相位，当该范围内的七组f i (θ)值全为负值时，取f i (θ)值绝对值最小的装配相位θ作为最终优化装配相位。
3.一种计算机可读存储介质，其特征在于，存储有计算机指令，所述计算机指令用于被处理器执行时实现如权利要求1或2所述的主从因素综合分析的转子堆叠装配优化方法。
