曲线拟合与分布拟合

打开实时脚本

此示例显示了如何执行曲线拟合和分布拟合，并讨论了每种方法何时适用。

在曲线拟合和分布拟合之间进行选择

曲线拟合和分布拟合是不同类型的数据分析。

当您想要将响应变量建模为预测变量的函数时，请使用曲线拟合。
当你想要对单个变量的概率分布建模时，使用分布拟合。

曲线拟合

在以下实验数据中，预测变量为时间，即服用药物后的时间。响应变量为浓缩的，血液中药物的浓度。假设只有响应数据浓缩的受实验误差的影响。

时间=[0.1 0.1 0.3 0.3 1.3 1.7 2.1 2.6 3.9 3.9．..5.1 5.6 6.2 6.4 7.7 8.1 8.2 8.9 9.0 9.5．..9.6 10.2 10.3 10.8 11.2 11.2 11.2 11.7 12.1 12.3．..12.3 13.1 13.2 13.4 13.7 14.0 14.3 15.4 16.1 16.1．..'; ';Conc = [0.01 0.08 0.13 0.16 0.55 0.90 1.11 1.62 1.79 1.59 .．..1.83 1.68 2.09 2.17 2.66 2.08 2.26 1.65 1.70 2.39．..2.08 2.02 1.65 1.96 1.91 1.30 1.62 1.57 1.32 1.56．..1.36 1.05 1.29 1.32 1.20 1.10 0.88 0.63 0.69 0.69．..0.49 0.53 0.42 0.48 0.41 0.27 0.36 0.33 0.17 0.20]';

假设你想建立血液浓度与时间的函数模型。情节浓缩的反对时间．

情节(时间、浓缩的“o”）;包含(“时间”）;ylabel (的血药浓度）;

图中包含一个轴对象。axis对象包含一个类型为line的对象。

假设浓缩的的函数为双参数威布尔曲线时间．威布尔曲线有其形式和参数

$y ＝ c （ x / 一个)^{（ b - 1)} e^{- （ x / 一个)^{b}} ，$

在哪里 $一个$ 是水平缩放， $b$ 是形状参数，并且 $c$ 是垂直缩放。

采用非线性最小二乘拟合威布尔模型。

modelFun = @ (p, x) p (3) . * (x / p (1)) ^ (p(2) 1)。* exp (- (x / p(1)) ^(2)页);startingVals = [10 2 5];nlModel = fitnlm(时间、浓缩的modelFun startingVals);

在数据上绘制威布尔曲线。

xgrid=linspace（0,20100）'；line（xgrid，predict（nlModel，xgrid），“颜色”，“r”）;

图中包含一个轴对象。轴对象包含两个类型为line的对象。

拟合的威布尔模型是有问题的。fitnlm假设实验误差是可加的，并且来自于具有常数方差的对称分布。然而，散点图显示，误差方差与曲线的高度成正比。此外，可加性对称误差暗示阴性血药浓度测量是可能的。

更现实的假设是乘法误差在对数尺度上是对称的。在该假设下，通过取两边的对数，将威布尔曲线拟合到数据上。使用非线性最小二乘法拟合曲线：

$日志（ y) ＝日志（ c) + （ b - 1) 日志（ x / 一个) - （ x / 一个)^{b} ．$

nlModel2=fitnlm（时间，对数（conc），@（p，x）对数（modelFun（p，x）），起始VAL）；

将新曲线添加到现有绘图中。

线(xgrid exp(预测(nlModel2 xgrid)),“颜色”(0 5 0),“线条样式”，'--'）;传奇({“原始数据”，“添加剂错误模式”，“乘法错误模式”}）;

图中包含一个轴对象。轴对象包含3个类型为line的对象。这些对象分别代表原始数据、加性误差模型、乘性误差模型。

模型对象nlModel2包含精度估计。最佳做法是检查模型的拟合优度。例如，在对数标度上绘制残差图，以检查乘法误差的恒定方差假设。

在本例中，使用乘法误差模型对模型预测几乎没有影响。有关模型类型影响更大的示例，请参阅将非线性模型转换为线性拟合的陷阱．

曲线拟合函数

Statistics and Machine Learning Toolbox™包含以下用于拟合模型的函数:fitnlm对于非线性最小二乘模型，fitglm对于广义线性模型，菲特尔格普对于高斯过程回归模型，和菲特斯姆支持向量金宝app机回归模型。
曲线拟合工具箱™ 提供简化曲线拟合任务的命令行和图形工具。例如，“工具箱”提供各种模型的起始系数值的自动选择，以及稳健和非参数拟合方法。
优化工具箱™ 具有执行复杂类型的曲线拟合分析的功能，例如分析带有系数约束的模型。
MATLAB®函数polyfit拟合多项式模型，以及MATLAB函数fminsearch在其他类型的曲线拟合中也很有用。

配电装置

假设您要对电气元件寿命的分布进行建模生活测量50个相同的电子元件失效的时间。

生活= [6.2 16.1 16.3 19.0 12.2 8.1 8.8 5.9 7.3 8.2．..16.1 12.8 9.8 11.3 5.1 10.8 6.7 1.2 8.3 2.3．..4.3 2.9 14.8 4.6 3.1 13.6 14.5 5.2 5.7 6.5．..5.3 6.4 3.5 11.4 9.3 12.4 18.3 15.9 4.0 10.4．..8.7 3.0 12.1 3.9 6.5 3.4 8.5 0.9 9.9 7.9]';

用直方图显示数据。

binWidth = 2;lastVal =装天花板(max ());binEdges = 0: binWidth: lastVal + 1;h =直方图(生活,binEdges);包含(“失败时间”）;ylabel (“频率”）;ylim ([0 10]);

图中包含一个轴对象。坐标轴对象包含一个直方图类型的对象。

由于生命周期数据通常遵循威布尔分布，一种方法可能是使用前面曲线拟合示例中的威布尔曲线来拟合直方图。要尝试这种方法，将直方图转换为一组点(x,y)，其中x是bin的中心，y是bin的高度，然后拟合这些点的曲线。

数量= histcounts(生活,binEdges);binCtrs = binEdges(1:end-1) + binWidth/2;h.FaceColor =[。9。9。9);持有在情节(binCtrs计数,“o”)；持有关

图中包含一个轴对象。轴对象包含两个直方图、直线类型的对象。

然而，将曲线拟合到直方图是有问题的，通常不建议这样做。

这个过程违背了最小二乘拟合的基本假设。箱数是非负的，这意味着测量误差不可能是对称的。此外，在尾部的仓数与在分布中心的仓数有不同的可变性。最后，垃圾箱计数有一个固定的总和，这意味着它们不是独立的测量。
如果将Weibull曲线拟合到钢筋高度，则必须约束曲线，因为直方图是经验概率密度函数（pdf）的缩放版本。
对于连续数据，拟合曲线到直方图而不是数据丢弃信息。
直方图中的条形高度依赖于箱边和箱宽的选择。

对于许多参数分布，极大似然是一种更好的估计参数的方法，因为它避免了这些问题。Weibull pdf几乎与Weibull曲线的形式相同:

$y ＝（ b / 一个) （ x / 一个)^{（ b - 1)} e^{- （ x / 一个)^{b}} ．$

然而 $b / 一个$ 替换刻度参数 $c$ 因为函数必须积分到1。要使用最大似然法将威布尔分布拟合到数据，请使用fitdist并具体说明“威布尔”作为发行版名。与最小二乘法不同，最大似然法发现的Weibull pdf最好地匹配缩放的直方图，而不最小化pdf和柱高之间的平方差的总和。

pd = fitdist(生活,“威布尔”）;

绘制数据的比例直方图并叠加拟合的pdf。

h=直方图（寿命、边缘、，“正常化”，“pdf”，“脸色”,[.9 .9 .9]); xlabel(“失败时间”）;ylabel (“概率密度”）;ylim (0.1 [0]);xgrid = linspace (20100 0) ';pdf = pdf (pd, xgrid);线(xgrid, pdf)

图中包含一个轴对象。轴对象包含两个直方图、直线类型的对象。

最好的做法是检查模型的拟合优度。

虽然通常不建议将曲线拟合到直方图，但该过程在某些情况下是合适的。例如，请参阅符合自定义分布．

分布拟合函数

统计和机器学习工具箱™ 包括函数fitdist用于将概率分布对象拟合到数据。它还包括专用配件功能（例如wblfit)，使用最大似然拟合参数分布，函数大中型企业用于拟合没有专用拟合函数的自定义分布，以及Ks密度用于拟合非参数分布模型到数据。
统计和机器学习工具箱还提供分布更健康该应用程序简化了分布拟合中的许多任务，例如生成可视化和诊断图。
最优化工具箱™中的函数使您能够适应复杂的分布，包括对参数有约束的分布。
MATLAB®函数fminsearch提供最大似然分布拟合。

另请参阅

曲线拟合与分布拟合

在曲线拟合和分布拟合之间进行选择

曲线拟合

曲线拟合函数

配电装置

分布拟合函数

另请参阅

相关的话题

统计和机器学习工具箱文档

金宝app

掌握机器学习:一步一步的指导与MATLAB