%% 拟合sigmoid曲线
sigmoid = @(params, x) params(1) ./ (1 + exp(-params(2) .* (x - params(3)))) + params(4);
%params(1) 是斜率参数,params(2) 是增长速率参数,params(3) 是 x 值的偏移参数,params(4) 是 y 值的偏移参数。
initialGuess = [1, 1, 70, 0];%代表参数的初始猜测值,根据数据和先验知识进行调整
%它们分别对应斜率、增长速率、x偏移、y偏移
paramsFit = lsqcurvefit(sigmoid, initialGuess, xData, yData);
a = paramsFit(1); % 斜率参数
b = paramsFit(2); % 增长速率参数
c = paramsFit(3); % x 偏移参数
d = paramsFit(4); % y 偏移参数 %这些参数估计值将用于绘制拟合的 sigmoid 曲线或进行其他分析
xFit = linspace(min(xData), max(xData), 100); % 创建用于绘制曲线的 x 值范围
yFit = sigmoid(paramsFit, xFit); % 计算对应的 y 值
plot(xData,yData,'o','LineWidth',2.5);
hold on;
plot(xFit,yFit,':','LineWidth',2);
hold on;
% plot(xData, yData, 'o', xFit, yFit, ':')
% set(gca,'xtick',0:180:1)
%set(gca,'ytick',0:1:0.01)
ylim([0 1])
grid on;
legend('数据点', '拟合曲线');
数据
xData=[0,128.1/2-0.2,128.1/2-0.1,128.1+0.2,128.1/2+0.1,128.1];
yData=[0.25,0.25-0.02,0.25-0.01,0.75+0.01,0.75+0.02,0.75];