Scribd
Upload
3 views

Objective_fun_matlab

Uploaded by

drbaskerphd
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
3 views

Objective_fun_matlab

Uploaded by

drbaskerphd
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

% Define parameters and weights

H_L = 1000; % Height on the left


H_R = 500; % Height on the right
R = 7; % Right component size
alpha = 0.5; % Weight for total height
beta = 0.3; % Weight for storage cost
gamma = 0.2; % Weight for imbalance penalty
delta = 0.1; % Weight for aging penalty
lambda = [0.05, 0.05, 0.05, 0.05]; % Weights for security constraints

% Define functions for Aging Penalty and Security Constraints


A = @(p) exp(-p); % Example aging penalty function
C1 = @(p) sin(p); % Example security constraint
C2 = @(p) cos(p);
C3 = @(p) p.^2;
C4 = @(p) sqrt(p);

% Define the storage cost function (choose one)


storage_cost = @(p, L) log(p .* L + 1); % Logarithmic
% storage_cost = @(p, L) p; % Linear
% storage_cost = @(p, L) 1; % Constant

% Define the objective function F(p, L)


F = @(p, L) abs(H_L - H_R) + ...
alpha * (H_L + H_R) + ...
beta * storage_cost(p, L) + ...
gamma * abs((1 - p) .* L - R) + ...
delta * A(p) + ...
lambda(1) * C1(p) + lambda(2) * C2(p) + lambda(3) * C3(p) +
lambda(4) * C4(p);

% Vary p and L to evaluate F(p, L)


p_values = linspace(0.01, 1, 50); % Varying p from 0.01 to 1
L_values = linspace(1, 20, 50); % Varying L from 1 to 20
[P, L_grid] = meshgrid(p_values, L_values); % Create grid for p and L

% Evaluate F(p, L) on the grid


F_values = arrayfun(F, P, L_grid);

% 3D Surface Plot
figure;
subplot(1, 2, 1); % Create a subplot for the 3D surface plot
surf(P, L_grid, F_values, 'EdgeColor', 'none');
xlabel('\bf p', 'FontSize', 14); % Bold label and larger font
ylabel('\bf L', 'FontSize', 14); % Bold label and larger font
zlabel('\bf F(p, L)', 'FontSize', 14); % Bold label and larger font
title('\bf 3D Surface Plot of F(p, L)', 'FontSize', 16); % Bold title and
larger font
set(gca, 'FontSize', 12); % Larger axis markings
colorbar;
grid on;

% Contour Plot
subplot(1, 2, 2); % Create a subplot for the contour plot
contourf(P, L_grid, F_values, 20, 'LineColor', 'none'); % Filled contour
plot
xlabel('\bf p', 'FontSize', 14); % Bold label and larger font
ylabel('\bf L', 'FontSize', 14); % Bold label and larger font
title('\bf Contour Plot of F(p, L)', 'FontSize', 16); % Bold title and
larger font
set(gca, 'FontSize', 12); % Larger axis markings
colorbar;
grid on;

You might also like