We consider a nonlinear stochastic differential equation driven by an \(\alpha \)-stable Lévy process (\(1<\alpha <2\)). We first prove existence and uniqueness of the invariant measure by the Bogoliubov-Krylov argument. Then we obtain some regularity results for the probability density of its invariant measure by establishing the a priori estimate of the corresponding stationary Fokker-Planck equation. Finally, by the a priori estimate of the Kolmogorov backward equation and the perturbation property of the Markov semigroup, we derive the response function and generalize the famous linear response theory in nonequilibrium statistical mechanics to non-Gaussian stochastic dynamic systems.