锁相环matlab代码

以下是一个简单的锁相环（PLL）的MATLAB代码，用于将输入信号与参考信号同步：

% 锁相环参数设置
f_ref = 1000;   % 参考信号频率
f_in = 100;     % 输入信号频率
Kp = 1;         % 比例增益
Ki = 0.01;      % 积分增益
Kd = 0.1;       % 微分增益
N = 1000;       % 采样点数

% 生成参考信号
t = linspace(0, 1, N);
ref_signal = sin(2*pi*f_ref*t);

% 生成输入信号
in_signal = sin(2*pi*f_in*t);

% 初始化PLL变量
theta = zeros(1, N);    % 相位
theta_dot = zeros(1, N);% 相位变化率
theta_ddot = zeros(1, N);% 相位变化率变化率
error = zeros(1, N);    % 相位误差
int_error = 0;          % 积分误差

% PLL循环
for i = 2:N
    % 计算相位误差
    error(i) = angle(ref_signal(i-1)*conj(exp(1j*theta(i-1)))) ;
    
    % 积分误差
    int_error = int_error + error(i);
    
    % 计算相位变化率和变化率变化率
    theta_dot(i) = Kp*error(i) + Ki*int_error + Kd*(error(i) - error(i-1));
    theta_ddot(i) = (theta_dot(i) - theta_dot(i-1))/(t(i) - t(i-1));
    
    % 更新相位
    theta(i) = theta(i-1) + theta_dot(i)*(t(i) - t(i-1)) + 0.5*theta_ddot(i)*(t(i) - t(i-1))^2;
end

% 绘制结果
figure;
subplot(2,1,1);
plot(t, in_signal, t, ref_signal);
title('Input and Reference Signals');
legend('Input Signal', 'Reference Signal');
xlabel('Time (s)');
ylabel('Amplitude');

subplot(2,1,2);
plot(t, theta);
title('Phase of Input Signal');
xlabel('Time (s)');
ylabel('Phase (rad)');

这段代码首先生成参考信号和输入信号，并初始化PLL变量。然后，在循环中计算相位误差、积分误差、相位变化率和变化率变化率，并使用欧拉法更新相位。最后，绘制输入信号的相位随时间的变化，以展示锁相环的同步效果。

请注意，这只是一个简单的示例代码，实际应用中需要根据具体的系统和信号特性进行参数调整和优化。