开关电容闭环运放

本文主要讨论基于开关电容的运放和积分器的建立过程。

0.Revision History

2025.4.24 1st manuscript

1.闭环系统

一个闭环系统如下图所示:

其闭环的传递函数可以表示为
$$\begin{array}{c}H(s)=\frac{v_o}{v_i}=\frac{a(s)}{1+a(s)f(s)}=\frac{a(s)}{1+T(s)}=\frac{1}{f(s)}\cdot \frac{T(s)}{1+T(s)}\end{array}$$
式$(1)$中环路增益可以表示为
$$T(s)=a(s)\cdot f(s)$$
闭环系统的稳定性分析在于环路，即$T(s)$的性质。根据Barkhausen稳定性判据，当
$$\begin{gathered} \angle T(s)=-180° \\ |T(s)|>1 \end{gathered}$$
时，环路不再稳定。从稳定性和建立速度的角度，通常要求环路的增益裕度GM为3-5，相位裕度为$60°$~$70°$。

2.开关电容电路

基于OTA的开关电容积分器和放大器如下图所示(Ref: Murmann’s Slide)：

开关电容积分器（放大器）在积分（放大）相位的单端等效电路如图所示(Ref: Nan Sun’s Slide)：

利用返回比（Return Ratio）分析方法，将电路中所有独立源置零，断开环路中的受控源，并添加激励信号，计算受控源的响应，返回比可以表示为
$$\begin{array}{c}ReturnRatio=-\frac{s_{return}}{s_{test}}\end{array}$$
输入源为独立源，输入接地；受控源为压控电流源，将其与输出断开；在断开的位置添加一个激励信号，由于受控源为电流源，因此添加电流激励$i_t$。闭环单级OTA的返回比计算小信号模型如下图所示

虚地点电压$v_x$可以表示为
$$\begin{array}{c}{v}_x=\frac{C_f}{C_f+C_s+C_{par}}\cdot v_o=\beta v_o=\beta \cdot [-i_t\cdot (R_L//\frac{1}{s{C}_{Ltot}})]\end{array}$$
其中反馈系数（feedback factor）$\beta$为
$$\begin{array}{c}\beta =\frac{C_f}{C_f+C_s+C_{par}}\end{array}$$
$C_{Ltot}$为从输出端看到的总的负载电容，表示为$C_f$串联($C_s$并联$C_{par}$)，然后并联$C_L$，表达式为
$$\begin{array}{c}{C}_{Ltot}=C_L+\frac{C_f\cdot (C_s+C_{par})}{C_f+C_s+C_{par}}=C_L+(1-\beta )\cdot C_f\end{array}$$
环路增益（返回比）可以表示为
$$\begin{array}{c}T(s)=-\frac{i_r}{i_t}=-\frac{G_m\cdot v_x}{i_t}=\beta \cdot \frac{G_m{R}_L}{1+s{R}_L{C}_{Ltot}}=\beta \cdot \frac{G_m{R}_L}{1+s{R}_L[C_L+(1-\beta )\cdot C_f]}\end{array}$$
由此可见，单级运放的环路传递函数中包含一个极点。当频率非常高时，环路增益可以简化为
$$\begin{array}{c}T(s){|}_{sishigh}\approx \beta \cdot \frac{G_m}{s{C}_{Ltot}}=\frac{{\omega }_c}{s}\end{array}$$
其中${\omega }_c$为拐点频率（同时也是T(s)的单位增益带宽${\omega }_{ugb}$）
$$\begin{array}{c}{\omega }_{ugb}=\frac{\beta \cdot G_m}{C_{Ltot}}\end{array}$$

3.闭环等效

由式$(1)$，闭环运放的传递函数可以表示为
$$H(s)=\frac{1}{f(s)}\cdot \frac{T(s)}{1+T(s)}$$
其中$\frac{1}{f(s)}$为环路增益无穷大时的理想闭环增益，计算原理图如下所示

当环路增益无穷大时，虚地点$v_x$严格为0，同时有
$$(v_i-v_x)\cdot s{C}_s=(v_x-v_0)\cdot s{C}_f$$
解得
$$\begin{array}{c}\frac{1}{f(s)}=\frac{v_o}{v_i}=-\frac{C_s}{C_f}\end{array}$$
闭环运放的DC增益为
$$\begin{array}{c}H(0)=\frac{1}{f(0)}\cdot \frac{T(0)}{1+T(0)}=-\frac{C_s}{C_f}\cdot \frac{\beta \cdot G_m{R}_L}{1+\beta \cdot G_m{R}_L}=-\frac{C_s}{C_f}\cdot \frac{1}{1+\frac{1}{\beta \cdot G_m{R}_L}}\end{array}$$
静态增益误差为
$$\begin{array}{c}|{ϵ}_0|=|\frac{H(0)-1/f(0)}{1/f(0)}|=|\frac{T(0)}{1+T(0)}-1|\approx \frac{1}{T(0)}=\frac{1}{\beta \cdot G_m{R}_L}\end{array}$$
由式(7)，随着频率增加，环路增益下降，闭环增益也随之下降，高频处的闭环增益可以表示为
$$\begin{array}{c}|H(j\omega )|=|\frac{1}{f}\cdot \frac{T(s)}{1+T(s)}|=\frac{C_s}{C_f}\cdot |\frac{{\omega }_c/s}{1+{\omega }_c/s}|=\frac{C_s}{C_f}\cdot \frac{{\omega }_c}{\sqrt{{\omega }^2+{\omega }_c^2}}\end{array}$$
当频率为$\omega ={\omega }_c({\omega }_{ugb})$时，闭环增益下降为$|H(0)|$的$1/\sqrt{2}$，即环路的单位增益带宽对应闭环系统的-3dB带宽，即
$$\begin{array}{c}{f}_{-3dB}=\frac{1}{2\pi }\cdot \frac{\beta G_m}{C_{Ltot}}\end{array}$$
单级OTA的增益带宽积（GBW）可以表示为
$$\begin{array}{c}GB{W}_{OTA}=G_m{R}_L\cdot \frac{1}{2\pi R_L{C}_{Ltot}}=\frac{1}{2\pi }\cdot \frac{G_m}{C_{Ltot}}\end{array}$$
则环路的-3dB带宽也可以表示为
$$\begin{array}{c}{f}_{-3dB}=\beta \cdot GB{W}_{OTA}\end{array}$$
根据以上分析，可以对开关电容闭环运放作单极点等效
$$\begin{array}{c}H(s)=\frac{H(0)}{1+\frac{s}{2\pi f_{-3dB}}}=\frac{H(0)}{1+\frac{s}{\beta G_m/C_{Ltot}}}=\frac{H(0)}{1+s\tau }\end{array}$$
上式中闭环运放建立时间常数为
$$\begin{array}{c}\tau =\frac{1}{\beta G_m}\cdot C_{Ltot}\end{array}$$
第一项为闭环运放输出端看过去的等效阻抗，第二项为输出端的等效负载电容。
以上推导中忽略了电容前馈的影响，在高频时，输入信号通过反馈电容前馈到输出端，当OTA中的跨导电流与前馈电流完全抵消时，OTA输出电压为0，即产生了一个零点。

在输出节点，根据KCL，零点的位置满足
$$v_x\cdot G_m=(v_x-0)\cdot s{C}_f$$
解得右半平面的零点为
$${\omega }_z=\frac{G_m}{C_f}$$
对比左半平面的极点为
$${\omega }_p={\omega }_c=\frac{\beta G_m}{C_{Ltot}}$$
当反馈电容$C_f$较小时，这个零点通常远远高于极点的频率。对于积分器的情况，由于反馈电容$C_f$通常较大，零点位置会远低于极点。考虑前馈电容导致的零点后的闭环运放传递函数改写为
$$\begin{array}{c}H(s)=H(0)\cdot \frac{1-s/{\omega }_z}{1+s/{\omega }_p}=-\frac{C_s}{C_f}\cdot [1-\frac{1}{T(0)}]\cdot \frac{1-\frac{s}{G_m/C_f}}{1+\frac{s}{\beta G_m/C_{Ltot}}}\end{array}$$

另一种输入结构

以上的闭环传递函数的推导中，输入位于采样电容的左侧；对于输入位于运放另一侧的情况（通常用于噪声计算），小信号等效模型如下图所示

当增益为无穷大时，虚地点电压$v_x$严格等于输入$v_i$，则有
$$v_i=\frac{C_f}{C_s+C_f+C_{par}}\cdot v_o$$
解得
$$\frac{1}{f(s)}=1+\frac{C_s+C_{par}}{C_f}=1/\beta$$
环路的计算与前文一致，基于此输入模式的运放闭环传递函数可以表示为
$$\begin{array}{c}H(s)=\frac{H(0)}{1+s\tau }=\frac{1}{\beta }\cdot [1-\frac{1}{T(0)}]\cdot \frac{1}{1+s\tau }\end{array}$$

4.线性建立

首先忽略前馈零点的影响，采用单极点闭环等效，即
$$\begin{array}{c}H(s)=\frac{H(0)}{1+\frac{s}{2\pi f_{-3dB}}}=\frac{H(0)}{1+s\tau }\end{array}$$
运放输入端加入幅度为$V_{step}$的阶跃信号后，输出端的时域表达式可以通过拉普拉斯反变换求解，即
$$\begin{array}{rl}{V}_{out}(t)& ={\mathcal{L}}^{-1}[H(s)\cdot \frac{V_{step}}{s}]=H(0)\cdot V_{step}\cdot {\mathcal{L}}^{-1}[\frac{1}{1+s\tau }\cdot \frac{1}{s}]=H(0)\cdot V_{step}\cdot (1-e^{-\frac{t}{\tau }})\\ & =-\frac{C_s}{C_f}\cdot V_{step}\cdot [1-\frac{1}{T(0)}]\cdot (1-e^{-\frac{t}{\tau }})\end{array}$$
其中
$${\mathcal{L}}^{-1}[\frac{1}{1+s\tau }\cdot \frac{1}{s}]={\mathcal{L}}^{-1}[\frac{1}{s}-\frac{\tau }{1+s\tau }]={\mathcal{L}}^{-1}[\frac{1}{s}-\frac{1}{s-(-1/\tau )}]=u(t)\cdot (1-e^{-\frac{t}{\tau }})$$
由时域表达式，OTA的有限直流增益引入静态建立误差，有限GBW引入动态建立误差，如下图所示

动态建立误差为
$${ϵ}_d=1-(1-e^{-\frac{t}{\tau }}){|}_{t=t_{settle}}=e^{-\frac{t_{settle}}{\tau }}$$
其中$t_{settle}$为建立时间，可以等价得表示为
$$t_{settle}=-\tau \cdot \ln {ϵ}_d$$
常见的建立时间要求如下表格所示

${ϵ}_d$	$t_{settle}/\tau$
10%	2.3
1%	4.6
0.1%	6.9
0.01%	9.2
$10^{-6}$	13.8

假设运放的阶跃建立过程完全线性，对于一个N-bit精度的系统，采样频率为$f_s$，时钟周期一半时间为建立时间，并且要求动态建立误差不超过1/2LSB，则需要满足
$$e^{-\frac{\frac{1}{2}\cdot t_s}{\tau }}=e^{-\frac{\frac{1}{2}\cdot \frac{1}{f_s}}{\frac{1}{2\pi \cdot f_{-3dB}}}}<\frac{1}{2}\cdot \frac{1}{2^N}$$
解得
$$\begin{array}{rl}& t_s>2\cdot (N+1)\cdot \ln 2\cdot \tau \\ & f_{-3dB}=\beta \cdot f_{GB{W}_{OTA}}>\frac{(N+1)}{\pi }\cdot \ln 2\cdot f_s\end{array}$$
其中ln2 = 0.7,对于24bit的建立精度，要求
$$\begin{array}{rl}& t_s>2\cdot 17.3\tau \\ & \beta \cdot f_{GB{W}_{OTA}}>5.7\cdot f_s\end{array}$$
每提高1bit的建立精度要求，对运放的带宽要求提高0.22$f_s$。

电容前馈的影响

当采样电容左侧开关闭合后，相当于在闭环运放输入端加入一个阶跃信号，在高频下，电阻视作断路，闭环运放等效为一个电容分压网络，如下图所示。

阶跃输入在OTA输入端分压产生的前馈量$V_{xstep}$完全由电容分压决定
$$\begin{array}{c}{V}_{xstep}=V_{istep}\cdot \frac{C_s}{C_s+C_{par}+\frac{C_f{C}_L}{C_f+C_L}}\end{array}$$
运放输出端产生的前馈量为
$$\begin{array}{rl}\frac{V_{ostep}}{V_{istep}}=\frac{V_{ostep}}{V_{xstep}}\cdot \frac{V_{xstep}}{V_{istep}}& =\frac{C_f}{C_f+C_L}\cdot \frac{C_s}{C_s+C_{par}+\frac{C_f{C}_L}{C_f+C_L}}=\frac{C_s}{C_f}\cdot \frac{\beta C_f}{C_L+(1-\beta )C_f}\\ & =\frac{C_s}{C_f}\cdot k_f\end{array}$$
上式中$k_f=\frac{\beta C_f}{C_L+(1-\beta )C_f}$为前馈电压系数，由上一小节，考虑电容前馈影响后的闭环运放传递函数为
$$H(s)=-\frac{C_s}{C_f}\cdot [1-\frac{1}{T(0)}]\cdot \frac{1-s/{\omega }_z}{1+s/{\omega }_p}$$
则输入为$V_{istep}$时，输出电压时域表达式为
$$\begin{array}{c}{V}_{out}(t)={\mathcal{L}}^{-1}[H(s)\cdot \frac{V_{istep}}{s}]=-V_{istep}\cdot \frac{C_s}{C_f}\cdot [1-\frac{1}{T(0)}]\cdot [1-e^{-\frac{t}{\tau }}\cdot (1+\frac{{\omega }_p}{{\omega }_z})]\end{array}$$
其中
$$1+\frac{{\omega }_p}{{\omega }_z}=1+\frac{\beta G_m}{C_{Ltot}}/\frac{G_m}{C_f}=\frac{C_L+C_f}{C_L+(1-\beta )C_f}=1+k_f>1$$
即前馈电压系数恰好等于极点与零点的比值，动态建立误差为
$${ϵ}_d=1-[1-e^{-\frac{t}{\tau }}\cdot (1+\frac{{\omega }_p}{{\omega }_z})]=e^{-\frac{t}{\tau }}\cdot (1+\frac{{\omega }_p}{{\omega }_z})=e^{-\frac{t}{\tau }}\cdot (1+k_f)$$
电容前馈导致动态建立误差提高。等价得可以得到建立时间与建立精度的关系为
$$t=-\tau \cdot \ln (\frac{{ϵ}_d}{1+k_f})=\tau \cdot [(N+1)\ln 2+\ln (1+k_f)]$$

5.压摆建立

对于全差分NMOS输入对运放，且输入管均工作在饱和区时，根据长沟道平方律关系
$$\begin{array}{c}{V}_{in1}-V_{in2}=V_{GS1}-V_{GS2}=\sqrt{\frac{2I_{D1}}{{\mu }_n{C}_{ox}\frac{W}{L}}}+V_{TH}-(\sqrt{\frac{2I_{D2}}{{\mu }_n{C}_{ox}\frac{W}{L}}}+V_{TH})\end{array}$$
输入管电流满足关系式$I_{D1}+I_{D2}=I_{SS}$，对上式两侧取平方
$$(V_{in1}-V_{in1}{)}^2=\frac{2}{{\mu }_n{C}_{ox}\frac{W}{L}}(I_{SS}-2\sqrt{I_{D1}{I}_{D2}})$$
同时有式子
$$4I_{D1}{I}_{D2}=(I_{D1}+I_{D2}{)}^2-(I_{D1}-I_{D2}{)}^2=I_{SS}^2-(I_{D1}-I_{D2}{)}^2$$
则上式可以解得
$$(I_{D1}-I_{D2}{)}^2=-\frac{1}{4}({\mu }_n{C}_{ox}\frac{W}{L}{)}^2\cdot (V_{in1}-V_{in2}{)}^4+I_{SS}\cdot {\mu }_n{C}_{ox}\frac{W}{L}\cdot (V_{in1}-V_{in2}{)}^2$$
$$I_{D1}-I_{D2}=\sqrt{{\mu }_n{C}_{ox}\frac{W}{L}{I}_{SS}}\cdot (V_{in1}-V_{in2})\cdot \sqrt{1-\frac{{\mu }_n{C}_{ox}\frac{W}{L}}{4I_{SS}}\cdot (V_{in1}-V_{in2}{)}^2}$$
差分输入为0时，两个输入管电流均为尾电流的一半，得到全差分运放的平衡态跨导
$$g_m=\frac{\partial I_{od}}{\partial V_{id}}{|}_{V_{id}\to 0}=\underset{V_{id}\to 0}{\lim }\frac{I_{D1}-I_{D2}}{V_{in1}-V_{in2}}=\sqrt{2{\mu }_n{C}_{ox}\frac{W}{L}{I}_{DS}}=\sqrt{{\mu }_n{C}_{ox}\frac{W}{L}{I}_{SS}}$$
此时输入管的过驱动电压为
$$V_{ov}=\frac{2I_D}{g_m}=\frac{I_{SS}}{g_m}$$
当差分输入电压为大信号时，使得只有一端的输入管开启，另一侧输入管截止，此时恰好满足$V_{GS2}=V_{TH}$，则差分输入为
$$\Delta V_{int1}=\sqrt{\frac{2I_{SS}}{{\mu }_n{C}_{ox}\frac{W}{L}}}=\sqrt{2}\cdot \frac{I_{SS}}{g_m}$$
即压摆发生时的输入电压只差为平衡态过驱动电压的$\sqrt{2}$倍。
当差分输入电压绝对值超过$\Delta V_{int1}$时，输出差分电流表达式不再成立。此时平衡态跨导为0，电路不具备压控电流源的能力，即超出工作范围。将平衡态跨导表达式带入，差分电流表达式可以重写为
$$\begin{array}{c}{I}_{od}=g_m{V}_{id}\cdot \sqrt{1-\frac{g_m^2}{4I_{SS}^2}\cdot V_{id}^2}\end{array}$$
$$\begin{array}{c}\frac{I_{od}}{I_{SS}}=\frac{V_{id}}{I_{SS}/g_m}\cdot \sqrt{1-\frac{1}{4}\cdot (\frac{V_{id}}{I_{SS}/g_m}{)}^2}\end{array}$$
以$\frac{V_{id}}{I_{SS}/g_m}$为x轴，$\frac{I_{od}}{I_{SS}}$为y轴，绘制曲线得

因此压摆建立的工作区间为
$$|V_{id}|>\sqrt{2}\cdot \frac{I_{SS}}{g_m}$$
输入管gm/id越大，越容易进入压摆区间。运放输入差分电压超过$\sqrt{2}\cdot \frac{I_{SS}}{g_m}$时，运放首先压摆建立时，运放输出电流由尾电流决定，当脱离压摆工作区间时，运放回归线性建立，建立过程中虚地点电压和输出差分电流如下图所示。

运放压摆率为
$$SR=\frac{\partial V_{od}}{\partial t}=\frac{I_{SS}}{C_{Ltot}}$$
输入端电压变化速度为$\beta \cdot SR$。当运放输入的电压阶跃为$V_{xstep}$($V_{xstep}$为负值)时，运放压摆的工作时间为
$$t_{slew}=\frac{\left|V_{xstep}\right|-\sqrt{2}\cdot \frac{I_{SS}}{g_m}}{\beta \cdot SR}=-\frac{V_{xstep}+2\sqrt{2}\cdot \frac{I_D}{g_m}}{\beta \cdot SR}$$
完成压摆建立后，运放输出电压为
$$\begin{array}{rl}{V}_{oslew}& =SR\cdot t_{slew}+v_{ostep}=-SR\cdot \left(\frac{V_{xstep}+\sqrt{2}\cdot \frac{I_{SS}}{g_m}}{\beta \cdot SR}\right)+\frac{C_s}{C_f}\cdot \frac{\beta C_f}{C_L+\left(1-\beta \right)C_f}\cdot V_{istep}\\ & =-\frac{1}{\beta }\left(V_{xstep}+\sqrt{2}\cdot \frac{I_{ss}}{g_m}\right)+\frac{C_s}{C_f}\cdot k_f\cdot V_{istep}\end{array}$$
运放输入从0完全线性建立到该输出电压的等效建立时间$t_{lin,eq}$可通过动态建立误差计算，即
$$dynamicerror=\frac{C_L+C_f}{C_L+\left(1-\beta \right)C_f}\cdot e^{-\frac{t_{lin,eq}}{\tau }}=\frac{-\frac{C_s}{C_f}\cdot V_{istep}-V_{oslew}}{-\frac{C_s}{C_f}\cdot V_{istep}}=1+\frac{SR\cdot t_{slew}+v_{ostep}}{\frac{C_s}{C_f}\cdot V_{istep}}$$
上式子中，
$$v_{ostep}=\frac{C_s}{C_f}\cdot \frac{\beta C_f}{C_L+\left(1-\beta \right)C_f}\cdot v_{istep}$$
上式可以改写为
$$\frac{C_L+C_f}{C_L+\left(1-\beta \right)C_f}\cdot e^{-\frac{t_{lin,eq}}{\tau }}=1+\frac{\beta C_f}{C_L+\left(1-\beta \right)C_f}+\frac{SR\cdot t_{slew}}{\frac{C_s}{C_f}\cdot V_{istep}}=(1+k_f)-\frac{SR\cdot t_{slew}}{\frac{C_s}{C_f}\cdot V_{istep}}$$
$$e^{-\frac{t_{lin,eq}}{\tau }}=1-\frac{C_L+\left(1-\beta \right)C_f}{C_f+C_L}\cdot \frac{SR\cdot t_{slew}}{\frac{C_s}{C_f}\cdot V_{istep}}=1-\left(1-\beta \frac{C_f}{C_f+C_L}\right)\cdot \frac{\frac{1}{\beta }\left(V_{xstep}+\sqrt{2}\cdot \frac{I_{ss}}{g_m}\right)}{\frac{C_s}{C_f}\cdot V_{istep}}$$
$$t_{lin,eq}=-\tau \cdot ln\left(1-\frac{C_L+\left(1-\beta \right)C_f}{C_f+C_L}\cdot \frac{SR\cdot t_{slew}}{\frac{C_s}{C_f}\cdot V_{istep}}\right)=-\tau \cdot ln\left(1-\left(1-\beta \frac{C_f}{C_f+C_L}\right)\cdot \frac{\frac{1}{\beta }\left(V_{xstep}+\sqrt{2}\cdot \frac{I_{ss}}{g_m}\right)}{\frac{C_s}{C_f}\cdot V_{istep}}\right)$$
运放脱离压摆后的线性建立时间为
$$\tau \cdot \left[\left(N+1\right)ln2+ln\left(1+k_f\right)\right]-t_{lin,eq}$$
总建立时间为
$$t_{slew}+\tau \cdot \left[\left(N+1\right)ln2+ln\left(1+k_f\right)\right]-t_{lin,eq}$$
注意到，当闭环放大器输入加入$V_{istep}$的阶跃后，OTA输入端电压变化为
$$V_{xstep}=\frac{C_s}{C_s+C_{in}+\frac{C_f{C}_L}{C_f+C_L}}\cdot V_{istep}<V_{istep}$$
当闭环增益（H(0)）足够大时，
$$\frac{V_{oswing}}{H(0)}>V_{istep}>V_{xstep}$$
当$\sqrt{2}\cdot \frac{I_{SS}}{g_m}>\frac{V_{oswing}}{H(0)}$，即满足
$$\frac{g_m}{I_D}=\frac{2}{V_{ov}}<2\sqrt{2}\cdot \frac{H(0)}{V_{oswing}}$$
$$V_{ov}>\frac{1}{\sqrt{2}H(0)}\cdot V_{oswing}$$
运放不会进入压摆的状态。

6.虚地点电压

假设建立过程中，OTA栅端始终不存在瞬态电流，则阶跃电压加入的瞬间，各个节点的瞬态压降由电容分压决定（认为运放尚未开始工作），电容网络如下图所示

根据KCL
$$\left(V_i-V_s\right)\cdot s{C}_s=V_x\cdot C_{par}+\left(V_x-V_o\right)\cdot C_f$$
$$V_x=\frac{V_i{C}_s+V_o{C}_f}{C_s+C_f+C_{par}}$$
仅考虑线性建立时
$$V_{out}\left(t\right)=-\frac{C_s}{C_f}\cdot V_{istep}\cdot \frac{T\left(0\right)}{1+T\left(0\right)}\cdot \left[1-\left(1+\frac{{\omega }_p}{{\omega }_z}\right)e^{-\frac{t}{\tau }}\right]$$
则虚地点的瞬态电压为
$$V_x\left(t\right)=\frac{C_s}{C_s+C_f+C_{par}}\cdot V_{istep}\cdot \left(1-\frac{T\left(0\right)}{1+T\left(0\right)}\left[1-\left(1+\frac{{\omega }_p}{{\omega }_z}\right)e^{-\frac{t}{\tau }}\right]\right)$$
忽略静态建立误差，则虚地点瞬态电压可以写为
$$V_x\left(t\right)=\frac{C_s}{C_s+C_f+C_{par}}\cdot V_{istep}\cdot \left(1+\frac{{\omega }_p}{{\omega }_z}\right)e^{-\frac{t}{\tau }}$$

Appendix

matlab建模，根据建立时间和建立精度要求迭代求解跨导

% Parameter analysis for SC-INT
% Consisting slew and linear settling
clear all
fs = 4e6;
tsettle = 0.9 * 1/2* 1/fs;
N_bit = 24;    % settling resolution
cs = 300e-15;    % equivalent sampling capacitor
cin1 = 1.8e-12;
cref = 2.4e-12;
cint1 = 12e-12;
cl = 2e-12;    % load cap
vistep = -16;

iter = 200;
gm = ones(1,iter+1)*1e-3;    % start from 1mS
gmoverid = 12;
cggovergm = 0.067e-15/(1e-6)    % gate cap w.r.t gm (fF/uS)

for i = 1:iter
    cpar = gm(i) * cggovergm + cin1 + cref - cs;
    cjunc = cpar/8 * 2;    % cdd of cascode FETs, approx 1/8 of Cgg for n/p MOSTs
    beta = cint1/(cs + cint1 + cpar)    % feedback ratio
    cltot = cl + (1-beta) * cint1;

	% input voltage step feedforward
	vxstep = vistep * cs / (cs + cpar + cint1 * cl / (cint1 + cl));    % virtual ground voltage
	vostep = vxstep * cint1 / (cint1 + cl);
	vo = zeros(1,ceil(tsettle * 1e9/0.1)) + vostep;    % calculate Vo according to settling time

	% slew settling
	Iss = 2 * gm(i) / gmoverid;
	SR = Iss / cltot;

	tau = cltot / (beta * gm(i));

	if abs(vxstep) <= sqrt(2) * Iss / gm(i)    % without slew
		tslew = 0;
		voslew = 0;
		tlineq = 0;
	else
	    tslew = (abs(vxstep) - sqrt(2) * Iss / gm(i)) / (beta * SR);
	    for j = 1:ceil(tslew*1e11)
	    	t = j / 1e11;
	    	vo(j) = vostep + t * SR;
	    end
	    voslew = SR * tslew + vostep;
	    tlineq = -tau * log(1 - (1 - beta * cint1 / (cint1 + cl)) * 1/beta * (vxstep + sqrt(2) * Iss / gm(i)) / (cs / cint1 * vistep));    
	end
	% linear settling
	fbw = 1/(2 * pi * tau);
	GBW = fbw / beta;
	tlin = (N_bit + 1) * log(2) * tau + tau * log((cl + cint1) / (cl + (1 - beta) * cint1)) - tlineq;
	for j = ceil(tslew * 1e11) + 1: ceil(tslew * 1e11)
		t = (j - ceil(tslew * 1e11)) / 1e11;
		vo(j) = -cs / cint1 * vistep * (1 - (1 + beta * cint1 / (cl + (1 - beta) * cint1)) * exp(-(tlineq + t) / tau));
	end
	ttot = tlin + tslew;
	gm(i + 1) = (1 + (ttot - tsettle) * 1e6) * gm(i);    % LMS based error feedback
end

vx = (vistep * cs + vo * cint1) / (cs + cint1 + cpar);    % KCL, virtual gnd voltage

figure(1)
plot(0.01 : 0.01 : ceil(tsettle * 1e11) / 100, [0 vo(1 : end - 1)], 'linewidth', 2); hold on
plot(0.01 : 0.01 : ceil(tsettle * 1e11) / 100, [0 vx(1 : end - 1)], 'linewidth', 2); hold on
plot([tslew tslew] * 1e9 + 1, [-5 2], 'linestyle', '--', 'color', 'g'); hold on
plot([0 tslew * 1e9 + 15], [0 0], 'linestyle', '--', 'color', 'm'); hold on
plot([0 tslew * 1e9 + 15], [cs / cint1 * (-vistep), cs / cint1 * (-vistep)], 'linestyle', '--', 'color', 'c'); hold off

xlabel('Settling Time /ns');
ylabel('Transient Voltage /V');
title('Transient Voltage of OTA output');
axis([0 tsettle * 1e9 + 15 -2.5 1.5]);
legend('Vout', 'Vx', 'slew time', 'loaction', 'southeast');

figure(2)
err_dynamic = 20 * log10(1 - vo ./ (-cs / cint1 * vistep));
plot(0.01 : 0.01 : ceil(tsettle * 1e11) / 100, err_dynamic); hold on
plot([tslew tslew] * 1e9 + 1, [-160 10], 'linestyle', '--', 'color', 'g'); hold on
line([0.01, ceil(tsettle * 1e11) / 100], [20 * log10(1/2 * 1/2^20), 20 * log10(1/2 * 1/2^20)]); hold off

xlabel('Settling Time /ns');
ylabel('Dynamic Settling Error /dB');
title('Dynamic Settling Error');
axis([0 tsettle*1e9+15 -160 10]);

瞬态建立波形如下图所示，红色曲线为虚地点电压变化，蓝色曲线为输出电压变化。

动态建立误差如下图所示，其中蓝色直线为20bit建立精度。