企业官网建设流程全解析

简介：单相锁相环（PLL）MATLAB/Simulink仿真资源，采用C语言S-Function实现数字PI控制器核心算法，面向电力电子、通信系统设计人员及高校相关专业学生。资源提供完整的Simulink仿真模型，涵盖鉴相器、数字PI控制器、压控振荡器等标准环节，可清晰观察相位误差收敛、锁定时间及稳态误差等动态行为，适用于载波同步、频率合成和时钟恢复等典型应用。压缩包共8个文件，其中.slx模型可直接运行，.c源文件便于修改算法，.mexw64动态库加速仿真，辅以.tlc代码生成支持文件，整体仅58KB，结构紧凑。已有2116人学习下载，可通过调整PI参数研究不同带宽与阻尼特性，也可学习C语言与Simulink混合建模方法。从工程角度看，该模型可作为数字锁相环算法验证的起点，帮助读者快速理解比例积分控制在相位同步中的作用，并进一步扩展到三相锁相环设计。 前段时间做并网逆变器项目，单相锁相环（PLL）是绕不开的坎。控制系统每毫秒都要知道电网电压的实时相位，而这套单相锁相环最终我用 MATLAB/Simulink 搭模型验证，又用 C 语言写了 S-Function 嵌回 Simulink 做代码级联合仿真，整个流程走完，收获很大。这篇内容就把我从方案选型、模型搭建、参数整定，到 S-Function 落地和问题排查的完整过程记录下来，写给做并网控制、电源开发、电机驱动，以及想把 Simulink 算法真正移植到嵌入式平台的朋友参考。

1. 单相锁相环：为什么说“锁相”不只是锁频率

1.1 并网控制里的核心矛盾

做并网设备的工程师都知道，逆变器要安全并网，控制器的 PWM 驱动信号必须和电网电压保持严格同步。这个同步说的不只是频率一致，而是相位实时一致。电网不是理想世界里的完美正弦波，它会有谐波、会有暂态跌落、频率也会在额定值附近小范围波动。锁相环要做的事，就是从一路被污染的正弦波里，实时估计出基波分量的相位、频率，甚至幅值。

三相系统里这件事儿相对好办，Clark 变换把三相静止坐标变成两相静止坐标，αβ 两路信号天然正交，角度直接算出来。单相系统就比较尴尬，物理上只有一路火线，你拿到的只有孤零零一路正弦信号，缺少另一路正交参考。这就好比人只有一只眼睛，能看见东西但判断不出纵深，必须靠其他办法把缺失的维度补出来。所有单相锁相环的难点，归根结底就是一件事：怎么构造出那路“不存在的正交信号”。

1.2 常见单相锁相方案怎么选

我实际梳理过几种主流方案，各有各的适用场景，关键看你的动态要求、谐波环境、计算资源和实现成本。

过零检测听起来最简单，实际并网场景里基本撑不住。谐波会让过零时刻来回抖动，相位输出像锯齿一样，而且每半个周期才更新一次，动态性能远远不够。T/4 延时法在工频 50Hz 固定时还好，一旦电网频率偏移，延时的角度就不再是 90°，锁相环输出会出现明显误差。希尔伯特变换精度高，但离散实现要挂一堆全通滤波器，DSP 上资源吃紧。最后我选了 SOGI 方案，综合看它兼顾了滤波能力、动态响应和嵌入式实现成本，工程上最靠谱。

1.3 为什么选 SO GI + S-Function 这条路线

锁相环本身有现成模块，比如 Simulink 的 PLL 库里面就有三相锁相环，但单相锁相环没有“开箱即用”的标准模块，而且很多工程场景需要你把算法移植到自己的控制板里，不能总依赖 MATLAB 环境。

我采用 S-Function 的核心原因有两条：第一，C 语言 S-Function 可以让我直接验证最终要在 DSP 或者 MCU 上运行的那一份 C 代码，仿真通过就意味着大部分算法逻辑已经确认；第二，S-Function 可以和 Simulink 里其他控制模块、被控对象模型联合仿真，比如接上一个并网逆变器主电路，看看锁相环在真实工况下的表现，比单独测一个算法模块有意义得多。SOGI 只是锁相环前端构造正交信号的手段，真正把它变成“能落地的代码”，S-Function 那条路最直接。

2. Simulink 参考模型搭建：让算法先转起来

2.1 一条清晰的信号流：从正弦波到相位角

在动笔写 C 代码之前，我习惯先在 Simulink 里把算法完整跑通，确认控制结构、参数整定和动态性能都没问题。整个单相锁相环的信号流是这样的：电网电压正弦波送进 SOGI，SOGI 输出两路正交信号 vα 和 vβ；vα、vβ 做 Park 变换，转换到 dq 旋转坐标系；q 轴分量 vq 作为相位误差信号送入 PI 调节器；PI 输出角频率增量，叠加上前馈角频率后积分，得到估计相位 θ；θ 再反馈回 Park 变换，形成闭环。

用一个生活化的类比，锁相环就像一个逐渐对准目标的自动云台：Park 变换负责测量“现在偏了多少度”（vq 就是偏差量），PI 调节器负责输出一个“转动的速度”，积分器负责把速度累加成角度。只要 vq 收敛到 0，就意味着云台已经稳稳对准了电网电压的相位。

2.2 SOGI 模块：正交信号怎么构造

SOGI 本质上是一个谐振在基波频率附近的二阶带通滤波结构。它的关键是让输出 vα 和输入同相，让另一路输出 vβ 严格滞后 90°，这样就有了一对正交基。用传递函数表示，从输入 vin 到 vα 和 vβ 的传递函数分别是：

Gα(s) = k·ω₀·s / (s² + k·ω₀·s + ω₀²)

Gβ(s) = k·ω₀² / (s² + k·ω₀·s + ω₀²)

这里 ω₀ 是谐振角频率，k 是阻尼系数。工程上 k 取 √2（约 1.414）比较常见，对应阻尼比 0.707，动态响应和滤波效果都比较均衡。k 越大，响应越快，但对谐波的抑制越差；k 越小，滤波越干净，但锁定过程会变得迟钝。我在模型里先用前端增益归一到 1 的理想情况，方便整定 PI 参数，后面再接实际电压幅值。

Simulink 里搭 SOGI 最直接的方法是用两个积分器构造状态空间。设两个状态变量 x1、x2，状态方程写成：

dx1/dt = x2

dx2/dt = -ω₀²·x1 - k·ω₀·x2 + vin

输出：

vα = k·ω₀·x2

vβ = k·ω₀²·x1

这个结构在 Simulink 里用两个 Integrator 模块和几个 Gain、Sum 模块就能搭出来，也可以在 S-Function 里实现。需要注意，x1、x2 的初始值要设为 0，否则启动瞬间会有短暂的自由振荡。我习惯再检查一下输出 vβ 是不是真滞后 90°——把两路信号用 Scope 拉出来比较，如果相位差明显不是 90°，多半是状态方程接反了。

2.3 Park 变换与 PI 参数整定

Park 变换是单相锁相环里比较绕的一步，本质是坐标旋转。把 vα、vβ 和估计角度 θ 代入旋转矩阵，得到 vd、vq。我这里用的变换关系是：

vd = vα·cosθ + vβ·sinθ

vq = -vα·sinθ + vβ·cosθ

当 θ 等于真实相位 θg 时，vd 等于输入幅值 Vm，vq 等于 0。所以 vq 就是天然的系统误差信号，把它送给 PI，PI 输出频率修正量，积分得到 θ，形成负反馈闭环。这里最容易搞错的是 vq 的符号方向，如果符号接反，锁相环会稳定在相差 180° 的状态，表现就是 θ 始终和真实相位差半个周期。

PI 参数整定我有现成的公式可以分享。锁相环线性化之后，vq ≈ Vm·Δθ，其中 Δθ 是相位误差。整个环路的开环传递函数近似为：

G_ol(s) = Vm·(Kp + Ki/s)·(1/s) = Vm·(Kp·s + Ki)/s²

闭环特征方程是 s² + Vm·Kp·s + Vm·Ki = 0，对照标准二阶系统 s² + 2ζωn·s + ωn² = 0，就能反推出：

Kp = 2ζωn / Vm

Ki = ωn² / Vm

我在模型里按标幺值处理，Vm 取 1，阻尼比 ζ 取 0.707，锁相环固有角频率 ωn 取 2π×20 rad/s（约 20Hz 带宽），算出来 Kp ≈ 177.7，Ki ≈ 15791。如果希望动态响应更快，可以把 ωn 提高到 2π×50，对应 Kp ≈ 444.2，Ki ≈ 98696。但带宽也不是越高越好，带宽太高会把谐波和噪声带进环路，输出角度就会抖。实际电网谐波大，我一般把锁相环带宽控制在 20～50Hz 之间再配合 SOGI 的滤波效果取折中。

3. C 语言 S-Function：从仿真模型到嵌入式代码

3.1 为什么非要用 S-Function

有些朋友问我，Simulink 模型跑得好好的，为什么要费劲写 C 语言的 S-Function？我的回答是：模型是“理想版本”，而最终要跑在控制板上的，是一段只能顺序执行的 C 代码。模型里连续积分器、理想 PI 的做法，到了嵌入式环境要变成离散递推、限幅、抗饱和处理，这些细节如果不提前验证，等到硬件联调再去查，代价会大得多。

S-Function 让我在 Simulink 环境里“模拟”一块真正的嵌入式芯片。它按固定的采样周期被反复调用，输入是当前时刻的采样电压，输出是算好的相位和频率，内部状态必须自己保存——这和单片机中断服务程序里的写法几乎一模一样。仿真通过，再把这一段代码搬到工程里，改一改接口就能跑，风险基本可控。

3.2 S-Function 的回调函数框架

C 语言 S-Function 的完整模板可以在 MATLAB 安装目录下找到，核心要了解几个回调函数：mdlInitializeSizes 定义输入输出端口、状态个数和参数个数；mdlInitializeSampleTimes 设置采样周期；mdlStart 初始化状态；mdlOutputs 计算输出；mdlUpdate 更新离散状态；mdlTerminate 做清理工作。

我用的锁相环 S-Function 有 1 个输入（电网电压采样值）、4 个离散状态（SOGI 两个状态、PI 积分器一个、相位角一个）、4 个输出（估计相位 θ、估计频率 f、vα、vβ）。采样周期设为 1e-4 秒，也就是 10kHz，和常见并网逆变器的控制频率保持一致。这样仿真里每个控制周期执行一次，和后续 DSP 上定时器中断调用的节奏完全相同。

有一个容易踩坑的地方：离散 S-Function 要求求解器选择定步长离散求解器，如果你用了变步长求解器，仿真会报错。我在 mdlInitializeSampleTimes 里设置采样时间后，Simulink 界面里 Solver 要选 discrete，Fixed-step size 设为 Ts。这点在联调阶段经常被忽略。

3.3 核心 C 代码实现与离散化细节

写 C 代码的时候，连续域的积分器要用离散递推公式替代。为了简单可靠，我用了前向欧拉离散，采样率 10kHz 时精度完全够用。SOGI 部分的 C 代码如下，状态数组 state 的前四个分量分别对应 x1、x2、PI 积分器输出、相位 θ：

#define PLL_PI 3.141592653589793 static real_T pll_step(real_T *state, const real_T *param, real_T vin, real_T *theta, real_T *freq, real_T *v_alpha, real_T *v_beta) { /* param 依次存放：Ts, k, w0, Kp, Ki, wff */ real_T Ts = param[0]; real_T k = param[1]; real_T w0 = param[2]; real_T Kp = param[3]; real_T Ki = param[4]; real_T wff = param[5]; real_T x1 = state[0]; real_T x2 = state[1]; real_T integ = state[2]; real_T theta_ = state[3]; /* 1. SOGI 状态更新 */ real_T dx1 = x2; real_T dx2 = -w0 * w0 * x1 - k * w0 * x2 + vin; x1 += Ts * dx1; x2 += Ts * dx2; /* 2. 由状态计算正交信号 */ real_T va = k * w0 * x2; real_T vb = k * w0 * w0 * x1; /* 3. Park 变换，得到 q 轴误差 */ real_T ct = cos(theta_); real_T st = sin(theta_); real_T vq = -va * st + vb * ct; /* 4. PI 积分器更新，带抗饱和限幅 */ integ += Ki * Ts * vq; if (integ > 40.0) integ = 40.0; if (integ < -40.0) integ = -40.0; /* 5. 频率叠加前馈，角度积分 */ real_T w = wff + Kp * vq + integ; theta_ += w * Ts; if (theta_ > 2.0 * PLL_PI) theta_ -= 2.0 * PLL_PI; state[0] = x1; state[1] = x2; state[2] = integ; state[3] = theta_; *theta = theta_; *freq = w / (2.0 * PLL_PI); *v_alpha = va; *v_beta = vb; return vq; }

这段代码有几个细节值得展开说。PI 积分器的限幅我放在 ±40 rad/s，也就是约 ±6.4Hz，这个限幅既能保证频率跳变时环路快速响应，又不会让积分器无限制饱和导致恢复迟缓。角度积分用了“超过 2π 就减去 2π”的包卷写法，避免浮点数无限增长。SOGI 的离散用前向欧拉，在采样率 10kHz 时极点离稳定边界很远，完全没问题；如果你把采样率降到 1kHz 以下，就要换用双线性变换离散，否则系统可能发散。

S-Function 的 mdlOutputs 里，我把当前状态直接换算成 θ、f、vα、vβ 输出；mdlUpdate 里才执行更新状态的操作。这样离散系统天然会有一拍滞后，即输出用的还是上一周期的状态，但 10kHz 采样率相对 50Hz 基波来说滞后量极小，工程上可以接受。如果要求更严格，可以在输出前先预测一步状态，不过会牺牲代码清晰度，我一般不做。

4. 仿真测试：四种工况验证锁相性能

4.1 测试场景与设置

模型和 S-Function 都跑起来之后，验证环节我设计了四种典型工况，覆盖稳态、动态、鲁棒性三个维度。测试一，输入 50Hz、幅值 1 的纯净正弦波，看稳态锁定情况；测试二，让频率在某个时刻从 49Hz 突变到 51Hz，看频率跟踪速度；测试三，让相位在 0.1 秒处跳变 30°，看相位纠偏能力；测试四，在基波上叠 5% 的三次谐波和 2% 的五次谐波，看谐波环境下锁相输出是否干净。

输入信号我直接用 Simulink 的 Sine Wave 模块生成，再叠加几个不同频率的正弦分量模拟谐波污染。S-Function 模块和前面搭的参考模型并联连接，同一个输入同时送给两套锁相环，方便直接对比。

4.2 波形怎么看、指标怎么算

锁相环性能好不好，我主要看三个量：vq 波形、θ 波形、f 波形。vq 收敛到 0 的速度和过程，直接反映锁相环的动态品质；θ 应该是随时间线性增长的锯齿波，叠在真实相位上，两者相减才看得出来误差；f 波形则直观展示频率估计的正确性和噪声水平。

以一个典型测试结果为例：输入标准 50Hz 正弦波，初始相位为零，PLL 在 0.01 秒内完成初始锁定，vq 从 1 快速收敛到 0，超调很小。频率突变测试里，49Hz 跳到 51Hz 后，约 40~60ms 内频率估计能回到新的稳定值，这跟 20Hz 环路带宽的设计预期基本一致。如果你看到 vq 收敛到 0 的时间太长，或者频率波形一直在振荡，就需要回头检查 PI 参数或者 SOGI 阻尼系数。

有一个容易忽略的问题：仿真开始时如果输入相位不是 0，而锁相环初始角度默认是 0，vq 起步就不是 0，锁相环要靠积分器慢慢转过去，锁定时间会变长。工程上可以在初始化时把 θ 设成 0，也可以利用过零检测粗略估计一个初始相位再启动锁相环。如果两者差距超过 180°，可能会先朝反方向追，直观表现就是锁定时间变长甚至来回摆动。

4.3 S-Function 与 Simulink 模型结果对比

我把 S-Function 模块和纯 Simulink 连续模型并联跑，结果很接近，主要在动态过程中有小幅差异。原因是参考模型用的是连续积分器，而 S-Function 用离散前向欧拉，天然有一拍滞后。稳态时两条曲线差别小到几乎重叠，只有相位输出之间有恒定的小相位差，大概在几度以内，这个量级的差距在工程上是允许的。

如果发现两条曲线差异明显，优先检查三件事：采样时间是否一致，离散化方法是否一致（一个连续一个离散当然会有偏差），以及 S-Function 里状态初始值和模型是否一致。我在调试时就遇到过状态初始值没设置的情况，SOGI 两个积分器从 0 开始，而参考模型另一个初值，结果启动阶段差很多。

5. 常见问题与排查技巧

5.1 编译和集成问题

我刚开始编译 C MEX S-Function 时，最常遇到的就是编译器没配置好。MATLAB 默认不认电脑上装的 VS，需要先跑一遍 mex -setup，把自己安装的编译器和 MATLAB 匹配起来。如果还没装编译器，装一个 VS Build Tools 的 C++ 桌面开发组件，一般就够了。

还有一种情况是模块名字对不上导致仿真报错。你在 Simulink 里放一个 S-Function 模块，模块参数里填的函数名必须和 C 文件里的 #define S_FUNCTION_NAME 一致。比如 C 文件里写的是 #define S_FUNCTION_NAME pll_sfun，那么模块参数就要填 pll_sfun，不能填文件名 pll_sfun.c。这个低级错误我见过不少人踩过。

另外提醒一下，更新了代码重新编译之前，必须先清除 MATLAB 工作空间里旧版本的 mex 文件，否则 mex 命令会提示文件被占用，或者一直在用旧版本跑。我习惯在编译命令前加 delete(['pll_sfun.', mexext]) 或者手动清一下工作区。

5.2 锁不住和锁太慢的问题

锁相环锁不住，我在调试中遇到的情况绝大多数是反馈符号问题。vq 的符号接反，锁相环会稳定在相位差 180° 的位置，波形看起来也“锁住”了，但实际上是锁反了。排查方法很简单：把输入信号和 θ 拉到同一个 Scope 里看，如果 θ 对应的过零时刻和输入差半个周期，说明符号反了。

锁得太慢则要从三方面看：PI 带宽是否太低、初始频率是否偏离太远、PI 积分限幅是否太窄。如果输入是 45Hz 的信号，而前馈频率按 50Hz 设置，PI 积分器要从 0 一点一点积出 5Hz 的修正量，如果限幅设成 ±2Hz，可能永远追不上。我一般在程序里把前馈频率做成一个可配置的输入参数，做的工程里工频稳定就写死 50Hz，做测试就附加一个外部变量可以动态修改。

5.3 稳态误差和谐波纹波问题

稳态相位误差不为零，优先检查 PI 积分器是否被限幅限制住了。另外，前馈角频率如果和实际电网频率不一致，稳态时 PI 需要持续输出一个固定修正量来维持锁定，如果积分器精度不足或限幅过窄，就会留下残差。解决方法是把前馈频率调准，或者适当放宽积分限幅。

输出角度有纹波，常见原因是环路带宽太高，把 SOGI 没有滤干净的谐波带进了反馈。我在测试四里叠加三次、五次谐波，如果锁相环带宽取 50Hz，θ 波形会有明显抖动；把带宽降回 20Hz 左右，抖动就小很多。谐波特别大的场合，还可以在 vq 后面再加一层低通滤波，或者用带陷波器的改进型 SOFI-PLL，但实现复杂度会上升，具体取舍看你实际谐波环境。

最终再分享一个我的习惯：锁相环的所有参数，包括 Ts、k、ω₀、Kp、Ki、前馈频率，都通过 S-Function 的入口传进来，不要在代码里写死。我当时图省事把几个系数定义为常量，后面调整参数每次都要重新编译，特别耽误事。改成参数接口之后，配合 MATLAB 脚本批量扫描参数，效率和安全感都上来了。你刚上手的话，建议一开始就留好参数接口，省得后面返工。

本文还有配套的精品资源，点击获取