GRID SEARCH FOR OPTIMISATION METHOD

------------------------------------------------------------------------My question is how to implant this procedure into my code-----------------------------------------------------------------------------

Hi, I am currently dealing with the paper "optimal signal processing FSCW radar data"

I need to find the best combination of t to get minimum of J

Here is the procedure to find the t:

1. I have 2 number of t

t=[t1 t2]

2. t will start with any value but greater than zero

3. t will insert to an algorithm to find Jminimum

4. t2 is held constant while the minimum of J is found along the t1 dimension.

5. t1 is held constant at the value previously found while the minimum of J is found along the t2 dimension.

6. 4and5 repeated until the time delay values are no longer updated.

------------------------------------------------------------------------My question is how to implant this procedure into my code-----------------------------------------------------------------------------

Below is my coding:

clear
clc
close all
load('C:\Users\User\Desktop\fyp\defect_free(H).mat');
load('C:\Users\User\Desktop\fyp\f_column.mat');
load('C:\Users\User\Desktop\fyp\fi_column.mat');
 load('C:\Users\User\Desktop\fyp\fr_column.mat');
fi=fi_column;
fr=fr_column;
f=f_column';
abs_H=abs(ifft(H));% getting magnitude of H from data in frequency domain
abs_H_temp=abs_H(1:2,:);% select from row 1 to row 2
 for z=1:100000; % number of iterations


t=[rand(1,1) rand(1,1)]; %generate set of t

%algorithm
for x=1:101;   % number of  frequency
    for y=1:n; %number of t
        Fr(x,y)=cos(2*pi*f(x)*t(y)); %real number of F
        Fi(x,y)=-i*sin(2*pi*f(x)*t(y)); % imaginary number of F
    end 
end
   Fr_T=Fr';
   Fi_T=Fi';
   A=(Fr_T*Fr)+(Fi_T*Fi);
   B=(Fr_T*fr)+(Fi_T*fi);
   a=lsqminnorm(A,B) ;    % solve tp get a set of a
J(z)=sum((a-abs_H_temp))^2; % to find the minimum of J from measured data
 end
J_min=min(J);