Manta rays have emerged as excellent bionic objects due to their efficient and flexible propulsion. A thorough understanding of the propulsion characteristics of manta ray robot pectoral fins is highly important for optimizing pectoral fin design and guidance strategies. Based on the silhouette and kinematics of the manta ray, we designed flexible pectoral fins and a 2DoF bionic mechanism that couples the flapping and pitching movements. Blade element theory was used to conduct a dynamic analysis of the bio-manta pectoral fins. The Strouhal number is used as a key dimensionless number that affects the hydrodynamic characteristics of a bio-manta robot during periodic flapping, as well as the flexibility of the pectoral fins and the angle of attack. We describe a series of hydrodynamic experiments conducted to analyze the single and interactive effects of the motion parameters of pectoral fin propulsion and regression analysis. The results show that when St is [0.2,0.4], the CT (thrust coefficient) of the flexible pectoral fin increases rapidly. When St increased beyond 0.4, the growth trend slows and remains essentially unchanged after 0.8. The CT reaches its peak, the propulsion efficiency is highest when St is between [0.4,0.6]. The angle of attack and the CT show a negative inverse quadratic trend, reaching a peak at 6.5°, and their influence is weaker than that of the frequency. The flexibility and CT showed a cubic trend, and the influence of flexibility was more significant when swimming at medium and low speeds. When swimming slowly, flexibility has a significant interactive effect on the frequency and angle of attack. However, these interactive effects reduce at high speeds. When the bio-manta robot starts at 0 speed, it should maneuver at high frequency and large amplitude while ensuring sufficient structural strength. After reaching a certain speed, the device maintains a large amplitude and low-frequency flapping operation to continuously drive propulsion.