I got FMINCON requires the following inputs to be of data type double: 'A' error,How it can be solved????

% Load Solar Cell data
load elec_solar_iv_data.mat
% Display the Solar Cell model
Model = 'elec_solar';
open_system(Model);

% List of parameters and initial values prior to optimization ParsListMain = {'Is', 'Iph', 'ec', 'Rs', 'Rp'}; InitGuessMain = [ 3e-7 3.8 1.5 .004 10 ]; ParsListTemp = {'TIPH1', 'EG', 'TXIS1'}; InitGuessTemp = [ .001 1.11 3 ];

%Plot Data Versus Solar Cell Output Using Initial Parameters
 load_system(Model);
set_param([Model '/Solar Cell'], 'prm', '3')
Pars = reshape([ParsListMain; cellstr(num2str(InitGuessMain'))'],1,[]);
set_param([Model '/Solar Cell'], Pars{:})
Pars = reshape([ParsListTemp; cellstr(num2str(InitGuessTemp'))'],1,[]);
set_param([Model '/Solar Cell'], Pars{:})

% Generate preliminary model curves and plot afminconinst data
num_lines = length(iv_data);
v_model = cell(1, num_lines);
i_model = cell(1, num_lines);
legend_info_data  = cell(1, num_lines);
legend_info_model = cell(1, num_lines);
for idx_data = 1:num_lines
    sim(Model);
    v_model{idx_data} = Vo.signals.values;
    i_model{idx_data} = Io.signals.values;
    legend_info_data{idx_data}  = [ 'Temp = '                         ...
        num2str(iv_data(idx_data).temperature) '\circC, Data'];
    legend_info_model{idx_data} = [ 'Temp = '                         ...
        num2str(iv_data(idx_data).temperature) '\circC, Model'];
end
plot([iv_data.v], [iv_data.i], 'd', [v_model{:}], [i_model{:}])
xlabel('Solar cell output voltage (V)');
ylabel('Solar cell output current (A)');
legend([legend_info_data legend_info_model], 'Location', 'Best');
title('Model with Initial Parameter Values');

%Optimize Main Tab Dialog Parameters at Room Temperature (Step 1)
% Find room temperature data index
idx_data = find([iv_data.temperature]==25);%#ok

% Optimize parameters in main dialog tab of Solar Cell
ParsList = ParsListMain;
OptParsMain = fmincon(@elec_solar_lse, InitGuessMain,              ...
    optimset('TolX', 1e-3));

% Update Solar Cell block with optimized parameters
Pars = reshape([ParsList; cellstr(num2str(OptParsMain'))'],1,[]);
set_param([Model '/Solar Cell'], Pars{:});
% Display optimized parameters
display(sprintf(['Optimized parameters for the solar cell main '      ...
    'dialog tab are:\n']));
display(sprintf('\t%5s = %s\n', Pars{:}));

%Optimize Parameters Controlling Temperature Dependence (Step 2)
 % Find index into data for non-room temperatures
idx_data = find([iv_data.temperature]~=25);

% Optimize parameters in temperature dialog tab of Solar Cell
ParsList = ParsListTemp;
OptParsTemp = fmincon(@elec_solar_lse, InitGuessTemp,              ...
    optimset('TolX', 1e-3));

% Update Solar Cell block with optimized temperature parameters
Pars = reshape([ParsList; cellstr(num2str(OptParsTemp'))'],1,[]);
set_param([Model '/Solar Cell'], Pars{:});
% Display optimized parameters
display(sprintf(['Optimized parameters for the solar cell '           ...
    'temperature dialog tab are:\n']));
display(sprintf('\t%5s = %s\n', Pars{:}));

%Display Optimized Curves
  for idx_data = 1:num_lines
   sim(Model);
   v_model{idx_data} = Vo.signals.values;
   i_model{idx_data} = Io.signals.values;
end
plot([iv_data.v], [iv_data.i], 'd', [v_model{:}], [i_model{:}])
xlabel('Solar cell output voltage (V)');
ylabel('Solar cell output current (A)');
legend([legend_info_data legend_info_model], 'Location', 'Best');
title('Model with Optimized Parameter Values');

Answer 1

[Steven Lord]

OptParsTemp = fmincon(@elec_solar_lse, InitGuessTemp,              ...
  optimset('TolX', 1e-3));

You can't just skip required arguments. If you want to pass the options structure into fmincon you MUST specify the fun, x0, A, b, Aeq, beq, lb, ub, and nonlcon inputs first. You've just specified fun and x0, leaving fmincon to treat your options structure as the A input.

Since you don't appear to have any constraints, you probably don't want to use fmincon. Instead, use one of the unconstrained optimization functions. Their signatures are much shorter than that of fmincon, since they don't accept constraints.

[Ashutosh Sharma]

Thanks Steven it must be helpful!!!

[Ashutosh Sharma]

But if i specify all the values then also it gives the same error.... I already use fminsearch and fminunc but i need to apply fmincon for my work....can you please describe me more???

[Steven Lord]

You have two calls to fmincon in your code. Did you modify them both?

[Ashutosh Sharma]

Dear i assign matrix A of 2*2 as i saw in other example..... Can you give me example how to apply this....

[muhammad ilyas khattak khattak]

@Steven Lord...steven, I hope you are well.... I just have a similar type of error.. I am trying to resolve... but I have not figured out/solve uptill now using the contents in this link/discussion shared by you/others....Following are both the code and error message , if you can help/guide, ...Many thanks in advance....

t = 1;

% FCVEC04January2022V1(DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3)

p = optimproblem('ObjectiveSense', 'minimize');

[NumRSUs,MuRSU,MaxNumtk,E,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3] = FCVECInitializeParameters08November2021V5();

% StoreInitialization = struct(NumRSUs,MuRSU,MaxNumtk,E,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

global StoreInitialization;

StoreInitialization = struct('NumRSUs_struct',NumRSUs,'MuRSU_struct',MuRSU,'MaxNumtk_struct',MaxNumtk,'E_struct',E,'DPrc_struct',DPrc,'pc_struct',pc,'Muc_struct',Muc,'Mui_struct',Mui,'Muk_struct',Muk,'B_struct',B,'pi_struct',pi,'pm_struct',pm,'W_struct',W,'Wfog_struct',Wfog,'t_struct',t,'T_struct',T,'cdk_struct',cdk,'j_struct',j,'i_struct',i,'TasksperVehicle_struct',TasksperVehicle,'NTotal_struct',NTotal,'Numv_struct',Numv,'Numtk_struct',Numtk,'TijPrime_struct',TijPrime,'loop1_struct',loop1,'Fprime_struct',Fprime,'PhikPrime_struct',PhikPrime,'D_struct',D,'F_struct',F,'Ravailk_struct',Ravailk,'Phik_struct',Phik,'theta1_struct',theta1,'matrixforAlgo3_struct',matrixforAlgo3);

% NumRSUs,MuRSU,MaxNumtk,E,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3;

D_Cloud = zeros(1,Numtk);

F_Cloud = zeros(1,Numtk);

theta1 = Algorithm1_24December2021V8(Numv,Numtk,Ravailk,i,j,F,theta1);

N_Total = Numv;

% non-orthogonal multiple access (NOMA) technique, the communication resource is divided into resource

L = 10*N_Total; % blocks (RBs), expressed as B1 = {1,2,.....L} for uplink and

M = 10*N_Total; % B2 = {1,2,.....,M} for downlink.

[Tctotal,Ttotal,tTijtoRij,ttransij,texeij,tkij] = FCVEC4AllOutputs05January2022V3(D_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

% Though the following used algorithm 3 that assigns tasks to appropriate

% fog-vehicles with OBCU using "matrixforAlgo3" variable....

% However, the above function "FCVEC4AllOutputs05January2022V3" pre-calculate probable time spent by tasks

for iloopagn12 = 1:1:Numtk

if (theta1(1,iloopagn12) == 1) % Only those tasks that are to be offloaded to cloud-server are segregated....

D_Cloud(1,iloopagn12) = D(1,iloopagn12);

end

end

[ElapsedTime,Ttotal,matrixforAlgo3,MatTkAllocVehsLM,summation2] = Algorithm3_04January2022V6(NumRSUs,MuRSU,L,M,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3,Ttotal,tTijtoRij,ttransij,texeij,Tctotal);

% c(1,:) = ttransij + texeij - cdk;

% [c,ceq] = FCVEC04January2022V1constraints(DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3);

% [Tctotal,Ttotal,tTijtoRij,ttransij,texeij,tkij] = FCVEC04January2022V1(DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3);

% Ttotal = FCVEC04January2022V1(D_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3);

% The objective is to minimize the time that is spent during tasks exeuction at OnBoard Cmoputation Units (OBCU)at fog-vehicles or at cloud... The above is final equation

% equation (13), taken from paper titled "Mobile Vehicles as Fog Nodes for Latency Optimization in Smart Cities"......

% Likd mathswork for optimvar careation..."https://www.mathworks.com/help/optim/ug/optimvar.html"

tic

theta = optimvar('theta',Numtk','Type','integer','LowerBound',0,'UpperBound',1);

ttransij1 = optimvar('ttransij1',Numtk,'Type','continuous','LowerBound',min(ttransij),'UpperBound',max(ttransij));

texeij1 = optimvar('texeij1',Numtk,'Type','continuous','LowerBound',min(texeij),'UpperBound',max(texeij));

rng(0) % For reproducibility

% OptimizationConstraints

constr = optimconstr(Numtk,1);

% Desctiption Link for Constraints on mathswork website hsing optimproblem builtinfucntion.....

% https://www.mathworks.com/help/optim/ug/optimproblem.html

X0.theta = theta1;

X0.ttransij1 = ttransij;

X0.texeij1 = texeij;

% p = optimproblem('Objective',FCVEC04January2022V1(D_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3),'Constraints',ttransij + texeij <= cdk,'ObjectiveSense', 'minimize');

% p.Constraints.constrt1 = ttransij + texeij <= cdk;

% Using vectorized objective and contraint functions....decription/tutorial help from mathworks link given below......

% https://www.mathworks.com/help/gads/vectorizing-the-objective-and-constraint-functions.html

% D_Updated =

for (loop1 = 1:1:Numtk)

if(matrixforAlgo3(loop1,4) == 0)

D_Cloud(1,loop1) = D(1,loop1);

F_Cloud(1,loop1) = F(1,loop1);

end

end

% p.Objective = FCVEC06January2022V2(ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

% fun = @(ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3)FCVEC06January2022V2(ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

% fun = @FCVEC06January2022V2;

% fun = @(ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3)FCVEC06January2022V2(ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

% fun = @(X0,StoreInitialization,ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3)FCVEC06January2022V2;

fun = @(x)FCVEC06January2022V2(x,ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

% fun = @(x0)FCVEC10January2022V3;

%ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3

% fun = @()FCVEC06January2022V2(ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

% ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3

% p.Constraints.cons1 = ttransij + texeij <= cdk;

% How to make constrains link of mathworks regarding making array of contraints using "optimconstr" builtin function.....

% https://www.mathworks.com/help/optim/ug/optimconstr.html

% for loop1 = 1:Numtk

% % p.Constraints.cons1 = ttransij + texeij <= cdk;

% constr(loop1,1) = ttransij1(loop1,1) + texeij1(loop1,1) <= cdk;

% end

% fmincon problem based optimiztion....

% https://www.mathworks.com/help/optim/ug/fmincon.html#d123e88857

A = [];

b = [];

Aeq = [];

beq = [];

% beq = constr;

lb = [];

ub = [];

% x0 = X0;

% x0 = theta1;

% function nlcons(Numtk,constr,ttransij1,texeij1,cdk)

% nonlcon =

% @(Numtk,constr,ttransij1,texeij1,cdk)nlcons(Numtk,constr,ttransij1,texeij1,cdk);\

% nonlcon = @(Numtk,constr,ttransij1,texeij1,cdk)nlcons;

x = theta;

ncon = @(x)nlcons(x,Numtk,constr,ttransij1,texeij1,cdk);

% Numtk,constr,ttransij1,texeij1,cdk

% Numtk,constr,ttransij1,texeij1,cdk

% p.nonlcon = nonlcon1;

% p.Constraints.constall = constr;

% [sol,fval] = solve(p,x0);

options = optimoptions(@fmincon,'Algorithm','interior-point',...

'SpecifyObjectiveGradient',true,'SpecifyConstraintGradient',true,...

'PlotFcn',{@optimplotx,@optimplotfval,@optimplotfirstorderopt});

% p.solver = 'fmincon';

% ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3

[Ttotal,theta,Tftotal,Tctotal] = FCVEC11January2022V4(X0,StoreInitialization,ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

TtotalHandle = @(theta,Tftotal,Tctotal)sum(theta(:,1).*Tftotal + (1 - theta(:,1)).*Tctotal','all'); % Equation (13)

% x = fmincon(fun(x0,StoreInitialization,ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3),x0,A,b,Aeq,beq,lb,ub,nonlcon(Numtk,constr,ttransij1,texeij1,cdk),options);

% x = fmincon(TtotalHandle(theta,Tftotal,Tctotal),X0,A,b,Aeq,beq,lb,ub,nonlcon(Numtk,constr,ttransij1,texeij1,cdk),options);

% x = fmincon(TtotalHandle(theta,Tftotal,Tctotal),X0,A,b,Aeq,beq,lb,ub,nonlcon(Numtk,constr,ttransij1,texeij1,cdk));

%Numtk,constr,ttransij1,texeij1,cdk

%ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3

x = fmincon(fun, X0, A, b, Aeq, beq, lb, ub, ncon, options);

ElapsedTime3 = toc;

Values4Graph = fval+ElapsedTime3;

%%%%%%%%%%%%%%%%%%%%% The FCVEC06January2022V2 function code is%%%%%%%%%

% function [i,j,texe2bepassed,ttrans2bepassed,Ttotal] = FCVEC12July2021(TasksperVehicle,Position,N_Total,t,theta,Phik)

% function [summoftexettrans] = FCVEC4Constraints14July2021(theta,Phik)

% function [Tctotal,Ttotal,tTijtoRij,ttransij,texeij,tkij] = FCVEC04January2022V1(DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3)

function Ttotal = FCVEC06January2022V2(x0,StoreInitialization,ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3)

%ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3

% DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta,matrixforAlgo3

% D,F,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,Fprime,t,T,cdk,theta,Phik,PhikPrime,j,i,TasksperVehicle,NTotal

% StoreInitialization.D_struct = D; % Darray = zeros(1,j);

% ttrans = zeros(NTotal,TasksperVehicle);

% texe = zeros(NTotal,TasksperVehicle);

% I = zeros(NTotal,TasksperVehicle);

W = 20*10^6; % Bandwidth of wireless channel = 20 MHz....

%%%%%%%%%%Task Offloading From MDs to RSU %%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%

% γ(i,j) = P(i)*σ2 + I(i,j); % Equation (1) % Suppose that task Tij is offloaded to R via a selected channel.

% Then the signal-to-interference-plus-noise-ratio (SINR) denoted as Rho at RSU from MDi can be defined as in Equation (1)....

% P(i) = 80*10^(-3); % P(i) = 80 mW = Pi;

Pi = 80*10^(-3); % Whereas "i" is the "ith" Mobile Device (MD)....

% Constant for all "i"

Prt = 120*10^(-3); % Whereas "i" is the "ith" Mobile Device (MD)....

% % Prt = 120 mW = Pr(t) = PR;

% Constant for all "t"

%%%%%% The following is website that specifies inter-conversion between dBm and mW....

% https://www.electronicdesign.com/technologies/communications/article/21796484/understanding-wireless-range-calculations

% P(dBm) = 10*log10(P(mW))

PidBm = 10*log10(Pi);

PrtdBm = 10*log10(Prt);

sigma2 = -70; % σ2 = the noise power = -70 dBm = sigma2

sigmafSquare = -80; % -80 dBm; "sigmafSquare" is the noise power in (VFC) vehicle fog computing...

% I(i,j) = -75; % Interference in dBm... same for all the nodes and vehicles...

Iij = -75; % Whereas "i" is the "ith" Mobile Device (MD)....and "j" is the "jth" task....

% Iij = I(i,j) = -75 dBm...since its constant for all values...

Ifij = -70; % Ifij = If(i,j) = -70 dBm...since its constant for all values... and specific for fog nodes interference...

% Rho(i,j) = P(i)*sigma.^2 + I(i,j); % Equation (1)

% Rho(i,j) = P(i)*sigma2 + I(i,j); % Equation (1)

% Rho(i,j) = PidBm/(sigma2 + Iij); % where P(i) denotes MDi’ transmission power and it is a constant.

Rhoij = PidBm/(sigma2 + Iij); % Iij is the interference caused by other MDs using the

% same channel as MDi. sigma2 or σ2 is the noise power. ...It is all the values being constant therefore,

% Rhoij , that is the interference is same for all the MDs....

% r(i,j) = W*log2(1 + Rho(i,j)); % Equation (2)

rij = W*log2(1 + Rhoij); % Then transmission rate (r) from MDi to R for an arbitrary wireless channel can be

% expressed as in Equation (2):

% tTijtoR(i,j) = D(i,j)/r(i,j); % Equation (3)

% tTijtoR(i,j) = D(i,j)/rij; % Accordingly, the time needed to offload Tij to R is Equation (3):

% tTijtoR(i,j) = D/rij; % Unlike other works which focus on energy optimization at MDs

tTijtoRij = D/rij; % by determining whether or not to offload the tasks, we in this

% paper aim to optimize the response latency of the tasks already

% offloaded to VFC, under an assumption that the tasks to be

% offloaded can be determined by MDs.

%%%%%%%%%%%%%%Task Scheduling Among Fog Vehicles%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%

% rtrans(i,j) = Wfog*log2(1 + Pr(t)/(sigmaf.^2+ If(i,j))); % Equation(4)

% rtrans(i,j) = Wfog*log2(1 + PrtdBm/(sigmafSquare + Ifij)); % Equation(4)

rtransij = Wfog*log2(1 + PrtdBm/(sigmafSquare + Ifij)); % Equation(4)

% rtrans(i,j) = rtransij...

% Similar to equation (2), the transmission rate from

% RSU to the fog vehicle via an arbitrary wireless channel can be

% expressed as in Equation (4)....

% Ifij is the interference caused by other task transmission using

% the same channel as Tij.

% where Pr(t) denotes RSU’s transmission power at tth time slot.

% Thus,ttrans ij can be expressed as in Equation (5):

% ttrans(i,j) = D(i,j)*rtrans(i,j); % Equation (5)

% ttrans(i,j) = D(i,j)/rtransij;

% ttrans(i,j) = D/rtransij;

% ttransij1 = D/rtransij; % Equation (5) % Similarly, the execution time texeij

% can be expressed as in subsequent Equation (6)..

% Equation (5) is substitted by

% importing pre-computed "ttransij1" from "main" function...

% ttransij1,texeij1

% ttrans2bepassed = ttrans(i,j);

ttrans2bepassed = ttransij1;

% texe(i,j) = F(i,j)/Muk; % Equation (6)

% texeij = F(i,j)/Muk; % Equation (6)

% texeij = F/Muk; % Where μk represents the process rate of kth vehicle to which

% Tij is offloaded. To calculate twqij , we need to know all the tasks

% allocated to the same vehicle prior to Tij. Let wqkij(t) denote

% the set of tasks allocated to the kth vehicle prior to Tij . Then

% the waiting time of Tij is the total execution time of tasks in wqkij(t),

% and thus obrains the equation (7):

% Equation (7) is substitted by

% importing pre-computed "texeij1" from "main" function...

texe2bepassed = texeij1;

% for loop3 = 1:T

if (size(texe2bepassed,2) < size(ttrans2bepassed,2))

IndxSub = size(ttrans2bepassed,2)-size(texe2bepassed,2);

xextzpd1(1,:) = wextend('1D','zpd',texe2bepassed(1,:),IndxSub,'r'); % Extending the "texe2bepassed" with zeros...

texe2bepassed = xextzpd1; % with length of extension equal to IndxSub...

end % to match it with "ttrans2bepassed".

if (size(texe2bepassed,2) > size(ttrans2bepassed,2))

IndxSub = size(texe2bepassed,2)-size(ttrans2bepassed,2);

xextzpd1(1,:) = wextend('1D','zpd',ttrans2bepassed(1,:),IndxSub,'r'); % Extending the "texe2bepassed" with zeros...

ttrans2bepassed = xextzpd1; % with length of extension equal to IndxSub...

end

% ...

% texeij = xextzpd1;

% texe2bepassed = xextzpd1;

summoftexettrans = texe2bepassed + ttrans2bepassed;

% if (T(i,j) wqkij(t))

% tkij(t) =sum(F(i,j)/Muk,'all'); % Equation (7)

% Fprime = randi([100 1000],1,j);

% PhikPrime = Phik(i,1:j);

% Fprime = PhikPrime.*Fprime;

% if (t == 1)

% tkij(t) = sum(F(t,:),'all')/Muk; % Equation (7)

% end

% if (t > 1)

% intermediaeSum = sum(Fprime(1:t-1,:),'All');

% tkij(t) = intermediaeSum/Muk; % Equation (7)

% end

% tkij(t) = Fprime/Muk; % Equation (7)

if (t == 1)

% if (size(Fprime,2) ~= size(PhikPrime,2))

% IndxSubAgn = size(PhikPrime,2)-size(Fprime,2);

% xextzpd1 = wextend('1D','zpd',Fprime(1,:),IndxSubAgn,'r'); % Extending the "texe2bepassed" with zeros...

% Fprime = xextzpd1; % with length of extension equal to IndxSub...

% end % to match it with "ttrans2bepassed"....

% Fprime = PhikPrime.*Fprime; % PhikPrime is previously assigned fog vehicles, "i" to tasks, "j"

% intermediaeSum = sum(Fprime(:,:),2);

% tkij = intermediaeSum/Muk; % Equation (7)

tkij = 0; % Equation (7)

end

if (t ~= 1)

% if (size(Fprime,2) ~= size(PhikPrime,2))

% IndxSubAgn = size(PhikPrime,2)-size(Fprime,2);

% xextzpd1 = wextend('1D','zpd',Fprime(1,:),IndxSubAgn,'r'); % Extending the "texe2bepassed" with zeros...

% Fprime = xextzpd1; % with length of extension equal to IndxSub...

% end

Fprime = PhikPrime.*Fprime; % PhikPrime is previously assigned fog vehicles, "i" to tasks, "j"

intermediaeSum = sum(Fprime(:,:),2);

tkij = intermediaeSum/Muk; % Equation (7)

end

% Given all the arrived tasks at the tth time slot, the total response

% latency Tf total(t) thus can be expressed as Equastion (8)

[Steven Lord]

That's a lot of unformatted code to read through. Could you please do one or both of the following:

1. Edit your post, select all your code, and press the first button in the Code section of the toolstrip to format it as MATLAB code.
2. Simplify your code down to no more than say 10-15 lines showing the fmincon call that throws the error and how the input arguments to that fmincon call are constructed. If possible run the code to just before the fmincon call and show us the output of whos at that point (which will tell us the size and type of those input argument.)

Including the full and exact text (all the text displayed in red and/or orange) of the error and/or warning messages (or stating that it is exactly the same as the error from the message that started this discussion) would also be useful. Sometimes small differences in the text of an error message can suggest what caused MATLAB to issue that error.

[muhammad ilyas khattak khattak]

@Steven Lord....Steven, Thanks for you tremendous effort and time .... I am lacking words to thank you.... I had to add to the code the functions that are used ....Therefore, it could not be limited to 10 or 15 lines only..... Your possible comments/guidance requested so that I can get the function minimized using defined scenario...Thanks..

[Steven Lord]

When I copied and pasted this into MATLAB Editor I received a Code Analyzer error that one of your functions was not ended with an end. If any of the functions in your file end with an end keyword then all of the functions in your file must end with end.

When I added the end where I think it should be located, I then received an error.

Output argument "Tctotal" (and possibly others) not assigned a value in the execution with
"example410972>FCVEC4AllOutputs12January2022V3" function.

Error in example410972 (line 17)
[Tctotal,Ttotal,tTijtoRij,ttransij,texeij,tkij] = FCVEC4AllOutputs12January2022V3(x0,StoreInitialization,ttransij,texeij,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

That function doesn't appear to define that variable at all.

I'm not sure I'm going to have the time to read through and debug the 350 lines of code you have here further. One suggestion I do have to increase the readability of this code is to reduce the number of input and output arguments your functions accept by having them return and accept struct arrays. Instead of something like:

function [NumRSUs,MuRSU,MaxNumtk,E,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3] = FCVECInitializeParameters08November2021V5()

which has 33 output arguments (three of which shadow built-in functions in MATLAB: pi, i, and j) I'd have it return a struct:

function params = FCVECInitializeParameters08November2021V5()

Then whenever you need (for example) Wfog you could access it as params.Wfog or (if you're going to access it many times) you could define a local variable for that specific parameter to save on typing.

Wfog = params.Wfog;

The same holds for packing the additional parameters you want to pass into your objective function into a struct, so instead of your fun being:

fun = @(x)FCVEC06January2022V2(x,ttransij1,texeij1,theta,D_Cloud,F_Cloud,DPrc,pc,Muc,Mui,Muk,B,pi,pm,W,Wfog,t,T,cdk,j,i,TasksperVehicle,NTotal,Numv,Numtk,TijPrime,loop1,Fprime,PhikPrime,D,F,Ravailk,Phik,theta1,matrixforAlgo3);

it could be something more like:

param2 = struct(...);
fun = @(x) FCVEC06January2022V2(x, param2);

[muhammad ilyas khattak khattak]

@Steven Lord🙏☺️☺️🙂...Thanks steven once again...I exatly did by encapsulating the parameters into a structure and then calling "fmincon"...Below are the final fmincon commands/statments.... Also in the attachment in pdf format (convertable/can be copied to .m format) , is the full code that contains all the helping functions/main functions in case it is needed to run.....If I get successful execution/obtain minimization of what I am doing , a great part of my work will be solved... Therefore, I am trying to contact again and again..... I am literarly lacking words to say thanks... THanks once again in advance, even if you spare a little time, it is much for me..