目录
3. ReLU (Rectified Linear Unit) 函数
(3)平移不变性(Translation Invariance)
(1)核心特点(也是卷积神经网络相比于全连接神经网络的特点)
一、神经元结构
1、生物神经元
神经细胞之间通过突触连接,突触通过复杂精巧的电化学过程传递信息。接收前面神经元的输入,汇总决策
传递。
2、人工神经元(MLP结构)
输入:向量 ,其中 D 是输入的维度。这些输入可以来自传感器或其他神经元的输出,它们代表了神经元接收到的信息或特征。
权重和偏置:每个输入信号 都有一个对应的权重
,这表示输入信号在神经元输出中的重要程度。权重向量
控制了每个输入信号对输出的影响。 此外,还有一个偏置 bbb,它是一个常数值,用来调整神经元的激活阈值。偏置的作用是确保模型的灵活性。
加权求和(加法器):所有输入信号 与相应的权重
相乘后求和,得到一个线性组合:
这个结果 z 是输入信号经过加权和偏置调整后的总和。
激活函数(激活层):计算得到的加权和 z 会输入到激活函数中。激活函数 用于引入非线性,使神经元能够学习到更复杂的映射关系。常见的激活函数包括Sigmoid、ReLU、Tanh等。在该图中,激活函数通过符号
来表示。经过激活函数处理后,输出
可以表示为:
这个输出 是神经元的最终输出,它可以被传递到下一层神经元或作为模型的最终预测结果。
输出:神经元的输出 将传递给下游神经元或最终用于分类或回归任务。通常,输出值
会根据实际任务的需求进行解释或处理。
3、激活函数
(1)激活函数的基本特性
- 连续并可导的非线性函数:一个合适的激活函数需要具备连续并且在其大部分区域可导的性质。这样才能确保在使用梯度下降法等优化算法时,神经网络能够有效地学习和更新参数。
- 简单性与效率:激活函数及其导函数应该尽可能简单,以提高网络计算的效率。复杂的激活函数会影响训练速度和模型的稳定性。
- 导函数的范围控制:导函数的值需要在一个合适的范围内,避免过大或过小,这有助于训练时梯度的稳定性。
(2)常见的激活函数
1. Sigmoid 函数
- 公式:
- 优点:常用于二分类问题的输出层。在早期的神经网络中较为普遍,逐渐被更简单的ReLU函数取代。
- 缺点:
- 在输入较大的正负值时,梯度几乎为零,容易导致梯度消失问题(vanishing gradient problem)。
- 输出范围在0到1之间,可能导致网络参数更新缓慢。
- 非零置中心化的输出会使得其后一层的神 经元的输入发生偏偏移(bias shift), 并进一步使得梯度下降的收敛速度变慢。
- 导数最大值为0.25:反向传播时,返回网络的误差将会在每一层收缩至少75%
2. Tanh 函数
- 公式:
- 优点:形状和Sigmoid函数相似,但输出范围是
,相比于Sigmoid可以避免一些偏置问题。
- 缺点:同样会遇到梯度消失的问题,特别是当输入的绝对值较大时。相比于Sigmoid缓解了梯度消失现象。
3. ReLU (Rectified Linear Unit) 函数
- 公式:
- 优点:
- 计算效率高,计算量小,因为没有复杂的数学运算,不涉及除法。
- 解决了梯度消失问题(当
时,梯度为1),收敛速度快。
- 一部分神经元的输出为 0 :网络稀疏, 减少了参数的相互依存关系,缓解过拟合。
- 缺点:
- 在
时,梯度为零,可能导致神经元“死亡”(即在某些情况下无法激活),从而影响模型的学习。
- 在
4. ReLU的变种
- Leaky ReLU:当
时,输出为
(其中
是一个小常数),它解决了ReLU函数死神经元问题。
- PreLU:类似于Leaky ReLU,但它引入了可学习的参数
,以使得每个神经元的负区域具有自适应的斜率。
- ELU (Exponential Linear Unit):当
时,使用指数函数
,避免了ReLU的死神经元问题,并能产生更好的训练效果。
- Softplus 函数:平滑的ReLU函数变种,它为负值提供了一个平滑的过渡,避免了ReLU的硬性零值。
(3)激活函数的选择
- ReLU:由于其计算速度快且有效,通常是首选激活函数,特别是在深度神经网络中。
- Leaky ReLU、ELU、MaxOut:如果ReLU的效果不佳,可以尝试这些ReLU的变种。
- Tanh:适合在LSTM等序列建模问题中使用。
- Sigmoid 和 Tanh:这些激活函数适用于RNN(LSTM、注意力机制等)结构中作为门控机制的激活函数,尤其在序列处理和控制模型中非常有用,但在其他网络中逐渐减少使用。
(4)激活函数的综合比较
- 计算效率:ReLU由于其简洁的计算形式,在计算效率上优于其他复杂的激活函数。
- 梯度消失问题:Sigmoid和Tanh容易在深层网络中遇到梯度消失问题,ReLU及其变种可以有效避免这一问题。
- 非零输出中心:Tanh的输出为-1到1之间,解决了Sigmoid函数输出非零中心问题,而ReLU则解决了梯度消失和学习效率的问题。
总体来说,ReLU和其变种因其高效性和相对简单的实现而广泛应用于现代神经网络中。
4、另一种人工神经元结构(KAN)
MLP结构和KAN结构对比
(1)KAN的理论基础
- KAN基于Kolmogorov-Arnold定理,该定理表明如果一个函数 f 是多元连续函数,它可以被表示成有限数量的单变换连续函数的两层嵌套组合形式。
- 公式表示为:
其中,
是外部函数,
是内部函数,表示多元输入到网络输出的嵌套结构。
(2)KAN的结构
-
外部函数
:外部函数是可学习的函数(比如 B样条函数),它们的作用是增强KAN的表示能力。
-
内部函数
:这些是与输入
相关的函数,通常是固定的(例如Sigmoid或线性函数),它们通过嵌套运算来处理输入数据。
-
KAN的计算模型:
- KAN在计算结构上与传统的多层感知机(MLP)有所不同:其节点上的激活函数是固定的,而边上的激活函数则是可学习的。
- KAN的每个层次都代表一个小型网络,表达了能量关系,更适合处理复杂的非线性映射。
(3)KAN的优势
- 与传统MLP对比:与传统的MLP网络(例如多层感知机)不同,KAN的模型是通过先变后加的方法来设计的,而不是先加后变(MLP结构)。
- 细粒度与多样化:KAN在学习过程中能处理更多细节,并通过可训练的内外部函数来实现更复杂的非线性转换。
(4)KAN的参数设计
- 可学习的激活函数:KAN的边(即每条连接)上的激活函数是可学习的,这意味着激活函数的形状和参数可以通过训练优化,以适应给定任务的需求。
- 可学习的B样条函数:B样条函数作为KAN的一部分,提供了更加细粒度的调控,并能通过其参数来细化模型的表现。
- 权重初始化:在KAN中,权重(如 w)会被初始化为适当的值(例如,使用Xavier初始化方法),以确保神经元的非线性转换能够有效训练。
(5)KAN的激活函数设计
- 内部函数
:内部函数设计为
,其中:
是一种Sigmoid线性单位(SiLU)。
是一个B样条函数,带有可训练的参数
来控制样条的曲线形状。
- 训练过程:这些内部函数的参数(如
)会在训练过程中进行优化,初始时这些参数是随机的,并且遵循正态分布。
(6)KAN的优势与应用
- KAN能够实现复杂的非线性映射,适用于对数据有高度要求的任务,如函数逼近和高维数据建模。
- 由于其模型设计和结构,KAN相对于传统的神经网络模型来说在处理某些特定问题时可能展现出更高的表现力和适应能力,特别是在多层次非线性转换的场景中。
- 表达能力更强,可用更紧凑的结构来刻画复杂的映射关系。
(7)与多层感知机(MLP)的对比
- 在传统的多层感知机中,每一层的激活函数通常是固定的(例如Sigmoid或ReLU),而KAN则将激活函数放置在边上,且这些激活函数是可训练的。
- KAN的模型设计使得它能够表达更加复杂的非线性映射,从而提升其表达能力。
总结来说,KAN通过将外部可学习激活函数和内部固定激活函数相结合,提供了一种更灵活、更强大的神经网络架构,能够更好地处理复杂的函数逼近任务,并具有比传统MLP更高的灵活性和表达能力。
《有关B样条的介绍见附录1》
二、全连接前馈神经网络
1、网络结构 和 符号体系
输入层:输入数据为 ,其维度为
(即有 D 个特征)。
第一个隐藏层:
输入层到第一隐藏层,输入 要乘上权重
(维度为
),
再加上偏置 (维度为
),
得到 (维度为
),
经过非线性变换 后得到 第一个隐藏层的输出
,
第二个隐藏层:
第一隐藏层到第二隐藏层, 要乘上权重
(维度为
)
再加上偏置 (维度为
),
得到 (维度为
),
经过非线性变换 后得到 第二个隐藏层的输出
,
. . .
. . .
最后一个隐藏层:
最后一个隐藏层的输出 ,记为
有时候要经过softmax层控制输出的维度。
综合一下:
对符号 的解释:
是从第
层到第
层的权重矩阵。这个矩阵的维度是
,其中:
是第
层的神经元个数。
是第
层的神经元个数。
矩阵 的每个元素:
- 矩阵中的每个元素
表示第
层的第
个神经元与第
层的第
个神经元之间的连接权重。
- 也就是说,矩阵
中的第
行、第
列的元素表示的是从第
层的第
个神经元传递到第
层的第
个神经元的权重。
矩阵结构:
- 矩阵
的结构可以展开为:
- 其中,
代表从第
层的第
个神经元到第
层的第
个神经元的连接权重。
2、模型训练:反向传播
(1)目标函数
模型参数为: 。
损失函数:损失函数 被定义为:
这是典型的交叉熵损失函数,用于分类问题。
目标函数:目标函数 用于寻找最小的网络参数
,其表达式为:
其中 λ 是正则化项的系数,W 是权重矩阵, 是权重的 Frobenius 范数。
(2)梯度下降
反向传播的目标是计算损失函数 J 相对于每一层的权重和偏置的梯度,使用链式法则进行计算。我们从输出层开始,逐层向后计算梯度。
损失函数关于参数 的梯度为
其中:
- 由于
(当
时,此公式为
),所以
容易计算,为:
需要使用反向传播算法计算。
(3)反向传播
即计算 ,记为
,被称为误差项。
3.1 输出层的误差项:
输出层的误差是损失函数对输出的梯度:
其中, 是损失函数关于输出的梯度,
是激活函数的导数,
表示Hadamard积(元素乘法)。
对于常见的损失函数(如均方误差),有:
3.2 计算隐藏层的梯度
对于第 层的误差
,我们有:
这里, 是上一层的误差,
是激活函数对第
层加权输入的导数 。
(4)权重更新
一旦得到了每一层的误差,接下来通过梯度下降法更新权重和偏置。
-
权重更新规则:
其中,学习率为 ,
是损失函数对权重的梯度。
梯度的计算为:
-
偏置更新规则:
梯度的计算为:
(5)梯度消失和梯度爆炸
反向传播的误差传递公式:
数值结果:
-
梯度爆炸:
-
梯度消失:
这部分展示了梯度爆炸和梯度消失问题,尤其在深层网络中,随着层数的增加,如果激活函数选择不适当,激活函数的导数不接近1,梯度传递过程可能导致梯度逐渐变大(爆炸)或逐渐变小(消失)。
- 梯度爆炸问题会导致训练过程对学习率敏感。
梯度消失问题会导致梯度趋近于0,无论学习率为多少,训练都没有进展,即只有顶层训练有效,网络变深也不会使得效果变好,即前边的隐藏层无论多少,都没有被训练。
三、为什么要有深度模型
1、两类模型比较
参数一定时,模型胖矮好还是瘦高好?
实验:
实验证明:瘦高更好。
如何解释呢?
2、瘦高模型表达能力更强的解释
(1)分层组合与层次结构
-
分层组合思想:
深层网络能够将一个复杂的函数分解为多个简单函数的嵌套组合。假设一个目标函数可以表示为多个子函数的组合:
每一层只需要学习相对简单的映射,逐层提取和转换特征,这种结构正好符合许多实际问题中数据的内在层次性。
-
高效表示:
通过深层结构,网络可以在每一层提取出中间表示,然后在下一层进行进一步组合。相比之下,浅层网络(胖矮模型)必须在一层内直接捕获输入到输出的复杂映射,这往往需要指数级更多的神经元才能达到同样的效果。
(2)参数利用率与“指数”优势
-
参数利用效率:
理论上,有研究表明对于某些特定函数类,深层网络可以用指数级更少的参数来表示,而浅层网络可能需要指数级的宽度才能逼近同样的函数。也就是说,深层结构具有“参数共享”或“复用”中间特征的优势,使得每个参数在构建复杂映射时能发挥更大的作用。 -
数学直观:
如果我们考虑一个由多层非线性激活函数构成的网络,假设每层用ReLU激活函数,深层网络可以形成许多线性区域。研究发现,随着深度的增加,网络能够划分出指数级别数量的线性区域,从而具有更高的函数拟合能力。这意味着深层网络在表达非线性、复杂函数方面具有理论上的优势。
(3)局部特征的逐层抽象
-
逐层抽象:
“瘦高”模型的每一层都能逐步抽象出更高级的特征,从原始输入到低级特征,再到中级和高级特征,最后得到用于决策的表示。这样的逐层抽象使得模型能够更精细地捕捉数据的结构信息。 -
对比浅层模型:
“胖矮”模型虽然单层拥有大量神经元,能够在单层内提取丰富信息,但由于层次结构较少,难以捕捉输入数据中的多层次、递归性的特征组合,从而在表达复杂关系上往往不如深层结构高效。
(4)理论支持
- 表达能力理论:
许多理论工作表明,对于特定类别的函数,深层网络在表达能力上具有“深度优势”(depth efficiency)。例如,一些文献证明,若使用深度为 LLL 的网络,其表示能力可以达到浅层网络所需参数数量的指数级减少。这为“瘦高”模型在理论上比“胖矮”模型更有效地表达复杂函数提供了数学依据。
四、卷积神将网络
1、对图像为什么要用卷积网络
(1)局部感受野与局部连接
-
局部相关性:
图像中相邻像素往往具有很强的相关性,例如边缘、纹理和颜色变化都是局部特征。卷积层通过局部感受野只关注输入图像的局部区域,能够高效捕捉这些低级特征,而不需要每个神经元与全图所有像素连接。 -
数学描述:
对于输入图像和卷积核
,卷积操作可表示为:
这种局部计算方式使得网络在提取特征时只需关注图像的局部信息。
(2)参数共享与降低模型复杂度
-
参数共享:
卷积核在整个图像上滑动应用,意味着同一组权重被用于不同的局部区域。这种权重共享机制显著减少了需要学习的参数数量,相比全连接层(每个像素都需要独立连接)能大幅降低模型复杂度。 -
好处:
降低参数数量不仅减少了内存需求和计算量,也降低了过拟合的风险,使得模型在数据量有限时更为稳定。
(3)平移不变性(Translation Invariance)
-
平移不变性原理:
由于卷积核在图像的各个位置进行相同的操作,卷积层对图像中物体的平移具有天然的不变性。也就是说,无论特征在图像的哪个位置出现,卷积层都能检测到相同的特征。 -
实际意义:
这使得卷积神经网络在处理实际图像时更具鲁棒性,能够有效识别出移动或位置略有偏移的目标。
(4)层次特征提取
-
多层结构:
CNN通常由多个卷积层和池化层堆叠而成,较低层可以捕捉到边缘、角点等简单特征,而较高层则可以组合低级特征构成更复杂的形状和语义信息(如物体的一部分或整体结构)。 -
特征层次化:
这种逐层抽象的能力使得卷积神经网络能够自动学习到从低级到高级的多层次特征表示,适应不同复杂度的任务需求。
(5)池化层的作用
-
降维与抗噪:
池化层(如最大池化、平均池化)对卷积层的输出进行下采样,降低了特征图的尺寸,减少计算量,同时对输入的微小变形和噪声具有一定的鲁棒性。 -
平滑特征:
通过池化层,可以进一步增强特征的不变性,使得后续的全连接层或分类器能够更好地利用这些抽象出的特征。
(6)总结
卷积神经网络之所以在图像处理上更为优秀,主要在于它利用了图像数据的局部性和空间结构,通过局部感受野和参数共享机制大幅降低模型复杂度,并通过多层结构逐步抽象出高层次特征。同时,平移不变性和池化层的应用使得CNN对于图像中的位移、旋转以及噪声具有较好的鲁棒性。正是这些特点,使得卷积神经网络成为图像分类、目标检测、图像分割等任务中的主流模型。
2、卷积层
(1)核心特点(也是卷积神经网络相比于全连接神经网络的特点)
局部链接:卷积神经网络(CNN)中,每个神经元仅与前一层的部分区域(即局部感受野)相连接,而不是与整张输入图像或上一层的所有神经元连接,从而减少了参数数量和计算量。
权值共享:在卷积层中,同一卷积核(滤波器)在整个输入图像上滑动并应用相同的权重,允许模型学习图像中的局部特征,并减少了需要学习的参数数量。
局部链接和权值共享使得卷积神经网络的复杂度相比于全连接神经网络的复杂度大大降低。
(2)边缘填充
卷积操作前通常进行边缘填充,以保证卷积操作不会使得数据缩减。
一般选择 0 填充,但我认为扩散填充更好,即填充的值和其里边的值一样。
(3)多输入通道
输入有几个通道就需要几个卷积核,输出结果是所有通道卷积结果的和:
有时输出通道也可以是多个,即多通道输入,多通道输出。
(4)一维卷积和三维卷积
最常见的卷积是二维卷积,但一维卷积和三维卷积也同样存在,
- 一维卷积:主要是对一维信号进行处理,如:文本、语音、时间序列
- 三维卷积:主要的三维信号有视频、医学图像等
(5)更多的花式卷积
5.1~5.4为在空间上的花式,5.5~5.7为在通道上的花式。
5.1 大小为1×1的卷积核
不识别空间模式,只是融合通道。
5.2 不同尺寸卷积核:inception v2
Inception系列网络(如Inception v2)采用了多尺度处理的思想,在同一层中并行使用不同尺寸的卷积核(例如1×1、3×3、甚至5×5)来捕捉不同尺度的信息。
- 多尺度特征提取: 大尺寸卷积核具有较大的感受野,能够捕捉全局或较大区域的特征;小尺寸卷积核则更适合提取细节信息。
- 1×1卷积的作用: 1×1卷积通常用于降维或升维,使得后续使用大尺寸卷积时能降低计算量,同时融合跨通道信息。
数学描述与优化:
- 如果直接用一个
卷积核,其参数量为
;
- Inception v2中常将
卷积拆分成两个连续的
卷积(或者使用更精细的分解方式),这样参数量大约为
既降低了参数量,也增加了非线性表达能力。
5.3 空洞卷积
5.4 可变形卷积核
标准卷积核:长方形或正方形
可变形卷积:卷积核形状可变,只看到感兴趣的图像区域
5.5 分组卷积
基本思想:
分组卷积将输入特征图的通道分成多个组,每一组单独进行卷积运算,再将结果拼接起来。
- 参数和计算量降低: 将通道分组后,每组卷积核只负责一部分通道,从而减少了参数数量和计算复杂度。
- 信息局部性: 这种方式有时可以迫使不同组学习到不同类型的特征,增强网络表达能力。
数学描述:
设输入通道数为 ,输出通道数为
,如果分为
组,则每组输入通道为
,输出通道为
。
标准卷积的参数量为
而分组卷积的参数量则为
5.6 深度可分离卷积
基本思想:
深度可分离卷积将标准卷积分解为两步:
- Depthwise Convolution(逐通道卷积): 对每个输入通道单独进行卷积,每个通道使用一个卷积核。
- Pointwise Convolution(1×1卷积): 将逐通道卷积的输出进行通道间的线性组合,从而融合各通道信息。
优势与数学描述:
- 参数和计算量大幅减少:
- 标准卷积的计算量为
- 深度可分离卷积的计算量为
- 标准卷积的计算量为
- 表达能力: 虽然参数较少,但经过逐通道和跨通道的两次卷积组合,仍能捕捉丰富特征。该结构在移动端网络(如MobileNet)中得到广泛应用。
5.7 通道洗牌
基本思想:
在使用分组卷积时,由于每个组内的通道之间信息相对独立,可能会导致跨组的信息交流不足。通道洗牌操作通过对通道进行重新排列,打乱各组内部的通道顺序,从而让后续分组卷积能够接收到来自不同组的混合信息。
实现方法与作用:
- 实现:
- 将输入特征图的通道维度按照分组数进行重塑(reshape),例如将形状从
变为
;
- 对分组后的通道维度进行转置;
- 再reshape回原始形状。
- 将输入特征图的通道维度按照分组数进行重塑(reshape),例如将形状从
- 作用:
通过这种方式,不同分组内的信息得以交换,使得后续层在执行分组卷积时能跨组提取特征,增强模型的表达能力。 - 应用:
通道洗牌被广泛应用于轻量级网络(如ShuffleNet),在保持较低计算量的同时提升了网络的性能。
3、池化层
池化层即降采样层,常用的池化操作有两种:
- 最大值池化:取一定区域的最大值
- 平均值池化:取一定区域的平均值
池化层的作用:
- 降低分辨率
- 提供一定程度的平移不变性
4、卷积网络的模型训练
(1)卷积层的训练
卷积层的反向传播类似全连接层的反向传播(只是有参数共享和 局部连接)
(2)池化层的训练
2.1池化层反向传播的核心思想
经过Pooling后,两层的feature map不一样大,为了保证梯度对位能传播下去,需要保证传递的loss(或者梯度) 总和不变。
2.2 池化层的反向传播:
反向传播的目标是计算损失函数对池化层输入的梯度,然后将其传递给上一层。在池化层的反向传播中,我们会根据前向传播时池化操作的类型(最大池化或平均池化)来更新梯度。
最大池化的反向传播:
- 前向传播:在一个池化区域内,最大池化会选择最大的值作为输出,并记录下该值的位置。
- 反向传播:反向传播时,梯度仅通过前向传播时记录的最大值位置传递。其他位置的梯度值为零,因为在这些位置上池化操作没有影响。
例如,如果池化窗口是2x2,且最大的值出现在位置(1,1),则在反向传播时,仅该位置的梯度会被传递到输入特征图,其他位置的梯度为0。保证池化前后的梯度(残差)之和保持不变。
平均池化的反向传播:
- 前向传播:平均池化计算区域内所有值的平均数作为输出。
- 反向传播:在反向传播中,平均池化的梯度被均匀分配到池化窗口中的每个位置。也就是说,所有位置的梯度是相等的,等于该池化区域的梯度值除以池化窗口的元素个数。
例如,对于一个2x2的池化窗口,假设在前向传播中输出的梯度为0.4,则在反向传播时,0.4的梯度会均匀分配到四个位置,每个位置的梯度为0.1。保证池化前后的梯度(残差)之和保持不变。不是把梯度复制 K 遍之后直接反向传播回 去。 这样会造成loss之和变为原来的 K 倍, 网络是会产生梯度爆炸。
2.3 梯度的传播:
无论是最大池化还是平均池化,池化层的反向传播都遵循以下的基本原则:
- 在 最大池化 中,只有前向传播时最大值的位置会接收到梯度信息。
- 在 平均池化 中,梯度被均匀分配到每个池化窗口的元素上。
2.4 池化层反向传播中的特殊注意事项:
- 梯度消失问题:池化层,尤其是最大池化层,由于只保留了最大值的梯度信息,可能会导致梯度在某些情况下非常小,影响模型的学习过程。
- 池化层参数更新:池化层本身通常没有可学习的参数(例如权重和偏置),因此它只负责将输入的梯度传递到上一层。
附录
1、贝塞尔曲线 与 B样条
(1)Langrange插值
如果知道三个点P0,P1,P2,那么如何确定一条曲线呢?
最简单的想法是,分别连接 PO,P1和 P1,P2,得到两个线段就可以确定一个曲线,但它不平滑,一阶导数不连续。二次曲线:三个点可以确定一条抛物线,最直接的做法是列三个方程,解方程组,得到二次曲线。类似地,可以递推到更多点。对于n个点,可以使用n-1阶次的函数来确定一个唯一的曲线。这就是多项式插值。
拉格朗日插值法:
- 简述:假设有 n 个点
,其中
互不相同。对于每一个 k (1 到 n 之间),定义 n−1 次多项式
如下:
-
Lagrange 基函数的性质:
对于
-
构建插值多项式:使用上述的基函数,构造 n−1 次的拉格朗日插值多项式
如下:
这个多项式
满足
,其中
。
Runge现象:
例子:
-
已知区间 [−1,1] 上函数
,取等距节点
,其中
。
-
通过拉格朗日插值法,构造
多项式:
其中,
为拉格朗日基函数。
-
从图中可以看出,在 x=0 附近,
能较好地通过
,但在某些地方,如区间 [−1,−0.8] 和 [0.8,1] 之间,
与
差异很大,这种现象称为 Runge现象。
-
由于插值法要求曲线一定要经过规定的一系列点,这样的限制很强,又由于点数很多,次数很高,会导致振荡的很厉害。
因此可采用拟合法,放松限制,让曲线不一定要必须经过某些点,常见的又最小二乘法,贝塞尔曲线等。
- 依据原有数据点,通过参数调整设置,使得生成曲线与原有点差距最小(最小二乘),因此曲线未必会经过原有数据点。
- 样条曲线(Spline curves):是给定一系列控制点而得到的一条曲线,曲线形状由这些点控制。一般分为插值样条和拟合样条。
(2)样条曲线 和 贝塞尔曲线
i、样条曲线:
一个简单且行之有效的方法是,把这些点作为限制点,然后在这些限制点中放置一条具有弹性的金属片,最后金属片绕过这些点后的最终状态即为所需曲线。而最终得到的形状曲线,就是样条曲线。这也是该名字的由来,其中金属片就是样条,形成的曲线就是样条曲线。
该想法虽然巧妙,但显然不具有推广性。因此问题就出来了,如何将其抽象出一个数学模型,从而在已知控制点条件下,仅仅通过数学公式从而获得平滑的样条。下面的贝塞尔曲线的限制比这个例子还要松,贝塞尔曲线不需要曲线过所有限制点。
ii、贝塞尔曲线:
三个点确定的贝塞尔曲线(二次贝塞尔曲线):
- 公式:
这是一个二次Bezier曲线的标准公式,其中 ,
,
是控制点,t 是参数,范围通常为 [0, 1]。
- 构建过程:
-
1. 定义新的控制点:
和
分别是线段
和
上的点,且 t 从0~1的过程中
和
从线段的一段到另一端。
-
2. 构建Bezier曲线:
-
3. 相关关系:
t 从0~1的变化过程使得这段贝塞尔曲线就生成了。这段内容介绍了如何通过线性插值方式逐步构建Bezier曲线,逐渐从控制点计算出曲线的具体点。
n+1 个点确定的贝塞尔曲线(n 次贝塞尔曲线):
假设有 n+1 个点,确定了 n 次的贝塞尔曲线,公式为:
另一种表示形式:
解释:
可以将 理解为基(当 n 确定时,这组基就是确定的),P 为这组基的线性组合系数。
这段内容描述了贝塞尔曲线的两种表示方式。第一种是通过控制点和对应的基函数来定义,第二种则通过权重 和控制点
来表示。
例子:三次基函数,也就是为了拟合四个点,共四个基。
iii、贝塞尔曲线的递推公式
当阶数很高时,贝塞尔曲线不易求解,所以可以由递推公式求解。
n 阶贝塞尔曲线 可以由前面 n 个点决定的 n−1 次贝塞尔曲线
线性组合推导而来,公式为:
这段内容描述了如何通过递归方式来构建贝塞尔曲线,每一步都利用前一阶贝塞尔曲线的线性组合来构造下一阶曲线。递归公式使得贝塞尔曲线的计算可以通过逐步逼近的方式进行。对于 n 阶贝塞尔曲线,可以通过逐步计算 (n−1) 阶的曲线,然后用一个简单的线性插值来结合控制点,得到更高阶的曲线。这不仅减少了复杂度,还方便了贝塞尔曲线的实现。
iv、贝塞尔曲线的性质与特点
- 贝塞尔曲线幂次=控制点个数-1,幂次越高,需要的计算也越复杂,效率低。
- 控制点与凸包性质:贝塞尔曲线总是位于其控制点构成的凸包内,这意味着曲线不会出现剧烈的波动,能够保证较好的稳定性和直观的几何解释。
- 端点插值性质:曲线在 t=0 时必定通过第一个控制点 P0,在 t=1 时通过最后一个控制点 Pn,这有利于曲线的拼接和连接。
- 切线性质:贝塞尔曲线在起点和终点处具有明确的切线方向:在 t=0 处,切线方向由向量 P1−P0 决定;在 t=1 处,切线方向由向量 Pn−Pn−1 决定。
- 仿射不变性:曲线对仿射变换(平移、旋转、缩放)具有不变性,即对控制点做任意仿射变换,曲线形状会随之相应变换,但基本性质不变。这为曲线在不同坐标系或图形环境中的应用提供了便利。
- 变差递减(Variation Diminishing)性质:贝塞尔曲线不会比其控制多边形“更加波动”,这意味着当与直线相交时,曲线的交点数不会超过控制多边形与直线交点数,从而保证了曲线的平滑性。
(3)B样条
对于复杂曲线,不要求用n-1次曲线,而是让次数低一点经常使用三次贝塞尔曲线一段一段地拼接成目标曲线,如Ps或Ai中使用钢笔工具画出物体轮廓所做的那样。
确保最终整体曲线一次光滑的条件是在连接如果使用这种方法,点两侧的斜率相等,即连接点和其两侧控制点共线。
i、基本概念
B样条曲线是利用一组控制点和一个节点向量(knot vector)定义的分段多项式曲线。它的主要特点在于:
- 局部控制:每个控制点只影响曲线的一部分,使得局部修改不影响整体形状。
- 高连续性:通过适当设置节点向量,可以控制曲线在各段之间的连续性(如
连续,其中 p 为曲线的阶数)。
- 灵活性:能够表示复杂形状,同时可推广到曲面(B样条曲面)和更高维度的数据拟合问题中。
ii、数学定义与构成要素
节点向量(Knot Vector)
节点向量是一系列非递减的实数序列,记为
- p 是曲线的阶数(通常阶数 p+1 表示多项式的次数,如三次曲线 p=3)。
- n+1 表示控制点的数量。
- 节点向量确定了曲线参数 t 的分段区间以及每个分段的多项式基函数的定义区间。常见的节点向量类型包括均匀节点向量、非均匀节点向量和开放节点向量(Open Uniform Knot Vector)。
B样条基函数
B样条的核心在于一组基函数 ,它们采用递归定义,通常称为 Cox-de Boor 递归公式:
-
零次基函数(阶数为 0,即 p=0)定义为:
-
递归定义(对于
):
其中如果分母为 0,则该项通常定义为 0。这一递归公式确保了每个
是一个分段多项式,且仅在区间
内有非零值。
B样条曲线的定义
给定控制点 和对应的节点向量 U,B样条曲线可定义为:
其中 为阶数为 p 的 B样条基函数。曲线的形状由控制点的位置和节点向量决定。
iii、主要性质
- 局部性:由于每个基函数
的非零区间仅局限于
,一个控制点
只影响曲线中对应的参数区间。这种局部控制特性使得修改某个控制点只对曲线的局部区域产生影响,而不会波及整体形状。
- 连续性控制:B样条曲线的连续性依赖于基函数的阶数 p 以及节点向量中节点的重复次数。一般来说,如果节点
重复次数为 k(
),那么在该节点处曲线的连续性为
。
- 高阶连续:合理选择节点向量可以使曲线在各段之间具有较高的连续性,例如三次 B样条(
)在没有重复节点的情况下具有
连续性。
- 凸包性:曲线始终位于其相应控制点的凸包内,这保证了曲线的数值稳定性和几何直观性。
- 仿射变换不变性:B样条曲线对仿射变换(平移、旋转、缩放)保持不变,即对控制点做任意仿射变换,曲线形状会相应变换,但基本性质不变。