We can express this in matrix form by altering the value of the homogeneous parameter, h:
[ x ] [1 0 0 0] [x] [ y ] = [0 1 0 0] [y] [ z ] [0 0 1 0] [z] [z/d] [0 0 1/d 0] [1]