Quasi-phasematching design to provide an approximation to a desired spectral amplitude response A(f) is provided. An initial phase response φ(f) corresponding to A(f) is generated. Preferably, d2φ(f)/df2 is proportional to A2(f). Alternatively, φ(f) can be a polynomial in f. A function h(x) is computed such that h(x) and H(f)=A(f)exp(iφ(f)) are a Fourier transform pair. A domain pattern function d(x) can be computed by binarizing h(x) (i.e., approximating h(x) with a constant-amplitude approximation). In some cases, the response provided by this d(x) is sufficiently close to A(f) that no further design work is necessary. In some embodiments of the invention, the need for binarization can be reduced or eliminated by providing amplitude modulation of the effective nonlinearity.