实现和CSS一样的easing动画？直接看MozillaChromium源码!
Posted Heavi的博客
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了实现和CSS一样的easing动画？直接看MozillaChromium源码!相关的知识，希望对你有一定的参考价值。
前言
在上一篇丝滑的贝塞尔曲线：从数学原理到应用介绍贝塞尔曲线实现动画时给自己留了一个坑，实现的动画效果和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、二分法介绍
**写在最后，如果大家有疑问可直接留言，一起探讨!感兴趣的可以点一波关注, 文章也在掘金上同步。
css3动画
ease：
1、ease：（逐渐变慢）默认值
2、linear：（匀速）
3、ease-in：(加速)
4、ease-out：（减速）
5、ease-in-out：（加速然后减速）
6、cubic-bezier
transition：all 1s ease 1s
代表的含义------执行变换的属性：transition-property,变换延续的时间：transition-duration,在延续时间段，变换的速率变化transition-timing-function,变换延迟时间transition-delay
以上是关于实现和CSS一样的easing动画？直接看MozillaChromium源码!的主要内容，如果未能解决你的问题，请参考以下文章