Help required in the code

1 view (last 30 days)
Raja Awais Liaqait
Raja Awais Liaqait on 10 Oct 2019
Hello everyone. I want to change the values of global variable
"global min_realvar ;
global max_realvar" .
For that I want these variables to be assigned as inpput. How I can do this.
Please check the code:
clc
clear all
close all
%initializing global variables used in the model
global nodesize ;
global nreal ;
global nbin ;
global nbits ;
global nobj ;
global ncon ;
global ngen ;
%loadin the data
binomial_rail=xlsread('locationbinomialrail_deneme.xlsx');
binomial_road=xlsread('locationbinomialroad_deneme.xlsx');
binomial_sea=xlsread('locationbinomialsea_deneme.xlsx');
distance_rail=xlsread('locationdistancerail_deneme.xlsx');
distance_road=xlsread('locationdistanceroad_deneme.xlsx');
distance_sea=xlsread('locationdistancesea_deneme.xlsx');
capacity_rail=xlsread('CapacityRail_deneme.xlsx');
capacity_road=xlsread('CapacityRoad_deneme.xlsx');
capacity_sea=xlsread('CapacitySea_deneme.xlsx');
dataread1=xlsread('locationbinomialrail_deneme.xlsx','B1:B32'); %
dataread2=xlsread('locationbinomialroad_deneme.xlsx','C1:C32');
dataread3=xlsread('locationbinomialsea_deneme.xlsx','D1:D32');
dataread4=xlsread('locationdistancerail_deneme.xlsx','E1:E32');
dataread5=xlsread('locationdistanceroad_deneme.xlsx','F1:F32');
dataread6=xlsread('locationdistancesea_deneme.xlsx','G1:G32');
dataread7=xlsread('CapacityRail_deneme.xlsx','H1:H32');
dataread8=xlsread('CapacityRoad_deneme.xlsx','I1:I32');
dataread9=xlsread('CapacitySea_deneme.xlsx','J1:J32');
load_input_data('data.in');
%initializing the local constraints
global min_realvar ;
global max_realvar ;
nodesize = 31 ;
nreal = 3 ;
ngen = 400 ; %take chunk from excel it sh0uld be 800 if even a single entry, if reduce then it can be same.
min_realvar = min_realvar(1:nreal);
max_realvar = max_realvar(1:nreal);
global min_binvar ;
global max_binvar ;
nbin = 2;
nbits = [3;3];
min_binvar = [0;0];
max_binvar = [-50;50]; %if want to add colums or row then range will be increse by -100 to 100
Beta=100;
%opening time: 0
%Closing time: 100 mins
% no external library other than math lib which is internal and tool box is
% maths tool box.
BPR_constant_a=0.15;
%Ship Speed
v1=18.5;
v2=44;
F1=(0.1727*(v1)^2-0.217*v1)/24;
F2=(0.1727*(v2)^2-0.217*v2)/24;
if(nreal > 0)
yijkt = nreal + 1 : nreal + nobj ;
elseif(nbin > 0)
yijkt = sum(nbits) + 1 : sum(nbits) + nobj ;
end
obj_func = @input_data ;
if(nreal > 0)
pareto_decision_maker = zeros(nodesize, nreal + nobj + ncon + 3);
mixed_pop = zeros(2 * nodesize, nreal + nobj + ncon + 3);
elseif(nbin > 0)
pareto_decision_maker = zeros(nodesize, sum(nbits) + nobj + ncon + 3);
mixed_pop = zeros(2 * nodesize, sum(nbits) + nobj + ncon + 3);
end
%plotting the data
tic;
TTijk = mass_balance(90); % it will be 120.
fuel_consumption = integrity_and_non_negativity_condition(TTijk, obj_func)
TTijk = transportation_times_of_the_modes(TTijk);
for i = 2:ngen
pareto_decision_maker = same_leaving_arrival_transportation_mode(TTijk, pareto_decision_maker);
pareto_decision_maker(:,yijkt) = 0;
pareto_decision_maker = integrity_and_non_negativity_condition(pareto_decision_maker, obj_func);
mixed_pop = transported_via_at_least_one_transportation_mode(TTijk, pareto_decision_maker);
TTijk = demand_of_nodes(mixed_pop);
end
pareto_decision_maker = capacity_level_for_transportation_mode(pareto_decision_maker);
pareto_decision_maker;
TTijk;

Answers (0)

Community Treasure Hunt

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

Start Hunting!