仿真分析报告)
并联有源谐波滤波器AHF仿真分析报告模型ActiveHarmonicFilter.slxMATLAB R2023b 自带示例Three-Phase Active Harmonic Filter分析日期2026-10-06工具链MATLAB R2023bSimulink / Simscape Electrical 专用电力系统元件库 纯 MATLAB 绘图数据来源原模型实机运行导出12 路信号 × 125001 采样点步长 2 µs目录概述电路拓扑与参数Simulink 模型截图仿真设置与数据导出工况时序稳态波形分析谐波与 THD 分析AHF 动态响应关键参数对比仿真结论文件清单与复现步骤分析过程中发现的问题与处理1. 概述1.1 装置功能并联有源谐波滤波器Shunt Active Harmonic Filter, AHF接在负载与电网之间的公共连接点PCC上通过实时检测负载电流中的谐波分量由电压型 IGBT 逆变器向 PCC注入等值反向的谐波电流使电网侧只看到基波电流从而把非线性负载产生的谐波就地消纳。与无源滤波器LC 调谐支路相比有源方案的特点是补偿特性不依赖电网阻抗不存在与电网谐振的风险能同时补偿多次谐波与无功且补偿容量可随负载变化需要自换相变流器有开关损耗容量受器件限制。1.2 本示例的控制思想示例采用参考电流跟踪方案PLL锁定电网相位生成与电网同相、幅值等于负载电流有效值的正弦参考电流i_ref参考电流与实测电网电流之差i_ref − i_s即为需要补偿的谐波误差滞环hysteresis电流控制器根据误差正负直接翻转 6 个 IGBT 的门极信号驱动逆变器输出该误差电流注入 PCC。控制目标是让电网电流尽可能等于纯正弦的i_ref——即让电源看不到负载的谐波。2. 电路拓扑与参数2.1 主电路结构位置环节关键参数左侧三相交流电源4.16 kV实测线电压有效值 4154 V 60 Hz→线路电感Ls电源内阻抗→PCC公共连接点有源滤波器并联接入点→负载支路两台二极管整流器 直流电流源↓并联支路有源谐波滤波器AHF负载由两台相移 30° 的二极管整流器构成Δ/Δ 接法整流器直流电流 2000 A → 3000 A 阶跃Δ/Y 接法整流器直流电流 3000 A在 t 10/60 s 经断路器投入两台整流器叠加即构成等效 12 脉波整流其 5 次、7 次谐波因相移而相互抵消。2.2 AHF 内部构成部件参数说明交流侧滤波电感Laf1 mH抑制开关纹波、决定电流跟踪速率直流母线电容C1、C2各 500 µF串联支撑直流母线IGBT 逆变桥两电平六脉冲由滞环控制器直接给出门极信号控制基准max_current4000 A决定误差信号的增益Gain1 1/max_current参考频率freq60 Hz供 PLL 与 RMS 环节使用数据文件ActiveHarmonicFilter.mat中另有AF.activefilter.voltage 440。经全模型参数引用扫描确认该字段未被任何模块引用属设计数据中的参考值不影响仿真。2.3 控制器内部结构三相误差分别经 3 个减法环节得到i_ref − i_meas再经 6 个 Relay阈值 ±eps等效过零比较生成 6 路门极信号g直接驱动逆变桥。2.4 测量与记录环节Measurements子系统汇集 6 路测量Vabc_source、Iabc_source、Vabc_AF、Iabc_AF、Vabc_load、Iabc_loadAHF 子系统内部另有Vdc、Idc共 8 路进入responses示波器。3. Simulink 模型截图原模型在 MATLAB R2023b 中打开后的实际界面print -s模型名渲染非示意图。3.1 模型总图图中可见三相电源 → 线路电感Ls→Shunt active harmonic filter青色AHF 本体→Load Bus→上下两路负载上Δ/Δ 整流器 负载输出Vdc1下经三相断路器的 Δ/Y 整流器 负载输出Vdc2。左下角为powergui显示Discrete 2e-06 s。3.2 powergui 求解器设置块面显示Discrete 2e-06 s.即电气网络的离散仿真步长为 2 µs与模型全局求解器设置变步长ode45 powergui 离散化一致。4. 仿真设置与数据导出4.1 求解器设置项目取值求解器ode45变步长powergui 离散化步长2 µsDiscrete 2e-06 s仿真时长15/60 0.25 s采样点数125001每路信号记录方式逐示波器挂接To WorkspaceStructure With Time4.2 数据导出方法原模型由powergui封装sim()在常规配置下只返回时间轴不产生任何工作区数据变量。本报告采用的方法是对 6 个示波器的每一个输入端口从其源端口引出一条支线接到新建的To Workspace模块共 14 路信号responses示波器有 8 个端口是信息量最大的一路。导出脚本ahf_export2.m。导出的信号清单变量通道含义ws_Vs_abc3电源相电压ws_responses_p13Vabc_source电源线电压ws_responses_p23Iabc_source电源电流补偿后ws_responses_p33Vabc_AF滤波器接入点电压ws_responses_p43Iabc_AF滤波器注入电流ws_responses_p53Vabc_load负载电压ws_responses_p63Iabc_load负载电流未补偿ws_responses_p71Vdc滤波器直流母线电压ws_responses_p81Idc滤波器直流母线电流ws_Iref_Isource_p12i_ref与i_sourcews_Vdc1/ws_Vdc21 / 1两路负载直流母线电压5. 工况时序时刻事件负载形态负载电流峰值t 0投入 Δ/Δ 整流器直流电流 2000 A6 脉波±2.2 kAt 5/60 83.3 ms直流电流阶跃 2000 A → 3000 A6 脉波±3.3 kAt 10/60 166.7 ms断路器闭合Δ/Y 整流器投入12 脉波±6.8 kA从全时段波形可以清楚看到三种工况下源电流蓝始终紧随负载电流橙的基波包络负载电流中的换流缺口6 脉波特征在源侧完全消失。6. 稳态波形分析6.1 补偿效果核心单个工频周期第 3 周期6 脉波 2000 A负载电流i_load橙平顶 ±2.2 kA每个周期两次换流缺口——典型的 6 脉波二极管整流电流源电流i_source蓝接近理想正弦换流缺口完全被补掉。6.2 三路电流关系滤波器注入电流i_AF绿恰好等于i_load与i_source之差且主要包含高次谐波分量其有效值远小于负载电流见下表这也是有源滤波方案容量利用率高的原因。工况负载电流有效值源电流有效值注入电流有效值注入容量占比6 脉波 2000 A1600 A1599 A351 A22.0 %6 脉波 3000 A2381 A2393 A480 A20.2 %12 脉波 3000 A4670 A4702 A223 A4.8 %12 脉波工况下负载自身的谐波已很低因此需要滤波器注入的电流显著减小——这正是无源多脉波整流 有源滤波组合的价值用便宜的变压器绕组抵消掉最主要的 5、7 次谐波剩余的 11、13 次再由 AHF 处理AHF 容量可以做得更小。6.3 滤波器直流母线窗口母线电压均值纹波峰峰值纹波率第 3 周期充电暂态12784 V5234 V40.9 %第 8 周期6 脉波 3000 A7632 V1082 V14.2 %第 13 周期12 脉波稳态7253 V621 V8.6 %母线电压在 t 27 ms 处出现16.4 kV 的充电过冲——模型未设直流电压闭环调节母线由 IGBT 桥的反并二极管自然整流充电能量平衡后稳定在约 7.25 kV对于 4154 V 线电压的系统峰值线电压为 5874 V母线必须高于该值才能注入电流7.25 kV 是合理的。6.4 三相电压与电流三相电压为对称正弦相电压有效值 2398 V线电压 4154 V三相源电流在运行到 12 脉波工况后幅值增大但波形依然保持正弦。7. 谐波与 THD 分析7.1 分析方法采用官方示例的 FFT 窗口口径一个基波周期60 HzW1 第 3 周期[2/60, 3/60]——6 脉波 / 2000 AW2 第 8 周期[7/60, 8/60]——6 脉波 / 3000 AW3 第 13 周期[12/60, 13/60]——12 脉波 / 3000 A各次谐波幅值由窗口内的直接 DFT 积分求得H(k·f₀) (2/T)·∫y(t)·e^(−j2πkf₀t)dtTHD √(ΣHₖ²)/H₁取至 50 次。7.2 频谱6 脉波工况负载电流呈现教科书式的 6 脉波谐波谱——只含6k±1 次5、7、11、13…幅值依次按 1/h 衰减谐波次数5 次7 次11 次13 次17 次19 次负载电流17.94 %11.52 %5.28 %3.57 %1.86 %1.37 %源电流0.03 %0.09 %0.12 %0.03 %——所有谐波在源侧被压到0.2 % 以下。7.3 频谱12 脉波工况Δ/Δ 与 Δ/Y 两台整流器相移 30° 后5 次与 7 次谐波相互抵消谐波谱跳变为12k±1 次谐波次数5 次7 次11 次13 次23 次25 次负载电流0.01 %0.00 %3.96 %2.43 %0.95 %0.81 %源电流0.17 %0.07 %0.40 %0.21 %——这正是多脉波整流的价值把能量从最低次最难滤除、幅值最大的 5、7 次挪到更高次使滤波任务大幅减轻。7.4 谐波结构对比7.5 THD 汇总与官方数值对照工况窗口负载电流 THD源电流 THD降幅6 脉波 2000 A第 3 周期22.41 %0.35 %98.5 %6 脉波 3000 A第 8 周期20.56 %0.94 %95.4 %12 脉波 3000 A第 13 周期4.80 %1.13 %76.6 %与官方文档给出的数值对照来源第 3 周期第 13 周期官方文档22.41 % → 0.69 %4.80 % → 0.93 %本报告实测22.41 % → 0.35 %4.80 % → 1.13 %负载侧 THD 完全吻合22.41 %、4.80 %误差 0.01 %源侧数值同量级差异来自 FFT 实现细节官方用 powergui 的 FFT 工具带插值本报告用原始采样直接积分。这说明分析口径与数据链路是正确的。8. AHF 动态响应上图给出全时段的i_ref与i_source对比负载从 2000 A 阶跃到 3000 At 83.3 ms、以及 12 脉波整流器投入t 166.7 ms两个扰动下实际源电流均在约 12 个工频周期内重新贴合参考电流与官方描述captures the new reference current within one cycle一致。下表为各参数工况下测得的扰动恢复时间以工频周期计判据为电流残差 RMS 降至初始值的 10 %参数工况恢复时间Laf 0.5 / 1 / 2 mH1.75 周期Laf 4 / 8 mH 3 周期未恢复C_dc 1001000 µF1.75 周期C_dc 2000 µF2.25 周期max_current 4000 A设计值1.75 周期max_current≤ 2000 A未恢复失稳9. 关键参数对比仿真对 5 个参数做了单参数扫描每次只改一个其余保持设计值共 23 次完整仿真。评价指标取三个官方窗口的源电流 THD、直流母线电压及其纹波、滤波器注入电流有效值、扰动恢复时间。9.1 灵敏度排序各指标在扫描范围内的最大相对变化以中位数为基准参数源 THD第3周期源 THD第13周期注入电流I_af母线纹波ΔV_dc控制器基准I_max3739 %260 %224 %159 %交流电感L_af2209 %171 %34 %77 %12 脉波投入时刻5067 %※16 %122 %123 %直流电容C_dc21 %18 %3 %328 %负载电流I_dc0 %26 %56 %65 %※ 该项数值偏大是因为该参数取2/60 s时第 3 周期窗口恰好落在投入暂态内属过渡过程而非稳态差异不能与其他项直接比较。结论本装置最关键的两个参数是控制器增益基准I_max与交流滤波电感L_af——它们直接决定电流跟踪的精度与稳定性直流电容只影响直流侧纹波对补偿效果几乎无影响。9.2 交流滤波电感L_af最关键设计参数L_af源 THD 第3周期源 THD 第8周期源 THD 第13周期母线电压注入电流0.5 mH0.50 %0.16 %0.95 %6945 V233 A1 mH设计值0.35 %0.94 %1.13 %7253 V223 A2 mH0.46 %3.46 %1.57 %8261 V212 A4 mH5.08 %7.46 %2.41 %9875 V190 A8 mH11.36 %11.59 %3.62 %11611 V159 A机理电感决定电流变化率di/dt (v_inv − v_pcc)/L。电感越大逆变器能提供的电流斜率越小滞环控制器来不及跟上参考误差 → 跟踪失败 → 谐波残余在源侧出现。同时跟踪能力下降使注入电流减小能量在直流母线累积母线电压被动抬升6.9 kV → 11.6 kV。设计建议L_af应取在刚好满足最大di/dt需求的下限附近。本模型 1 mH 为最优增大到 4 mH 以上补偿效果急剧恶化THD 上升一个数量级。9.3 控制器基准I_max稳定性最敏感Gain1 1/max_current因此该值越小则环路增益越大。I_max环路增益比源 THD 第3周期源 THD 第13周期注入电流母线纹波状态500 A8×45.85 %33.53 %2066 A16387 V失稳1000 A4×5.15 %12.53 %1023 A11324 V振荡2000 A2×0.50 %12.70 %2500 A9989 V过补偿4000 A设计值1×0.35 %1.13 %223 A621 V正常8000 A0.5×1.22 %0.98 %213 A524 V增益偏低机理该增益把误差信号放大后作为参考电流增益过大时注入电流远超实际所需2000 A 时注入 2500 A超过负载谐波电流本身形成过补偿与振荡增益过小时8000 A响应变慢暂态窗口 THD 上升。设计建议设计值 4000 A 恰在最佳点。环路增益放大 2 倍以上即进入失稳区该参数没有向更小方向的安全裕度整定时须特别谨慎。9.4 直流母线电容C_dcC_dc母线纹波峰峰值纹波率源 THD 第3周期源 THD 第13周期100 µF2230 V31.8 %0.33 %1.06 %250 µF1057 V14.8 %0.36 %1.10 %500 µF设计值621 V8.6 %0.35 %1.13 %1000 µF284 V3.9 %0.32 %1.14 %2000 µF195 V2.8 %0.29 %1.26 %机理直流纹波与电容近似成反比100 µF → 2000 µF纹波 2230 V → 195 V降 11.4 倍而容量增 20 倍因纹波同时受直流电压水平影响。交流侧补偿效果几乎不受电容影响THD 在 0.29 %–0.36 % 之间这是直流侧储能与交流侧补偿解耦的体现。设计建议若只关注交流补偿质量电容可以取得较小取大容量的意义在于降低直流纹波、提高抗扰动能力2000 µF 时恢复时间反而从 1.75 升到 2.25 周期说明过大电容会让直流侧动态变慢。9.5 负载电流水平与平衡度改变 Δ/Δ 整流器的直流电流Δ/Y 侧固定 3000 A观察 12 脉波工况Δ/Δ 侧I_dc负载 THD第13周期源 THD第13周期注入电流1000 A12.36 %1.45 %380 A2000 A7.07 %1.24 %271 A3000 A平衡4.80 %1.13 %223 A4000 A5.22 %1.21 %292 A重要发现12 脉波的谐波抵消效果严格依赖两台整流器电流平衡。当 Δ/Δ 侧只有 1000 A 而 Δ/Y 侧 3000 A 时相移抵消不再完全负载 THD 从 4.80 % 恶化到 12.36 %5、7 次谐波重新出现。设计多脉波整流系统时两套绕组的负载电流必须匹配。值得注意的是即便负载 THD 恶化到 12.36 %源侧仍只有 1.45 % —— AHF 对不平衡工况同样有效。9.6 12 脉波整流器投入时刻断路器动作时刻负载 THD第13周期源 THD第13周期注入电流2/60 s4.80 %0.98 %213 A5/60 s4.80 %0.96 %214 A10/60 s设计值4.80 %1.13 %223 A20/60 s超出仿真时长等于不投入20.56 %1.07 %479 A关键结论当 12 脉波整流器始终不投入时负载电流 THD 高达20.56 %而源电流 THD 仍只有1.07 %—— 这说明AHF 单凭自身即可完成约 95 % 的谐波治理多脉波整流并非必需只是能显著降低 AHF 的补偿容量注入电流从 479 A 降到 223 A减少 53 %。这给出了两种工程路线的量化权衡方案负载侧 THD源侧 THDAHF 注入电流代价仅 AHF20.56 %1.07 %479 AAHF 容量大、开关损耗高多脉波 AHF4.80 %1.13 %223 A需要专用移相变压器10. 结论补偿效果本示例的并联 AHF 把源电流 THD 从负载侧的 22.41 % 降到0.35 %6 脉波工况降幅 98.5 %12 脉波工况下从 4.80 % 降到 1.13 %。负载侧 THD 与官方文档完全吻合。谐波选择性6 脉波负载的谐波集中在 5、7、11、13 次5 次高达 17.94 %加入 Δ/Y 整流器构成 12 脉波后5、7 次被抵消至 0.01 % 以下谐波谱跳变到 11、13 次。AHF 对两种工况均有效源侧所有谐波都被压到 0.4 % 以下。动态性能负载阶跃与整流器投切两种扰动下AHF 在1.75 个工频周期内重新跟踪参考电流。参数敏感度按影响从大到小控制器增益基准I_max 交流电感L_af 负载平衡度 直流电容C_dc。I_max是稳定性开关环路增益放大 2 倍即进入过补偿/失稳区设计值 4000 A 已在最佳点L_af是跟踪能力的开关超过 2 mH 后 THD 急剧恶化4 mH 时 5.08 %8 mH 时 11.36 %设计值 1 mH 为最优C_dc只影响直流纹波100 µF → 2000 µF 纹波降 11 倍对交流补偿几乎无影响。工程权衡AHF 单独使用即可达到约 95 % 的谐波治理率引入 12 脉波整流可把 AHF 所需注入电流减少 53 %479 A → 223 A但代价是移相变压器与严格的负载平衡要求两侧电流不匹配时负载 THD 从 4.80 % 恶化到 12.36 %。11. 文件清单与复现步骤11.1 目录结构ActiveHarmonicFilterExample/ ├── 有源谐波滤波器仿真分析报告.md 本报告 ├── ahf_export2.m 数据导出脚本挂接 To Workspace导出 14 路信号 ├── ahf_plot.m 结果绘图脚本波形/频谱/THD纯 MATLAB ├── ahf_sweep_run.m 参数对比仿真脚本脚本形式需在 base 工作区执行 ├── ahf_plot_sweep.m 参数扫描绘图脚本 ├── extra_stats.m 附加统计谐波明细、电压等级、母线纹波 ├── shoot_model.m / shoot_deep.m / shoot_more.m / shoot_powergui.m 模型界面截图脚本 ├── archive_261006/ 一次性侦察脚本留档 ├── ahf_scope.mat 基准工况全量数据14 路信号 ├── sweep_laf.mat / sweep_cap.mat / sweep_gain.mat ├── sweep_step.mat / sweep_breaker.mat 五组扫描结果 ├── slprj/ Simulink 编译缓存可再生 └── figs/ 全部图片22 张与汇总表 │ ├── —— 模型界面截图6 张—— ├── m1_model.png 模型总图 ├── m2_ahf.png AHF 子系统内部 ├── m3_hysteresis.png 滞环电流控制器 ├── m4_measurements.png 测量子系统 ├── m5_inverter.png IGBT 逆变桥 ├── m6_powergui.png powergui 求解器设置Discrete 2e-06 s │ ├── —— 仿真结果10 张—— ├── w1_补偿效果_负载vs源电流.png 核心负载电流 vs 源电流 ├── w2_三路电流.png 负载/源/滤波器注入电流 ├── w3_全时段_负载与源电流.png 三个阶段全览 ├── w4_直流母线电压与电流.png Vdc、Idc 全时段 ├── w5_频谱_第3周期_6脉波.png 6 脉波谐波谱 ├── w6_频谱_第13周期_12脉波.png 12 脉波谐波谱 ├── w7_THD对比柱状图.png 三工况 THD 对比 ├── w8_参考电流跟踪.png 动态响应 ├── w9_谐波结构对比.png 6 脉波 vs 12 脉波 ├── w10_三相电压电流.png 三相电压电流 │ ├── —— 参数扫描6 张—— ├── s1_电感Laf扫描.png 交流电感影响 ├── s2_直流电容扫描.png 直流电容影响 ├── s3_控制器增益扫描.png 控制器增益影响 ├── s4_负载电流扫描.png 负载水平与平衡度 ├── s5_投入时刻扫描.png 12 脉波投入时刻 ├── s6_灵敏度矩阵.png 参数灵敏度热力图 │ ├── ahf_stats.csv 三工况统计表 └── sweep_summary.csv 全部 23 次扫描运行的指标明细11.2 复现步骤第一步导出仿真数据MATLAB R2023bcd(本示例目录)ahf_export2% 自动挂接 To Workspace 并运行仿真生成 ahf_scope.mat第二步绘制结果图ahf_plot% 生成 figs/w1..w10 与 ahf_stats.csvextra_stats% 打印谐波明细等附加统计第三步模型界面截图shoot_model shoot_deep shoot_more% 生成 figs/m1..m6第四步参数对比仿真逐组运行每组约 1 分钟SWEEP_NAMElaf;ahf_sweep_run SWEEP_NAMEcap;ahf_sweep_run SWEEP_NAMEgain;ahf_sweep_run SWEEP_NAMEstep;ahf_sweep_run SWEEP_NAMEbreaker;ahf_sweep_run ahf_plot_sweep% 生成 figs/s1..s6 与 sweep_summary.csv⚠️ahf_sweep_run必须以脚本形式在 base 工作区执行先赋SWEEP_NAME再运行不能写成函数——原因见第 12 章第 2 条。12. 分析过程中发现的问题与处理本章记录本次分析中踩到的坑供后续复用。12.1 powergui 模型的数据导出模型的电气网络由powergui求解sim()在默认配置下只返回时间轴ReturnWorkspaceOutputsoff时返回tout不产生任何数据变量。早期尝试过但全部失败的路径包括尝试结果模型级LoggingToFile只生成含SimulationMetadata的空 mat 文件线对象DataLogging/SignalLoggingR2023b 已移除这些线参数Scope 块自身的DataLoggingDataLoggingVariableName参数可设置成功但变量不落盘powergui 内置 FFT 工具saveon不生成文件从 Scope 图形对象读取数据ScopeLineAnimator无可读数据属性从 Scope 截图数字化反解波形曲线跳跃、中文变方框THD 误差超过 100 %不可用可行方案为每个示波器的每一个输入端口从其源端口引支线接到To WorkspaceSaveFormat Structure With Time。关键点是必须逐端口处理——responses示波器有 8 个输入端口只处理第一个会丢掉 7 路最有价值的信号。12.2 从函数内调用sim时 To Workspace 变量不落盘这是本次最耗时的问题。现象把参数扫描写成function后5 次不同参数运行的指标完全相同。排查结论当sim()从函数内部发出时To Workspace模块的变量不会写入 base 工作区因此evalin(base, ...)读到的一直是上一次留下的陈旧数据看起来就像参数不起作用。处置把扫描脚本改写为脚本文件ahf_sweep_run.m在 base 工作区执行并先用SWEEP_NAME传入配置。脚本内的局部函数有独立作用域故所有度量计算直接内联在脚本主体中。12.3save_system的另存为会重命名内存中的模型为做参数扫描曾用save_system(mdl, ahf_work.slx)生成工作副本。但该调用具有Save As 语义——它把内存中的模型重命名为ahf_work导致ActiveHarmonicFilter从内存中消失。此后所有load_system(ActiveHarmonicFilter)都从磁盘重新加载原始模型之前挂接的 14 个To Workspace块全部丢失仿真自然也就不再输出数据。处置不再用save_system生成副本参数变更通过set_param(块, 参数, 值)set_param(mdl,SimulationCommand,update)直接施加于当前已加载且已挂接记录块的模型。12.4 遮罩间接引用导致参数变更不被检测AHF 子系统的块参数写作表达式platform.activefilter.Laf其中platform是遮罩工作区变量遮罩参数platform AF。在 base 工作区重新赋值AF结构体后Simulink 不会检测到该变化仿真仍沿用已编译模型表现为改参数无效果。同时经全模型扫描确认AF.activefilter.voltage字段未被任何模块引用改动它必然无效——这也是早期扫描结果完全相同的原因之一。处置直接修改块本身的对话参数set_param(LAF,Inductance,8e-3)并在每次修改后强制update触发重编译。扫描结束后再还原为原表达式。12.5 中文字体名在中文 Windows 上是本地化名MATLAB 绘图默认字体用SimSun会找不到字体在中文版 Windows 上listfonts返回的是本地化名称宋体。本报告全部图件使用set(ax,FontName,宋体)与set(0,defaultAxesFontName,宋体)设置中文与 LaTeX 数学符号i_{load}等均正常渲染。12.6 谐波谱图必须用离散竖线而非柱状图在对数纵轴上使用bar绘制谐波谱时柱体从 0 起画而 0 在对数轴上不存在MATLAB 会把柱体拉到坐标轴底部导致所有柱体连成一片色块完全无法读数。处置改为手工绘制离散竖线stem 标记点这也是电力电子谐波谱的标准表现形式。报告完