最近在做LCL三相并网逆变器的Simulink仿真控制方案走的是dq变换锁相环PLL全状态反馈LQ线性二次型控制这条路。做这个项目的初衷很直接LCL滤波器在并网逆变器里有天然的谐振尖峰传统PI双环不加有源阻尼很难稳得住而状态空间方法可以把LCL的完整模型装进控制器把谐振阻尼和电流跟踪一起解决。整套仿真搭建下来踩了不少坑也把参数整定的直觉摸清楚了。这篇就做个完整记录希望给正在做并网逆变器控制仿真、或者准备把LQR从课本挪到工程里的同学一些可以照抄的细节。1. 项目整体思路为什么LCLLQR是个好组合LCL三相并网逆变器是光伏、储能、充电桩这些并网设备里最核心的一环。这几年做相关方向的工程师和研究生越来越多但大多数人习惯一上来就套PI双闭环。真把LCL滤波器摆上去才会发现问题单靠网侧电流反馈系统的稳定裕度被谐振尖峰吃掉波形很容易在某个频率上抖起来。LCL滤波器本身是三阶系统两个电感加一个电容好处是高频谐波衰减能力比单L强很多坏处是幅频特性在谐振频率附近冒出一个很高的尖峰。这个尖峰如果不处理任何一点谐波激励或者控制扰动都会把它点着。传统做法是在电容支路串无源阻尼电阻简单有效但损耗大大功率场合不划算现代做法是有源阻尼也就是靠控制算法把谐振峰“摁”住。我选LQR全状态反馈核心原因是它把LCL在dq坐标系下的六阶模型一次性放进状态方程里求解所有内部状态逆变器侧电流、电容电压、网侧电流都参与反馈。谐振阻尼天然包含在最优控制律里不需要额外设计虚拟电阻回路。相比PI双环电容电流反馈的有源阻尼方案LQR的思路更“整体”代价是需要精确的模型矩阵和相对复杂的调参。这套方案适合谁一是刚接触并网逆变器的研究生需要一份能复现的仿真架构二是工程师想对比经典PI控制和现代状态空间方法在并网场景下的差异。后面的内容会从原理讲到Simulink搭建细节参数、矩阵、模块连线都尽量给出可直接抄的值。1.1 先拆一盘棋dq变换、锁相环、全状态反馈各自做什么这套系统拆成四个部分分开看都不难难在怎么拼起来。锁相环管的是角度。并网逆变器必须跟电网相位同步不能自己随便生成正弦波输出必须实时跟踪电网电压矢量的相位。我用的SRF-PLL同步旋转坐标系锁相环原理很简单把三相电压变换到dq坐标系让q轴电压趋近于0这时d轴跟电网电压矢量重合这个角度就是控制要用的θ。dq变换管的是坐标。三相交流量在静止abc坐标系里是正弦量控制器追正弦很别扭。但转进同步旋转坐标系后50Hz的正弦变成直流量可以像控制直流电机一样做闭环。前提是Park变换的θ必须来自PLL角度不对dq轴解出来的电流就不是直流量控制器看到的全是“交流干扰”。全状态反馈管的是所有变量的闭环。我不只看网侧电流还把逆变器侧电流和电容电压一起拉进反馈LQR帮我把这些状态量乘上一个反馈增益矩阵K算出逆变器应输出的电压指令。这样能一次性把LCL的谐振模态压制住比只反馈网侧电流的单环控制器稳得多。LQ线性二次型管的是反馈增益怎么算。状态反馈的关键问题是K怎么取。LQR的思路很直接定义代价函数J∫(xᵀQxuᵀRu)dtQ加权状态误差R加权控制动作求解Riccati方程得到最优K。在MATLAB里就是一行Klqr(A,B,Q,R)但Q、R怎么选需要物理直觉这是后文第3节的重点。1.2 方案选型背后的几个判断有人会问PI双闭环成熟又简单为什么折腾LQR我用一个真实感受回答PI双闭环在LCL上不是不能做只是需要额外的有源阻尼策略。常见做法是电容电流反馈虚拟电阻法本质上也是状态反馈只是只反馈了部分状态。LQR则是系统的完整状态反馈不需要单独设计阻尼策略代价是状态矩阵的建模和调试难度上来了。从控制理论角度看LQR本质是线性状态调节器直接作用于系统内部动态特性。LCL是典型的欠阻尼系统谐振峰在频域看就是一个明显的模态LQR能把这个模态的极点直接配置到期望位置这是它天然的优势。另一个原因是我个人的验证需求想知道基于模型的控制方法在并网逆变器上能做到什么程度尤其是有功/无功阶跃时动态响应怎么变化。设计目标也比较明确并网电流THD小于5%有功功率10kW无功0kVar动态响应时间控制在50ms以内超调尽量不超过20%。后面所有参数整定都围绕这几个指标展开。1.3 硬件参数与仿真案例设定这里给出整套仿真的参数基准后面所有矩阵计算和模块搭建都基于这份设定。电网三相380V/50Hz相电压有效值220V峰值约311V直流侧电压750V逆变器侧电感L13mH网侧电感L21mH滤波电容C10uF。两个电感的比例31网侧电感取小一点能减少电网阻抗对控制的影响。电容选10uF算出来的谐振频率约1.84kHz在开关频率10kHz的五分之一附近滤波和动态响应的权衡比较合理。开关频率10kHz控制采样周期Ts1/fsw1e-4s。这个参数组合对应额定有功10kW下的网侧电流有效值约15.2Ad轴电流峰值指令大约在21A量级。这个参数也决定了我后面LQR权重里电流项的初始量级物理量级匹配是调Q矩阵的第一步。2. 从模型到算法把核心环节在Simulink里逐一落地2.1 锁相环SRF-PLL的搭建与调试锁相环在Simulink里有两个选择用Simscape Electrical的现成PLL模块或者自己搭。我的建议是至少自己搭一次因为能真正理解内部机制也能处理现成模块的一些坑比如初始角度、仿真步长引起的角度抖动。SRF-PLL闭环结构三相电压va、vb、vc经过Clarke变换到αβ再经过Park变换依赖当前θ得到vd、vq把vq作为误差信号送进PI调节器PI输出叠加到额定角频率ω_nom2π×50上积分后得到θθ再反馈回Park变换。当vq0说明θ超前了PI输出降低Δω让积分慢一点vq0则相反。收敛后vq稳定在0θ就锁在了电网电压矢量方向上。Simulink实现细节有几个容易翻车的地方。Clarke变换和Park变换可以用Fcn模块或者Matlab Function写推荐后者的状态明确方便调试。PI调节器的输出需要限幅一般限制频率偏移在±5Hz以内防止启动阶段角度猛冲。积分器的初始角度可以设0仿真开始后PLL会自动拉着θ往电网实际相位靠收敛时间一般在几十毫秒以内这个期间电流指令先保持0等PLL锁定后再闭合电流环。有个细节是关于θ的归一化。积分器内部的角度是持续累加的仿真时间长了数值可能变得很大影响后续cos/sin的精度。我习惯在送给PWM和dq变换模块之前做mod(θ,2π)归一化但积分器内部保持累计值这样角度不会跳变。如果用Simulink的wrap-to-2pi模块注意它会让角度在0到2π之间跳变但cos/sin是周期函数跳变不会引起波形阶跃只有你把这个跳变直接当“相位信号”去做比较时才会出问题。实测下来连续累计角度输出端归一化是最稳的组合。锁相环的PI参数可以按经典方法设计。开环传递函数近似为Kp×(sKi/Kp)/s × 1/s × Vg其中Vg是电网电压d轴幅值约311V。我按自然频率100Hz、阻尼比0.707设计得到Kp在0.5到2之间、Ki在50到200这个量级。带宽不用太高太高容易被电网谐波干扰。实测中先固定Ki再调Kp看dq电压波形收敛速度和vq抖动幅度抓到手感后就固定下来不用频繁动。2.2 dq变换与坐标对齐的细节dq变换的坐标对齐是整套控制的地基。我用等幅值Clarke变换d轴跟电网电压矢量对齐。这样d轴电压分量等于电网相电压峰值311Vq轴电压分量理论为0控制上直白好观察。具体变换公式我贴一下方便直接抄进Matlab Function。Clarke变换等幅值写法vα (2/3)×(va - 0.5×vb - 0.5×vc)vβ (2/3)×(0.866×vb - 0.866×vc)Park变换我用的标准形式vd vα×cosθ vβ×sinθvq -vα×sinθ vβ×cosθ这里符号约定一定要和PLL闭环极性对上。有的教材用另一套符号其实都能锁住只要PI输出的正负极性跟vq的误差方向匹配。我在调试时走过弯路PLL锁不住频率一直飘后来发现是Park变换里sin项符号反了导致vq误差方向反了PI把一个正反馈闭环包了进去。这个如果不去追溯原理光调PI参数死活调不出来。三相电流的dq变换跟电压完全一样可以用同一个变换模块复制一份输入换成三相电流。所有状态量转成dq后在旋转坐标系下都是直流量控制带宽的概念才成立。还有个容易忽略的点是dq变换输出的d、q分量是峰值表达不是有效值表达。比如额定10kW工况下等幅值变换出来的i2d峰值约21A有效值约15.2A。后续做功率计算、指令换算时千万别拿峰值去套有效值公式单位搞混了LQR参考值会差√3倍波形对不上。2.3 状态空间模型与LQR增益计算LQR是模型驱动的控制方法状态方程写错一个符号后面全白搭。这里给出LCL在dq旋转坐标系下的完整状态方程。状态向量取x[i1d, i1q, v_cd, v_cq, i2d, i2q]ᵀ其中i1是逆变器侧电流v_c是滤波电容电压i2是网侧电流。控制输入u[v_id, v_iq]ᵀ是逆变器桥臂中点电压相电压的dq分量忽略PWM谐波后的平均值。扰动项是电网电压的dq分量w[v_gd, v_gq]ᵀ以电网电压定向时v_gd311Vv_gq0。d轴状态方程L1 di1d/dt v_id - v_cd ωL1 i1qL2 di2d/dt v_cd - v_gd ωL2 i2qC dv_cd/dt i1d - i2d ωC v_cqq轴状态方程L1 di1q/dt v_iq - v_cq - ωL1 i1dL2 di2q/dt v_cq - v_gq - ωL2 i2dC dv_cq/dt i1q - i2q - ωC v_cd注意这里的交叉耦合项ωLi和ωCv是因为在同步旋转坐标系下观察三相交流量天然会看到对方轴的贡献。很多人在建模时会把这几个项漏掉导致LQR模型与实际系统不匹配控制器在仿真里跑出来根本不对。耦合项符号也很重要d轴方程里是正ωq轴方程里是负ω这个不能凭感觉写要从Park变换的微分推导出来。代入参数L13e-3、L21e-3、C10e-6、ω2π×50得到A矩阵A [0, 314.16, -333.33, 0, 0, 0; -314.16, 0, 0, -333.33, 0, 0; 100000, 0, 0, 314.16, -100000, 0; 0, 100000, -314.16, 0, 0, -100000; 0, 0, 1000, 0, 0, 314.16; 0, 0, 0, 1000, -314.16, 0]B矩阵B [333.33, 0; 0, 333.33; 0, 0; 0, 0; 0, 0; 0, 0]输出矩阵C选网侧电流dq分量C [0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 1]有了这个连续状态空间模型接下来做离散化和LQR求解。在MATLAB脚本里按下面的方式操作% 参数定义 L1 3e-3; L2 1e-3; C 10e-6; fg 50; fsw 10e3; Ts 1/fsw; w 2*pi*fg; % 连续状态矩阵 A [0, w, -1/L1, 0, 0, 0; -w, 0, 0, -1/L1, 0, 0; 1/C, 0, 0, w, -1/C, 0; 0, 1/C, -w, 0, 0, -1/C; 0, 0, 1/L2, 0, 0, w; 0, 0, 0, 1/L2, -w, 0]; B [1/L1, 0; 0, 1/L1; 0, 0; 0, 0; 0, 0; 0, 0]; % 离散化控制周期Ts sys_c ss(A, B, eye(6), zeros(6,2)); sys_d c2d(sys_c, Ts, zoh); Ad sys_d.A; Bd sys_d.B; % 输出矩阵网侧电流d/q Cd [0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 1]; % 增广积分状态消除稳态误差 A_aug [Ad, zeros(6,2); -Cd, eye(2)]; B_aug [Bd; zeros(2,2)]; % LQR权重矩阵 Q diag([10, 10, 50, 50, 1000, 1000, 3000, 3000]); R diag([0.001, 0.001]); % 求解离散LQR增益 K dlqr(A_aug, B_aug, Q, R);这里解释几个设计决策。增广积分状态是为了消除稳态误差LQR本质是个状态调节器没有积分项时模型参数稍有偏差或者存在电网电压扰动网侧电流就会偏离指令。把参考电流误差的积分量作为两个新增状态放进去系统就从“调节器”升级成了“伺服跟踪器”阶跃参考信号下稳态误差理论上为零。Q矩阵的量纲是按物理量级定的。i2d和i2q是最终要严格控制的网侧电流权重最高1000。电容电压v_cd、v_cq是维持LCL滤波特性的中间变量权重50主要用来阻尼谐振。i1d、i1q是逆变器侧电流权重10它起到限制内环电流和间接抑制过流的作用。最后两个积分状态权重3000保证积分动作足够积极但又不会把系统推成振荡。R矩阵取0.001对应控制电压的代价因为逆变器输出电压最大也只有几百伏控制能量相对充裕R取小一些能让系统响应更快。如果想让响应更快可以增大Q整体数值比如把i2的权重提到3000到5000如果波形出现高频振铃说明状态权重相对R太大把R从0.001提到0.01或者降低Q即可。后面调试部分我会再细说。3. Simulink主电路与控制环路的完整搭法3.1 主电路搭建与仿真环境配置打开Simulink第一步把主电路搭出来。电路拓扑从直流侧到电网侧依序是直流电压源750V、Universal Bridge三相逆变桥、L1电感、C电容带阻尼电阻可选项、L2电感、三相电网。直流电源用DC Voltage Source就行电压参数750V。Universal Bridge选三相桥臂功率开关器件先用Ideal Semiconductor仿真速度快适合验证控制逻辑。等控制方案验证完毕再换成IGBT带导通损耗模型做更接近实测的仿真这是两个阶段的事千万别一上来就用详细器件的开关暂态模型仿真速度会拖到让你怀疑人生。LCL滤波器参数在桥臂输出和电网之间串联。L13mHC10uFL21mH。电容支路可以串一个0.1Ω到2Ω的小电阻看你要不要保留无源阻尼。由于LQR本身自带阻尼我在最终方案里把阻尼电阻去掉了这样能看出状态反馈真实压制谐振的能力。但是第一次调试时建议保留一个1Ω阻尼电阻万一控制器没调好它能拉一把不让波形彻底发散等LQR稳住后再逐步减小到0。电网模型用Three-Phase Source相电压220V RMS频率50Hz。在网侧L2和电网之间可以串一个很小的线路阻抗比如1mH0.1Ω模拟真实电网内阻。这个阻抗会影响PLL观测到的电压幅值和相位但控制策略本身有鲁棒性不会因为这种程度的变化而失稳。仿真参数设置上求解器我推荐ode23tb变步长。最大步长设1e-5秒或者更小否则PWM载波和开关纹波可能被仿真步长“漏过去”波形看起来毛刺很小但实际是失真。仿真时长先跑0.5秒足够看到从启动到稳态再到阶跃响应的全过程。如果只验证PLL可以短到0.2秒。测量模块用Three-Phase V-I Measurement分别测电网侧电压和网侧电流再用Current Measurement测逆变器侧电流、Voltage Measurement测电容电压。这些测量值输出的是三相瞬时波形统一送进三相测量总线然后分别接到dq变换模块。注意测量模块的质量可能会让信号总线比较乱建议用Goto/From标签管理信号避免一堆信号线在模型里交叉。PWM生成我用的是自搭三角波比较法不依赖库里的PWM Generator模块。三角波频率10kHz幅值从-1到1。控制侧经逆dq变换得到三相调制波指令也归一化到-1到1两者比较后输出0/1脉冲一对桥臂上下管互补。这种做法的好处是完全掌握PWM机制要加死区也方便。实际搭建时用三角波发生器Repeating Sequence频率设1e4比较器输出进Universal Bridge的桥臂控制端口。3.2 控制环路的信号流程与关键连接控制电路是整个仿真的核心信号流我按顺序理一遍。从测量总线出来的三相电流分别进三套abc_to_dq变换模块逆变器侧电流i1转到dq得到i1d、i1q电容电压vc转到dq得到v_cd、v_cq网侧电流i2转到dq得到i2d、i2q。三者用的θ都来自同一个PLL模块的输出。参考指令i2d_ref是有功电流指令。额定10kW功率下等幅值dq变换的d轴电流峰值约21A计算过程P3/2×vd×idvd≈311Vid≈2×10000/(3×311)≈21.4A。i2q_ref设0对应单位功率因数并网。调试时可以再加一个阶跃源0.1秒时把i2d_ref从20A切到30A观察动态响应。误差信号e_di2d_ref-i2d、e_qi2q_ref-i2q分别送入积分器得到两个积分状态e_intd、e_intq。积分器输出需要限幅我一般限制在±20左右防止启动瞬间误差大导致积分饱和控制器输出一直顶在限幅上好几毫秒不回来。状态反馈的接线顺序必须跟K矩阵的列序完全一致。K_aug的列序是[i1d, i1q, v_cd, v_cq, i2d, i2q, e_intd, e_intq]。把这8路信号按这个顺序排列然后进入一个Mux再接上一个Gain模块Gain的增益矩阵填K_aug。Gain模块输出的两路就是待补偿前的控制量u_d和u_q。最后一步是前馈和归一化。u_dq叠加电网电压前馈量v_gd这里约311V再把结果除以直流电压的一半Vdc/2375得到调制波幅值。为什么是Vdc/2因为三相桥双极性SPWM下相电压峰值可以用Vdc/2×调制比来线性近似反过来说明调制比m2×v_ref/Vdc。不除以这个系数调制波幅度可能超过1PWM会进入过调制区波形出现削顶畸变。调制波dq量经过逆Park变换到αβu_alphaud×cosθ-uq×sinθu_betaud×sinθuq×cosθ再从αβ重构三相va_refu_alphavb_ref-0.5×u_alpha0.866×u_betavc_ref-0.5×u_alpha-0.866×u_beta。这三路调制波输入PWM比较器。控制回路中还需要注意信号类型。Simulink里总线信号和向量信号混合使用时容易出维度错误我习惯把所有状态量用向量Mux汇总Gain模块直接处理向量。其他模块输出标量再用Mux合成。遇到维度不匹配的错误提示多半是哪个模块的Output端口维度和Mux预期不一致去Signal Attributes里检查尺寸就好。3.3 LQR参数整定的实操方法与收敛技巧Q、R的整定是一个“量级匹配→粗调→微调”的过程。我个人的经验是先用Bryson规则确定量级Q的对角元素i_i取1/(最大允许偏差的平方)R对角元素j_j取1/(最大控制量的平方)。比如允许网侧电流峰值偏差20A那Q[5,5]1/4000.0025按这个数量级再放大到1000因为直接按Bryson算出来的权重太小响应偏慢。最后我选用的Q矩阵是多次仿真后的结果既满足动态响应又留了稳定裕度。粗调阶段看阶跃波形。把i2d_ref从20A阶跃到30A网侧电流d轴应该在5到10毫秒内跟上且无振荡。如果电流爬得太慢说明Q相对R太小提高i2对应的权重或者降低R如果电流有高频振铃说明状态权重偏大或者积分权重太高降低积分权重换稳定性。微调阶段看解耦效果。设置iq_ref为一个阶跃比如从0阶跃到5A观察id轴是否被拉偏。如果id出现明显扰动说明dq交叉耦合补偿没有做好或者PLL角度偏差导致坐标系没有真正对齐。LQR是全状态反馈理论上耦合都在模型里处理了如果还能看到明显耦合优先去检查A矩阵里ωL和ωC耦合项的符号和量纲。还有一个关键技巧是逐步“加狠”LQR。第一次仿真先用保守的Q所有电流权重在100以下R在0.01先让系统稳定下来波形至少不发散。稳定后每轮更新一个参数先把i2的权重翻倍看响应是否变快再微调积分权重看稳态误差清零时间最后动R看能不能进一步压超调。这样单变量逐个试比一次性把所有权重调到最终值更容易定位问题。在Simulink里修改Q矩阵不用重建模型只要在MATLAB命令行窗口改Q、重算K、然后把Gain模块的矩阵参数替换掉再点仿真即可。这个循环做熟了一套参数整定半小时内就能完成。如果参数改完波形剧烈变化用“回到上一组参数”能帮你快速缩小问题范围。4. 调试中踩过的坑与排查实录4.1 PLL失锁与角度对齐问题PLL失锁是最先遇到的坑现象是dq电压波形在启动后持续震荡或者vq收敛不到0。我排查的顺序是三层。第一层先断开PLL的输出对控制的影响手动给θωt一个理想斜坡看dq变换出来是不是平滑的直流量。如果是说明dq变换模块本身没问题。第二层把PLL闭环接上但只接PLL不接电流环单独观察vq收敛速度和θ波形。第三层才接完整控制。PLL参数不对会让锁定过程有超调甚至失锁。Ki太小则锁定速度很慢50ms以上还没平稳Ki太大会在锁定后出现持续振荡。调试时先把Ki设小比如50确认能锁住再逐步增大到响应可以接受。如果电网电压里谐波分量大可以在vq误差信号后加一个截止频率500Hz的低通滤波滤掉谐波干扰。仿真里用理想电压源没这个问题但实机信号采集后必须考虑。角度对齐还有一个隐含的坑PLL锁定的是电网电压矢量角但dq变换的初始θ跟电网电压初始相角不一定同步。仿真从t0开始电网电压源在t0时a相电压为0默认但PLL内部积分器初始值可能设成0这时θ0对应的是a相电压峰值时刻并不是t0时的状态。PLL会自动追赶相位差表现为启动后θ有一个快速调整过程dq电压会出现短暂的低频摆动。要消除这个初始摆动可以把PLL积分器初值设为电网电压在t0时刻的相角比如设成0对应a相过零但更省事的做法是仿真前0.05秒不闭合电流环等PLL锁定后再给电流指令。4.2 LCL谐振引起的振荡与权重对策把LQR接上后最典型的坏现象是并网电流出现高频振荡频率跟谐振频率1.84kHz接近。这个问题的本质是谐振模态没有被充分阻尼。LQR反馈矩阵K包含了所有状态的信息理论上能压住谐振但如果Q矩阵里电容电压权重太小或者控制输出R权重太大LQR不舍得用力“踩”谐振模态振荡就浮出来。我的对策是分两步。第一步增大电容电压v_cd、v_cq在Q矩阵里的权重从50拉到200甚至500让控制器对电容电压的偏差更敏感反馈增益中对电容电压那一列的输出会更大等效于把谐振峰旁边的阻尼做得更厚。第二步是适当降低R比如从0.01降到0.001让控制器有更多“动作空间”去抑制振荡。这两个动作交替微调通常几轮就能把振荡压下去。还有一个容易被忽略的点PWM的离散化效应。桥臂每10kHz才更新一次输出电压这个零阶保持特性在谐振频率1.84kHz处会带来额外相位延迟。如果LQR是按连续系统设计然后直接丢进离散仿真这个延迟没被建模实际稳定裕度会低于设计值。更工程的做法是按第2.3节的方式先离散化再解dlqr把零阶保持器放进状态矩阵里。我一开始直接用连续lqr凑忽略了离散化结果同样的Q/R在连续模型里很稳接进Simulink离散PWM后反而振起来后来换成dlqr才解决。4.3 稳态误差、代数环与离散化细节稳态误差的根源通常是积分增广没加或者积分权重太小。如果i2d稳态值与指令差了0.5A以上先检查积分器的误差输入是不是真的接到了i2d_ref-i2d再检查积分增益是否够。积分权重3000在我的系统里稳态误差在一个仿真周期内就消到0.01A以下。但积分权重不能无限加大太大会导致参考阶跃时积分输出冲过头造成明显的超调和回复振荡。代数环是Simulink仿真里最容易卡住的报错。全状态反馈测量值经过Gain生成调制波调制波经过PWM影响主电路主电路再影响测量值这个环路在数值求解时如果没有一个采样周期彻底断开就会形成代数环。Simulink能解但是求解速度大幅下降严重时直接报错。解决方案是在反馈通道的测量信号上插入一个Memory模块或Unit Delay。我更推荐Unit Delay因为控制周期明确也与离散化的DSP实现一致。插入一拍的延迟意味着控制量落后真实状态一个采样周期对1kHz带宽的控制响应影响很小可以接受。离散化细节还包括积分器的选择。既然控制回路按离散化实现所有积分器、滤波器都换成Discrete-Time Integrator和Discrete Filter采样周期严格设为Ts1e-4。特别注意PLL里的积分器也要离散化否则混合连续/离散模型容易在变步长求解器里产生误差。我的PLL内部PI、积分器以及电流环的积分增广全部用离散模块这样整个控制侧是离散系统主电路保持连续仿真速度和稳定性都比较好。状态反馈信号顺序错位也是个高发问题。有一次我把i2d、i2q接到反馈输入端时和e_intd、e_intq的顺序弄反K矩阵的乘数关系完全乱了波形剧烈振荡。排查办法是在Gain模块前接一个Display把8路信号逐一对照状态向量的顺序。别嫌麻烦这个步骤的纠错效率最高。最后一个小建议仿真波形用Scope同时观察i2d、i2q和调制波m_d、m_q。调制波如果长时间贴着±1限幅说明控制量饱和LQR的线性假设就不成立了这时系统会出现奇怪的次同步振荡。解决方法是降低电流指令幅度或者在Controller里加入抗饱和逻辑把积分器在限幅期间冻结更新。我在实际调试中最深的体会是LQR在Simulink里的搭建难点不在数学而在“离散化信号顺序数值归一化”这些工程细节。模型矩阵推导对了一遍就能稳定跑通但工程细节需要反反复复打磨。这套仿真跑稳定之后后续扩展方向也很自然把PLL换成带谐波抑制的增强型PLL在LQR优化里加入输出限幅的抗饱和设计或者把状态反馈中的测量模型加进低通滤波和延时环节逼近实机传感器特性。这些都是在当前框架上小步修改的事。
