MATLAB Weibull Reliability Fit Function for Simulink

function [y] = reliability_fit(days, failure_rates) % This function provides a fit for reliability data using the Weibull distribution % Inputs: % days: a vector of time intervals for which the failure rates have been measured % failure_rates: a vector of corresponding failure rates % Output: % y: a vector of predicted reliability values based on the fit

% Define the Weibull distribution function weibull = @(x, m, n) (m/n)*(x/n)^(m-1)*exp(-(x/n)^m);

% Define the objective function to minimize the sum of squared errors objective = @(p) sum((failure_rates - weibull(days, p(1), p(2))).^2);

% Set initial parameter values and bounds for optimization p0 = [1, 1]; lb = [0, 0];

% Perform optimization using fmincon options = optimoptions('fmincon', 'Algorithm', 'interior-point', 'Display', 'off'); p = fmincon(objective, p0, [], [], [], [], lb, [], [], options);

% Generate predicted reliability values based on the optimized parameters y = weibull(days, p(1), p(2));

end