随着天体力学和天体物理学的不断发展,人们对宇宙了解的深入,引力波的探测成为了几十年来的热点问题。引力波到达地面极其微弱,探测起来非常困难,这就需要高精度的理论模板作为支撑。连续性引力波为探测的首选,因此,由黑洞或者中子星组成的致密双星系统成为了最为理想的引力波探测波源。本文的主要工作也正是在以此为背景进行展开的。本文主要是对后牛顿哈密顿形式的双黑洞系统,含有牛顿项、轨道部分精确到3PN,自旋—轨道部分精确到3.5PN,此时,系统含有7个独立的孤立积分,由于自旋变量非正则的,经典哈密顿理论不再适用,当对自旋变量作正则化处理后,系统由12维变成10维,并且含有5个独立的孤立积分,有总能量积分、三维角动量、轨道角动量大小积分,当系统含有自旋—自旋耦合项精确到2PN,系统是不可积的,通过数值模拟可以验证。自旋—自旋哈密顿量占总哈密顿量的比重的临界值,也可以作为区分系统有序和混沌。最后,对系统的构型图,截面图和引力波波形,进行数值模拟,研究发现有序轨道的引力波是有序的,呈现周期性变化,而混沌轨道的引力波波形是混沌的,表现出不可预知性,明显杂乱无章。这些工作也为引力波探测做了充分的理论准备,提供了参考。