MATLAB Simulink 韦伯尔分布拟合函数

以下是一个可以在 Simulink 中运行的 MATLAB 函数，该函数可以对输入的天数和失效率数据进行韦伯尔分布拟合：

function [params, goodness] = weibull_fit(days, failure_rates)
    % 韦伯尔分布拟合
    
    % 确保输入数据的长度相同
    assert(length(days) == length(failure_rates), '输入数据长度不同！');
    
    % 定义韦伯尔分布函数
    weibull_pdf = @(x, A, B) (B / A) * (x / A).^(B-1) .* exp(-(x / A).^B);
    
    % 初始参数猜测
    A0 = mean(days);
    B0 = 1 / var(log(failure_rates));
    
    % 最小二乘拟合
    [params, goodness] = fit(days, failure_rates, weibull_pdf, 'StartPoint', [A0, B0]);
end

该函数首先确保输入数据的长度相同，然后定义了韦伯尔分布函数。接着，函数使用平均天数和失效率的对数方差作为初始参数猜测，然后使用 MATLAB 的 fit 函数进行最小二乘拟合。函数返回拟合的参数和拟合优度。在 Simulink 中使用该函数时，将天数和失效率作为输入，将拟合参数和拟合优度作为输出即可。