In this paper, the bias‐dependent current–voltage (I–V) characteristics and their high‐order derivatives of GaAs pseudomorphic high electron mobility transistors (pHEMTs) have been modeled over a wide temperature range. To simulate these characteristics at different temperatures, the model is developed considering the dependence on the ambient temperature. It is the first time that the temperature‐dependent high‐order derivatives of I–V characteristics of pHEMT are predicted, which can guarantee their accuracy under different bias conditions. The artificial neural networks are employed with the temperature as one of the input variables. The validity of this model has been demonstrated by comparing the measured and modeled Ids and its derivatives (gm, gm2 and gm3, derived from the I–V characteristics numerically) of a GaAs pHEMT at different temperature range (250–400 K, with step of 50 K). The results show that the proposed model has a better agreement of high‐order derivatives than the popularly used Angelov model, especially for the third‐order derivative. Copyright © 2016 John Wiley & Sons, Ltd.