Respiratory sinus arrhythmia (RSA) represents a physiological phenomenon of cardiopulmonary interaction. It is known as a measure of efficiency of the circulation system, as well as a biomarker of cardiac vagal and well-being. In this article, RSA is modeled as modulation of heart rate by respiration in an interactive cardiopulmonary system with the most effective system state of resonance. By mathematically modeling of this modulation, we propose a quantitative measurement for RSA referred to as "Cardiopulmonary Resonance Function (CRF) and Cardiopulmonary Resonance Indices (CRI)," which are derived by disentanglement of the RR-intervals series into respiratory-modulation component, R-HRV, and the rest, NR-HRV using spectral G-causality. Evaluation of CRI performance in quantifying RSA has been conducted in the scenarios of paced breathing and in the different sleep stages. The preliminary experimental results have shown superior representation ability of CRF and CRI compared to Heart Rate Variability (HRV) and Cardiopulmonary Coupling index (CPC).