请用MATLAB语言写一个可以运行的锁频环代码

以下是一个简单的锁频环代码示例，包含了一个二阶锁频环和一个简单的正弦波信号作为输入：

% 锁频环参数
Kp = 1; % 比例增益
Ki = 1; % 积分增益
Kd = 1; % 微分增益
N = 2; % 阶数
fref = 10; % 参考频率
fs = 100; % 采样频率
T = 1/fs; % 采样周期
t = 0:T:1-T; % 时间向量

% 生成输入信号
fin = 3; % 输入信号频率
Vin = sin(2*pi*fin*t); % 输入信号

% 初始化锁频环状态
theta = 0; % 相位
omega = 0; % 频率偏差
Vout = 0; % 输出信号

% 锁频环迭代
for i = 1:length(t)
    % 计算频率偏差
    omega = 2*pi*fin - theta;
    
    % 计算PID控制器输出
    Vpid = Kp*omega + Ki*sum(omega*T) + Kd*(omega - Vout)/T;
    
    % 更新相位
    theta = theta + (2*pi*fref + Vpid)*T;
    
    % 生成输出信号
    Vout = sin(theta*N);
    
    % 绘制结果
    plot(t(1:i), Vin(1:i), t(1:i), Vout(1:i))
    xlabel('时间 (s)')
    ylabel('幅度')
    legend('输入信号', '输出信号')
    drawnow
end

该代码实现了一个简单的二阶锁频环，可以通过调整比例、积分和微分增益来优化锁定效果。输入信号可以通过修改频率和幅度来模拟不同的实际场景。