分岐限定アルゴリズムの再帰関数を実装しようとしています。アルゴリズムを再帰的に呼び出すたびに、再帰呼び出しで lb と ub の値を変更しています。エラーが表示されていますUB must be a real valued nx by 1 column vertor error in MATLAB
。私のコードは以下に添付されています:
function [cost] = BB(c, A, b, lb, ub, ctype, vartype, s, xmin, fmin, fid, NumberofNodes,min_cost, val);
[xmin_last, fmin_last, status] = glpk (c, A, b, lb, ub, ctype, vartype, s) %GLPK function
frac_value_1= find(xmin_last > 0.999);
xmin(frac_value_1)=1;
frac_value_2= find(xmin_last < 0.0001);
xmin(frac_value_2)=0;
frac_value=intersect(find (xmin_last > 0) , find( xmin_last < 1)) ; %positions of fractional variables
size(frac_value,1);
if (isempty(frac_value))
fprintf(fid,'no fractional value');
fprintf(fid,'\n\n');
cost = fmin_last
return
else
round_value=frac_value(1)
one_value_cost=find ( xmin_last== 1 ) %variables with value 1
zero_value_cost=find ( xmin_last==0) %variables with value 0
ub(zero_value_cost)=0 %changing ub for 0 values
%val=0
[cost_0] = BB(c, A, b, lb, ub(round_value)=0, ctype, vartype, s, xmin_last(round_value)=0, fmin_last, fid, NumberofNodes,min_cost,val=0);
%val=1
[cost_1] = BB(c, A, b, lb(round_value)=1, ub, ctype, vartype, s, xmin_last(round_value)=1, fmin_last, fid, NumberofNodes,min_cost, val=1);
end
min_cost = min(cost_0, cost_1)
endfunction;
と のエラーが表示されub(round_value)=0
ていlb(round_value)=1
ます。どんな助けでも大歓迎です