平滑M样条

Question

我目前正在使用d=4阶Mspline近似模型中所需的基线函数：

r(t)=γ1 M1(t|4,x)+...γp Mp(t|4,x)

为避免局部波动，我想通过惩罚基线函数的曲率来惩罚可能性。

l(O)-∫r''(u)^2du

我的问题是我不知道如何计算R中的以下部分：∫r''(u)^2du

我发现∫r''(u)^2du = γT Ωγ

Ω=∫M''(u)M''(u)du

M''是M样条的二阶导数矩阵。