前言
在上一篇丝滑的贝塞尔曲线：从数学原理到应用介绍贝塞尔曲线实现动画时给自己留了一个坑，实现的动画效果和CSS的transition-timing-function: cubic-bezier差别较大，如下图所示，红色为Linear、绿色为CSS的cubic-beizer、蓝色为自己实现的cbezier。本着有坑必填的原则，直接把Mozilla、Chromium的cubic-bezier实现源码给翻出来。
为什么和CSS效果不一致？假设贝塞尔曲线为cubic-bezier(0.25, 0.1, 0.25, 1.0)，X轴表示时间，Y轴表示元素样式属性变化量，设贝塞尔曲线为Q(t), 我直接把动画时间当做贝塞尔曲线的t，用Qy(X)来计算Y值。而正确的流程是已知X值，根据Qx(t)=X对方程求根，其根即为t值，再根据t值求Y = Qy(t)。
整个流程的难点在于如何对Qx(t)=X求根t值，已知贝塞尔曲线方程为
Q(t) = X1 * (1 - t)³ + X2 * 3t(1 -t)² + X3 * 3t²(1 - t) + X4 * t³
由于曲线X轴范围为[0, 1]，因此X1=0、X4=1, 从而Q(t)可写为
Q(t) = X2 * 3t(1 -t)² + X3 * 3t²(1 - t) + t³
将方程各项展开再合并成如下形式，主要是为了更方便的介绍浏览器源码实现。
Q(t) = ( ( (1.0 - 3.0 * X2 + 3.0 * X1) * t + (3.0 * X2 - 6.0 * X1) ) * t + (3.0 * X1) ) * t
浏览器是如何实现的
浏览器实现流程：使用参数cubic-beizer(X1, Y1, X2, Y2)生成贝塞尔曲线方程，根据传入的X值对方程求解得到t值，再将t传入Y方向的曲线方程计算Y值。难点在于如何求解方程的根，而求根一般可通过Newton-Raphson method牛顿法、Bisection method二分法。在深入浏览器源码之前，先介绍这两种求根方法，Mozilla和Chromium都基于这两种方法求解贝塞尔曲线的根。
Newton-Raphson method牛顿法
牛顿法是求解数学方程的一种常用方法。它是一种迭代算法，通过不断逼近函数的根来寻找解。牛顿法可以用于求解实函数的根，也可以用于求解实函数的最小值或最大值。
牛顿法的基本思想是：选择一个初始值 \(x_0\)，然后通过递归计算出一个数列 \(x_1\), \(x_2\), ...，\(x_n\)，使得每一项\(x_{n+1}\) 都是函数 f(x) 在 \(x_n\) 处的切线与x轴的交点。也就是说，我们用当前的切线来近似替代原函数，从而得到更接近真实根的近似解。假设Xn处的函数值为\(f(x_n)\)、斜率为\(f'(x_n)\)，\(X_{n+1}\)为\(f'(x_n)\)与X轴的交点，从而有:
其中 f'(x) 是函数 f(x) 在 x 处的导数。每一次迭代中，我们通过计算 \(f(x_n)\) 和 \(f'(x_n)\) 来求出 \(x_{n+1}\) 的值，并将其作为下一次迭代的初始值。这个过程会不断重复，直到得到一个满足预定精度要求的解。
牛顿法的缺陷是，当斜率很小时，与X轴平行，得不到解。
Bisection method二分法
Bisection method是一种求解方程的数值方法，可以用来找到一个函数f(x)=0的解。这种方法利用了连续函数的介值定理（Intermediate Value Theorem），即在一个区间[a,b]内，如果f(a)和f(b)的符号不同，那么必定存在至少一个数c属于[a,b]，使得f(c)=0。
算法的基本思想是将一个区间[a,b]不断分成两半，然后找到包含根的那一半，并继续在该子区间内进行分割。每次分割之后，都会得到一个新的子区间，其长度是原来的一半。通过不断重复这个过程，可以越来越逼近根的位置。
上图中f(x)递增，并且f(a) < 0, f(b) > 0。我们继续获取了a和b的中点x0。由上图可知f(x0) > 0，所以我们可以把x0看成是新的b。于是我们继续寻找a和x0的中点，重复上述过程，由于我们最大的误差就是区间的长度，所以当我们区间的长度缩减到足够小，那么就说明我们已经找到了一个足够近似的解。
Mozilla实现
由于@greweb已经将Mozilla的缓动动画通过JS实现，Mozilla的实现直接通过JS代码介绍。
主体结构和我们介绍的流程一致，定义bezier方法，参数分别为P1、P2坐标，返回结果为函数BezierEasing, 其内部会根据传入的X值，调用getTForX(x)求根得到t，最后调用calcBezier(t, mY1, mY2)计算Y轴的曲线值。
export function bezier (mX1, mY1, mX2, mY2) {
// ...其他逻辑
return function BezierEasing (x) {
// Because JavaScript number are imprecise, we should guarantee the extremes are right.
if (x === 0 || x === 1) {
return x;
}
return calcBezier(getTForX(x), mY1, mY2);
};
};
calcBezier函数比较简单，就是根据贝塞尔曲线方程以及传入的t值计算结果，由于要计算一阶导数，将Q(t)拆成了三部分，1.0 - 3.0 * X2 + 3.0 * X1对应A函数，3.0 * X2 - 6.0 * X1对应B函数，3.0 * X1对应C函数。
// Q(t) = ( ( (1.0 - 3.0 * X2 + 3.0 * X1) * t + (3.0 * X2 - 6.0 * X1) ) * t + (3.0 * X1) ) * t
function A (aA1, aA2) { return 1.0 - 3.0 * aA2 + 3.0 * aA1; }
function B (aA1, aA2) { return 3.0 * aA2 - 6.0 * aA1; }
function C (aA1) { return 3.0 * aA1; }
// at为t，aA1表示x1或者x2,aA2表示y1或者y2
// 其形式和p0 * Math.pow(1 - t, 3) + p1 * 3 * t * Math.pow(1 - t, 2) + p2 * 3 * Math.pow(t, 2) * (1 - t) + p3 * Math.pow(t, 3)一致
function calcBezier (aT, aA1, aA2) { return ((A(aA1, aA2) * aT + B(aA1, aA2)) * aT + C(aA1)) * aT; }
难点在于getTForX函数，其作用是根据X值计算t，涉及到曲线方程求根，而求根一般可通过Newton-Raphson method牛顿法、Bisection method二分法, 也就是我上面介绍的两种方法。
而Mozilla在求根时会根据斜率选择不同的方法，为了提升计算性能，在求根之前先采样，使用sampleValues缓存X轴方向当t=0,0.1，...,0.9, 1时对应的曲线值。
// kSplineTableSize为11，kSampleStepSize为1.0 / (11 - 1.0) = 0.1;
// sampleValues存储样本值的目的是提升性能，不用每次都计算。
var sampleValues = float32ArraySupported ? new Float32Array(kSplineTableSize) : new Array(kSplineTableSize);
// i从0到10，sampleValues长度为11
for (var i = 0; i < kSplineTableSize; ++i) {
// i * kSampleStepSize的范围0到1(10 * 0.1);
// sampleValues[0] = calcBezier(0, mX1, mX2);
// sampleValues[1] = calcBezier(0.1, mX1, mX2);
// ...
// sampleValues[9] = calcBezier(0.9, mX1, mX2);
// sampleValues[10] = calcBezier(1, mX1, mX2);
sampleValues[i] = calcBezier(i * kSampleStepSize, mX1, mX2);
}
对于getTForX函数的实现，首先遍历sampleValues数组，找到小于目标值aX的最大值，为了减少求根的遍历次数，为了使guessForT尽量接近目标根，使用(aX - sampleValues[currentSample]) / (sampleValues[currentSample + 1] - sampleValues[currentSample])计算出差值在两个样本区间的百分比dist，那么intervalStart + dist * KSampleStepSize即为根据样本值能找到最接近目标根guessForT的初始值。
// 已知X值，根据X值求解T值
var NEWTON_MIN_SLOPE = 0.001;
function getTForX (aX) {
var intervalStart = 0.0;
var currentSample = 1;
// lastSample为10
var lastSample = kSplineTableSize - 1;
// sampleValues[i]表示i从0以0.1为step，每一步对应的曲线的X坐标值，直到X坐标值小于等于aX
// 假如aX=0.4，则sampleValues[currentSample]<=aX为止
for (; currentSample !== lastSample && sampleValues[currentSample] <= aX; ++currentSample) {
// intervalStart为到aX经过的step步骤
intervalStart += kSampleStepSize; // kSampleStepSize为0.1
}
//TODO:currentSample为什么要减1？sampleValues[currentSample]大于了ax，所以要--，使得sampleValues[currentSample]<=ax
--currentSample;
// Interpolate to provide an initial guess for t
// ax-sampleValues[currentSample]为两者之间的差值，而(sampleValues[currentSample + 1] - sampleValues[currentSample])一个步骤之间的总差值。
var dist = (aX - sampleValues[currentSample]) / (sampleValues[currentSample + 1] - sampleValues[currentSample]);
// guessForT为预计的初始T值，很粗糙的一个值，接下来会基于该值求根(t值)。
var guessForT = intervalStart + dist * kSampleStepSize;
// 预测的T值对应位置的斜率
var initialSlope = getSlope(guessForT, mX1, mX2);
// 当斜率大于0.05729°时，使用newtonRaphsonIterate算法预测T值。0.05729是一个很小的斜率
if (initialSlope >= NEWTON_MIN_SLOPE) {
return newtonRaphsonIterate(aX, guessForT, mX1, mX2);
} else if (initialSlope === 0.0) { // 当斜率为0，则直接返回
return guessForT;
} else { // 当斜率小于0.05729并且不等于0时，使用binarySubdivide
// 求得的根t，位于intervalStart和intervalStart + kSampleStepSize之间, mX1、mX2分别对应p1、p2的X坐标
return binarySubdivide(aX, intervalStart, intervalStart + kSampleStepSize, mX1, mX2);
}
}
使用getSlope函数计算guessForT的斜率initialSlope，NEWTON_MIN_SLOPE为测试的一个斜率临界(0.001)，当大于等于该斜率时使用牛顿法球根，否则使用二分法球根。牛顿法实现：
var NEWTON_ITERATIONS = 4;
function newtonRaphsonIterate (aX, aGuessT, mX1, mX2) {
// NEWTON_ITERATIONS为4， 只进行了4次迭代， 根据精度和性能之间做了平衡。
for (var i = 0; i < NEWTON_ITERATIONS; ++i) {
// 计算t值对应位置的斜率
var currentSlope = getSlope(aGuessT, mX1, mX2);
if (currentSlope === 0.0) {
return aGuessT;
}
// 假设f(t) = 0，求解方程的根。其f(t)=calcBezier(t) - ax
// 牛顿-拉佛森方法: Xn-1 = Xn - f(t) / f'(t)，应用到求贝塞尔曲线的根：Tn = Tn+1 - (calcBezier(t) - ax) / getSlope(t)
var currentX = calcBezier(aGuessT, mX1, mX2) - aX;
aGuessT -= currentX / currentSlope;
}
// 这里只迭代了4次，求得近似值
return aGuessT;
}
牛顿法求根会不断根据上一个值\(x_n\)迭代下一个值\(x_{n+1}\)，而newtonRaphsonIterate函数仅仅迭代了4次，是根据精度和性能之间的平衡考量。根据牛顿法公式:
这里的f(t) = calcBezier(t) - aX, f'(t)即计算在t位置的斜率用getSlope(aGuessT, mX1, mX2)表示，所以aGuessT的迭代值可写为aGuessT = aGuessT - f(t)/f'(t)。当斜率为0时，aGuessT即为最终值；否则直到迭代4次结束。
当斜率小于NEWTON_MIN_SLOPE时，降级使用Bisection method二分法,要计算的t值肯定在[intervalStart、intervalStart + kSampleStepSize]之间，并且f(intervalStart)和f(intervalStart + kSampleStepSize)的乘积小于0，满足二分法条件。
binarySubdivide(aX, intervalStart, intervalStart + kSampleStepSize, mX1, mX2);
二分法会不断缩小[a,b]范围，直到计算的值和目标值aX差值小于SUBDIVISION_PRECISION,也即小于0.0000001。 或者大于最大迭代次数SUBDIVISION_MAX_ITERATIONS也会终止。
var SUBDIVISION_PRECISION = 0.0000001;
var SUBDIVISION_MAX_ITERATIONS = 10;
// 二分球根法：
// 求得的根t，位于aA和aB之间, mX1、mX2分别对应p1、p2的X坐标
// https://zhuanlan.zhihu.com/p/112845185
function binarySubdivide (aX, aA, aB, mX1, mX2) {
var currentX, currentT, i = 0;
do {
currentT = aA + (aB - aA) / 2.0;
// 假设f(t) = 0，求解方程的根。其f(t)=calcBezier(t) - ax
currentX = calcBezier(currentT, mX1, mX2) - aX;
if (currentX > 0.0) {
aB = currentT;
} else {
aA = currentT;
}
// 如果currentX小于等于最小精度(SUBDIVISION_PRECISION)或者超过迭代次数SUBDIVISION_MAX_ITERATIONS，则终止
} while (Math.abs(currentX) > SUBDIVISION_PRECISION && ++i < SUBDIVISION_MAX_ITERATIONS);
return currentT;
}
aA、aB为初始范围，每次迭代计算中点currentT = aA + (aB - aA) / 2.0, 然后使用其值计算currentX(calcBezier(currentT, mX1, mX2))。当currentX大于aX时，说明currentT还大于目标t，则将aB赋值为currentT缩小右边界；否则将aA赋值为currentT缩小左边界。
至此，Mozilla的实现介绍完毕。
Chromium实现
实现和Mozilla类似，核心部分SolveCurveX也是对已知的X和曲线方程求解t。先迭代样本数组，找到第一个大于等于X的样本值，设置最接近的值为t2, 在区间[t0, t1]中。
在选择牛顿法和二分法时，Chromium和Mozilla有些区别，Chromium先使用牛顿法迭代kMaxNewtonIterations4次直到计算值小于精度newton_epsilon，或者一阶导(斜率)SampleCurveDerivativeX(t2)小于最小斜率kBezierEpsilon也会中止。
如果牛顿法未计算出满足精度要求的根，则继续使用二分法判断。和Mozilla的区别在于没限制迭代次数，一直根据t0<t1条件迭代，直到找到满足精度的t2。r
double CubicBezier::SolveCurveX(double x, double epsilon) const {
double t0;
double t1;
double t2 = x;
double x2;
double d2;
int i;
// Linear interpolation of spline curve for initial guess.
// 迭代样本值，找到和x最接近的初始t，这里为t2，在t0和t1之间。
double delta_t = 1.0 / (CUBIC_BEZIER_SPLINE_SAMPLES - 1);
for (i = 1; i < CUBIC_BEZIER_SPLINE_SAMPLES; i++) {
if (x <= spline_samples_[i]) {
t1 = delta_t * i;
t0 = t1 - delta_t;
// 根据x差值在所在区间的百分占比计算t在阶段间的百分占比，t2为接近x根的初始值
t2 = t0 + (t1 - t0) * (x - spline_samples_[i - 1]) /
(spline_samples_[i] - spline_samples_[i - 1]);
break;
}
}
// Perform a few iterations of Newton's method -- normally very fast.
// See https://en.wikipedia.org/wiki/Newton%27s_method.
// 牛顿法的精度
double newton_epsilon = std::min(kBezierEpsilon, epsilon);
for (i = 0; i < kMaxNewtonIterations; i++) {
// x2为t2所在x值和目标x差值
x2 = SampleCurveX(t2) - x;
// 如果差值小于精度则任务t2即为方程根
if (fabs(x2) < newton_epsilon)
return t2;
// d2为t2位置的斜率
d2 = SampleCurveDerivativeX(t2);
// 如果斜率小于最小精度kBezierEpsilon,任务与x轴平行，无法继续计算
if (fabs(d2) < kBezierEpsilon)
break;
// 否则使用牛顿法继续迭代
t2 = t2 - x2 / d2;
}
if (fabs(x2) < epsilon)
return t2;
// Fall back to the bisection method for reliability.
// 降级使用二分法
while (t0 < t1) {
x2 = SampleCurveX(t2);
// 如果x2小于最小精度，则任务为要求的根
if (fabs(x2 - x) < epsilon)
return t2;
// 下面的流程和二分法流程一致
if (x > x2)
t0 = t2;
else
t1 = t2;
t2 = (t1 + t0) * .5;
}
// Failure.
return t2;
}
效果对比
使用封装的动画函数animate执行动画，查看动画效果。easing为bezier-easing(0.25, 0.1, 0.25, 1.0)，和CSS指定贝塞尔曲线一致。
animate(bezieRef2.current, {
duration: 2000,
easing: 'bezier-easing(0.25, 0.1, 0.25, 1.0)',
styles: [{ left: !ended ? '0%' : '100%' }, { left: !ended ? '100%' : '0%' } ]
})
animate函数内部调用resolveEasing解析传入的字符串bezier-easing(0.25, 0.1, 0.25, 1.0),并生成贝塞尔曲线函数，生成的曲线函数形式为function BezierEasing (x): number，根据传入的X(时间百分比)求解t值，再根据t计算Y值(元素位置)。
/**
* 元素动画
* @param el DOM元素
* @param props 动画属性
*/
function animate(el, props) {
const duration = props.duration;
const easingFunc = resolveEasing(props.easing);
const styleFuncs = resolveStyles(el, props.styles);
const start = Date.now();
const animationHandle = () => {
const timeRatio = (Date.now() - start) / duration;
if (timeRatio <= 1) {
const percent = easingFunc(timeRatio);
for (const key in styleFuncs) {
const elementStyle = el.style;
elementStyle[key] = styleFuncs[key](percent);
}
requestAnimationFrame(animationHandle);
}
};
animationHandle();
}
原始的CSS实现缓动动画，非常简单，直接设置样式即可：
.fade-in-cbezier {
transition: all 2s;
transition-timing-function: cubic-bezier(0.25, 0.1, 0.25, 1.0);
}
设置相同的曲线参数，自己实现的贝塞尔缓动动画和CSS对比，其效果几乎完全一致，图中绿色的为CSS，黄色为自己实现的动画效果。
CSS 默认提供了几个固定的贝塞尔曲线linear、ease-in、ease-out、ease-in-out。
- linear 匀速移动，cubic-bezier(0.0, 0.0, 1.0, 1.0)
- ease-in 先慢后快 cubic-bezier(0.42, 0, 1.0, 1.0)
- ease-out 先快后慢 cubic-bezier(0, 0, 0.58, 1.0)
- ease-in-out 慢速、提速、减速三个阶段 cubic-bezier(0.42, 0, 0.58, 1.0)
可通过在线预览查看动贝塞尔曲线。
运行效果：
动画库入门
标星44.8k的animejs，其核心部分也使用了Mozilla的贝塞尔曲线动画，另外一个重要的问题是如何转换单位，如将left设置为50%、20em或者50vw，只有将单位统一才能计算中间差值。
animejs将元素的初始值使用convertPxToUnit函数转换为目标单位，以convertPxToUnit(el, '50', 'em') )为例，convertPxToUnit为元素el的父元素附加一个临时子要素tempEl，目标单位为em，设置tempEl的宽度为100em，然后得到因子factor=100em/offestWidth，即目标长度/像素，那么转换后的值为factor * parseFloat ('50')。
/**
* 将像素值转换为目标单位值
* @param {*} el
* @param {*} value
* @param {*} unit
* @returns
*/
function convertPxToUnit(el: Element, value: string, unit: string) {
const baseline = 100;
// 创建一个和el类型一样的要素
const tempEl = document.createElement(el.tagName);
// 获取要素的parent
const parentEl = (el.parentNode && (el.parentNode !== document)) ? el.parentNode : document.body;
parentEl.appendChild(tempEl);
tempEl.style.position = 'absolute';
// 设置基线为100个目标单位
tempEl.style.width = baseline + unit;
// 宽度因子，目标长度/一个像素
const factor = baseline / tempEl.offsetWidth;
parentEl.removeChild(tempEl);
// parseFloat会将最后的单位忽略得到数值
const convertedUnit = factor * parseFloat(value);
return convertedUnit;
}
在执行动画时，可同时设置多个属性，如opacity、height、width、left、top等，resolveStyles函数遍历每个属性并为其生成一个获取中间值的函数，如设置{width: '30%'}, 那么目标值destValue为30。如何获取width的当前值并转换单位为%？先调用getComputedStyle(el).getPropertyValue(key.toLowerCase())获取其像素值，再通过前面实现的convertPxToUnit将像素值转换为%。有了startValue和destValue，中间值即可通过startValue + total * percent计算。
export type StyleProperties = {
[x: string]: string | number;
}
function resolveStyles(el: Element, styles: StyleProperties) {
const keys = Object.keys(styles);
const styleFuncs: { [x: string]: (t: number) => number | string } = {};
for (const key of keys) {
const value = styles[key] + '';
const unit = <string>getUnit(value);
const destValue = parseFloat(value);
const styleValue = getComputedStyle(el).getPropertyValue(key.toLowerCase());
const startValue = unit ? convertPxToUnit(el, styleValue, unit) : parseFloat(styleValue);
const total = destValue - startValue;
styleFuncs[key] = (percent: number) => {
const curVal = startValue + total * percent;
return unit ? curVal + unit : curVal;
}
}
return styleFuncs;
}
对外提供的API格式为animate(el: HTMLElement, props: AnimationProps),props包括duration、styles、easing三个参数，函数会根据(Date.now() - start) / duration计算动画执行进度，然后遍历styleFuncs设置每个样式属性的中间值。
export type AnimationProps = {
duration: number;
styles: StyleProperties;
easing: string;
}
/**
* 元素动画
* @param el DOM元素
* @param props 动画属性
*/
export function animate(el: HTMLElement, props: AnimationProps): { pause: () => void } {
const duration = props.duration;
const easingFunc = resolveEasing(props.easing);
const styleFuncs = resolveStyles(el, props.styles);
const cAniInstance = {
paused: false,
}
const start = Date.now();
const animationHandle = () => {
if (cAniInstance.paused) {
return;
}
const timeRatio = (Date.now() - start) / duration;
if (timeRatio <= 1) {
const percent = easingFunc(timeRatio);
for (const key in styleFuncs) {
const elementStyle = el.style as any;
elementStyle[key] = styleFuncs[key](percent);
}
requestAnimationFrame(animationHandle);
}
}
requestAnimationFrame(animationHandle);
return {
pause: () => {
cAniInstance.paused = true;
}
}
}
使用比较简单，实现一个弹窗效果：
animate(ref1.current, {
duration: 500,
easing: 'cbezier(0.25, 0.1, 0.25, 1.0)',
styles: {
opacity: '1',
width: '300px',
height: '200px',
left: '30%',
top: '30%',
}
})
实现鼠标跟踪，由于gif图像有丢帧，效果不是很明显，Demo已上传到cbezier。
const mouseMove = throttle((e) => {
const x = e.offsetX, y = e.offsetY;
if (aniPlay) {
aniPlay.pause();
}
// console.log(ref1.current.style.left, ref1.current.style.top);
aniPlay = animate(ref1.current, {
duration: 15,
// easing: 'ease-out',
easing: 'cbezier(0.25, 0.1, 0.25, 1.0)',
styles: {
left: x + 'px',
top: y + 'px',
}
})
}, 15);
parentRef.current.addEventListener('mousemove', mouseMove);
下一步计划： 实现完整的动画库hiani.js。以上仅仅是一个入门版动画，要实现和animejs一样的效果，还需要支持rotate、scale、translate、skew等transform样式以及颜色变化，下一步将完善动画库支持更多样式属性、提供play、pause、reset等完整的函数控制动画。
参考
1、Bezier Curve as Easing Function
2、Bezier Curve based easing functions – from concept to implementation
3、mozilla贝塞尔实现
4、Chromium贝塞尔实现
5、Newton Raphson method
6、Newton Raphson formula
7、牛顿法和二分法的区别
8、二分法介绍
**写在最后，如果大家有疑问可直接留言，一起探讨!感兴趣的可以点一波关注, 文章也在掘金上同步。
标签：el,const,动画,Mozilla,mX2,mX1,源码,easing,1.0 From： https://www.cnblogs.com/w-wanglei/p/hiani.html