1、A 题：级联型 H 桥变换器的阶梯波特定消谐技术研究摘要我们的模型主要包括四部分：1）基于经验公式与波形拟合相结合的初值设定模型，用于获取对不同单元级联个数 n 能使非线性消谐方程尽快收敛的初值；2）等步长的检验搜索模型，用于确定（1） 、 （2） 、 （3）问中特定单元个数下使非线性消谐方程有解的调制比幅值 m 的取值范围；3）基于步长搜索模型中开关角范围的包含关系分析模型，用以分析出 m 特定时随 n 增加所解得开关角的包含关系以及输出电压波形质量改善效果；4）通过设定控制角改变单元输出正负的功率均衡求解模型使各 H 桥变换器单元在一定周期内实现功率均衡。在问题一中，主要思路就是初值设定使
2、阶梯波尽可能拟合正弦波；基于此想法建立了一种基于波形拟合的初值设定模型，在n和m 为已知的情况下，把阶梯波的转折上边沿点和正弦曲线拟合，根据在正弦曲线上的阶梯拐点的纵坐标运用反三角函数确定初始的开关角；针对不能取反三角的 m范围对最后一个初值点采用前点切线取角的方法，得到完整的开关角初值。解方程采用matlab中fsolve函数迭代求解。但是该模型对 m的范围要求较高，超出范围便不适用；采用一种综合简便性和精确性考虑较优的经验公式，建立了改进的补充经验公式的初值设定模型，得到较为完整的初值设定方法。对 n=3， m=0.8和m=0.5 以及 n=5,15 时不同 m 值的情况设定初值后求解结果
3、表明初值与结果差值很小，而且均能较快收敛。针对问题二，对特定的模块数 n ， m只有在一定范围内开关角求解方程才会有解；因为要直接解非线性超越方程中 m的范围难以实现；所以建立了一种等步长的检验搜索模型，采用对 m的值在一定精度下进行搜索使方程有解来确定其范围的方法，在给定步长下变化 m的值，代入方程看是否有解，逐步搜索出 m在该精度下的范围。在 n=3,5,15 的情况下搜索得出m的范围，在 m范围内根据公式计算THD的结果表明THD随着m变化波动范围很小，说明求出的范围较为合理。在问题三求解上，沿用步长搜索方法建立了基于步长搜索模型中开关角范围的包含关系分析模型。基于问题二的 m 范围，在
4、 n=3,5,15 都适用的 m 范围内搜索，过程中对确定 m 值分别求解不同 n 的开关角来分析解的相互包含关系，并计算 THD 来反映波形质量；最终结合所有 m 值的分析结果来综合评估得到不同 n 的解的相互包含关系，以及 n 增加的波形质量总体改善特性。取 m=0.85 的开关角解的相互包含关系结果表明了模型反映了解具有包含关系，波形质量的改善特性也反映了随 n 增加质量更优的结果。对于问题四，为了实现功率均衡，对每一模块单元设定三个控制角，使单元的输出为先正后负再正，综合功率相等与谐波消除列出方程，增加一定数目功率相等约束方程,同时减少一定数目消谐波方程，解方程求解控制角。对应控制角的
5、值，控制各模块开关 S1i，S 2i，S 3i，S 4i 。综上所述，模型在初值设定下开关角的求解值较为精确；使开关角有解的m 范围求解较为可靠；分析解的包含关系中 THD 随着 n 增大而减小，反映对谐波的消除效果较好；并且在单元功率均衡优化开关控制策略上能给出可行的方案。模型具有一般性，可推广到实际工程的消谐技术中。关键词： 非线性消谐方程，初值设定，调制比幅值范围，功率平衡I 问题重述与分析1.1 问题背景在电力系统中实现高电压、多电平输出的研究中，解决脉冲宽度调制（PWM ）较广泛的方法主要是 SPWM(Sinusoidal PWM)法，其中运用三电平级联 H 桥变换器的特定谐波消除脉
6、宽调制技术（Selected Harmonic Elimination Pulse Width Modulation，SHEPWM）是实现 SPWM 法的一种方案。关键是通过选择特定的开关时刻，在满足期望的输出基波电压 vac 的同时，来消除选定的低次谐波，进而改善输出电压的波形质量，根据题意分析和已有的研究成果，发现初值设定对于求解结果的准确性和收敛性有影响。1.2 问题的分析与拆解1.2.1 问题剖析在 n 个单元串联的 H 桥变换器系统中，对于第 i 个 H 桥变换器单元，控制开关可得到的输出电压 vaci 主要有 Vdci，0，- Vdci；当 H 桥变换器单元直流侧独立电压 Vdci
7、 都为 Vdc 时，可输出（ 2n+1）电平数的阶梯型电压 vac，通过对该波形进行傅里叶级数分解，以导通角 表示 n 个单元分别导通时序，那么 vac 的第is 个奇数次谐波的幅值的傅里叶级数可以改为式（ 1）：121,354cos()()cos()dcs nsVs （1）其中： 1202n 。因此消谐技术的关键就是试找出一组 i（i=1, n） ，使输出电压的基波分量幅值为V 1m，且不含有低次谐波。具体方法如下：定义调制比幅值m=V 1m/(nVdc)，那么根据上述消除谐波的要求，令低频率的展开方程部分为0，便可写出消谐关于 i（i=1, n）的非线性代数方程组式（2） 。 11212c
8、os()s()cos()45c50s(7)s()s(7).mndcnV (2)另外，对于特定的级联模块数n，调制比m在一定范围内取值才能使非线性方程组有解。并且需要一组初值 i0（i=1, n）使方程组求解收敛。如需校核阶梯波特定消谐技术的输出电压v ac波形质量，可通过总电压谐波畸变率（Total Harmonics Distortion，THD）来描述，如式（3） 。(3)%10%10214372523,16 VVTHDkii 可见低次谐波若占有输出的部分越少，THD 的值就越小，消谐得到的波形质量越高。因此，综上分析，式（2）为模型求解约束条件（根据此约束条件进行求解），式（3）数值尽量
9、小为约束求解下的目标函数。1.2.2 问题拆解综合题意指定及剖析，原问题的五问可拆解分为以下几个子问题模型的建立和求解：（1）初值设定模型：如何建立合理的初值确定模型使式（2）的非线性超越方程求解尽快收敛，并满足一定精度要求；（2）迭代求解模型：给定 n， m 时，如何对方程（ 2）中 i（i=1, n）进行快速求解：针对问题 1 至问题 3 中 i 的角度求取；（3）有解判定及搜索模型：给定 n 时如何迭代收搜使方程（2）有解的 m范围：针对问题 1 至问题 3 中 m 范围的求取；（4）解的延伸与判定模型：给定 m 情况下，n 从小变大时 i（i=1, n）解是如何拓展延伸的，具备何种规律
10、，以及如何编程建立基于 THD 的波形质量判定程序：针对前三问中 THD 求取及第四问；（5）功率均衡求解模型：如何对各单元的通断角度进行综合匹配设定，使得各模块单元功率相等：针对第五问。1.3 总体建模求解思路根据以上对问题的剖析和拆解，我们的总体建模求解思路如下：（1）初值设定模型：阶梯型波上拐角点与正弦曲线重合；（2）迭代求解模型：对不同的 n，通过初值设定模型求取初值，通过牛顿迭代求解；（3）有解判定及搜索模型：取给定步长下变化 m 的值，求出对应初值，再进行迭代求解，确定 m 的取值范围；（4）解的延伸与判定模型：通过 m 范围选取合适的 m 值，求出不同模块数 n 所对应开关角进行
11、讨论分析；（5）功率均衡求解模型：所有单元在 1/4 周期内先正后负再正。II 基本假设1、假设每个模块单元只有三个控制角度。2、假设每个模块单元完全相同。III 符号说明M 调制比幅值N 串联的电平数i0 开关角的初值i 开关角的解Vdc H桥变换器单元直流侧独立电压V1m 电压的基波分量幅值搜索 m时的步长THD 总电压谐波畸变率开关角初值与求解结果误差Vs 第s个奇数次谐波幅值Vci 单元i的输出电压基波幅值S 各模块单元的变换功率方差IV 初值设定模型及迭代求解模型4.1 迭代求解模型方程组式（2）为非线性超越方程组，解决该类方程组比较常用的一种数值方法是牛顿迭代法，Matlab 内置
12、的基于最小二乘的 fsolve 函数，可采用高斯-牛顿迭代法搜索，加快对求解的收敛速度，因次我们考虑使用 fsolve 函数进行求解。4.2 初值设定模型采用迭代求解模型的缺点是对初值的选取有一定的要求，初值选取对方程组求解收敛性和结果的精度有直接影响 1。所以解决问题一的建模关键在于建立一个能给出可靠初值的模型方法。 方程组式（2）为非线性超越方程组，解决该类方程组比较常用的一种数值方法是迭代法，我们考虑使用 Matlab 软件内置的 fsolve 函数进行求解，但是该方法对于初值的选取有一定的要求，初值选取对方程组求解收敛性和结果的精度有直接影响。所以解决问题一的建模便转换为建立一个能给出
13、可靠初值的模型方法。4.2.1 模型建立的准备给出 n=3,m=0.8;在运用模型给出初值 、 、 的情况下采用 Matlab10230中 fsolve 函数迭代求解方程对应的变量值为 、 、 。4.2.2 模型建立基于波形拟合的初值设定模型对于 n 已知的情况，H 桥变换器单元输出电压叠加消谐得出的阶梯形电压如图 1 所示。V d c2 V d c3 V d cn V d c 123n- n V d c- 3 V d c20图 1 叠加消谐得出的阶梯形电压实际采用交流波形叠加通过开关角消谐得到阶梯波，使其更接近正弦波。基于此原理，当阶梯波的阶梯高度 和 n 已知时，采用阶梯曲线与正弦曲线拟d
14、cV合的方法，视纵坐标为 （k=1,2.n）的阶梯拐角点 刚好都在正dckdcVk，弦线上时为较好拟合；则可以得到较可靠的初值 便为 的值。i0k由调制比公式 m=V1m/(nVdc)可得，(4)nmdc1取V 1m=1，可以得到初值，= （k=1,2.n-1） (5)i0k1-s由于可能m值不同可能出现如下图出现的两种拟合情况(a)m1(b)m1:正弦曲线能包括阶梯波形，所以= (6)n0m1si-(2) m nVdc，对第n个初值只能采用近似取值的方法，采用(n-1)点出的导数作为斜率，取 至 终点为区间的斜直线的中点，其纵坐标对应的 值作为01-）（2/ 第n个初值为= (7)n00)1
15、(0)1(0)1(2/c/ nnnos最终可以得到n个基本的初值，并代入方程求得开关角的解。4.2.3 模型的结果和分析n=3，m=0.8 的初值分别为：=0.4298 =0.9851 =1.1470102030代入方程运用fsolve函数求解结果为：=0.5103 =0.9501 =1.1255123回代到原方程组得：(8)12345123*3cos()()cos()1.890.0.855 -7(7)()().61从结果可知，此时求出的开关角的值非常精确，对于此情况得出结果的迭代收敛，该初值选取的模型比较可靠。 4.3 初值设定模型的改进和求解4.3.1 初始模型的局限性对n取某一特定值，通
16、过分析可知当(n-1)/(n*m)1时， 没有实nmk1-si数解，模型计算初值不再适用；即m的范围为： ，在此范围内改变mn)(对应不同的n的解的情况如表(1)。表1 不同m和n情况下开关角的求解m=0.9 m=1.04初值0结果i误差 初值0结果i误差 0.2241 0.2401 -0.0364 0.1935 0.1766 -0.0607 0.4606 0.4703 -0.0253 0.3948 0.2277 -0.03370.7297 0.7596 -0.0012 0.6150 0.4779 -0.00431.0949 0.9933 -0.1467 0.8766 0.7161 0.330
17、0n=5(m0.8)1.2039 1.1587 -0.1441 1.2925 1.0666 0.6476m=0.94 m=1.05初值0结果i误差 初值0结果i误差 0.0710 0.0378 -0.3612 0.0606 0.0321 -0.4269 0.1423 0.1230 -0.0336 0.1215 0.0973 -0.05670.2144 0.2215 -0.0021 0.1828 0.1928 0.03090.2876 0.2222 -0.0452 0.2449 0.1889 0.01210.3625 0.3319 0.0660 0.3079 0.2759 -0.00510.43
18、95 0.4707 0.0276 0.3722 0.3190 0.00520.5195 0.4801 0.0075 0.4381 0.3820 -0.00080.6033 0.6203 0.0120 0.5062 0.4424 -0.00860.6923 0.6826 -0.0502 0.5769 0.5098 0.00550.7884 0.7462 -0.0202 0.6511 0.5845 0.00810.8949 0.9021 -0.0117 0.7297 0.6702 0.0009 1.0180 1.0064 0.0092 0.8143 0.7759 -0.00251.1732 1.1
19、384 0.0238 0.9074 0.9105 0.00091.4516 1.4063 0.0196 1.0131 1.0067 -0.0135n=15(m0.93)1.4587 1.5858 0.0342 1.1411 1.1316 -0.0220在m的范围满足约束时初值选取较为可靠，但是经过验证可以发现在m小于该约束时模型不适用，但取适当的初值方程仍然有解，模型具有局限性。4.3.2 模型的改进补充经验公式的初值设定模型在拟合曲线的初值确定模型的基础上，为扩大模型的适用范围，使得模型在 时仍然能得出可靠的初始值，通过查阅相关文献选取了一种经验n)1(m公式来改进模型；从文献 2的初值确定
20、思路来看，该经验公式在实际工程应用广泛，该公式是以给定适当 的初值得到的解为基准，不断给 m一个很小的0增量，逐步迭代得到所求的m的初值，经验表明用牛顿迭代法可以很快收敛于新的解，所以采用经验公式改进原模型。当 时，据经验公式可以得到 时，初值 计算方法(i 为奇数且n)1( 0i3iN-3)：(9)1020(1)(2)0(1)0.6().6()3)iinnni给迭代初值时，不能让任意两个开关角相等，否则在解线性方程组时 会发生奇异，故采取以下修正方法：(10)0(1)6()0.1.iini将此计算方法添入基于波形拟合的初值确定模型，得到改进后补充经验公式的初值确定模型。4.3.3 改进模型的
21、结果和分析对于n=3的情况，取m=0.5，此时m超出原模型范围(m2/3),此时采用改进模型赋值得：=0.5236 =0.7854 =1.0472102030代入求解得到开关角的结果为：=0.7116 =1.1489 =1.5595123此时(11)1234123*3cos()()cos()1.780.555 29(7)()()-6.1从结果可知，此时求出的开关角的值非常精确，初值的结果迭代收敛，非常可靠。所以改进后的初值确定模型适用范围大，迭代收敛的效果较好；对初值确定体现了优化效果。V 有解判定及搜索模型5.1 模型建立的准备首先对给定 n 的值，分别取 n 为 3,5,15 来对 m 的
22、范围进行研究；当 m 的值确定时，采用问题一的初值设定模型给出初值，求解开关角沿用 fsolve 函数的迭代方法。5.2 模型的建立等步长检验搜索的有解判定模型对于仅给出了 n 的情况，m 此时是变量，我们先考虑 m 最大可能取到的范围。由公式(12)1231 121114cos()().cos()cos().4/()().cos()d mndcmdcnnVnV知 m0。所以 0m4/。由于方程组（2）中的各项均任意次可微，说明此方程组的可解范围也应连续。当 m 在(0, 4/)时，采用逐个 m值点的搜索思路，选取等步长依次增加 m 的值进行搜索，用求解方法找出使方程有解的 m 值，若设步长规
23、定为 ，则有(i=0,1,2.n) (13)1ii其中 m0=0，m n=(4/ - )；把确定的 n，m i 代入方程 (2)求解即可确定 mi 是否符合有解的条件。最后把 mi 对应的 通过图像描绘就可以得到使得方程有解的 mi的范围。具体方法如下：m 从 0 开始，设置 =0.01，m i 的取值依次为 0，0.01，0.021.27。对于每一个 m 值，根据补充经验公式的初值设定模型，分别求取相应的初值，然后通过 fsolve 函数进行开关角 的求解，如果能够求解出符合要求的 值，则认为m 符合要求。根据此时求解出的 值，计算相应的 Vs ，求出对应的 THD。5.3 模型的结果与分析绘制 m 值与求解出的开关角关系图，在图上，若对应一个 m 的取值满足关系 1202n ，即开关角之间大小不能相等；实际上两个开关角大小相差大于 0.05/n 时我们认为不相等，且最大的开关角小于 1.57，最小的开关角大于 0。则认为该 m 值满足要求。当 n=3，m 与开关角 的关系如图：(a) 全局情况(b)m 值左极限放大情况
