关于锁相环的z域传输函数画幅频响应一直画不对的问题

linxin2018

Halumolu的二阶锁相环z域，我按他论文里的公式和参数写了一个matlab代码，一直没法画出他的原图，是我哪里写的有问题呢？

Halumolu的公式：

其中：



他的结果：

我的代码：
KVCO=1e6;
Icp=130e-6;
C1=451e-12;
C2=14.5e-12;
R=10e3;
wref=2e6;
T=1/wref;
w=1:10:1e8;
f=w/2/pi;
K=KVCO*Icp*R*C1/(wref*(C1+C2));
wp3=1/(R*(C1*C2/(C1+C2)));
a=exp(-wp3*T);
z=exp(1j*w*T);
coe1=(C1*(1-a)/(C1+C2)+2*pi/R/C1/wref);
coe2=(C1*(1-a)/(C1+C2)+2*pi/R/C1/wref*a);
LoopGain=K*(z.^2*coe1-z*coe2)./(z.*z.*z-z.*z*(2+a)+z*(1+2*a)-a);
semilogx(f,20*log10(abs(LoopGain)));


w_i他没写清楚是什么，我觉得是输入频率或者角频率


附上他的论文：
Analysis_of_charge-pump_phase-locked_loops.pdf (927.49 KB , 下载次数: 21 )

2025-7-11 03:01 上传

点击文件名下载附件

大佬们帮忙看看，谢谢啦