用matlab写一个可以在simulink中运行的函数该函数可以对输入的天数和失效率数据进行y = mnxn^m-1exp-xn^m拟合输出mn可以将函数转换后进行线性拟合

function [m, n] = fit_data(days, failure_rate)

% 将天数和失效率数据转换为列向量 days = days(:); failure_rate = failure_rate(:);

% 对失效率数据取自然对数 ln_failure_rate = log(failure_rate);

% 进行线性拟合 P = polyfit(days, ln_failure_rate, 1);

% 计算m和n m = 1/P(1); n = exp(-P(2)/P(1));

end