Using the formula of Wigner function in coherent representation, we have obtained the analytical expression for Wigner function of N00N state. Based on phase space method, we study the quantum interference with N00N state as input. We derive the analytical expression of conditional probability related with the input parameter N and phase parameter φ and analyze it numerically. It is shown that, when φ is 0 or π, the output is just N00N state. It is also shown that, for 2002 state as input, the output must be 2002 state, which is independent of phase parameters. Moreover, as the number of input photon N increases, the phase probability distributions remain to have one, two, three and four peaks and get narrower. All these results can offer theoretical reference for experiments.