In this paper, we consider an insurance risk model with mixed premium income, in which both constant premium income and stochastic premium income are considered. We assume that the stochastic premium income process follows a compound Poisson process and the premium sizes are exponentially distributed. A new method for estimating the expected discounted penalty function by Fourier-cosine series expansion is proposed. We show that the estimation is easily computed, and it has a fast convergence rate. Some numerical examples are also provided to show the good properties of the estimation when the sample size is finite.