Frequency-space domain high-order modeling based on an average-derivative optimal method