发明名称 一种适用于滑模变结构控制逆变器的新型分岔图绘制方法
摘要 本发明公开了一种适用于滑模变结构控制逆变器的新型分岔图绘制方法,包括如下步骤:1、建立起滑模变结构控制离散迭代映射模型;2、分析现有分岔图绘制方法,确定其判定稳定域和分岔行为不准确的原因,以此寻找解决方法;3、研究状态突变现象的机理;4、根据状态突变机理,推导出变采样点绘制逆变器分岔图的方法,并绘制相应的图形。本发明提出一种变采样点分岔图绘制方法,经验证该方法绘制的分岔图不仅能准确地确定系统稳定工作域,而且能用于分析系统的分岔行为,最后通过电路实验验证了理论分析的正确性,研究结果为研究滑模变结构控制逆变器的非线性行为及其分岔图的绘制提供参考。
申请公布号 CN106126833A 申请公布日期 2016.11.16
申请号 CN201610497495.2 申请日期 2016.06.28
申请人 西南石油大学 发明人 代云中;蒋林;吴显松;吴俊;任海军;赵鹏程;樊新建
分类号 G06F17/50(2006.01)I;H02J3/38(2006.01)I;H02J3/00(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 成都正华专利代理事务所(普通合伙) 51229 代理人 李林合;李蕊
主权项 一种适用于滑模变结构控制逆变器的新型分岔图绘制方法,其特征在于,步骤一、离散迭代映射建模:ε为滑模控制系数;i为并网电流,i<sub>ref</sub>为并网参考电流,其中i<sub>ref</sub>(t)=I<sub>ref</sub>×sin(2πf<sub>1</sub>t),I<sub>ref</sub>为i<sub>ref</sub>的幅值;为了消除SMC控制不连续性带来的抖振现象,在SMC控制器中引入了比例环节,比例系数k<sub>p</sub>,得到调制信号(u<sub>control</sub>)为:u<sub>control</sub>=ε×sign(i<sub>ref</sub>‑i)+k<sub>p</sub>×(i<sub>ref</sub>‑i)                (1);将u<sub>control</sub>和单极性三角波载波u<sub>carrier</sub>进行比较,产生PWM驱动信号,并生成S<sub>1</sub>~S<sub>4</sub>的PWM驱动信号图形;u<sub>control</sub>与零信号比较,产生S<sub>5</sub>和S<sub>6</sub>的PWM驱动信号;其中i<sub>ref</sub>(t)=I<sub>ref</sub>*sin(2πf<sub>1</sub>t),I<sub>ref</sub>为i<sub>ref</sub>的幅值,f<sub>1</sub>为i<sub>ref</sub>的频率;将n时刻的参考电流i<sub>ref</sub>(n)与n时刻的负载电流i(n)比较后,经SMC控制得到系统第n时刻的占空比d<sub>n</sub>,其计算方法如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>n</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><mi>&epsiv;</mi><mo>&times;</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><mi>i</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><mi>k</mi><mi>p</mi></msub><mo>&times;</mo><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><mi>i</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><mi>i</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>&GreaterEqual;</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>&epsiv;</mi><mo>&times;</mo><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><msub><mi>i</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><mi>k</mi><mi>p</mi></msub><mo>&times;</mo><mrow><mo>(</mo><mi>i</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><msub><mi>i</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><mi>i</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>&lt;</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001033262830000011.GIF" wi="1558" he="158" /></maths>设逆变器两桥臂中点A,B间的电压为U<sub>AB</sub>;当U<sub>AB</sub>输出为±E和0时分别定义为1和0电平,且定义当i从A流向B时为正;根据开关管的导通状况可将其分为多种开关模态,再根据负载电感电流i方向不同分为多种工作状态,并形成各自的状态方程;系统离散迭代映射建模采用频闪映射思想,将i在n+1时刻值i<sub>n+1</sub>用n时刻值i<sub>n</sub>来表示;在不同的时间段内频闪映射模型不同;L和R为负载电感、负载电阻,令a=E/R,b=L/R,u<sub>grid</sub>=U<sub>m</sub>×sin(ωt),U<sub>m</sub>为u<sub>grid</sub>的幅值,T<sub>s</sub>=1/f<sub>s</sub>为载波周期;当i<sub>ref</sub>&gt;0时,U<sub>AB</sub>为+E和0,以T<sub>s</sub>作为频闪采样周期,由一种模态的状态方程可得i的频闪采样模型为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>i</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><msub><mi>i</mi><mi>n</mi></msub><mo>-</mo><mn>0.5</mn><mi>a</mi><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>(</mo><msub><mi>d</mi><mi>n</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo><mo>&times;</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mi>m</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&omega;nT</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mn>2</mn><mi>R</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001033262830000012.GIF" wi="1548" he="111" /></maths>当i<sub>ref</sub>&lt;0,U<sub>AB</sub>为‑E和0,由另一种模态的状态方程得到i的频闪采样模型为<maths num="0003"><math><![CDATA[<mrow><msub><mi>i</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><msub><mi>i</mi><mi>n</mi></msub><mo>-</mo><mn>0.5</mn><mi>a</mi><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>(</mo><msub><mi>d</mi><mi>n</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo><mo>&times;</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mi>m</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&omega;nT</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>b</mi></mfrac></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mn>2</mn><mi>R</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001033262830000013.GIF" wi="1551" he="118" /></maths>步骤二、滑模变结构控制逆变器状态突变位置产生的原理:根据步骤一种得到的结构,滑模变结构控制逆变器状态突变位置产生的原理为:误差信号e=i<sub>ref</sub>‑i,当系统误差信号e在第n个开关周期出现e(n)=i<sub>ref</sub>(n)‑i(n)&lt;0的情况时,u<sub>control</sub>(n)=‑ε+k<sub>p</sub>×e(n)&lt;0,u<sub>control</sub>(n)&lt;u<sub>carrier</sub>(n),因此占空比d<sub>n</sub>将突变为0,电感电流将急剧下降,导致在n+1个开关周期出现e(n+1)=i<sub>ref</sub>(n+1)‑i(n+1)&gt;0,u<sub>control</sub>(n)=‑ε+k<sub>p</sub>×e(n)&gt;0,u<sub>carrier</sub>(n+1)与u<sub>control</sub>(n+1)必有两个交点,第n+1个开关周期的占空d<sub>n+1</sub>&gt;0;因此在电感电流正半周期下降段,脉宽调制PWM控制信号的占空比出现了d<sub>n</sub>&lt;d<sub>n+1</sub>非单调下降的情况,导致系统出现了状态突变的不稳定现象;步骤三、突变位置的确定:在每个正弦周期采样N个点,其中N=f<sub>s</sub>/f<sub>1</sub>,f<sub>s</sub>为采样率,f<sub>1</sub>为i<sub>ref</sub>的频率;将步骤二中得到的占空比代入,选择从N<sub>0</sub>时(N<sub>0</sub>为大于π/2的任意选定值)的采样点开始取电流波形下降段M个开关周期连续两个开关周期的占空比d<sub>n</sub>和d<sub>n+1</sub>作差再除以两者差的绝对值,把计算出的M个数相加得到P<sub>M</sub>可表示如下:<maths num="0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>M</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><msub><mi>N</mi><mn>0</mn></msub></mrow><mrow><msub><mi>N</mi><mn>0</mn></msub><mo>+</mo><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><mfrac><mrow><msub><mi>d</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>d</mi><mi>n</mi></msub></mrow><mrow><mo>|</mo><msub><mi>d</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>d</mi><mi>n</mi></msub><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001033262830000021.GIF" wi="1126" he="150" /></maths>若在M个开关周期内占空比一直保持单调变化则P<sub>M</sub>=M时,系统即处于稳定工作状态;若P<sub>M+1</sub>&lt;M+1,系统处于分岔不稳定状态,连续计算出前M和M+1个开关周期的P<sub>M</sub>和P<sub>M+1</sub>的值,若P<sub>M</sub>=M且P<sub>M+1</sub>&lt;M+1,则确定状态突变周期为第N<sub>0</sub>+M个开关周期;步骤四、分岔图绘制方法:根据步骤三得到的突变点数据绘制分岔图,在合理范围内,取分岔参数ε为任意初始值代入离散迭代方程中,当P<sub>M</sub>=M时,系统处于稳定工作状态,可选择任意采样点为多个电源周期的采样位置绘制分岔图;当P<sub>M</sub>=M且P<sub>M+1</sub>&lt;M+1,则表示系统出现了分岔不稳定现象,通过状态突变周期判据可判断状态突变周期为N<sub>0</sub>+M,则以第N<sub>0</sub>+M个采样点为多个电源周期的固定采样位置绘制分岔图。
地址 610500 四川省成都市新都区新都大道8号
您可能感兴趣的专利