零基础使用MATLAB求解三阶电荷泵锁相环CPPLL相关系统参数

madson2022

已于 2023-06-09 15:01:55 修改

阅读量1k

点赞数

文章标签： matlab 开发语言

于 2023-06-09 15:00:39 首次发布

本文链接：https://blog.csdn.net/madson2012/article/details/131125973

版权

参考资料：

4. 三阶电荷泵PLL稳定性分析_哔哩哔哩_bilibili

一、三阶锁相环模型

每个模块的传递函数如下：

由每个模块传递函数，得到环路增益传递函数，于是便可以得到零点、极点。

根据传递函数，也可知，这是一个三阶系统，三个极点，一个零点

二、理论计算

那么，现在是已知一个系统的传递函数：

零、极点表达式为：

目标函数为相位裕度：

现在的问题转化为：已知G(s)，求R1/C1/C2，使PM在增益交点处满足要求（其中N、Ip、Kvco，这三个量为已知设计量）。

第一步：根据系统要求（GBW/相位裕度），求零点、极点。

对目标函数(相位裕度PM)求导，使取极大值时的频率 $\omega$ 等于GBW，且极大值为相位裕度要求值（如60度）。这样就可以求得零点、极点（用GBW、相位裕度表示）。

第二步：根据零、极点，求C1/C2/R1

最终结果为：

三、MATLAB代码

%-----程序功能---------%
%根据给写的系统需要参数、以及相关模块的设计参数，求出R1，C1，C2

clear
%系统需求参数
wu=2*pi*0.8e6;% 带宽，单位：rad/s，通常小于输入参考时钟的1/10
pm=60/180*pi;%相位裕度，单位：rad

%设计参数
N=10;%最大分频比
Icp=44e-6;%电荷泵充放电电流，单位：A
Ios=0;%电荷泵失调电流，单位:A
Kpd=Icp/2/pi; %鉴频鉴相器增益
Kvco=250e6*2*pi;%压控振荡器增益，单位：rad/s/V

%计算相关数值
wp=wu*cos(pm)/(1-sin(pm));%计算出极点，单位：rad/s
wz=wu^2/wp;%计算出零点，单位：rad/s
C2=1/N*Kpd*Kvco/wu^2*sqrt(wz/wp) %单位：F
C1=(wp/wz-1)*C2 %单位：F
R1=1/C1/wz %单位：欧姆
%% 

%绘制波特图
G=1/N*Kpd*Kvco*tf([R1*C1 1],[R1*C1*C2 C1+C2 0 0]);%PLL环路增益传输函数
p=bodeoptions;
p.FreqUnits='Hz';
bode(G,p)
grid on;
%% 

%绘制闭环系统的单位阶跃响应
Gf=feedback(G*N,1/N); %设置闭环系统
figure
step(Gf,20e-6); 
grid on;

仿真结果图：