二阶锁相环matlab代码

以下是一个简单的二阶锁相环（PLL）的 MATLAB 代码示例：

clc;
clear all;

% 模拟输入信号
fs = 100e3; % 采样率
t = 0:1/fs:0.1; % 时间范围
f0 = 1e3; % 输入信号频率
x = sin(2*pi*f0*t); % 输入信号

% 设计二阶锁相环参数
fc = 10e3; % VCO 中心频率
Kp = 1; % 相位检测器增益
Ki = 1; % 积分器增益
Kv = 1; % VCO 灵敏度
zeta = 0.707; % 阻尼比
wn = 2*pi*fc; % 自然频率
K1 = wn^2/Kv; % 一次环节增益
K2 = 2*zeta*wn/Kv; % 二次环节增益

% 模拟 PLL
theta = zeros(size(t)); % 相位误差
theta_dot = zeros(size(t)); % 相位误差导数
vco = zeros(size(t)); % VCO 输出
phi = zeros(size(t)); % 相位检测器输出
int = 0; % 积分器状态

for i = 2:length(t)
    % 相位检测器
    phi(i) = x(i)*sin(theta(i-1));
    % 积分器
    int = int + Ki*phi(i)/fs;
    % 二次环节
    theta_dot(i) = theta_dot(i-1) + K1*vco(i-1)/fs - K2*theta(i-1)/fs - K2*int/fs;
    % 一次环节
    theta(i) = theta(i-1) + theta_dot(i)/fs;
    % VCO 输出
    vco(i) = sin(2*pi*(fc + Kv*theta(i))*t(i));
end

% 绘制图形
figure;
subplot(2,1,1);
plot(t,x,t,vco);
title('输入信号和 VCO 输出');
xlabel('时间 (s)');
ylabel('幅度');
legend('输入信号','VCO 输出');
subplot(2,1,2);
plot(t,theta);
title('相位误差');
xlabel('时间 (s)');
ylabel('相位误差 (rad)');

该代码生成了一个包含输入信号和 VCO 输出的波形图以及相位误差的时间图。可以通过调整 PLL 参数来观察它们对输出的影响。