A mathematical model used to study cardiopulmonary interaction was developed based on our previous multi-element nonlinear cardiovascular model, and the impacts of respiratory movement on cardiovascular hemodynamics were studied by computer simulation. To verify the model, the comparisons of simulation results, which are the high frequency power (HF, 0.15-0.40Hz) of heart rate variability (HRV) and blood pressure variability (BPV) spectrum, with the experimental data reported in references were performed. It indicated that the simulation results agreed well with the experimental data, and the extended cardiovascular model could reflect the cardiopulmonary interaction. In addition, the changes in HRV caused by combined variations of arterial baroreflex sensitivity (ABS) and cardiopulmonary baroreflex sensitivity (CBS) were investigated through this model. The simulation results showed that the lowered ABS and CBS were the main cause of the decrease in respiratory sinus arrhythmia (RSA) with human body aging.