三次谐波失真的计算

论文中的一个基于源极跟随器的滤波器结构，我在计算谐波失真HD3遇到了困难。我的思路是先解出传输函数A，vo=Avi，a1=A,然后对A的表达式求三次导可得到a3，所以HD3=1/4*a3/a1*u^2。可是怎么计算传输函数呢?（低频不考虑电容,管子工作在亚阈值区）

各路大神看看啊