In this paper, we first design the corresponding integration algorithm and matlab program according to the Gauss–Legendre integration principle. Then, we select the Lorenz system, the Duffing system, the hidden attractor chaotic system and the Multi-wing hidden chaotic attractor system for chaotic dynamics analysis. We apply the Gauss–Legendre integral and the Runge–Kutta algorithm to the solution of dissipative chaotic systems for the first time and analyze and compare the differences between the two algorithms. Then, we propose for the first time a chaotic basin of the attraction estimation method based on the Gauss–Legendre integral and Lyapunov exponent and the decision criterion of this method. This method can better obtain the region of chaotic basin of attraction and can better distinguish the attractor and pseudo-attractor, which provides a new way for chaotic system analysis. Finally, we use FPGA technology to realize four corresponding chaotic systems based on the Gauss–Legendre integration algorithm.