We have carried out a theoretical study of double-δ-doped InAlAs/InGaAs/InP high electron mobility transistor (HEMT) by means of the finite differential method. The electronic states in the quantum well of the HEMT are calculated self-consistently. Instead of boundary conditions, initial conditions are used to solve the Poisson equation. The concentration of two-dimensional electron gas (2DEG) and its distribution in the HEMT have been obtained. By changing the doping density of upper and lower impurity layers we find that the 2DEG concentration confined in the channel is greatly affected by these two doping layers. But the electrons depleted by the Schottky contact are hardly affected by the lower impurity layer. It is only related to the doping density of upper impurity layer. This means that we can deal with the doping concentrations of the two impurity layers and optimize them separately. Considering the sheet concentration and the mobility of the electrons in the channel, the optimized doping densities are found to be 5 × 1012and 3 × 1012 cm-2 for the upper and lower impurity layers, respectively, in the double-δ-doped InAlAs/InGaAs/InP HEMTs.