A method for automatically tuning fully cross-coupled multivariable PID controllers from decentralized relay feedback is developed with new techniques for process frequency-response matrix estimation and multivariable decoupling design.