k = size(n_1);
t = k* rt ;
alt_20hz= lat_01(1,1);
off_nadir_angle_wf_20hz_ku = F_off_nadir_angle_wf_ocean_01_ku(1,1) ;
c = 299792458 ;
Re = 6371000 ;
antenna_w = rad2deg(1.28) ;
rt = 3.125 * 10.^-6 ;
sigma_p = 0.513 * rt ;
lamda = (sin(antenna_w).^ 2 )* (1 / ( 2 * log(2)));
a = (4 * c ) / (lamda .* alt_20hz * (1+ (alt_20hz / Re)));
b_offnadirangel = cos( 2 .* sqrt(off_nadir_angle_wf_20hz_ku)) - (( sin(2 * off_nadir_angle_wf_20hz_ku)) / lamda ) ;
c_offnadirangel = b_offnadirangel .* a ;
a_offnadirangel = exp(-4* sin(off_nadir_angle_wf_20hz_ku))./ (lamda);
MaxIter =600 ;
f1 = n_1 ;
fun = @(P)sum(f1 - (a_offnadirangel .* P(1).^(( 1 + erf((t - P(2) - c_offnadirangel .* P(3).^2 )/...
( sqrt(2) .* P(3))) ) / 2) .* exp(-1 .*(c_offnadirangel .* ( t - P(2) - (0.5 .* c_offnadirangel .* P(3).^2)))) + Tn)).^2
options = optimset;
options = optimset(options,'Display', 'off');
options = optimset(options,'MaxIter', 600);
options = optimset(options,'MaxFunEvals', 600);
options = optimset(options,'TolFun', 10^-10);
options = optimset(options,'TolX', 10^-10);
x0 =[0,0,0];
MaxIter = 600 ;
Tol_fun = 10^-10 ;
while c <= MaxIter
P = fminsearch(fun,x0,options);
if abs(fun) < Tol_fun
break;
else
stop_gate = stop_gate +1
n_1 = waveform_norm(edge_foot(1):stop_gate(1))' ;
P = fminsearch(fun, x0, options);
end
c=c+1 ;
end
Unable to perform assignment because the size of the left side is 1-by-1 and the size of the right side is 1-by-2.
Error in fminsearch (line 200)
fv(:,1) = funfcn(x,varargin{:});
