关于锁相环的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他没写清楚是什么，我觉得是输入频率或者角频率
[size=13.3333px]

[size=13.3333px]附上他的论文：
Analysis_of_charge-pump_phase-locked_loops.pdf(927.49 KB, 下载次数: 0 ) [size=0.83em]1 分钟前 上传
点击文件名下载附件
下载积分: 资产 -2 信元, 下载支出 2 信元


[size=13.3333px]

大佬们帮忙看看，谢谢啦

linxin2018

顶一下

laojun001

就这样一点一点学下去，坚持3-5年你就感觉不一样了，会有点底气了，尤其对底层技术有那么一点自信了

xuwenwei

新人回答，你想画出波特图，你直接abs计算绝对值是可以表示z函数的幅值的吗？

linxin2018

xuwenwei 发表于 2025-7-15 17:26
新人回答，你想画出波特图，你直接abs计算绝对值是可以表示z函数的幅值的吗？ ...

抱歉，我不是很懂这个，请问正确的做法是什么呢

tanborui123

你把w转置为列向量再跑一下