You can specify the structure of a transfer function to estimate (the number of poles and zeros), and can specify that say all numerator coefficients should be positive. This will the absence of non-minimum phase zeros.
See this doc page . Search for the section "specify estimated transfer function model structure and coefficient constraints"
