Effacer les filtres
Effacer les filtres

GRID SEARCH FOR OPTIMISATION METHOD

11 vues (au cours des 30 derniers jours)
jason law
jason law le 5 Mai 2020
------------------------------------------------------------------------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);

Réponses (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by