query about optimizing all pass filter parameters

I have a function H(f) = 10 –(α*L/10)exp –j[ D*λ0^2*L*ω^2/ (2*c)] .
And θap(phase of all pass filter in frequency domain) is a function of T, r, ω0 and f. Expression of θap is known.
I have to optimize all pass filter parameters T, r, ω0 such that minimum mean square condition is achieved

MMSE( T, r, ω0 )= ∑ { arg( 1/ H(f) ) – θap( T, r, ω0, f) }^2

[Note: Summation is across f.]

All the input values α, L, D, λ0, ω ,c are known except T, r, ω0 .
I have only bound constraints for r (0
Can anybody tell me which optimization function will be appropriate for optimizing these values.
I’ll really be grateful if anybody could provide me Matlab code for this. Start a New Thread