Answered step by step
Verified Expert Solution
Link Copied!

Question

1 Approved Answer

I got all the tests failed with wrong answers: For example I know Ws = 2 . 1 5 5 7 , the equations for

I got all the tests failed with wrong answers: For example I know Ws =2.1557, the equations for p_state, c_util, p_denom, Lq. Please help me to correct the codes and make all teh tests passed. Here is the code: function [Ws, Wq, c_util, p_drop, p_state]= MMCQ(lambda, mu, c, Nwait)
rho = lambda /(c * mu); % Calculate the utilization factor
c_bar = rho *(1+(c * rho)^c / factorial(c)/(1- rho));
c_util = c_bar / c; % Calculate the utilization
% Calculate p0(probability of no customers in the system)
p0_denom = sum((c * rho).^(0:(c -1))/ factorial(0:(c -1)));
p0_denom = p0_denom +((c * rho)^c /(factorial(c)*(1- rho)))*(1-(rho^Nwait /(factorial(Nwait)*(1- rho))));
p0=1/ p0_denom;
p_state = zeros(1, Nwait +1);
% Calculate pN (probability of having N customers in the system)
for i =0:Nwait
if i <= c
p_state(i +1)=(rho^i / factorial(i))* p0;
else
p_state(i +1)=(rho^i /(factorial(c)* c^(i - c)))* p0;
end
end
pN = p_state(Nwait +1); % probability of having N customers in the system
% Calculate lambda_loss and lambda_eff
lambda_loss = lambda * pN;
lambda_eff = lambda - lambda_loss;
% Calculate Lq (average number of customers in the queue)
Lq = sum((0:Nwait).* p_state);
% Calculate Ls (average number of customers in the system)
Ls = Lq + lambda_eff / mu;
% Calculate p_drop (probability of a customer being dropped)
p_drop = pN /(1+ pN);
% Calculate Ws and Wq
Ws = Ls / lambda_eff;
Wq = Lq / lambda_eff;
end% Run this test case to check your code
c =2; % Number of servers
lambda =2; % Average Arrival Rate per minute
mu =0.1; % Service rate
Nwait =15;
% Call the MMCQ function with the given parameters
%[Ws, Wq, c_util, p_drop, p_state]= MMCQ(lambda, mu, c, Nwait);
[Ws, Wq, c_util, p_drop, p_state]= MMCQ(1,0.5,4,5);
% Define a relative error function
rel_error = @(x, y) abs(x - y)/ y;
% Validate the results and print "PASS" or "FAIL" for each metric
if rel_error(Ws,0.1557)<0.005
fprintf('Ws =%6.3f PASS
', Ws);
else
fprintf('Ws =%6.3f FAIL
', Ws);
end
if rel_error(Wq,0.1557)<0.005
fprintf('Wq =%6.3f PASS
', Wq);
else
fprintf('Wq =%6.3f FAIL
', Wq);
end
if rel_error(c_util, 0.4986)<0.005
fprintf('c_util =%6.3f PASS
', c_util);
else
fprintf('c_util =%6.3f FAIL
', c_util);
end
if rel_error(p_drop, 0.00272)<0.005
fprintf('p_drop =%8.5f PASS
', p_drop);
else
fprintf('p_drop =%8.5f FAIL
', p_drop);
end
fprintf('p_state: %s
', num2str(p_state));
fprintf('sum(p_state): %f
', sum(p_state));
% Extra Tests: p_state
% p_state should have Nwait+1 probabilities
% disp(['The probability of being in state ', num2str(Nwait),' is ', num2str(p_state(Nwait +1))]);
disp(['The probability of being in state ',' is ', num2str(p_state)]);
% Check that the probabilities sum to 1
if abs(sum(p_state)-1)<1e-5
disp('The probabilities sum to 1. Test PASSED.');
else
disp('The probabilities do not sum to 1. Test FAILED.');
end

Step by Step Solution

There are 3 Steps involved in it

Step: 1

blur-text-image

Get Instant Access with AI-Powered Solutions

See step-by-step solutions with expert insights and AI powered tools for academic success

Step: 2

blur-text-image

Step: 3

blur-text-image

Ace Your Homework with AI

Get the answers you need in no time with our AI-driven, step-by-step assistance

Get Started

Students also viewed these Databases questions