发明名称 一种兆瓦级三相变流器神经网络内模控制方法
摘要 本发明公开了一种兆瓦级三相变流器神经网络内模控制方法,将神经网络运用到内模控制中,针对兆瓦级变流器在旋转αβ坐标系下的特点,提出了一种基于神经网络内模原理的控制方案,将平衡三相变流器模型经3S/2S坐标变换之后,等效为两个相互独立的单相变流器。因此,只需两个完全相同的单相控制器分别对α相和β相进行控制,就可以实现对整个兆瓦级变流器的控制。本发明的方法能够克服系统参数摄动或外界扰动对系统造成的影响,使系统具有良好的鲁棒性和抗扰性。
申请公布号 CN103368443B 申请公布日期 2015.12.02
申请号 CN201310290454.2 申请日期 2013.07.10
申请人 温州大学 发明人 戴瑜兴;郜克存;郑崇伟;全惠敏;曾国强
分类号 H02M7/5395(2006.01)I 主分类号 H02M7/5395(2006.01)I
代理机构 长沙正奇专利事务所有限责任公司 43113 代理人 马强
主权项 一种兆瓦级三相变流器神经网络内模控制方法,包括兆瓦级三相变流器,所述兆瓦级三相变流器包括IGBT逆变模块,所述IGBT逆变模块通过LC滤波电路接负载;所述IGBT逆变模块包括直流侧电容和三个并联的桥臂,所述直流侧电容与所述桥臂并联,所述桥臂由两个IGBT串联组成,其特征在于,该方法为:1)计算兆瓦级变流器在静止abc坐标下的数学模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><msup><mi>A</mi><mi>T</mi></msup></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>A</mi></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><mi>u</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>I</mi><mo>&CenterDot;</mo></mover><mi>l</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>/</mo><mi>C</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mn>3</mn></msub></mrow></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mi>r</mi><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>A</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>u</mi><mi>u</mi></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>l</mi></msub></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mi>C</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>A</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>U</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mrow><mi>o</mi><mi>p</mi><mi>h</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>,</mo></mrow>]]></math><img file="FDA0000794472210000011.GIF" wi="855" he="374" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000794472210000012.GIF" wi="373" he="206" /></maths>I<sub>3</sub>为三阶单位矩阵,u<sub>u</sub>=[u<sub>ab</sub> u<sub>bc</sub> u<sub>ca</sub>]<sup>T</sup>为LC滤波电路电容电压,I<sub>l</sub>=[I<sub>a</sub>I<sub>b</sub>I<sub>c</sub>]<sup>T</sup>为LC滤波电路电感电流,I<sub>oph</sub>=[I<sub>oa</sub> I<sub>ob</sub> I<sub>oc</sub>]<sup>T</sup>为LC滤波电路输出电流,r为LC滤波电路电感损耗,C为LC滤波电路滤波电容值,L为LC滤波电路电感值,U<sub>i</sub>为IGBT逆变模块三相输出电压,i=A、B、C;2)对所述数学模型进行坐标变换,得到:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = '' close = ''><mtable><mtr><mtd><mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mi>C</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>T</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mo>(</mo><mrow><msub><mi>U</mi><mi>d</mi></msub><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>AT</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><mrow><mi>&alpha;</mi><mi>&beta;</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mrow><mi>&alpha;</mi><mi>&beta;</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>/</mo><mi>C</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>T</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>AT</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mi>r</mi><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>AT</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mi>C</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>T</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mo>(</mo><mrow><msub><mi>U</mi><mi>d</mi></msub><mo>/</mo><mi>L</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>AT</mi><mrow><mn>2</mn><mi>S</mi><mo>/</mo><mn>3</mn><mi>S</mi></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>S</mi><mrow><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mrow><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000794472210000013.GIF" wi="1313" he="271" /></maths>其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mn>3</mn><mi>S</mi><mo>/</mo><mn>2</mn><mi>S</mi></mrow></msub><mo>=</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></mtd><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msqrt><mn>3</mn></msqrt><mo>/</mo><mn>2</mn></mrow></mtd><mtd><mrow><mo>-</mo><msqrt><mn>3</mn></msqrt><mo>/</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000794472210000014.GIF" wi="646" he="172" /></maths><img file="FDA0000794472210000015.GIF" wi="81" he="75" />为αβ坐标系下IGBT逆变模块三相输出电压;<img file="FDA0000794472210000016.GIF" wi="71" he="88" />为αβ坐标系下IGBT逆变模块三相输出电流;S<sub>α,β</sub>为αβ坐标系下IGBT逆变模块开关状态;I<sub>α,β</sub>为αβ坐标系下LC滤波电路输出电流;U<sub>d</sub>为直流母线电压;3)利用步骤2)的公式,得到兆瓦级变流器在αβ坐标系下数学模型:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mn>1</mn><mrow><mn>3</mn><mi>C</mi></mrow></mfrac><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>i</mi><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mfrac><mn>1</mn><mrow><mn>3</mn><mi>C</mi></mrow></mfrac><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>I</mi><mrow><mi>o</mi><mi>&alpha;</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mrow><mi>o</mi><mi>&beta;</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000794472210000021.GIF" wi="806" he="227" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><mo>-</mo><mi>r</mi><mo>/</mo><mi>L</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>r</mi><mo>/</mo><mi>L</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>i</mi><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>u</mi><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>U</mi><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000794472210000022.GIF" wi="994" he="175" /></maths>[u<sub>α</sub>u<sub>β</sub>]<sup>T</sup>=T<sub>3S/2S</sub>·[u<sub>ab</sub> u<sub>bc</sub> u<sub>ca</sub>]<sup>T</sup>,其中,[I<sub>oα</sub>I<sub>oβ</sub>]<sup>T</sup>=T<sub>3S/2S</sub>·[I<sub>oa</sub> I<sub>ob</sub> I<sub>oc</sub>]<sup>T</sup>,[i<sub>α</sub>i<sub>β</sub>]<sup>T</sup>=T<sub>3S/2S</sub>·[I<sub>a</sub> I<sub>b</sub> I<sub>c</sub>]<sup>T</sup>,i<sub>α</sub>、i<sub>β</sub>为αβ坐标系下兆瓦级变流器输出电流、U<sub>α</sub>、U<sub>β</sub>为αβ坐标系下IGBT逆变模块输出电压、u<sub>α</sub>、u<sub>β</sub>为αβ坐标系下兆瓦级变流器输出电压;4)利用αβ坐标系下α相兆瓦级变流器输出电压建立神经网络内模控制器模型:y<sub>α</sub>(k)=h<sub>1</sub>u<sub>α</sub>(k‑1)+h<sub>2</sub>u<sub>α</sub>(k‑2)‑g<sub>1</sub>y<sub>α</sub>(k‑1)‑g<sub>2</sub>y<sub>α</sub>(k‑2),其中,y<sub>α</sub>(k)为α相兆瓦级变流器输出电压的离散值,h<sub>1</sub>、h<sub>2</sub>、g<sub>1</sub>、g<sub>2</sub>均为d轴变流器传递函数系数;u<sub>α</sub>(k)为神经网络逆模型;<sub>k</sub>表示第<sub>k</sub>时刻;5)将上述y<sub>α</sub>(k)、u<sub>α</sub>(k)作为神经网络内部模型输入,得到神经网络内部模型表达式:y<sub>α</sub>(k+1)=f[y<sub>α</sub>(k),…,y<sub>α</sub>(k‑n+1),u<sub>α</sub>(k),…,u<sub>α</sub>(k‑m+1)]+d<sub>α</sub>(k),其中d<sub>α</sub>(k)为d轴干扰噪声;f(x)为神经网络隐层输入和隐层输出之间的函数映射关系;<img file="FDA0000794472210000023.GIF" wi="376" he="125" /><img file="FDA0000794472210000024.GIF" wi="1180" he="120" /><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>v</mi></munderover><msub><mi>w</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>v</mi><mo>=</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000794472210000031.GIF" wi="865" he="143" /></maths>w<sub>ij</sub>为网络输入节点与隐节点间的网络权值,θ<sub>i</sub>(k)为输入层的阈值,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>u</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>j</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mn>0</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>m</mi><mo>-</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>j</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>m</mi><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000794472210000032.GIF" wi="892" he="156" /></maths>m为u<sub>α</sub>(k)的阶次;n为y<sub>α</sub>(k+1)的阶次;6)假设输出无偏差跟踪输入,则y<sub>α</sub>(k)≈r<sub>α</sub>(k),…,y<sub>α</sub>(k‑n+1)≈r<sub>α</sub>(k‑n+1),得到神经网络逆模型表达式:u<sub>α</sub>(k)=f<sup>‑1</sup>[r<sub>α</sub>(k),…,r<sub>α</sub>(k‑n+1),u<sub>α</sub>(k‑1),…,u<sub>α</sub>(k‑m+1),e<sub>f</sub>(k)],其中,e<sub>f</sub>(k)为e<sub>f</sub>(s)的离散值,e<sub>f</sub>(s)=G<sub>f</sub>(s)e<sub>m</sub>(s)=G<sub>f</sub>(s)[y<sub>α</sub>(s)‑y<sub>m</sub>(s)],G<sub>f</sub>(s)为反馈函数,y<sub>α</sub>(s)为y<sub>α</sub>(k)的连续函数,y<sub>m</sub>(s)为神经网络输出层输出y<sub>m</sub>(k)的连续函数,<img file="FDA0000794472210000033.GIF" wi="470" he="141" />v<sub>i</sub>(k)为神经网络隐层节点与输出节点间的网络权值,z<sub>i</sub>(k)为神经网络隐层输出,r<sub>α</sub>(k)为当前时刻α相兆瓦级变流器参考输入电压;w为神经网络隐层输出个数;7)定义J(k)为神经网络学习目标函数,则神经网络训练性能指标函数J为:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>J</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><msup><mrow><mo>&lsqb;</mo><msub><mi>r</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>u</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mn>2</mn></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>e</mi><mi>c</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&epsiv;</mi><mn>2</mn></msub><mo>,</mo></mrow>]]></math><img file="FDA0000794472210000034.GIF" wi="884" he="138" /></maths>其中,ε<sub>2</sub>&gt;0;N为神经网络输入或输出样本的长度;8)求神经网络训练性能指标函数J的最小值,确定u<sub>α</sub>(k),从而求得y<sub>α</sub>(k),使y<sub>α</sub>(k)很好地跟踪输入信号r<sub>α</sub>(k)。
地址 325035 浙江省温州市高教园区