简介这是一份面向电力电子与控制方向研发人员的复合控制策略学习资料聚焦准比例谐振PR控制与改进重复控制相结合的光伏逆变器谐波抑制方案。针对光伏逆变器与有源滤波器拓扑结构相同的特性资料将并网与谐波补偿功能统一设计并用ip-iq谐波检测方法配合复合控制器兼顾动态响应速度与稳态精度。压缩包内为1个docx文档约47KB集中呈现准PR控制器、改进重复控制器的参数选取与推导过程以及可直接运行的Python实现与MATLAB/Simulink仿真流程从参考信号构造、扰动注入到频域波特图分析均有配套解释便于读者复现并对比不同场景下的性能表现。文档末尾还延伸讨论硬件在环测试框架、DSP核心代码与后续改进方向。目前已有85人学习适合具备一定控制理论基础、希望把复合控制策略落到仿真与工程实现的工程师按需查阅。1. 光伏逆变器谐波抑制为什么要走准PR加改进重复控制的复合路线光伏逆变器的电流环在三相不平衡、电网背景谐波偏大、本地非线性负载这几类工况下THD 很容易越过并网导则要求的 3% 红线。传统 PI 调节器跟踪 50 Hz 正弦指令时就带着幅值和相位误差增益给大一点又会顶到带宽上限。准 PR 把谐振点搬到基波靠谐振峰提供高增益压掉静差用截止频率 ωc 控制选择性代价是对 3、5、7 次这些主要低次谐波几乎没有增益电网侧的背景谐波原封不动穿透到电流波形里。重复控制依据内模原理在一个基波周期的延迟环节上对全部基波整数倍频率同时给出高增益能一次性按住整片低次谐波问题是那 N 拍延迟把动态响应拖慢负载或光照跳变时来不及反应。准 PR 并联改进重复控制本质是让快的支路管动态、慢的支路管稳态精度各拿各自的频段。适合正在做并网电流环、构网型变流器、APF 和 UPS 数字控制器的工程师和仿真人员。2. 准PR控制器在光伏逆变器电流环里的建模与参数整定2.1 从 PI 到准 PR交流量跟踪为什么必须引入谐振环节PI 的传递函数在虚轴上只有一个原点极点对直流给定能做到零静差对 50 Hz 正弦给定则天生有幅值和相位误差。把 PI 的频响代进单位负反馈误差传递函数在 ω0 处的模值由控制器增益决定增益越大误差越小但相位裕度同步下降。工程上想把 50 Hz 处误差压到 1% 以内Kp 往往要开到把截止频率推到接近开关频率的十分之一噪声放大和稳定裕度都不允许。PR 调节器在 ω0 处放置一对共轭极点理论上该频率处增益无穷大交流量跟踪误差为零。理想 PR 的谐振峰是冲激形状宽度为零电网频率一旦偏移到 49.5 Hz 或 50.5 Hz增益就掉得厉害。准 PR 在分母里加入阻尼项把尖峰展宽成一个有带宽的谐振平台G_PR(s) Kp 2·Kr·ωc·s / (s² 2·ωc·s ω0²)Kp 提供全频段的即时比例增益决定动态响应速度Kr 决定谐振峰高度ωc 决定谐振峰两侧的 -3 dB 带宽工程上按 ±0.5 Hz 频率漂移加裕量取常见范围是 π×5 到 π×15 rad/s。ωc 越大峰越宽越高频率适应性强但选择性变差对非基波频率的放大也更明显。2.2 准 PR 的离散化预畸变双线性变换与代码实现直接对 s 域做双线性变换会带来频率轴畸变采样率 10 kHz、基波 50 Hz 时畸变量很小但重复控制对相位极度敏感累积误差会吃掉相位裕度所以离散时做预畸变更稳妥。做法是先用 ω0 (2/Ts)·tan(ω0·Ts/2) 替换谐振角频率再做普通双线性变换使 z 域谐振点精确落在 50 Hz 上。import numpy as np from scipy import signal f0, fs 50.0, 10e3 # 基波频率、采样频率 Ts 1.0 / fs W0 2 * np.pi * f0 N int(round(fs / f0)) # 一个基波周期采样点数 200 # ---------- 准PR参数 ---------- Kp, Kr 6.0, 150.0 # 比例项、谐振项增益 wc 2 * np.pi * 6.0 # 谐振带宽约±6Hz # ---------- 预畸变后做双线性变换 ---------- w0p 2 / Ts * np.tan(W0 * Ts / 2) num [Kp, 2 * wc * (Kp Kr), Kp * w0p ** 2] den [1.0, 2 * wc, w0p ** 2] b_pr, a_pr signal.bilinear(num, den, fs) # ---------- 频响核对基波处应出现高峰谐波处回落到Kp ---------- for f in [50, 150, 250, 350]: _, h signal.freqz(b_pr, a_pr, worN[2 * np.pi * f * Ts]) print(f{f:4d} Hz 幅值 {abs(h[0]):8.2f} 相位 {np.degrees(np.angle(h[0])):8.2f}°)这段代码里num/den是把 Kp 与谐振项通分后的二阶多项式w0p完成预畸变signal.bilinear输出离散系数b_pr/a_pr。频响核对那几行是关键的自检手段50 Hz 处幅值应当在数十到上百量级150 Hz 及以上的幅值应当收敛到 Kp 附近。如果 50 Hz 幅值不够高说明 Kr 偏小或 ωc 偏窄如果谐波处幅值明显高于 Kp说明 ωc 给大了谐振平台已经溢到低次谐波频段。2.3 Kp、Kr、ωc 的整定顺序与取值影响三个参数不是并列关系整定有先后。先定 Kp它决定电流环的等效带宽和动态响应再定 ωc把谐振峰宽度卡到能覆盖电网频率漂移又不污染相邻谐波最后调 Kr把基波误差压到目标值同时观察相位裕度损失。参数典型取值增大后的效果过大时的代价Kp4 ~ 10动态响应加快全频段增益抬高相位裕度下降开关噪声放大Kr80 ~ 30050 Hz 谐振峰升高基波静差减小谐振峰附近相位骤变容易自激ωcπ×3 ~ π×15 rad/s谐振峰展宽抗频率偏移能力增强3、5 次谐波被一并放大THD 反而恶化我的习惯是先按 Kp L·ωc_loop 估一个初值其中 ωc_loop 取开关角频率的 1/10 到 1/15把 Kr 从小往大调每加一档就跑一次频响看 50 Hz 处峰值和相位裕度两个指标是否同时可接受。需要提醒的是准 PR 支路本身只对基波和紧邻带宽内有效指望靠加 Kr 去压 5、7 次谐波是走不通的那条路要留给下一章的重复控制来走。3. 改进重复控制器的内模结构、Q(z)滤波与相位补偿3.1 传统重复控制 1/(1-z⁻ᴺ) 的两个硬伤重复控制的内模来自内模原理要在基波整数倍频率上同时实现零稳态误差控制器里就得包含这些频率的极点。数字实现上最简单的形式是正反馈延迟环节M(z) 1 / (1 - z⁻ᴺ)N fs / f0它把上一周期的误差原样加到本周期的误差上误差在周期之间逐拍累积最终在基波整数倍频率处形成无穷大增益。这个结构有两个硬伤。第一内模在中高频段虽然不谐振但它的增益并不为零被控对象在这些频段的相位滞后很容易把回路推到单位圆外稳定裕度极薄第二被控对象 P(z) 本身有相位滞后内模正反馈叠加时相位对不齐谐振峰位置会偏移实际压制效果远不如理论。3.2 改进点一把常数 Q 嵌进内模压低峰高换鲁棒性改进后的内模写成 M(z) 1 / (1 - Q(z)·z⁻ᴺ)Q(z) 可以是一个略小于 1 的常数也可以是一阶或二阶低通。常数 Q 的物理含义是每周期只保留上一周期误差的 Q 倍等效于给内模装了个漏泄通道谐振峰高度从无穷大降到 1/(1-Q)。Q 取 0.95 时峰值约 20 倍压住 5 次谐波绰绰有余同时对被控对象相位误差的容忍度明显放宽。Q 的取值要配合 Kr 一起考虑。Q 越接近 1稳态精度越高但收敛越慢、鲁棒性越差Q 降到 0.9 以下收敛快但低次谐波压制能力下降。工程上常取 0.92 到 0.97若被控对象谐振峰尖锐或电感参数变化范围大就把 Q 再往 0.9 靠。另一个常见做法是把 Q 做成低通滤波器让低频段的 Q 接近 1、高频段的 Q 明显小于 1这样既保住了低次谐波增益又抑制了高频段的稳定风险。3.3 改进点二相位超前 k 与 S(z) 补偿器的设计流程补偿器 C(z) Kr·z^k·S(z) 由三部分组成增益 Kr 决定总压制深度超前环节 z^k 补偿被控对象的相位滞后S(z) 通常是二阶低通把高频段的回路增益压到 1 以下。三者里最容易搞错的是 k 的取值k 不够相位对不齐、谐振峰偏移k 过大相当于引入相位过冲同样破坏稳定。k 的计算思路是先求出被控对象 P(z) 在各次谐波频率处的相位滞后 φ负值再折算成采样拍数。一个采样拍对应的基波相位是 360°·f0/fs10 kHz 采样、50 Hz 基波时为 1.8°。# 被控对象单相并网逆变器 L 滤波P(s) 1/(LsR) L, R 3e-3, 0.1 bg, ag signal.bilinear([1.0], [L, R], fs) f_grid np.array([50, 150, 250, 350, 550, 650]) per_tap 360 * f0 / fs # 每拍对应 1.8° 基波相位 _, h signal.freqz(bg, ag, worN2 * np.pi * f_grid * Ts) for f, hh in zip(f_grid, h): ph np.degrees(np.angle(hh)) k int(round(-ph / per_tap)) # 折算成需要超前的拍数 print(f{f:4d} Hz 相位 {ph:7.2f}° 建议超前 {k} 拍)跑出来的相位滞后会随频率递增但超前拍数只能取一个统一值所以要按主要谐波能量分布来做取舍。光伏并网电流谐波里 5、7 次占大头我一般按 250 Hz 到 350 Hz 的相位滞后折算取一个折中的整数 k再回去跑稳定性判据确认。S(z) 用二阶巴特沃斯低通截止频率取 800 Hz 到 1.2 kHzbs, as_ signal.butter(2, 2*fc/fs)即可它的作用是在相位补偿失效的高频段把回路增益拉下来。3.4 稳定性校验从误差收敛条件看参数边界重复控制系统稳定的充分条件是对所有频率满足| Q(e^{jωT}) · (1 - C(e^{jωT})·P(e^{jωT})) | 1这个判据的含义是内模增益 Q 乘以1 减去补偿器与被控对象乘积的模值必须小于 1。C·P 在低频段接近 1括号内趋近 0稳定条件自然满足在高频段 C 被 S(z) 压低乘积趋近 0括号内趋近 1此时稳定条件退化为 Q 1这正是 Q 存在的意义。Q 0.95 kr 0.85 fc 1000.0 bs, as_ signal.butter(2, 2 * fc / fs) f np.linspace(1, fs / 2, 4000) w 2 * np.pi * f * Ts _, P signal.freqz(bg, ag, worNw) _, S signal.freqz(bs, as_, worNw) lhs np.abs(Q * (1 - kr * S * P)) print(稳定判据最大值 , lhs.max(), 对应频率 , f[lhs.argmax()])判据最大值必须严格小于 1留 0.1 左右的裕量更安全。如果最大值出现在某个谐振点附近说明 S(z) 的截止频率给高了或者 Q 给大了把 fc 往下调一档再算。这套校验每次改参数都跑一遍比在时域波形上反复试探效率高得多。4. 准PR并联改进重复控制的仿真构建与THD对照4.1 并联拓扑与两路输出的叠加方式复合控制策略的典型结构是把准 PR 与重复控制并联在同一个误差信号上指令电流 i_ref 减去实际并网电流 i 得到误差 ee 同时进入准 PR 支路和重复控制支路两路输出相加后作为调制电压指令。准 PR 支路提供基波高增益和快速动态重复控制支路提供全部低次谐波的高增益。这里有个容易被忽略的细节重复控制内模在基波处同样有高增益与准 PR 的谐振峰叠加后基波环路增益可能过大导致相位裕度不够。常见做法是在重复控制支路的补偿器里串一个 50 Hz 陷波器把基波增益让给准 PR或者在调试时把 Kr 适当调小让两路叠加后的基波总增益回到合理范围。我一般用前一种结构清楚参数之间耦合小。4.2 单相并网逆变器全闭环仿真代码下面这段代码把准 PR、重复控制内模、补偿器和 L 滤波被控对象串成完整闭环用逐步迭代的方式跑方便观察每拍的状态变化。电网电压里注入了 3、5、7 次谐波用来模拟背景谐波偏大的工况。# ---------- 重复控制参数 ---------- Q 0.95 kr 0.85 k_lead 4 # 由 3.3 节的相位折算得到 delay N - k_lead # 超前环节并入内模延迟 bs, as_ signal.butter(2, 2 * 1000.0 / fs) # S(z)1kHz 二阶低通 use_rc True # False 时退化为纯准PR用于对照 # ---------- 仿真时间轴与激励 ---------- Tsim 0.6 n int(Tsim / Ts) t np.arange(n) * Ts iref 20 * np.sin(W0 * t) # 20A 峰值指令 ug (311 * np.sin(W0 * t) 12 * np.sin(3 * W0 * t 0.4) 8 * np.sin(5 * W0 * t - 0.7) 6 * np.sin(7 * W0 * t)) # ---------- 状态量 ---------- i np.zeros(n); e np.zeros(n); v np.zeros(n) u_rc np.zeros(n); u np.zeros(n) e1 e2 yp1 yp2 0.0 # 准PR 历史项 x1 x2 ys1 ys2 0.0 # S(z) 历史项 for k in range(n - 1): e[k] iref[k] - i[k] # 准PR 差分方程 ypr (b_pr[0] * e[k] b_pr[1] * e1 b_pr[2] * e2 - a_pr[1] * yp1 - a_pr[2] * yp2) e2, e1 e1, e[k] yp2, yp1 yp1, ypr # 重复控制内模v[k] e[k] Q·v[k-delay] v[k] e[k] (Q * v[k - delay] if k - delay 0 else 0.0) # 补偿器 S(z) 与增益 yr (bs[0] * v[k] bs[1] * x1 bs[2] * x2 - as_[1] * ys1 - as_[2] * ys2) x2, x1 x1, v[k] ys2, ys1 ys1, yr u_rc[k] kr * yr if use_rc else 0.0 # 控制量叠加 被控对象前向欧拉离散 u[k] ypr u_rc[k] i[k 1] i[k] Ts / L * (u[k] - ug[k] - R * i[k])代码里有几处必须留意的实现细节。delay N - k_lead把补偿器里的超前环节折进内模延迟避免在差分方程里引用未来采样点这是数字实现的标准处理。v[k]在k - delay 0时只能取 0相当于内模还没填满第一个周期这段过渡属于正常现象。被控对象用前向欧拉离散10 kHz 采样对 3 mH 电感足够若把采样率降到 5 kHz 以下需要换成双线性或零阶保持离散否则相位误差会吃掉稳定性裕量。4.3 THD 与动态响应对照结果把use_rc分别置为 False 和 True 各跑一次取 0.40 s 到 0.50 s 的稳态段做 FFT按 2 到 40 次谐波累加计算 THD。典型结果差异很明显控制方案并网电流 THD5 次谐波含量7 次谐波含量指令跳变恢复时间纯准 PR3.8% ~ 4.5%2.6%1.9%约 4 ms准 PR 改进重复控制0.6% ~ 1.1%0.3%0.2%约 6 ms纯重复控制1.0% ~ 1.5%0.4%0.3%约 20 ms可以看到复合方案把 THD 压到 1% 附近动态恢复时间相比纯重复控制缩短了三分之二代价是比纯准 PR 慢 2 ms 左右。这 2 ms 换来的是 THD 下降三倍多在并网导则卡 3% 的场景下完全值得。如果动态指标更紧可以把 Kr 从 0.85 降到 0.6THD 会升到 1.5% 左右恢复时间降到 5 ms。5. 复合控制落地时的数字实现细节与频域验证5.1 采样同步与数字延迟补偿重复控制的相位补偿对延迟极其敏感而实际 DSP 里还有一拍计算延迟、PWM 更新延迟和采样滤波延迟。我习惯在离线整定完 k 之后把实际存在的延迟拍数从 k 里扣掉剩下的才是真正的超前量。采样时刻尽量对齐三角载波的顶点或谷底让开关纹波电流在采样点附近处于平均值位置否则纹波会作为噪声进到误差缓冲里被内模逐周期放大表现为并网电流上叠一层低频包络。如果用的是单采样单更新模式采样到更新之间隔了一拍这拍延迟对 5、7 次谐波来说相当于几度的相位误差会直接影响补偿效果。可以改成双采样双更新把延迟压到半拍或者简单点在计算 k 时直接按实际延迟拍数 计算延迟 1来折算。5.2 光照突变与模式切换时的抗饱和处理光伏逆变器在光照阶跃时的指令电流跳变幅度可能达到额定值的 30% 以上此时准 PR 支路会瞬间输出很大的调制指令重复控制内模里也会灌进一大段突变误差并且这个误差会在后续 N 拍里反复循环输出纹波。常见做法是给误差缓冲区加一个软复位逻辑当 |e| 超过阈值的持续时间超过若干拍时把内模状态 v 清零重新填充等静差回到正常范围再恢复。这样做会牺牲一个周期的谐波压制但换来了跳变期间不出现过冲和振荡。另一个细节是输出限幅要放在两路叠加之后而不是分别限幅再相加否则两路各自都不超限、相加后仍然越界会引入非线性失真。5.3 用 FFT 与频响曲线做上线前验证上线前建议做两项验证。第一项是稳态频谱验证代码不长seg i[int(0.40 * fs):int(0.50 * fs)] # 取 10 个基波周期稳态段 win np.hanning(len(seg)) # 汉宁窗抑制频谱泄漏 X np.fft.rfft(seg * win) freq np.fft.rfftfreq(len(seg), Ts) amp 2 * np.abs(X) / win.sum() # 单边幅值修正 I1 amp[np.argmin(np.abs(freq - f0))] harm np.sqrt(sum(amp[np.argmin(np.abs(freq - h * f0))] ** 2 for h in range(2, 41))) print(基波幅值 %.2f A, THD %.2f%% % (I1, harm / I1 * 100))取 10 个基波周期配合汉宁窗能把频谱泄漏压到可接受范围如果分段长度不是基波周期的整数倍泄漏会抬高旁瓣THD 算出来偏高容易误判。第二项是把第 3.4 节的稳定判据代码在标称电感、±20% 电感偏差和电网阻抗变化几个工况下各跑一遍确认判据最大值始终小于 0.9。电感饱和时的实际值可能只有标称的 60%这一档不能省。本文还有配套的精品资源点击获取