How to iterate to find stable controller/feedback gains in closed loop system and performance metrics

Background: I have the identified transfer function of a plant for an experimental setup that I have. It is in the form of a vector containing the real and imaginary parts for the frequency response. I also have the transfer function for the feedback and controller that is controlling this plant. The controller has both feedforward and feedback gains. I could use the identification toolbox to fit a system to the frequency response data, however for now I prefer to not do so since I want to stay with measured data as much as possible.

Question: 1) I'm wondering if there is a built in function to matlab to iterate through all the combinations of gains in the controller and feedback to first figure out when the system is stable. Then iterate to minimize steady state error. Then iterate to minimize settling time. Etc..

2) If this would be easier using a fitted plant. Please let me know. I have an 8 order plant that seems to fit ok, but I want to explore other options before using this.

Objective: Right now my system is stable, but the gains were found through testing, I pretty much want to come up with a table that shows the settling time, overshoot, and steady state error for various gain inputs.