多响应纵向变量广泛产生于各个前沿领域的定量化研究中,其协方差分析比单一响应纵向变量复杂许多,而目前可用的方法却不多见。本文提出一种将协方差矩阵和精度矩阵按照响应变量维数进行分块的正则化估计方法,首先将Cholesky分解方法推广到分块矩阵情形以确保估计中协方差矩阵的正定性,分解的因子矩阵也可分别对应于自回归模型的线性预测系数和预测误差方差。其次将带惩罚的负对数似然函数极小化问题化为关于该回归模型参数矩阵的双凸优化问题,建立交替凸搜索(Alternate Convex Search)算法加以解决,其中运用了Glasso方法和分块的坐标递减算法进行交替的凸优化,并证明了所提算法的收敛性。最后对提出的算法进行了模拟数据实验,通过对估计值的分析和与其他方法的比较验证了所提方法的有效性,尤其对于具有特殊相依结构的多响应纵向变量,验证了本文估计方法对其结构特点的体现。