The solubility of edaravone in four mixed solvents including (ethyl acetate + n-propanol), (acetonitrile + methanol), (acetonitrile + ethanol), and (acetonitrile + n-propanol) was determined by using the isothermal saturation method at temperatures ranging from (273.15 to 313.15) K under 101.3 kPa. At a certain composition in the four mixed solvents, the solubility data of edaravone increased with the rising temperature. At a certain temperature, the solubility data decreased monotonically with the increasing mass fraction of acetonitrile in the mixture of (acetonitrile + n-propanol). Whereas in (ethyl acetate + n-propanol), (acetonitrile + methanol), and (acetonitrile + ethanol), it increased first and then decreased with the increasing mass fraction of ethyl acetate or acetonitrile. The composition dependence of solubility has a maximum at around w(ethyl acetate) = 0.2 in (ethyl acetate + n-propanol) or w(acetonitrile) = 0.2 in (acetonitrile + ethanol) and w(acetonitrile) = 0.4 for mixtures of (acetonitrile + methanol). The results obtained in the selected four mixtures were correlated by the Jouyban-Acree model and CNIBS/R-K model. The largest relative average deviation and root-mean-square deviation values were found in (acetonitrile + methanol) calculated by the Jouyban-Acree model and were no more than 4.71 x 10(-2) and 22.32 x 10(-4) , respectively. Statistical analysis of experimental and calculated values shows that the residual of CNIBS/R-K model is smaller, and the linear correlation coefficient is closer to 1. Therefore, the CNIBS/R-K model is more suitable to correlate the solubility values of edaravone in the four mixed solvents including (ethyl acetate + npropanol), (acetonitrile + methanol), (acetonitrile + ethanol), and (acetonitrile + n-propanol). The experimental data and model parameters obtained will play an important role in the preparation and purification of edaravone.