Using the thermal entangled state representation η|, we examine the master equation (ME) describing phase-sensitive reservoirs. We present the analytical expression of solution to the ME, i.e., the Kraus operator-sum representation of density operator ρ is given, and its normalization is also proved by using the IWOP technique. Further, by converting the characteristic function χ(λ) into an overlap between two "pure states" in enlarged Fock space, i.e., χ(λ) = η =−λ |ρ|η =0 , we consider time evolution of distribution functions, such as Wigner, Q-and P-function. As applications, the photon-count distribution and the evolution of Wigner function of photon-added coherent state are examined in phase-sensitive reservoirs. It is shown that the Wigner function has a negative value when κt 1 2 ln(1 + μ ∞ ) is satisfied, where μ ∞ depends on the squeezing parameter |M| 2 of environment, and increases as the increase of |M|.