Wavelets have been applied to many areas including image processing and numerical PDEs. In this talk, we discuss wavelet-based methods for numerically solving the elliptic interface problem and the Helmholtz equation. The elliptic interface problem has piecewise smooth diffusion coefficients and source terms which are discontinuous across interface curves. Hence, the exact solution $u$ has low regularity $u\in H^{1+\epsilon} (\Omega)$ with $\epsilon<1/2$. The solutions of Helmholtz equations have high frequencies leading to ill-conditioned indefinite linear systems for large wavenumbers. These impose challenges to traditional methods. All current known methods require complicated treatment along the interface and lack uniformly bounded condition numbers. First, we shall present some basics on wavelets in Sobolev spaces for the purpose of numerical PDEs. For 2D elliptic interface problems, then we discuss high-order wavelet methods for Helmholtz equations with or without interfaces. Finally, we present a second-order wavelet method with proved convergence and uniformly bounded condition numbers. This talk is based on joint work with M. Michelle.