CN101242101A - 电力系统负荷谐波稳态模型的建立方法及仿真方法 - Google Patents

电力系统负荷谐波稳态模型的建立方法及仿真方法 Download PDF

Info

Publication number
CN101242101A
CN101242101A CNA2007100373596A CN200710037359A CN101242101A CN 101242101 A CN101242101 A CN 101242101A CN A2007100373596 A CNA2007100373596 A CN A2007100373596A CN 200710037359 A CN200710037359 A CN 200710037359A CN 101242101 A CN101242101 A CN 101242101A
Authority
CN
China
Prior art keywords
cos
harmonic
model
sin
load
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CNA2007100373596A
Other languages
English (en)
Other versions
CN101242101B (zh
Inventor
程浩忠
胡丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Shanghai Maglev Transportation Engineering Technology Research Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Maglev Transportation Engineering Technology Research Center filed Critical Shanghai Maglev Transportation Engineering Technology Research Center
Priority to CN2007100373596A priority Critical patent/CN101242101B/zh
Publication of CN101242101A publication Critical patent/CN101242101A/zh
Application granted granted Critical
Publication of CN101242101B publication Critical patent/CN101242101B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种电力系统负荷谐波稳态模型的建立方法及仿真方法,基于Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)h=1,2,…,N所表达的黑箱模型,将由幅值和相角表示的输入数据和输出数据变换为其复数直角坐标形式,把由新的输入输出决定的新产生的黑箱模型作为模型训练和仿真的核心,并添加相应步骤使新的建模方法在接口上和原方法兼容,可减小模型复杂程度,提高建模过程效率。

Description

电力系统负荷谐波稳态模型的建立方法及仿真方法
技术领域
本发明涉及一种建模方法及其配套的仿真方法,特别涉及一种电力系统负荷谐波稳态黑箱模型的建立方法及其配套的仿真方法。
背景技术
供电网络中理想的电压和电流是一种正弦波形信号,但实际的电压和电流信号往往受各种因素制约无法完全理想地呈现正弦波型。电力系统负荷稳态谐波模型可以表示如下:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N  (1)
式(1)中Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的电压向量;C为负荷的特征参数集合向量,不同的负荷包含不同的特征参数,一般电力电子设备负荷的特征参数是其电路结构参数,含有电弧和铁磁非线性设备的谐波负荷的特征参数为表征其伏安特性的各参数;h从1到N,代表从基波到N次谐波。
式(1)表述的谐波稳态模型只提供了模型的形式,至于输入输出之间的函数关系Fh(h=1,2,…,N)需要对应具体负荷具体求取。实质上,这是一种“黑箱”模型,一般与函数拟合的方法配合使用建立模型。“黑箱”模型通常是指对内部结构和行为不清楚的系统,依据可控因素的输入所引起的可观测因素变化的各种实验数据来确定系统运行规律,从而建立的系统模型。
程浩忠和占勇等人2005年8月发表在《中国电机工程学报》第25卷第16期上的论文《基于广义生长-剪枝径向基函数神经网络的谐波源建模》公开了一种运用广义生长-剪枝径向基函数神经网络建立式(1)所述模型的方法,并通过对脉宽调制电压源逆变可调速传动系统的谐波源建模研究,证明了式(1)所述模型的实用性。
广义生长-剪枝径向基函数神经网络是神经网络的一种。简单的说,人工神经网络可以在不清楚非线性系统内部结构的前提下,通过对实际系统一系列已经测得的输入输出数据对的训练,充分地逼近实际系统输入输出的非线性关系。总之,神经网络是一种非线性系统黑箱建模的有效手段。广义生长-剪枝径向基函数神经网络的建模手段是一种改进的神经网络,该方法在生长和剪枝策略中引入神经元“权”的概念。在训练中,发现某个神经元的“权”比学习精度小,这个神经元将被去掉,且只需要检查与最新输入数据距离最近的神经元的“权”。如果新的输入数据并不需要增加新的神经元,那么只有距离最近的那个神经元的参数被调整。这样,计算量将减少,学习的速度也提高了。也就是说,广义生长-剪枝径向基函数神经网络具有根据实际待拟合系统自适应地改变自身结构和连接方式的能力。
如果在建模过程中所要拟合的非线性函数复杂度越高,一般所需计算量越大。广义生长-剪枝径向基函数神经网络是一种可以随建模系统复杂程度自动改变自身结构的人工神经网络,如果待拟合系统复杂程度较高,最后在给定模型精度要求下训练成熟的神经网络就比较庞大,这样会导致神经网络的训练时间较长,还会导致利用训练成熟的神经网络模型进行仿真的计算量较大。另外,如果使用不能调节自身结构的神经网络建立式(1)所指模型,在同样的初始神经网络、同样的训练数据和同样的训练计算量的条件下,结果是较复杂的待拟合系统模型的训练精度会较低。总之,建模过程的复杂程度和待拟合的非线性系统的复杂度相关,非线性系统越复杂,建模过程的复杂程度越高。
建模过程的效率可以定义为一定的计算量所能达到的模型精度。也就是说,建立相同精度的模型,建模过程中计算量越小则建模效率越高;同样也可以说,消耗相同的计算量,所建模型精度越高则建模效率越高。
式(1)所表述的稳态谐波特性“黑箱”模型的建模过程如图1所示。图中,模型输入和输出数据框标识了模型的输入输出变量。这种建模过程中使用基波和各次谐波电压向量的幅值和相角作为模型输入变量,同时使用负荷吸收基波和各次谐波电流向量的幅值和相角作为模型的输出向量。
图1中标出了现有“黑箱”建模和仿真过程数据形式及其流向。所谓模型仿真是利用建立的模型对现实中系统进行模拟的过程。标有“输入/输出端模型训练”的箭头代表“黑箱”建模过程的数据流向,“黑箱”建模过程实际上是通过分析和逼近未知系统输入/输出实际数据的关系得到对实际系统的模型描述,所以其流向都是从输入/输出数据变量指向待建模型。标有“输入/输出端模型仿真”的箭头代表利用已经建立的模型进行仿真的数据流向,通过把输入数据代入所建模型可以计算出模型的输出数据。
式(1)所表示的模型是频域稳态谐波特性模型的一般形式,其输入输出变量实质上都是复数。因此,所要拟合的非线性函数也是一种复变函数。而SVC(静止无功补偿器)、TCSC(可控串补)、PMW-VSC(电压源型脉宽调制变流器)、STATCOM(静止同步补偿器)、和UPFC(统一潮流控制器)等典型非线性电力电子设备稳态频域模型中经常涉及的运算是复数加减法(更多的是累加)。为了简化讨论,可以列举一个最简单的4输入、2输出复数加减法系统作为分析示例,具体形式如下:
O∠θo=I1∠θi1±I2∠θi2    (2)
假设式(2)在工程上上式可以代表某一实际系统,此系统的输出变量为O和θo,输入变量是I1、θi1、I2和θi2。根据复数加法公式,可得:
O ∠ θ o
= I 1 ∠ θ i 1 ± I 2 ∠ θ i 2 = ( I 1 cos θ i 1 + j I 1 sin θ i 1 ) ± ( I 2 cos θ i 2 + j I 2 sin θ i 2 )
= ( I 1 cos θ i 1 ± I 2 cos θ i 2 ) + j ( I 1 sin θ i 1 ± I 2 sin θ i 2 )
= ( I 1 cos θ i 1 ± I 2 cos θ i 2 ) 2 + ( I 1 sin θ i 1 ± I 2 sin θ i 2 ) 2 ∠ ( arctan I 1 sin θ i 1 ± I 2 sin θ i 2 I 1 cos θ i 1 ± I 2 cos θ i 2 ( + π ) )
由此,输入输出变量之间函数关系为:
O = ( I 1 cos θ i 1 ± I 2 cos θ i 2 ) 2 + ( I 1 sin θ i 1 ± I 2 sin θ i 2 ) 2 - - - ( 3 )
θ o = ( arctan ( I 1 sin θ i 1 ± I 2 sin θ i 2 ) ( I 1 cos θ i 1 ± I 2 cos θ i 2 ) ( + π ) ) - - - ( 4 )
可见:即使是最简单的复数加减法,采用式(1)所述模型时,“黑箱”模型内部需要拟合式(3)和式(4)两个非线性函数,其中涉及正弦、余弦、平方、开方和反正切这些复杂度很高的函数。
综上所述,如果直接使用式(1)表示的模型,在建立电力负荷谐波稳态“黑箱”模型时,由于需要拟合更多复数加减法,模型需要负担拟合更多的复杂非线性函数,最终会造成模型复杂程度较大,建模过程效率较低的问题。
发明内容
本发明要解决的技术问题是提供一种电力系统负荷谐波稳态模型的建立方法,可以提高电力负荷谐波稳态模型建模效率。
为解决上述技术问题,本发明电力系统负荷谐波稳态模型的建立方法,建立如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N  (1)
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据。
为解决上述技术问题,本发明电力系统负荷谐波稳态模型的建立方法提供的另一种技术方案为,建立如式(1)所表示的电力系统负荷谐波稳态黑箱模型,包括预处理输出数据的步骤,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi和Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,...,INcosθNi]作为训练及建立模型的输出数据。
为解决上述技术问题,本发明电力系统负荷谐波稳态模型的建立方法还提供了一种技术方案,即,建立如式(1)所表示的电力系统负荷谐波稳态黑箱模型,包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据;还包括预处理输出数据的步骤,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,...,INcosθNi]作为训练及建立模型的输出数据。
本发明在电力系统负荷谐波稳态黑箱模型建立和训练过程中,将以幅值和相角形式表示的输入数据或输出数据转化为复数直角坐标形式作为训练及建立模型的输入数据或输出数据,可减小模型复杂程度,提高建模过程效率。
本发明所要解决的另一问题是为上述建模方法提供一种配套的电力系统负荷谐波稳态模型的仿真方法,使按上述方法新产生的黑箱模型的输入输出形式兼容原有模型的输出输入接口。
为解决上述问题,本发明电力系统负荷谐波稳态模型的仿真方法,对如下式所表示的电力系统负荷谐波稳态黑箱模型进行仿真:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,θ2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
如果模型建立过程包括预处理输入数据的步骤,即,将所述供电电压中的基波和各次谐波的电压向量V2∠θ1u,V1∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据;则本发明所提供的仿真方法包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为在上述建模过程中新产生的黑箱模型模型的输入数据。
如果模型建立过程包括预处理输出数据的步骤,即,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)  转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,...,INcosθNi]作为训练及建立模型的输出数据;则本发明所提供的仿真方法包括输出数据的反变换的步骤,在按上述建模过程新产生的黑箱模型的输出数据输出前对其预处理,将表达形式为[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]的数据转化为[I1,I2,…IN;θ1i,θ2i,…,θNi]形式的数据作为输出数据输出,其中,ReIh=Ihcosθhi,ImIh=Ihsinθhi  h=1,2,…,N。
如果模型建立过程包括预处理输入数据的步骤,即,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据;还包括预处理输出数据的步骤,即,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,...,INcosθNi]作为训练及建立模型的输出数据;则本发明所提供的仿真方法包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为在上述建模过程中新产生的黑箱模型的输入数据;并包括输出数据的反变换的步骤,在新产生的黑箱模型的输出数据输出前进行处理,将表达形式为[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]的数据转化为[I1,I2,…IN;θ1i,θ2i,…,θNi]形式的数据作为输出数据输出,其中,ReIh=Ihcosθhi,ImIh=Ihsinθhi  h=1,2,…,N。
本发明电力系统负荷谐波稳态模型的仿真方法,在仿真过程中采用与建模过程相配套的步骤,使新产生的黑箱模型的输出数据兼容原有模型的输出和输入接口。
附图说明
图1是现有技术中稳态谐波特性建模方法过程示意图;
图2是本发明模型训练数据流流向示意图;
图3是模型仿真数据流流向示意图;
图4是本发明模型训练的流程示意图;
图5是模型仿真的流程示意图;
图6是作为本发明实施例数据来源的三相桥式六脉冲全控晶闸管整流电路的仿真电路图;
图7是对图6所示电路进行仿真所得到电流和电压的第一种波形示意图;
图8是对图6所示电路进行仿真所得到电流和电压的第二种波形示意图;
图9是本发明实施例中对比现有技术方案和本发明提出的技术方案两者效果的训练流程示意图。
具体实施方式
下面结合附图对本发明作进一步详细的说明。
本发明针对现有模型的缺点,在建立模型时增加输入数据和输出数据的预处理环节,下面分述之。
其一,输入数据的预处理。
对于式(1)中的输入变量,即[V1,V2,…,VN;θ1u,θ2u,…,θNu,C],不直接使用它们作为输入训练及建立“黑箱”模型,而是在此之前加入一个输入数据预处理的步骤,将以幅值和相角形式表示的输入数据中Vh,θhu(h=1,2,…,N)转化为复数的直角坐标表示形式,即:Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),转化按下式进行:
Vh∠θhu=Vh·cosθhu+j·Vh·sinθhu(h=1,2,…,N)。
获得Vh·cosθhu和Vh·sinθhu(h=1,2,…,N)后,最终形成的输入变量形式为[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C],用这组新的输入向量作为“黑箱”模型的输入进行训练。使用这种新的输入向量,可以减少“黑箱”模型建模过程中正弦函数和余弦函数的拟合负担。
其二,输出数据的预处理。
对于式(1)中的输出变量,I1∠θ1i,I2∠θ2i,…,IN∠θNi,即[I1,I2,…IN;θ1i,θ2i,…,θNi],不直接使用它们作为输出训练及建立“黑箱”模型,在此之前加入一个输出数据预处理的步骤,把原输入数据中Ih,θhi(h=1,2,…,N)表示方法转化为复数在直角坐标系下的表示方法,即:Ih·cosθhi和Ih·sinθhi(h=1,2,…,N),转化按下式进行:
Ih∠θhi=Ihcosθhi+jIhsinθhi  h=1,2,…,N。
取上式等式右边的实部和虚部作为输出变量,首先定义符号:ReIh和ImIh为:ReIh=Ihcosθhi,ImIh=Ihsinθhi  h=1,2,…,N,最终形成新的输出向量形式:[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]。用这组新的输出向量作为“黑箱”模型的输出进行训练。使用这种新的输出向量,可以减少“黑箱”模型建模过程中的平方、开方和反正切函数的拟合负担。
上述对输入数据的预处理和对输出数据的预处理既可以分开使用,也可以同时使用,由此形成三种不同的技术方案,即:方案一,对输入数据进行预处理,而不对输出数据进行预处理;方案二,对输出数据进行预处理,而不对输入数据进行预处理;方案三,对输入数据进行预处理,同时也对输出数据进行预处理。上述三种方案基于同一构思,为简明起见,下面仅以技术方案三为例描述按本发明所提供的方法进行模型训练和仿真的过程。
使用本发明提供的方法建立电力系统负荷谐波稳态模型时,模型训练数据流和模型仿真数据流可分别用图2和图3表示。图2在图1所示现有技术中的模型的基础上结合本发明所提供的方法把现有技术中的“原有模型”内部分为“输入数据的预处理”、“新产生的黑箱模型”、“输出数据的预处理”三部分;图3在图1所示现有技术中的模型的基础上结合本发明把现有技术中的模型内部分为“输入数据的预处理”、“新产生的黑箱模型”、  以及“仿真时数据的反变换”三部分。
其中“新产生的黑箱模型”像原模型一样也是一个“黑箱”模型,需要结合某种“黑箱”建模手段才能完成建模过程,如按照前述《基于广义生长-剪枝径向基函数神经网络的谐波源建模》一文中描述的广义生长-剪枝径向基函数神经网络手段。而“新产生的黑箱模型”的输入数据和输出数据就是按照本发明所提供的方法预处理后的输入数据和输出数据,即,输入数据为:[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…VN·cosθNu,VN·sinθNu,C],输出数据为[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]。
输入数据的预处理和输出数据的预处理在训练黑箱模型时需要用到。图2和图3中出现一个“新产生的黑箱模型”,是由于输入、输出变量的预处理过程使“原有模型”的输入、输出变量发生改变,在新的输入、输出变量的条件下自然地产生出另一个有别于“原有模型”的新模型。本发明的根本目的就是使这个“新产生的黑箱模型”所需拟合的非线性关系比“原有模型”简单,因此可以通过间接地训练、仿真“新产生的黑箱模型”,代替直接操作“原有模型”,从而达到提高建模过程效率的结果。
仿真过程要与建模过程相对应,如果在训练时同时包括“输入数据的预处理”和“输出数据的预处理”步骤,则在仿真时必须同时包括“输入数据的预处理”和“仿真时的反变换”步骤;如果在训练时只包括“输入数据的预处理”步骤,则在仿真时必须只包括“输入数据的预处理”步骤;如果在训练时只包括“输出数据的预处理”步骤,则在仿真时必须只包括“仿真时的反变换”步骤。只有符合以上三点,才能在输入输出接口上保持和原有模型相同。
图3中“输入数据的预处理”和“仿真时反变换”这两部分,是运用已经建立的“新产生的黑箱模型”进行仿真时所需的数据处理步骤。其中“输入数据的预处理”部分数据处理方式和模型训练时相同,因此图中使用同一名称的模块表示。而“仿真时反变换”这部分,是输出数据预处理步骤的反变换。即把[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]变换为[I1,I2,…IN;θ1i,θ2i,…,θNi]的过程,具体处理过程如下:
&theta; hi = &pi; 2 Re I h = 0 , Im I h > 0 sgn ( Re I h ) &CenterDot; arctan ( Im I h Re I h ) Re I h &NotEqual; 0 - &pi; 2 Re I h = 0 , Im I h < 0 h = 1,2 , &CenterDot; &CenterDot; &CenterDot; , N
I h = ReI h cos &theta; hi &theta; hi = 0 ImI h sin &theta; hi &theta; hi &NotEqual; 0
其中函数sgn()是符号函数,其取值如下:
sgn ( x ) = 1 x > 0 1 x = 0 - 1 x < 0
经上述处理后,可以在模型仿真时由“新产生的黑箱模型”的输出数据[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]得到“原有模型”的输出数据[I1,I2,…IN;θ1i,θ2i,…,θNi]。“仿真时反变换”是进行模型仿真时需要的数据处理步骤,图3表示了模型仿真时数据流的流向。
综上所述,本发明的实现可用图4表示的模型训练流程图和图5表示的模型仿真流程图表示。
本发明所提供的建模方法可以提高电力负荷谐波稳态模型建模过程的效率。
现有技术中电力负荷谐波稳态模型建模过程即使在最简单的4输入、2输出情况下,模型内部的复数加减法系统也要拟合正弦、余弦、平方、开方和反正切这些复杂度很高的函数,导致建模时效率比较低下。
对于同一简单示例,如果通过使用本发明所提供的新的建模方法,在建模过程中进行输入数据和输出数据的预处理过程,则输入数据由幅值相角形式变换为[I1cosθi1,I1sinθi1,I1cosθi2,I2sinθi2],输出数据变换为[Ocosθo,Osinθo]。这时建立下式所示的模型:O∠θo=I1∠θi1±I2∠θi2
根据复数加法公式,可得:
Ocosθo±jOsinθo
=I1∠θi1±I2∠θi2=(I1cosθi1+jI1sinθi1)±(I2cosθi2+jI2sinθi2)
=(I1cosθi1±I2cosθi2)+j(I1sinθi1±I2sinθi2)
由此得出,输入数据和输出数据之间的函数关系为:
Ocosθo=I1cosθi1±I2cosθi2
Osinθo=I1sinθi1±I2sinθi2
可见,经过输入输出数据预处理,“新产生的黑箱模型”输入输出变量之间需要拟合的关系只是单纯的加法和减法。比原有建模方法需要拟合正弦、余弦、平方、开方和反正切等运算函数要简单的多,所以新的建模方法的建模效率较高。
以上只是讨论最简单的复数加减法模型的建立过程,已可发现本发明所提供的方法可有效提高建模效率。如果考虑典型非线性电力电子设备稳态频域模型中大量的复数加减法(更多的是以累加的形式出现),则使用本发明所提供的方法可以较大地提高建模的效率。建模效率的提高具体表现为可以减少模型训练所需的计算量,或者可以在相同训练计算量的前提下提高模型的精度。
下面以图6所示的三相桥式六脉冲全控晶闸管整流电路为例来对本发明的建模效率进行说明,输入数据和输出数据来源于用Simulink对图6所示电路的仿真。整流电路是一种典型的非线性元件电路,是电力电子器件在实际电网中最普遍的应用之一。图6所示整流电路的作用是通过对6个晶闸管导通角和导通次序的控制,把三相交流电变换为单相直流电,或进行反向的变换。
图6中,整流器导通角固定在30度。由于仿真电路是三相对称的,本实施例只取整流器A相电流和A相电压的实时值ia(t)和ua(t)进行分析。具体的分析过程是,不断对前一周期(20ms)时间内的信号使用傅立叶变换,分析其1-30次谐波(包括基波)的幅值和相角。为尽量使所有暂态过程衰减,仿真时间设置为0-1秒(50个周期)。
通过对结果的分析,每次仿真所得第1-30次谐波分量(包括基波)的傅立叶分析结果表现为两种类型。第一种,经过一段时间,暂态分量消失后,电压和电流幅值保持平稳,或有微量波动,如图6所示,可以认为,此次谐波分量已经到达稳态。第二种,谐波分量的幅值很小,相角不断随机波动,如图7所示,可以认为,此次谐波稳态分量为零。
通过对100次仿真结果的分析综合,在1-30次谐波(包括基波)范围内,此电路包括1、3、5、7、11、13、17、19、23、25、29,共11次谐波稳态分量(包括基波)。因此,确定此电路根据式(1)的“黑箱”模型的输入为:[V1,V3,V5,V7,…,V29;θ1u,θ3u,θ5u,θ7u,…,θ29u],共22个变量;输出为:[I1,I3,I5,I7,…,I29;θ1i,θ3i,θ5i,θ7i,…,θ29i],共22个变量。
图6中的基波及3、5、7次谐波发生器,可以产生基波线电压110kV,3次谐波电压占基波的3~30%,5次谐波电压占基波的2~20%,7次谐波电压占基波的1~10%。其中3、5、7次谐波含量的范围是由随机变量实现的,每次仿真时其含量总是在其范围内随机产生。图6所示电路总共仿真了100次,保证可以采集100组不同的电流ia(t)和电压ua(t)波形。通过它们进行傅里叶分析,可以得到100组不同的输入[V1,V3,V5,V7,…,V29;θ1u,θ3u,θ5u,θ7u,…,θ29u]和与之对应的输出[I1,I3,I5,I7,…,I29;θ1i,θ3i,θ5i,θ7i,…,θ29i]。这100组输入数据和输出数据就是图6中三相桥式6脉冲全控晶闸管整流电路稳态谐波特性“黑箱”模型的训练数据,输入数据和输出数据的基波和各次谐波幅值和相角数值见附表一和附表二。
用上述100组输入输出数据,可以训练一个传统的稳态谐波特性模型,也就是图1中的“现有模型”。根据本发明所提供的方法对原始输入数据和输出数据分别进行预处理,可以得到100组新的输入数据[V1·cosθ1u,V1·sinθ1u,V3·cosθ3u,V3·sinθ3u,…V29·cosθ29u,V29·sinθ29u]和100组输出数据[I1·cosθ1i,I1·sinθ1i,I3·cosθ3i,I3·sinθ3i,…I29·cosθ29i,V29·sinθ29i],用这些数据可以训练一个经简化的稳态谐波特性模型,也就是图2中的“新产生的黑箱模型”。
最后,分别使用现有技术和本发明所提供的方法训练同一初始状态和相同结构的神经网络模型,训练使用同一种训练函数,即带动量因子的负梯度法;这两个神经网络都被训练3000个周期,所有100组输入输出数据顺序地参与训练1遍称为1个周期。在此之后,观察现有技术和本发明所提供的方法对训练数据的均方差误差(MeanSquared Error),较大的误差反应了网络训练难度较大,也就是网络的复杂度较大,建模效率较低。
这些神经网络都是22个输入节点、45个隐层节点和22个输出节点的三层反向传播网络。而且,这种训练总共进行了100次,每次神经网络的初始状态都随机产生,通过对100次结果的对比分析,可以得到比较可靠的判断。
仿真采用软件Matlab进行。图9表示以现有技术和通过本发明所提供的方法对神经网络的训练流程。图9中所得MSE(i)1、MSE(i)2是分别基于现有建模方法和本发明建模方法对相同结构、相同初始条件的神经网络使用带动量因子的负梯度法各训练3000周期后,得到的两个神经网络的均方误差。R(i)是一个比率,可以代表新的建模方法所得网络的均方差比现有建模方法的均方差减少的百分比。最终得到100组MSE(i)1、MSE(i)2和R(i),其值请详见附表二。从附表二中可以看出,R(i)的最小值为6.0102%,最大值为47.411%,平均值为30.221%。这说明本实施例中100个随机初值的神经网络建模过程中,使用本发明提出的建模方法比现有建模方法在同样多的训练周期内得到均方根误差均有减少,至少减少6.0102%,平均可以减少30.221%。
附表三是基于现有建模方法和按本发明中只对输入数据进行预处理的建模方法对相同结构、相同初始条件的神经网络使用带动量因子的负梯度法各训练3000周期后,得到的两个神经网络的均方误差表。从表中可以看出,使用本发明提出的建模方法比现有建模方法在同样多的训练周期内得到均方根误差均有减少,平均可以减少4.2689%。
附表四是基于现有建模方法和按本发明中只对输出数据进行预处理的建模方法对相同结构、相同初始条件的神经网络使用带动量因子的负梯度法各训练3000周期后,得到的两个神经网络的均方误差表。从表中可以看出,使用本发明提出的建模方法比现有建模方法在同样多的训练周期内得到均方根误差均有减少,平均可以减少91.649%。
这些数据说明本发明提出的建模方法可以较大程度地提高电力负荷谐波稳态特性模型的建模效率。
附表一是图6中三相桥式6脉冲全控晶闸管整流电路稳态谐波特性“黑箱”模型的训练输入数据。
附表二是是图6中三相桥式6脉冲全控晶闸管整流电路稳态谐波特性“黑箱”模型的训练输出数据。
附表三是基于现有建模方法和按本发明中对输入数据和输出数据同时进行预处理的建模方法对相同结构、相同初始条件的神经网络使用带动量因子的负梯度法各训练3000周期后,得到的两个神经网络的均方误差表。
附表四是基于现有建模方法和按本发明中只对输入数据进行预处理的建模方法对相同结构、相同初始条件的神经网络使用带动量因子的负梯度法各训练3000周期后,得到的两个神经网络的均方误差表。
附表五是基于现有建模方法和按本发明中只对输出数据进行预处理的建模方法对相同结构、相同初始条件的神经网络使用带动量因子的负梯度法各训练3000周期后,得到的两个神经网络的均方误差表。
附表一(电压单位V,角度单位rad)
  V1   θ1u   V3   θ3u   V5   θ5u   V7   θ7u   V11   θ11u   V13   θ13u
  7.1843E+04   -0.087   2.5713E+04   1.453   2.0842E+04   -2.513   2.0401E+04   3.021   7.5116E+03   -2.808   5.0163E+03   2.147
  7.2110E+04   -0.084   1.3732E+04   0.115   2.1503E+04   -2.756   1.8680E+04   2.910   8.8605E+03   -2.865   6.7717E+03   2.077
  7.1729E+04   -0.082   2.5059E+04   -1.645   1.4824E+04   -2.357   2.0998E+04   2.747   8.7037E+03   -3.014   6.4192E+03   1.891
  7.1692E+04   -0.082   1.2639E+04   -0.669   1.3870E+04   -2.232   1.9566E+04   2.633   9.2416E+03   -3.037   7.0825E+03   1.862
  7.1193E+04   -0.103   6.0592E+03   1.274   1.8912E+04   -1.927   1.5395E+04   -2.901   7.3170E+03   -2.410   4.5095E+03   2.564
  7.0516E+04   -0.126   3.0259E+03   -1.580   2.1188E+04   -1.095   9.9928E+03   -2.253   6.0067E+03   -2.074   3.5841E+03   2.620
  7.0993E+04   -0.105   2.3197E+04   -2.984   1.9516E+04   -1.814   1.7595E+04   -2.636   6.0018E+03   -2.485   3.1073E+03   2.299
  7.1941E+04   -0.089   1.9224E+04   2.386   2.4142E+04   -2.541   1.8825E+04   3.098   7.2663E+03   -2.696   4.7641E+03   2.294
  7.1033E+04   -0.108   1.0094E+04   1.191   1.9647E+04   -1.761   1.4631E+04   -2.734   6.6634E+03   -2.316   3.7310E+03   2.614
  7.0573E+04   -0.121   6.3472E+03   -1.898   2.1676E+04   -1.312   1.2079E+04   -2.076   5.7576E+03   -2.282   3.6040E+03   2.324
  7.1365E+04   -0.100   1.4723E+04   -0.632   2.7856E+04   -2.021   1.6380E+04   -2.720   3.5245E+03   -2.680   1.0589E+03   1.279
  7.1472E+04   -0.094   1.0979E+04   1.822   1.9207E+04   -2.206   1.8970E+04   -3.084   7.1459E+03   -2.657   4.4342E+03   2.294
  7.1013E+04   -0.108   2.3022E+04   -2.714   2.2169E+04   -1.745   1.4967E+04   -2.605   5.3547E+03   -2.396   2.4847E+03   2.311
  7.0324E+04   -0.126   1.9463E+04   -2.377   2.4278E+04   -0.904   8.4925E+03   -2.323   6.2469E+03   -2.295   4.5403E+03   2.381
  7.2111E+04   -0.070   6.8761E+03   -0.129   1.0568E+04   -2.587   1.8702E+04   2.257   1.0629E+04   2.885   9.0105E+03   1.433
  7.2488E+04   -0.059   5.9948E+03   0.073   7.7702E+03   2.374   1.5384E+04   2.190   1.1981E+04   2.493   1.0875E+04   0.966
  7.0659E+04   -0.124   9.5653E+03   2.943   1.6611E+04   -1.314   9.3662E+03   -2.302   8.2830E+03   -1.876   5.3810E+03   -3.128
  7.2083E+04   -0.070   1.5196E+04   2.098   8.3944E+03   -2.770   1.6574E+04   2.294   1.1244E+04   2.847   9.7687E+03   1.387
  7.1773E+04   -0.081   1.5550E+04   -2.259   1.4233E+04   -2.405   1.9751E+04   2.699   9.2441E+03   -3.026   7.1097E+03   1.876
  7.2086E+04   -0.057   1.3873E+04   -2.715   4.5379E+03   -1.126   1.3155E+04   1.864   1.2317E+04   2.334   1.1158E+04   0.775
  7.0182E+04   -0.133   1.2680E+04   1.925   2.2863E+04   -0.666   7.1755E+03   -2.139   74037E+03   -2.006   5.3095E+03   2.829
  7.1574E+04   -0.095   2.6709E+04   -1.327   1.9161E+04   -2.305   1.5634E+04   -3.129   8.4707E+03   -2.573   6.0411E+03   2.411
  7.2156E+04   -0.080   1.0437E+04   -0.253   1.9029E+04   -2.835   1.9706E+04   2.772   9.2951E+03   -3.015   7.3281E+03   1.897
  7.2646E+04   -0.058   1.3360E+04   -0.419   8.1383E+03   2.796   1.8489E+04   2.008   1.1446E+04   2.474   1.0260E+04   0.946
  7.1353E+04   -0.103   5.9180E+03   1.304   2.4683E+04   -1.996   1.4837E+04   -2.814   5.3454E+03   -2.386   2.3530E+03   2.533
  7.1208E+04   -0.096   1.3646E+04   0.276   1.3768E404   -1.944   1.4078E+04   3.086   9.6519E+03   -2.596   7.3124E+03   2.372
  7.1066E+04   -0.101   1.9265E+04   0.583   1.6900E+04   -1.860   1.7649E+04   -2.856   7.2671E+03   -2.500   4.4219E+03   2.428
  7.0433E+04   -0.115   3.1097E+03   1.183   1.5662E+04   -0.911   1.0366E+04   -2.851   9.8319E+03   -2.158   7.1397E+03   2.854
  7.1388E+04   -0.092   2.0103E+04   -1.931   1.5463E+04   -2.102   1.7454E+04   3.020   8.6680E+03   -2.711   6.2029E+03   2.239
  7.1043E+04   -0.109   6.4085E+03   -2.037   2.6875E+04   -1.694   1.2931E+04   -2.570   3.6877E+03   -2.605   1.8544E+03   1.424
  7.0513E+04   -0.121   5.6161E+03   2.828   2.5216E+04   -1.128   9.7243E+03   -2.483   4.9904E+03   -2.478   3.5386E+03   1.978
  7.0480E+04   -0.122   2.3671E+04   1.459   2.5889E+04   -1.064   9.4541E+03   -2.559   5.1141E+03   -2.530   3.9042E+03   1.966
  7.2615E+04   -0.057   3.8984E+03   0.493   6.6833E+03   2.873   1.7994E+04   1.947   1.1539E+04   2.433   1.0366E+04   0.898
  7.1850E+04   -0.08   6.8309E+03   -0.037   1.3034E+04   -2.681   1.6353E+04   2.726   1.0478E+04   -3.033   8.6799E+03   1.865
  7.0739E+04   -0.116   1.2225E+04   -2.565   1.5452E+04   -1.364   9.3816E+03   -2.732   9.2764E+03   -2.084   6.6039E+03   2.949
  7.0737E+04   -0.117   2.5386E+04   1.660   1.8760E+04   -1.400   1.1918E+04   -2.529   6.9562E+03   -2.090   4.0029E+03   2.836
  7.0241E+04   -0.130   2.6110E+04   -2.105   2.2780E+04   -0.678   1.0304E+04   -2.527   6.3914E+03   -2.098   4.3210E+03   2.613
  7.1224E+04   -0.107   1.3123E+04   -0.690   2.7212E+04   -1.826   1.3162E+04   -2.691   3.8362E+03   -2.443   1.1993E+03   1.787
  7.0649E+04   -0.120   6.5224E+03   1.197   2.2307E+04   -1.309   1.1344E+04   -2.312   5.4020E+03   -2.257   3.0808E+03   2.316
  7.1952E+04   -0.08   1.3889E+04   2.875   1.6566E+04   -2.629   2.0835E+04   2.758   8.8951E+03   -3.027   6.7519E+03   1.880
  7.1264E+04   -0.099   9.8704E+03   0.308   2.6108E+04   -1.992   1.7514E+04   -2.685   3.8775E+03   -2.744   1.4621E+03   1.386
  7.1021E+04   -0.110   1.2405E+04   1.249   2.5893E+04   -1.664   1.2696E+04   -2.592   4.0445E+03   -2.479   1.7807E+03   1.734
  7.0336E+04   -0.120   1.2890E+04   -1.544   2.5379E+04   -1.057   8.5507E+03   -2.125   6.4450E+03   -2.477   4.8306E+03   2.193
  7.0630E+04   -0.124   1.7730E+04   -1.688   1.8897E+04   -1.283   1.1171E+04   -2.201   6.6979E+03   -1.989   3.8819E+03   2.869
  V1   θ1u   V3   θ3u   V5   θ5u   V7   θ7u   V11   θ11u   V13   θ13u
  7.2538E+04   -0.070   8.3979E+03   -2.834   1.7244E+04   2.980   1.8864E+04   2.487   1.0730E+04   2.932   9.2917E+03   1.493
  7.1824E+04   -0.084   1.7915E+04   -1.891   1.6076E+04   -2.540   1.9026E+04   2.820   9.1664E+03   -2.936   7.0180E+03   1.985
  7.1540E+04   -0.095   1.1702E+04   2.672   2.3108E+04   -2.228   1.8136E+04   -2.979   6.0986E+03   -2.580   3.2196E+03   2.388
  7.0894E+04   -0.107   2.1521E+04   3.074   1.9416E+04   -1.752   1.7751E+04   -2.515   5.8695E+03   -2.492   3.0802E+03   2.216
  7.1023E+04   -0.106   1.9749E+04   3.047   1.8420E+04   -1.783   1.5372E+04   -2.761   7.0443E+03   -2.356   4.1449E+03   2.588
  7.0276E+04   -0.134   1.6464E+04   -1.108   2.1855E+04   -0.797   7.2046E+03   -1.843   7.1818E+03   -1.957   5.0848E+03   2.867
  7.1248E+04   -0.104   1.1435E+04   -2.837   1.9639E+04   -1.951   1.3946E+04   -2.930   7.5616E+03   -2.353   4.8085E+03   2.647
  7.1014E+04   -0.109   1.9702E+04   -0.236   2.7602E+04   -1.651   1.2237E+04   -2.565   3.4796E+03   -2.673   2.0533E+03   1.284
  7.0308E+04   -0.116   2.2255E+04   0.525   2.7121E+04   -1.049   8.3663E+03   -2.431   6.2921E+03   -2.624   4.7336E+03   2.005
  7.2186E+04   -0.079   2.3907E+04   -3.057   1.9742E+04   -2.820   2.1262E+04   2.772   8.7748E+03   -3.035   6.7042E+03   1.877
  7.2473E+04   -0.071   1.6108E+04   1.836   1.6482E+04   3.042   1.9028E+04   2.491   1.0644E+04   2.948   9.1575E+03   1.512
  7.1447E+04   -0.097   1.1361E+04   3.132   2.0981E+04   -2.165   1.7776E+04   -2.990   6.6942E+03   -2.562   3.8902E+03   2.400
  7.0702E+04   -0.113   2.4468E+04   -1.540   2.7591E+04   -1.403   1.0812E+04   -2.330   4.4468E+03   -2.890   3.8180E+03   1.472
  7.0584E+04   -0.105   9.2926E+03   -2.348   1.3814E+04   -1.062   1.2647E+04   -3.046   1.0140E+04   -2.431   7.6620E+03   2.548
  7.1641E+04   -0.100   1.2649E+04   2.981   2.8101E+04   -2.181   1.4616E+04   -2.893   5.1528E+03   -2.368   2.2154E+03   2.701
  7.0322E+04   -0.126   1.7105E+04   1.013   2.4578E+04   -0.906   8.2924E+03   -2.303   6.4058E+03   -2.327   4.7667E+03   2.368
  7.1881E+04   -0.063   2.2378E+04   -2.450   3.5995E+03   -1.487   1.3471E+04   2.130   1.2182E+04   2.577   1.0894E+04   1.062
  7.2168E+04   -0.082   6.7264E+03   -1.054   2.2273E+04   -2.764   2.1061E+04   2.884   8.1409E+03   -2.930   5.9175E+03   2.009
  7.1260E+04   -0.099   2.3800E+04   -1.449   2.2215E+04   -2.027   1.8772E+04   -2.790   5.2583E+03   -2.625   2.3400E+03   2.114
  7.1412E+04   -0.094   8.6865E+03   2.218   1.6748E+04   -2.151   1.7152E+04   3.110   8.3663E+03   -2.646   5.8392E+03   2.315
  7.1253E+04   -0.103   1.7512E+04   0.447   2.0453E+04   -1.956   1.4551E+04   -2.895   7.0689E+03   -2.364   4.2341E+03   2.624
  7.1767E+04   -0.088   8.6700E+03   -2.593   1.8819E+04   -2.482   1.8118E+04   3.000   8.5429E+03   -2.764   6.2212E+03   2.190
  7.1618E+04   -0.089   1.0984E+04   2.525   1.5269E+04   -2.387   1.6352E+04   2.941   9.5271E+03   -2.787   7.3714E+03   2.155
  7.1312E+04   -0.084   1.4881E+04   -1.746   9.7762E+03   -1.763   1.2212E+04   2.585   1.1338E+04   -3.008   9.5251E+03   1.887
  7.2193E+04   -0.071   3.0648E+03   -2.113   1.1401E+04   3.038   1.5826E+04   2.512   1.1402E+04   2.925   1.0006E+04   1.479
  7.0565E+04   -0.124   2.1724E+04   -0.086   2.2121E+04   -1.163   1.0344E+04   -2.289   5.6144E+03   -2.186   3.3995E+03   2.424
  7.1546E+04   -0.092   2.2767E+04   1.045   1.9464E+04   -2.271   1.9001E+04   -3.126   7.3101E+03   -2.679   4.6496E+03   2.278
  7.1386E+04   -0.093   1.3869E+04   2.798   1.5300E+04   -2.104   1.6892E+04   3.017   8.8555E+03   -2.699   6.4264E+03   2.254
  7.1797E+04   -0.099   2.3358E+04   -1.512   2.7415E+04   -2.315   1.3248E+04   -2.974   6.9095E+03   -2.335   4.2718E+03   2.735
  7.0970E+04   -0.109   2.2488E+04   -0.144   2.0090E+04   -1.725   1.5279E+04   -2.617   6.1215E+03   -2.347   3.2028E+03   2.489
  7.2462E+04   -0.060   1.9988E+04   1.433   6.6749E+03   -2.848   1.7798E+04   1.902   1.1418E+04   2.503   1.0159E+04   0.980
  7.2390E+04   -0.068   3.8857E+03   -1.538   1.3332E+04   2.821   1.5844E+04   2.476   1.1523E+04   2.847   1.0233E+04   1.388
  7.0677E+04   -0.118   5.8759E+03   -0.370   2.5990E+04   -1.280   1.0215E+04   -2.554   4.3385E+03   2.554   3.0086E+03   1.751
  7.2037E+04   -0.076   4.6594E+03   -0.938   1.2647E+04   -2.914   1.6691E+04   2.607   1.0811E+04   3.100   9.1840E+03   1.687
  7.0672E+04   -0.117   1.1739E+04   0.462   1.5312E+04   -1.268   8.6326E+03   -2726   9.5152E+03   -2.051   6.8475E+03   2.988
  7.2333E+04   -0.075   1.5329E+04   0.736   1.7775E+04   -3.017   2.0239E+04   2.615   9.8559E+03   3.097   8.1120E+03   1.691
  7.0710E+04   -0.118   1.7207E+04   -0.318   2.0024E+04   -1.312   1.1697E+04   -2.614   6.4182E+03   -2.113   35143E+03   2.743
  7.1950E+04   -0.074   8.3687E+03   -0.422   1.2686E+04   -2.305   1.9323E+04   2.299   1.0024E+04   3.009   8.1846E+03   1.582
  7.0356E+04   -0.129   1.4791E+04   1.823   2.3009E+04   -0.928   8.3382E+03   -2.032   6.4422E+03   -2186   4.6633E+03   2.531
  7.1745E+04   -0.080   5.8612E+03   0.600   1.3473E+04   -2.188   1.8381E+04   2.472   9.7582E+03   -3.102   7.7552E+03   1.785
  7.1572E+04   0.099   3.1177E+03   -1.137   2.3231E+04   -2.216   1.5125E+04   -2.985   7.0170E+03   -2.432   4.3031E+03   2.585
  7.2418E+04   -0.062   9.0780E+03   -1.548   7.6074E+03   2.923   1.6946E+04   2.146   1.1567E+04   2.601   1.0327E+04   1.096
  7.2297E+04   -0.061   1.8691E+04   1.152   4.7546E+03   3.088   1.5708E+04   2.091   1.1810E+04   2.549   1.0581E+04   1.032
  7.0314E+04   -0.132   1.9338E+04   -2.022   2.2706E+04   -0.677   9.5116E+03   -2.465   6.5862E+03   -2.059   4.7052E+03   2.694
  7.0908E+04   -0.089   1.0123E+04   -0.466   9.9267E+03   -1.086   1.3429E+04   2.859   1.1497E+04   -2.953   9.6301E+03   1.947
  7.1033E+04   -0.108   8.1983E+03   -3.040   2.3394E+04   -1.755   1.4842E+04   -2.593   4.8871E+03   -2.448   2.1289E+03   2.122
  V1   θ1u   V3   θ3u   V5   θ5u   V7   θ7u   V11   θ11u   V13   θ13u
  7.2133E+04   -0.071   4.4827E+03   1.216   1.0817E+04   -2.758   1.8645E+04   2.345   1.0659E+04   2.914   9.0520E+03   1.468
  7.0427E+04   -0.129   1.8171E+04   -0.085   2.1526E+04   -1.016   9.1951E+03   -1.933   6.3471E+03   -2.079   4.2902E+03   2.638
  7.2259E+04   -0.059   5.4351E+03   -1.179   4.1131E+03   2.037   1.4358E+04   2.236   1.2247E+04   2.461   1.1119E+04   0.926
  7.1257E+04   -0.101   2.0217E+04   -2.189   2.7127E+04   -1.949   1.6412E+04   -2.654   3.5686E+03   -2.749   1.4642E+03   1.191
  7.1050E+04   -0.101   1.2894E+04   1.337   1.4429E+04   -1.667   7.8852E+03   2.829   1.0760E+04   -2.468   8.6135E+03   2.523
  7.0654E+04   -0.121   1.4005E+04   -0.541   1.9784E+04   -1.350   1.2186E+04   -2.191   6.2173E+03   -2.109   3.5464E+03   2.650
  7.1450E+04   -0.098   1.4532E+04   -0.046   1.9456E+04   -2.174   1.5215E+04   -3.049   8.0441E+03   -2.491   5.4609E+03   2.503
  7.0424E+04   -0.126   1.2905E+04   -2.133   2.2806E+04   -0.934   1.0047E+04   -2.493   5.7241E+03   -2.184   3.6806E+03   2.431
  7.0873E+04   -0.111   1.6415E+04   -0.192   1.7009E+04   -1.557   1.2684E+04   -2.790   8.0193E+03   -2.220   5.1969E+03   2.774
  7.0667E+04   -0.119   7.3415E+03   -1.785   1.8686E+04   -1.265   1.1524E+04   -2.606   7.0798E+03   -2.068   4.1583E+03   2.859
附表一(续)
  V17   θ17u   V19   θ19u   V23   θ23u   V25   θ25u   V29   θ29
  4.5314E+02   -3.023   1.4568E+03   -2.207   3.5228E+03   -1.633   3.8192E+03   -2.977   3.0142E+03   -2.528
  2.4319E+03   2.683   8.5825E+02   1.938   2.3666E+03   -2.020   3.1685E+03   2.991   3.4958E+03   -2.825
  1.9257E+03   2.289   1.6927E+02   2.010   2.6571E+03   -1.984   3.3616E+03   2.907   3.5537E+03   -3.005
  2.6996E+03   2.242   8.4739E+02   0.992   2.1043E+03   -2.051   2.9790E+03   2.843   3.5796E+03   -3.070
  7.9680E+02   1.255   2.1189E+03   -0.837   3.8324E+03   -0.374   3.7211E+03   -1.685   2.4518E+03   -1.176
  3.8214E+03   2.265   3.8805E+03   0.852   3.2396E+03   1.335   2.0759E+03   -0.086   1.8403E+03   -0.607
  2.6082E+03   1.030   3.5196E+03   -0.531   4.0247E+03   -0.080   3.1881E+03   -1.408   1.2249E+03   -1.364
  5.7994E+02   -2.160   1.7569E+03   -2.163   3.6117E+03   -1.514   3.8321E+03   -2.827   2.8063E+03   -2.317
  1.7045E+03   1.188   2.8705E+03   -0.479   3.9754E+03   0.021   3.4737E+03   -1.274   1.7246E+03   -0.869
  4.1419E+03   2.045   4.0541E+03   0.638   3.1081E+03   1.045   1.9074E+03   -0.465   2.0454E+03   -1.050
  3.9529E+03   0.372   4.4301E+03   -0.993   3.2756E+03   -0.517   1.8530E+03   -1.882   1.1870E+03   2.507
  4.7913E+02   0.473   2.0509E+03   -1.482   3.8207E+03   -1.042   3.8228E+03   -2.383   2.5945E+03   -1.925
  3.0004E+03   1.056   3.7989E+03   -0.398   3.7894E+03   0.114   2.7432E+03   -1.201   8.3586E+02   -1.695
  4.5223E+03   2.369   3.7529E+03   0.948   2.6772E+03   1.135   2.0822E+03   -0.556   2.7246E+03   -0.605
  5.3574E+03   1.682   3.6948E+03   0.254   6.1073E+02   1.040   8.5158E+02   1.578   2.6098E+03   2.118
  8.2994E+03   1.053   7.0221E+03   -0.473   4.2353E+03   -0.379   3.0159E+03   -1.884   6.3745E+02   -1.546
  1.7994E+03   2.831   2.1269E+03   0.751   3.7566E+03   1.149   3.5466E+03   -0.112   2.5367E+03   0.481
  6.4221E+03   1.608   4.8381E+03   0.156   1.6795E+03   0.451   4.5424E+02   -0.533   1.8042E+03   2.044
  2.7504E+03   2.274   9.1420E+02   1.057   2.0576E+03   -2.075   2.9441E+03   2.835   3.5763E+03   -3.068
  8.6581E+03   0.798   7.3236E+03   -0.762   4.5784E+03   -0.757   3.3049E+03   -2.308   8.7983E+02   -2.398
  4.2004E+03   2.876   3.7245E+03   1.425   3.0399E+03   1.671   2.4774E+03   0.155   2.5136E+03   0.232
  1.3626E+03   2.977   5.0456E+02   -1.776   3.0228E+03   -1.047   3.5628E+03   -2.360   3.3327E+03   -1.832
  31223E+03   2.377   1.4093E+03   1.239   1.8146E+03   -2.389   2.7397E+03   2.643   3.4828E+03   3.075
  7.4276E+03   1.033   6.0818E+03   -0.490   3.1734E+03   -0.359   1.9592E+03   -1.811   5.8565E+02   0.372
  2.3777E+03   0.465   3.5425E+03   -0.913   3.8760E+03   -0.350   3.0772E+03   -1.637   7.5420E+02   -1.381
  2.8168E+03   2.793   8.9462E+02   1.277   2.2174E+03   -0.818   3.0111E+03   -2.218   3.6340E+03   -1.733
  1.2566E+03   1.254   2.3598E+03   -0.765   3.9509E+03   -0.386   3.7292E+03   -1.723   2.4078E+03   -1.291
  2.7819E+03   3.129   1.4780E+03   1.124   3.0634E+03   0.616   3.4041E+03   -0.767   3.5713E+03   -0.214
  1.5795E+03   2.579   3.7113E+02   -1.222   2.9718E+03   -1.172   3.5177E+03   -2.543   3.4567E+03   -2.086
  4.2251E+03   1.004   4.3365E+03   -0.348   2.8004E+03   0.105   1.3357E+03   -1.423   1.8462E+03   -2.792
  4.5517E+03   1.949   3.9572E+03   0.583   2.4186E+03   0.769   1.5584E+03   -1.092   2.6459E+03   -1.281
  4.6819E+03   1.967   3.8294E+03   0.581   2.3087E+03   0.698   1.7761E+03   -1.183   2.8029E+03   -1.259
  75916E+03   0.968   6.2557E+03   -0.564   3.3747E+03   -0.463   2.1562E+03   -1.937   4.2757E+02   0.045
  V17   θ17u   V19   θ19u   V23   θ23u   V25   θ25u   V29   θ29
  4.5314E+02   -3.023   1.4568E+03   -2.207   3.5228E+03   -1.633   3.8192E+03   -2.977   3.0142E+03   -2.528
  4.7448E+03   2.242   2.9762E+03   0.882   3.6898E+02   -2.713   1.4924E+03   2.722   3.0582E+03   -3.116
  2.0629E+03   -3.067   1.0965E+03   0.810   3.1388E+03   0.551   3.4982E+03   -0.766   3.3493E+03   -0.154
  2.1904E+03   1.897   2.9998E+03   0.227   3.9319E+03   0.743   3.2757E+03   -0.531   1.6286E+03   -0.221
  4.2816E+03   2.560   3.8566E+03   1.157   2.9343E+03   1.469   2.0409E+03   -0.116   2.4836E+03   -0.242
  3.7283E+03   0.736   4.2442E+03   -0.602   3.1984E+03   -0.062   1.8426E+03   -1.404   1.1571E+03   -3.129
  3.9151E+03   1.873   4.0174E+03   0.472   3.1777E+03   0.941   1.9175E+03   -0.503   1.7805E+03   -1.263
  2.3725E+03   2.343   6.4756E+02   1.454   2.3391E+03   -2.182   3.1557E+03   2.753   3.5536E+03   3.141
  3.9275E+03   0.474   4.4078E+03   -0.923   3.4034E+03   -0.481   2.0061E+03   -1.867   1.1513E+03   2.755
  3.9496E+03   1.068   4.2440E+03   -0.289   3.0519E+03   0.207   1.6440E+03   -1.217   1.5362E+03   -2.612
  4.6192E+03   2.170   3.8137E+03   0.719   2.6932E+03   0.824   2.1703E+03   -0.896   2.7635E+03   -0.962
  2.9283E+03   2.214   3.4607E+03   0.682   3.7571E+03   1.219   2.8695E+03   -0.065   1.3860E+03   -0.124
  5.8533E+03   1.772   4.3184E+03   0.355   1.3299E+03   0.989   6.7423E+02   0.723   2.1845E+03   2.041
  2.6352E+03   2.439   8.2681E+02   1.336   2.1523E+03   -1.936   3.0129E+03   3.004   3.5664E+03   -2.861
  1.5175E+03   -0.003   2.9607E+03   -1.425   4.0257E+03   -0.914   3.5943E+03   -2.233   1.6967E+03   -1.771
  3.0617E+03   1.190   3.7698E+03   -0.342   3.9521E+03   0.099   2.9587E+03   -1.243   1.1424E+03   -1.499
  1.4315E+03   1.303   2.5777E+03   -0.531   3.9608E+03   -0.074   3.6044E+03   -1.381   2.0924E+03   -0.932
  4.3098E+03   2.894   3.6311E+03   1.458   2.9562E+03   1.773   2.3273E+03   0.239   2.6087E+03   0.255
  4.8314E+02   1.671   1.8121E+03   -0.853   3.6990E+03   -0.350   3.7179E+03   -1.646   2.6229E+03   -1.090
  4.4174E+03   1.055   4.3389E+03   -0.289   2.5334E+03   0.118   1.0943E+03   -1.583   2.1364E+03   -2.742
  4.6869E+03   1.943   3.8296E+03   0.477   2.6847E+03   0.534   2.1415E+03   -1.221   2.8188E+03   -1.336
  2.4075E+03   2.415   8.3998E+02   1.646   2.3697E+03   -2.369   3.1754E+03   2.608   3.5200E+03   3.012
  5.6480E+03   1.798   4.0811E+03   0.388   1.1029E+03   1.120   7.3144E+02   1.110   2.3510E+03   2.123
  9.6781E+02   0.255   2.5035E+03   -1.320   3.9364E+03   -0.831   3.7265E+03   -2.153   2.1608E+03   -1.674
  4.9255E+03   1.538   3.9814E+03   0.159   2.0151E+03   0.140   1.6200E+03   -1.962   2.9115E+03   -2.000
  3.2570E+03   2.881   1.4870E+03   1.144   2.4805E+03   0.002   3.0999E+03   -1.477   3.7186E+03   -1.008
  2.3551E+03   -0.007   3.5729E+03   -1.267   3.8461E+03   -0.660   3.1060E+03   -1.931   6.9503E+02   -1.306
  4.4669E+03   2.348   3.8119E+03   0.912   2.6862E+03   1.080   2.2126E+03   -0.597   2.6635E+03   -0.644
  8.1114E+03   1.172   6.6567E+03   -0.343   3.7208E+03   -0.258   2.3968E+03   -1.765   1.3142E+02   2.707
  1.5565E+03   2.779   8.4873E+02   2.982   2.9853E+03   -2.120   3.5863E+03   2.842   3.3729E+03   -3.008
  2.7544E+03   0.467   3.7686E+03   -1.005   3.9780E+03   -0.552   3.0736E+03   -1.890   8.1613E+02   -1.979
  1.1625E+03   2.640   7.3330E+02   -1.386   3.2019E+03   -1.048   3.6341E+03   -2.400   3.3176E+03   -1.918
  8.3370E+02   0.941   2.2871E+03   -0.866   3.8597E+03   -0.342   3.6735E+03   -1.636   2.2456E+03   -1.110
  1.6575E+03   2.763   3.6498E+02   -2.858   2.8393E+03   -1.544   3.4755E+03   -2.873   3.4417E+03   -2.403
  2.9696E+03   2.617   1.0942E+03   1.350   1.9292E+03   -1.508   2.8464E+03   -2.854   3.5661E+03   -2.398
  5.7640E+03   2.235   3.9377E+03   0.842   6.5744E+02   0.815   7.2354E+02   -3.024   2.6806E+03   -2.826
  6.7604E+03   1.731   5.2361E+03   0.291   2.1056E+03   0.611   8.5004E+02   -0.610   1.4808E+03   2.164
  3.9887E+03   2.102   3.9873E+03   0.710   3.0555E+03   1.156   1.8826E+03   -0.336   1.9625E+03   -0.931
  1.9540E+02   0.647   1.8283E+03   -1.622   3.7288E+03   -1.158   3.8257E+03   -2.499   2.7509E+03   -2.034
  1.8260E+03   2.630   1.7482E+02   -0.648   2.8140E+03   -1.152   3.4211E+03   -2.524   3.5189E+03   -2.061
  8.3081E+02   -0.856   2.2410E+03   -1.524   3.7363E+03   -0.790   3.7113E+03   -2.041   2.2605E+03   -1.382
  2.4219E+03   1.195   3.3920E+03   -0.354   3.9604E+03   0.142   3.1857E+03   -1.162   1.2237E+03   -1.034
  7.1998E+03   1.077   5.7884E+03   -0.440   2.8149E+03   -0.298   1.5757E+03   -1.720   8.4403E+02   0.814
  7.1863E+03   1.613   5.7496E+03   0.159   2.7061E+03   0.439   1.4578E+03   -0.899   9.8660E+02   1.752
  4.4988E+03   1.665   4.0925E+03   0.314   2.3357E+03   0.581   1.3172E+03   1.331   2.4975E+03   -1.741
  5.5160E+03   2.012   3.8407E+03   0.618   6.7653E+02   1.346   7.5253E+02   2.066   2.5624E+03   2.709
  V17   θ17u   V19   θ19u   V23   θ23u   V25   θ25u   V29   θ29
  4.5314E+02   -3.023   1.4568E+03   -2.207   3.5228E+03   -1.633   3.8192E+03   -2.977   3.0142E+03   -2.528
  2.3182E+03   -2.973   1.1274E+03   1.086   3.0446E+03   0.663   3.4346E+03   -0.663   3.4336E+03   -0.045
  4.1863E+03   2.070   2.5042E+03   0.750   9.9385E+02   3.102   1.9411E+03   2.129   3.1920E+03   2.574
  2.6138E+03   1.782   3.3757E+03   0.239   3.8261E+03   0.775   2.9802E+03   -0.508   1.2458E+03   -0.496
  4.1895E+03   1.890   2.4283E+03   0.518   8.3794E+02   -2.975   1.9145E+03   2.129   3.2396E+03   2.482
  4.4246E+03   2.563   3.7677E+03   1.158   2.7234E+03   1.347   2.1341E+03   -0.301   2.6502E+03   -0.289
  3.5508E+03   2.146   1.7226E+03   0.818   1.4084E+03   -2.296   2.4321E+03   2.637   3.4671E+03   3.003
  4.9186E+02   -0.181   2.1410E+03   -1.354   3.7735E+03   -0.749   3.7357E+03   -2.035   2.3736E+03   -1.464
  7.4013E+03   1.227   6.0035E+03   -0.275   3.0285E+03   -0.111   1.7794E+03   -1.547   6.6715E+02   1.027
  7.7756E+03   1.139   6.3843E+03   -0.375   3.4564E+03   -0.259   2.1899E+03   -1.736   2.7075E+02   0.694
  4.1938E+03   2.596   3.9317E+03   1.183   2.8660E+03   1.563   2.2641E+03   -0.017   2.4161E+03   -0.204
  5.7370E+03   2.284   3.8840E+03   0.857   8.7450E+02   0.348   1.1284E+03   -2.434   2.8965E+03   -2.467
  3.3664E+03   0.997   4.0138E+03   -0.421   3.6138E+03   0.085   2.4143E+03   -1.249   8.9837E+02   -2.266
  5.4152E+03   1.728   3.7586E+03   0.304   6.6286E+02   1.050   8.0345E+02   1.607   2.5721E+03   2.193
  4.2053E+03   2.563   3.8854E+03   1.168   2.9050E+03   1.516   1.9968E+03   -0.053   2.3599E+03   -0.257
  8.6312E+03   0.997   7.3282E+03   -0.538   4.5898E+03   -0.476   3.3424E+03   -1.998   9.1000E+02   -1.936
  4.1237E+03   0.536   4.4690E+03   -0.840   3.1569E+03   -0.396   1.6928E+03   -1.816   1.4464E+03   2.755
  4.3294E+03   3.038   2.3699E+03   1.713   1.0495E+03   -0.358   2.0937E+03   -1.900   3.4401E+03   -1.389
  3.3268E+03   2.034   3.7484E+03   0.562   3.6264E+03   1.078   2.5801E+03   -0.248   1.3689E+03   -0.652
  7.0740E+02   2.906   1.1116E+03   -1.322   3.3872E+03   -0.786   3.7012E+03   -2.097   3.0810E+03   -1.557
  4.1816E+03   2.253   3.9545E+03   0.866   2.8721E+03   1.243   1.8134E+03   -0.351   2.2782E+03   -0.727
  1.2176E+03   2.262   1.8408E+03   -0.142   3.7165E+03   0.215   3.6694E+03   -1.087   2.7120E+03   -0.544
  2.2082E+03   2.000   2.9915E+03   0.316   3.9039E+03   0.821   3.2928E+03   -0.456   1.6653E+03   -0.138
附表二(电流单位A,角度单位rad)
  I1   θ1i   I3   θ3i   I5   θ5i   I7   θ7i   I11   θ11i   I13
  1.0354E+01   -0.676   2.5245E-03   -2.179   1.8048E+00   -0.249   1.1156E+00   -1.619   4.3561E-01   -1.171   2.4780E-01
  1.0193E+01   -0.681   2.5112E-03   -2.200   1.8264E+00   -0.276   1.1653E+00   -1.658   5.1716E-01   -1.229   3.3369E-01
  1.0335E+01   -0.694   6.1086E-04   -2.271   1.8505E+00   -0.342   1.1608E+00   -1.743   5.0348E-01   -1.381   3.1509E-01
  1.0351E+01   -0.696   7.2075E-04   -2.241   1.8723E+00   -0.354   1.1854E+00   -1.755   5.3530E-01   -1.405   3.4765E-01
  1.0913E+01   -0.635   2.2781E-04   1.431   1.8789E+00   -0.053   1.1438E+00   -1.338   4.1916E-01   -0.784   2.2050E-01
  1.1631E+01   -0.579   5.2168E-03   1.344   1.9034E+00   0.211   1.0931E+00   -0.983   3.3836E 01   -0.473   1.7911E-01
  1.1029E+01   -0.632   2.5926E-03   1.140   1.8399E+00   -0.041   1.0805E+00   -1.334   3.3950E-01   -0.868   1.5194E-01
  1.0347E+01   -0.667   2.9320E-03   -2.152   1.7913E+00   -0.204   1.1071E+00   -1.554   4.2314E-01   -1.057   2.3621E-01
  1.1068E+01   -0.623   1.2993E-03   1.172   1.8747E+00   0.006   1.1210E+00   -1.257   3.7949E-01   -0.694   1.8225E-01
  11502E+01   -0.592   5.1186E-04   0.966   1.8605E+00   0.144   1.0669E+00   -1.087   3.2302E-01   -0.685   1.8149E-01
  1.0768E+01   -0.641   2.2274E-03   -2.132   1.6909E+00   -0.089   9.4055E-01   -1.402   1.9631E-01   -1.040   5.1175E-02
  1.0635E+01   -0.659   8.4287E-04   -2.167   1.8338E+00   -0.170   1.1169E+00   -1.506   4.1131E-01   -1.024   2.1755E-01
  1.1081E+01   -0.621   1.3778E-03   1.218   1.8174E+00   0.010   1.0520E+00   -1.259   3.0217E-01   -0.777   1.2189E-01
  1.1712E+01   -0.584   1.2033E-02   1.389   1.8870E+00   0.171   1.0623E+00   -1.055   3.5616E-01   -0.707   2.2929E-01
  9.9949E+00   -0.729   9.6167E-04   -2.338   1.8641E+00   -0.518   1.2146E+00   -1.986   6.1737E 01   -1.766   4.4170E-01
  9.6629E+00   -0764   4.5181E-04   -2.503   1.8542E+00   -0.688   1.2400E+00   -2.237   6.9477E-01   -2.159   5.3240E-01
  1.1530E+01   -0.584   4.9580E-03   1.311   2.0028E+00   0.202   1.2282E+00   -0.979   4.7239E-01   -0.264   2.6318E-01
  9.9961E+00   -0.732   3.3427E-05   1.269   1.8831E+00   -0.533   1.2394E+00   -2.009   6.5249E-01   -1.806   4.7846E-01
  1.0311E+01   -0.695   3.8458E-04   -2.259   1.8651E+00   -0.348   1.1833E+00   -1.751   5.3588E-01   -1.394   3.4886E-01
  9.8287E+00   -0.779   2.4203E-04   0.872   1.8959E+00   -0.767   1.2634E+00   -2.333   7.1201E-01   -2.322   5.4623E-01
  I1   θ1i   I3   θ3i   I5   θ5i   I7   θ7i   I11   θ11i   I13
  1.1917E+01   -0.571   3.0634E-03   3.043   1.9857E+00   0.247   1.1695E+00   -0.939   4.2477E-01   -0.389   2.6822E-01
  1.0615E+01   -0.654   4.9694E-04   -2.156   1.8766E+00   -0.141   1.1784E+00   -1.462   4.9071E-01   -0.942   2.9636E-01
  1.0108E+01   -0.694   1.9611E-03   -2.222   1.8308E+00   -0.343   1.1747E+00   -1.750   5.4174E-01   -1.379   3.6051E-01
  9.5759E+00   -0.765   3.4820E-04   -2.465   1.8259E+00   -0.700   1.2108E+00   -2.248   6.6503E-01   -2.177   5.0231E-01
  1.0846E+01   -0.631   1.0556E-03   -2.169   1.7875E+00   -0.031   1.0467E+00   -1.309   3.0537E-01   -0.750   1.1488E-01
  1.0803E+01   -0.655   4.3324E-04   1.142   1.9490E+00   -0.150   1.2392E+00   -1.470   5.5531E-01   -0.972   3.5789E-01
  1.0943E+01   -0.641   7.6996E-04   1.144   1.8827E+00   -0.083   1.1397E+00   -1.386   4.1392E-01   -0.878   2.1605E-01
  1.1484E+01   -0.613   5.5040E-03   1.223   2.0609E+00   0.060   1.2897E+00   -1.171   5.6020E-01   -0.543   3.4825E-01
  1.0650E+01   -0.665   1.5607E-04   1.294   1.8932E+00   -0.201   1.1839E+00   -1.544   4.9900E-01   -1.083   3.0366E-01
  1.1074E+01   -0.619   1.6021E-03   1.058   1.7272E+00   0.014   9.4785E-01   -1.260   2.0468E-01   -0.991   9.2660E-02
  1.1522E+01   -0.593   3.3000E-03   1.311   1.8056E+00   0.130   9.9078E-01   -1.112   2.7893E-01   -0.893   1.8142E-01
  1.1552E+01   -0.593   4.8522E-03   1.284   1.8031E+00   0.122   9.8332E-01   -1.128   2.9018E-01   -0.953   1.9945E-01
  9.5798E+00   -0.769   2.9684E-04   -2.495   1.8300E+00   -0.719   1.2142E+00   -2.273   6.7001E-01   -2.218   5.0753E-01
  1.0272E+01   -0.696   3.5741E-04   -2.312   1.8970E+00   -0.349   1.2346E+00   -1.755   6.0850E-01   -1.402   4.2536E-01
  1.1357E+01   -0.606   2.5017E-03   1.244   2.0169E+00   0.092   1.2635E+00   -1.127   5.3108E-01   -0.464   3.2217E-01
  1.1376E+01   -0.599   5.7749E-03   1.282   1.9296E+00   0.122   1.1496E+00   -1.094   3.9514E-01   -0.479   1.9575E-01
  1.1827E+01   -0.574   5.9156E-03   1.464   1.9326E+00   0.230   1.0968E+00   -0.963   3.6159E-01   -0.502   2.1831E-01
  1.0965E+01   -0.622   1.1817E-04   -2.278   1.7336E+00   0.006   9.6907E-01   -1.260   2.1594E-01   -0.812   5.8408E-02
  1.1456E+01   -0.593   6.0732E-03   1.281   1.8530E+00   0.141   1.0520E+00   -1.084   3.0373E-01   -0.655   1.5442E-01
  1.0211E+01   -0.695   2.0479E-03   -2.232   1.8358E+00   -0.348   1.1631E+00   -1.756   5.1666E-01   -1.392   3.3232E-01
  1.0808E+01   -0.643   2.2896E-03   -2.143   1.7090E+00   -0.101   9.5640E-01   -1.424   2.1491E-01   -1.111   7.1405E-02
  1.1113E+01   -0.615   1.3317E-03   1.121   1.7550E+00   0.035   9.7614E-01   -1.225   2.2583E-01   -0.861   8.8553E-02
  1.1602E+01   -0.600   6.9124E-03   1.358   1.8700E+00   0.090   1.0660E+00   -1.171   3.6746E-01   -0.885   2.4263E-01
  1.1544E+01   -0.584   3.0254E-03   1.338   1.9345E+00   0.198   1.1441E+00   -0.995   3.7944E-01   -0.379   1.9109E-01
  9.7810E+00   -0.724   1.3878E-03   -2.309   1.8283E+00   -0.490   1.2055E+00   -1.961   6.2520E-01   -1.716   4.5559E-01
  1.0318E+01   -0.687   7.0839E-04   -2.220   1.8604E+00   -0.307   1.1826E+00   -1.696   5.3161E-01   -1.304   3.4445E-01
  1.0625E+01   -0.652   1.0256E-03   -1.935   1.7883E+00   -0.134   1.0682E+00   -1.455   3.4950E-01   -0.943   1.5852E-01
  1.1102E+01   -0.628   2.5470E-03   1.287   1.8403E+00   -0.023   1.0746E+00   -1.314   3.3134E-01   -0.879   1.5220E-01
  1.1044E+01   -0.628   1.4860E-03   1.174   1.8870E+00   -0.016   1.1370E+00   -1.289   4.0123E-01   -0.735   2.0249E-01
  1.1886E+01   -0.565   1.0182E-02   1.488   1.9639E+00   0.276   1.1429E+00   -0.901   4.1036E-01   -0.363   2.5496E-01
  1.0908E+01   -0.632   2.2082E-04   1.339   1.8864E+00   -0.034   1.1568E+00   -1.308   4.3468E-01   -0.727   2.3529E-01
  1.1098E+01   -0.617   1.9381E-03   1.030   1.7130E+00   0.019   9.2865E-01   -1.255   1.9211E-01   -1.066   1.0301E-01
  1.1537E+01   -0.611   1.1068E-02   1.272   1.8548E+00   0.030   1.0455E+00   -1.254   3.5890E-01   -1.040   2.3870E-01
  1.0082E+01   -0.697   2.5074E-03   -2.232   1.8099E+00   -0.353   1.1503E+00   -1.766   5.1172E-01   -1.398   3.3041E-01
  9.8213E+00   -0.723   1.4656E-03   -2.309   1.8324E+00   -0.483   1.2051E+00   -1.951   6.2007E-01   -1.700   4.4916E-01
  1.0692E+01   -0.650   8.0244E-04   -2.259   1.8229E+00   -0.125   1.1008E+00   -1.442   3.8475E-01   -0.929   1.9069E-01
  1.1300E+01   -0.611   4.4000E-03   1.012   1.7169E+00   0.030   9.1748E-01   -1.269   2.4960E-01   -1.317   1.9222E-01
  1.1240E+01   -0.639   2.4490E-03   1.185   2.0374E+00   -0.070   1.2882E+00   -1.355   5.8084E-01   -0.814   3.7478E-01
  1.0660E+01   -0.636   2.5740E-03   -2.120   1.7511E+00   -0.053   1.0296E+00   -1.336   2.9810E-01   -0.726   1.1054E-01
  1.1714E+01   -0.588   5.3426E-03   1.746   1.8959E+00   0.153   1.0808E+00   -1.082   3.6805E-01   -0.720   2.3997E-01
  1.0010E+01   -0.757   5.4347E-04   0.861   1.9183E+00   -0.656   1.2744E+00   -2.179   7.0427E-01   -2.080   5.3312E-01
  1.0128E+01   -0.688   2.7759E-03   -2.218   1.7927E+00   -0.308   1.1278E+00   -1.703   4.7502E-01   -1.291   2.9225E-01
  1.0806E+01   -0.646   2.9891E-04   -2.265   1.7786E+00   -0.110   1.0337E+00   -1.430   2.9765E-01   -0.995   1.1382E-01
  1.0669E+01   -0.659   5.2641E-05   2.988   1.8838E+00   -0.170   1.1734E+00   -1.502   4.8159E-01   -1.016   2.8593E-01
  1.0900E+01   -0.632   6.7789E-04   1.197   1.8664E+00   -0.035   1.1323E+00   -1.311   4.0536E-01   -0.738   2.0690E-01
  1.0421E+01   -0.672   6.0160E-04   -2.176   1.8519E+00   -0.229   1.1658E+00   -1.587   4.9575E-01   -1.132   3.0570E-01
  I1   θ1i   I3   θ3i   I5   θ5i   I7   θ7i   I11   θ11i   I13
  1.0497E+01   -0.673   3.4084E-04   -2.311   1.8986E+00   -0.238   1.2130E+00   -1.598   5.5206E-01   -1.157   3.6142E-01
  1.0570E+01   -0.693   1.0977E-03   1.008   1.9733E+00   -0.341   1.2895E+00   -1.733   6.5448E-01   -1.383   4.6597E-01
  9.9640E+00   -0.725   8.0574E-04   -2.363   1.8793E+00   -0.494   1.2460E+00   -1.962   6.6215E-01   -1.726   4.9018E-01
  1.1568E+01   -0.585   3.7943E-03   1.334   1.8729E+00   0.180   1.0653E+00   -1.029   3.1646E-01   -0.584   1.7018E-01
  1.0579E+01   -0.662   1.3711E-03   -2.182   1.8312E+00   -0.184   1.1210E+00   -1.525   4.2126E-01   -1.045   2.2837E-01
  1.0657E+01   -0.665   3.7411E-04   1.184   1.9006E+00   -0.197   1.1931E+00   -1.537   5.1012E-01   -1.072   3.1472E-01
  1.0589E+01   -0.635   6.6605E-04   -2.072   1.8086E+00   -0.048   1.1105E+00   -1.327   4.0250E-01   -0.702   2.1087E-01
  1.1113E+01   -0.621   6.9752E-04   1.263   1.8564E+00   0.014   1.0949E+00   -1.251   3.4702E-01   -0.726   1.5662E-01
  9.6830E+00   -0.763   5.4650E-04   -2.452   1.8425E+00   -0.690   1.2183E+00   -2.227   6.6298E-01   -2.148   4.9748E-01
  9.8303E+00   -0.732   3.2389E-04   -2.370   1.8615E+00   -0.527   1.2397E+00   -2.012   6.6955E-01   -1.803   5.0112E-01
  1.1402E+01   -0.599   1.8163E-03   1.253   1.7716E+00   0.103   9.6358E-01   -1.145   2.4226E-01   -0.958   1.5297E-01
  1.0111E+01   -0.709   8.5902E-04   -2.291   1.8830E+00   -0.416   1.2348E+00   -1.851   6.2799E-01   -1.551   4.5007E-01
  1.1412E+01   -0.604   3.6202E-03   1.254   2.0341E+00   0.107   1.2772E+00   -1.106   5.4448E-01   -0.432   3.3401E-01
  9.9497E+00   -0.710   1.7877E-03   -2.280   1.8268E+00   -0.419   1.1845E+00   -1.858   5.7501E-01   -1.550   3.9863E-01
  1.1407E+01   -0.598   3.1907E-03   1.277   1.9107E+00   0.127   1.1231E+00   -1.087   3.6387E-01   -0.499   1.7206E-01
  1.0126E+01   -0.718   8.4246E-04   -2.337   1.8652E+00   -0.464   1.2001E+00   -1.906   5.8230E-01   -1.641   4.0150E-01
  1.1753E+01   -0.576   9.7100E-03   1.626   1.9030E+00   0.213   1.0854E+00   -0.998   3.6700E-01   -0.584   2.3560E-01
  1.0300E+01   -0.702   5.3200E-04   -2.290   1.8821E+00   -0.385   1.2034E+00   -1.796   5.6599E-01   -1.470   3.8040E-01
  1.0679E+01   -0.641   1.0994E-04   -2.345   1.8306E+00   -0.078   1.1183E+00   -1.372   4.0606E-01   -0.800   2.1140E-01
  9.7378E+00   -0.754   2.5654E-04   -2.387   1.8539E+00   -0.642   1.2302E+00   -2.166   6.7181E-01   -2.050   5.0549E-01
  9.7823E+00   -0.759   7.3265E-05   1.047   1.8693E+00   -0.667   1.2422E+00   -2.199   6.8465E-01   -2.104   5.1798E-01
  1.1863E+01   -0.571   1.3328E-02   1.469   1.9430E+00   0.243   1.1033E+00   -0.947   3.7773E-01   -0.477   2.3304E-01
  1.0831E+01   -0.688   2.3095E-03   -2.222   2.0227E+00   -0.311   1.3175E+00   -1.693   6.6313E-01   -1.326   4.7194E-01
  1.1061E+01   -0.622   1.2144E-03   1.154   1.7917E+00   0.004   1.0242E+00   -1.269   2.7438E-01   -0.829   1.0463E-01
  9.9922E+00   -0.726   1.0396E-03   -2.330   1.8636E+00   -0.503   1.2163E+00   -1.968   6.1917E-01   -1.737   4.4377E-01
  1.1732E+01   -0.573   2.8212E-03   1.563   1.9110E+00   0.237   1.1014E+00   -0.956   3.5968E-01   -0.478   2.1523E-01
  9.7723E+00   -0.767   1.3711E-04   0.849   1.8813E+00   -0.704   1.2581E+00   -2.257   7.0880E-01   -2.195   5.4442E-01
  1.0838E+01   -0.639   9.6378E-04   -1.948   1.6973E+00   -0.083   9.3940E-01   -1.397   1.9800E-01   -1.120   7.1369E-02
  1.0960E+01   -0.644   2.1853E-03   1.142   2.0084E+00   -0.097   1.3004E+00   -1.386   6.2085E-01   -0.847   4.2104E 01
  1.1485E+01   -0590   3.3052E-03   1.329   1.9015E+00   0.164   1.1111E+00   -1.047   3.5167E-01   -0.502   1.7520E-01
  1.0721E+01   -0.646   1.5113E-04   -2.182   1.8767E+00   -0.102   1.1661E+00   -1.406   4.6419E-01   -0.862   2.6767E 01
  1.1675E+01   -0.581   5.4829E-03   1.411   1.8868E+00   0.197   1.0629E+00   -1.007   3.2316E-01   -0.585   1.8560E-01
  1.1197E+01   -0.617   2.8258E-03   1.199   1.9480E+00   0.037   1.1931E+00   -1.209   4.5817E-01   -0.602   2.5348E-01
  1.1440E+01   -0.597   2.7047E-03   1.280   1.9441E+00   0.134   1.1603E+00   -1.076   4.0230E-01   -0.453   2.0336E-01
附表二(续)
  I17   θ17i   I19   θ19i   I23   θ23i   I25   θ25i   I29   θ29i
  1.8453E-02   -1.397   4.7481E-02   -0.634   9.7620E-02   -0.034   9.7184E-02   -1.380   6.6908E-02   -0.932
  9.3989E-02   -1.998   3.1565E-02   -2.771   6.4907E-02   -0.405   7.9718E-02   -1.699   7.8542E-02   -1.221
  7.2852E-02   -2.376   6.5037E-03   -2.750   7.3465E-02   -0.385   8.5854E-02   -1.780   7.8608E-02   -1.407
  1.0201E-01   -2.428   2.9498E-02   2.610   5.8031E-02   -0.452   7.5686E-02   -1.845   7.9244E-02   -1.477
  2.9477E-02   2.799   7.2498E-02   0.770   1.0560E-01   1.210   9.5647E-02   -0.095   5.2767E-02   0.405
  1.4048E-01   -2.453   1.3440E-01   2.431   8.5677E-02   2.896   5.4165E-02   1.435   3.8895E-02   0.902
  9.7467E-02   2.597   1.2127E-01   1.066   1.0922E-01   1.499   8.2240E-02   0.165   2.4655E-02   0.181
  2.2411E-02   -0.611   5.7332E-02   -0.591   1.0048E-01   0.093   9.7511E-02   -1.228   6.3460E-02   -0.720
  6.3645E-02   2.749   9.8633E-02   1.124   1.0870E-01   1.600   8.9572E-02   0.304   3.6391E-02   0.708
  1.5261E-01   -2.669   1.4046E-01   2.218   8.1628E-02   2.615   5.0543E-02   1.053   4.3034E-02   0.455
  I17   θ17i   I19   θ19i   I23   θ23i   I25   θ25i   I29   θ29i
  1.4786E-01   1.954   1.5080E-01   0.602   8.9156E-02   1.076   4.8116E-02   -0.316   2.5877E-02   -2.242
  1.7444E-02   1.992   6.8880E-02   0.115   1.0600E-01   0.550   9.8200E-02   -0.792   5.7308E-02   -0.329
  1.1219E-01   2.630   1.3032E-01   1.199   1.0260E-01   1.695   7.0125E-02   0.367   1.6873E-02   -0.205
  1.6294E-01   -2.305   1.2969E-01   2.562   7.0988E-02   2.657   5.6691E-02   0.946   5.6351E-02   0.990
  2.0300E-01   -2.991   1.2517E-01   1.864   1.7559E-02   2.566   2.1258E-02   3.135   5.7391E-02   -2.557
  3.1294E-01   2.663   2.3598E-01   1.132   1.1912E-01   1.216   7.7571E-02   -0.287   1.4895E-02   -0.019
  6.5676E-02   -1.891   7.6120E-02   2.348   1.0227E-01   2.719   9.2495E-02   1.464   5.3453E-02   2.047
  2.4265E-01   -3.067   1.6290E-01   1.762   4.7645E-02   2.033   1.2228E-02   1.064   3.9455E-02   -2.635
  1.0416E-01   -2.397   3.1898E-02   2.679   5.6674E-02   -0.473   7.4751E-02   -1.853   7.9072E-02   -1.473
  3.2454E-01   2.403   2.4558E-01   0.838   1.2715E-01   0.836   8.4087E-02   -0.720   1.9509E-02   -0.825
  1.5091E-01   -1.820   1.2668E-01   3.012   8.2065E-02   -3.059   6.6725E-02   1.691   5.2885E-02   1.801
  5.1645E-02   -1.696   1.6488E-02   -0.202   8.3890E-02   0.555   9.0961E-02   -0.761   7.4074E-02   -0.237
  1.1983E-01   -2.299   5.0145E-02   2.854   4.9317E-02   -0.775   6.8713E-02   -2.052   7.7679E-02   -1.606
  2.8107E-01   2.644   2.0467E-01   1.118   8.9980E-02   1.231   5.0927E-02   -0.206   1.1563E-02   1.961
  8.9295E-02   2.046   1.1986E-01   0.686   1.0688E-01   1.241   7.9064E-02   -0.053   1.6041E-02   0.215
  1.0418E-01   -1.886   2.9773E-02   2.825   6.1224E-02   0.767   7.7975E-02   -0.620   7.9219E-02   -0.147
  4.6443E-02   2.793   8.1610E-02   0.837   1.0863E-01   1.192   9.6309E-02   -0.141   5.1398E-02   0.286
  1.0104E-01   -1.582   5.2772E-02   2.648   8.4360E-02   2.198   9.0086E-02   0.829   7.6009E-02   1.362
  5.8187E-02   -2.086   1.3019E-02   0.376   8.2170E-02   0.418   9.0595E-02   -0.948   7.5739E-02   -0.496
  1.5758E-01   2.586   1.4825E-01   1.241   7.5059E-02   1.694   3.4587E-02   0.111   4.0076E-02   -1.263
  1.6593E-01   -2.766   1.3694E-01   2.152   6.1856E-02   2.314   4.3657E-02   0.404   5.5807E-02   0.243
  1.6851E-01   -2.706   1.3291E-01   2.201   6.0409E-02   2.214   4.8987E-02   0.319   5.7821E-02   0.338
  2.8692E-01   2.578   2.1042E-01   1.043   9.5376E-02   1.130   5.5801E-02   -0.335   8.3400E-03   1.576
  1.7961E-01   -2.434   1.0061E-01   2.492   9.1985E-03   -1.114   3.7689E-02   -1.978   6.7764E-02   -1.516
  7.5412E-02   -1.486   3.9578E-02   2.367   8.6664E-02   2.140   9.1320E-02   0.830   7.2304E-02   1.432
  8.1167E-02   -2.828   1.0472E-01   1.822   1.0649E-01   2.311   8.4465E-02   1.032   3.3266E-02   1.329
  1.5575E-01   -2.160   1.3360E-01   2.724   7.7282E-02   3.010   5.5592E-02   1.388   5.2546E-02   1.283
  1.3964E-01   2.323   1.4406E-01   0.994   8.7349E-02   1.532   4.7009E-02   0.163   2.5369E-02   -1.586
  1.4442E-01   -2.842   1.3855E-01   2.053   8.3937E-02   2.506   4.9597E-02   1.007   3.8191E-02   0.237
  9.0720E-02   -2.326   2.4030E-02   3.049   6.4447E-02   -0.575   7.9765E-02   -1.935   7.9176E-02   -1.540
  1.4666E-01   2.050   1.5071E-01   0.670   9.2125E-02   1.110   5.2252E-02   -0.301   2.4824E-02   -2.017
  1.4716E-01   2.648   1.4473E-01   1.302   8.1821E-02   1.791   4.1865E-02   0.327   3.4020E-02   -1.087
  1.6693E-01   -2.502   1.3168E-01   2.335   7.1917E-02   2.352   5.8728E-02   0.622   5.7608E-02   0.643
  1.0828E-01   -2.504   1.2005E-01   2.269   1.0137E-01   2.787   7.4292E-02   1.486   2.8566E-02   1.389
  2.2313E-01   -2.900   1.4673E-01   1.969   3.8093E-02   2.536   1.7833E-02   2.222   4.7530E-02   -2.627
  1.0012E-01   -2.235   2.8947E-02   2.943   5.9262E-02   -0.326   7.6692E-02   -1.682   7.9167E-02   -1.259
  5.7200E-02   1.556   9.9324E-02   0.171   1.1095E-01   0.674   9.1984E-02   -0.642   3.6668E-02   -0.179
  1.1376E-01   2.764   1.3044E-01   1.256   1.0625E-01   1.679   7.6352E-02   0.326   2.2931E-02   0.000
  5.3246E-02   2.854   8.9055E-02   1.075   1.0851E-01   1.506   9.3331E 02   0.200   4.4348E-02   0.651
  1.5452E-01   -1.777   1.2602E-01   3.076   7.9099E-02   -2.986   6.2816E-02   1.743   5.4094E-02   1.863
  1.7407E-02   -3.059   6.1640E-02   0.755   1.0215E-01   1.238   9.5522E-02   -0.054   5.6860E-02   0.490
  1.6441E-01   2.636   1.4830E-01   1.296   6.7366E-02   1.705   2.8590E-02   -0.074   4.6432E-02   -1.203
  1.6953E-01   -2.730   1.3256E-01   2.095   7.1316E-02   2.057   5.7806E-02   0.283   5.8273E-02   0.259
  9.2909E-02   -2.260   3.1489E-02   -3.054   6.4925E-02   -0.759   7.9555E-02   -2.085   7.8695E-02   -1.669
  2.1522E-01   -2.874   1.3880E-01   2.002   3.1618E-02   2.659   1.8703E-02   2.608   5.1315E-02   -2.547
  3.6096E-02   1.813   8.4362E-02   0.281   1.0880E-01   0.760   9.5877E-02   -0.559   4.7194E-02   -0.082
  I17   θ17i   I19   θ19i   I23   θ23i   I25   θ25i   I29   θ29i
  1.7976E-01   3.138   1.3750E-01   1.761   5.2038E-02   1.662   4.5054E-02   -0.449   6.1646E-02   -0.420
  1.1899E-01   -1.810   5.1069E-02   2.669   6.8344E-02   1.569   8.1888E-02   0.116   8.0474E-02   0.562
  8.7805E-02   1.586   1.1906E-01   0.332   1.0655E-01   0.941   7.9186E-02   -0.329   1.5453E-02   0.312
  1.6073E-01   -2.313   1.3026E-01   2.541   7.2473E-02   2.616   5.9299E-02   0.927   5.6336E-02   0.976
  3.0392E-01   2.776   2.2317E-01   1.257   1.0328E-01   1.333   6.0720E-02   -0.179   3.3453E-03   -1.952
  6.0762E-02   -1.908   2.8946E-02   -1.810   8.2490E-02   -0.511   9.0139E-02   -1.846   7.5884E-02   -1.407
  1.0343E-01   2.042   1.2829E-01   0.592   1.0906E-01   1.038   7.9172E-02   -0.308   1.6379E-02   -0.379
  4.2486E-02   -2.022   2.5036E-02   0.217   8.8882E-02   0.543   9.3357E-02   -0.803   7.3047E-02   -0.329
  3.1314E-02   2.496   7.7981E-02   0.741   1.0651E-01   1.246   9.4318E-02   -0.046   4.8560E-02   0.474
  6.3417E-02   -1.910   1.2354E-02   -1.378   7.8472E-02   0.059   8.8327E-02   -1.277   7.6459E-02   -0.808
  1.1222E-01   -2.059   3.7384E-02   2.956   5.3184E-02   0.093   7.2547E-02   -1.260   7.9048E-02   -0.806
  2.1528E-01   -2.446   1.3150E-01   2.434   1.8046E-02   2.392   1.9763E-02   -1.406   5.9110E-02   -1.239
  2.5584E-01   -2.943   1.7654E-01   1.898   5.9935E-02   2.198   2.2737E-02   0.994   3.1810E-02   -2.508
  1.4700E-01   -2.610   1.3752E-01   2.291   8.0865E-02   2.722   4.9261E-02   1.187   4.1845E-02   0.589
  6.3476E-03   2.066   6.1276E-02   -0.027   1.0354E-01   0.435   9.8324E-02   -0.907   6.1021E-02   -0.438
  6.7520E-02   -2.039   6.4942E-03   0.939   7.7995E-02   0.434   8.8189E-02   -0.931   7.7075E-02   -0.477
  2.9915E-02   0.738   7.4031E-02   0.067   1.0407E-01   0.819   9.4252E-02   -0.438   5.1162E-02   0.208
  9.0336E-02   2.764   1.1648E-01   1.244   1.0770E-01   1.723   8.2093E-02   0.411   2.4843E-02   0.536
  2.7229E-01   2.687   1.9506E-01   1.167   7.9791E-02   1.293   4.1087E-02   -0.117   1.7506E-02   2.436
  2.7201E-01   -3.061   1.9345E-01   1.767   7.6774E-02   2.028   3.8284E-02   0.709   2.0773E-02   -2.918
  1.6530E-01   -3.042   1.4057E-01   1.894   6.0621E-02   2.144   3.6379E-02   0.179   5.3419E-02   -0.201
  2.0889E-01   -2.663   1.2979E-01   2.227   1.9565E-02   2.886   1.8755E-02   -2.664   5.6485E-02   -1.971
  8.4679E-02   -1.393   4.0578E-02   2.624   8.3936E-02   2.250   8.9819E-02   0.934   7.3769E-02   1.536
  1.6038E-01   -2.604   8.6535E-02   2.369   2.5933E-02   -1.594   4.8240E-02   -2.582   7.0713E-02   -2.101
  9.7206E-02   -2.935   1.1682E-01   1.833   1.0347E-01   2.348   7.6742E-02   1.057   2.5176E-02   1.027
  1.5900E-01   -2.784   8.2775E-02   2.131   2.2193E-02   -1.373   4.8094E-02   -2.569   7.1590E-02   -2.202
  1.6011E-01   -2.131   1.2914E-01   2.744   7.2093E-02   2.890   5.7969E-02   1.213   5.5427E-02   1.275
  1.3424E-01   -2.528   5.8758E-02   2.431   3.8635E-02   -0.692   6.1683E-02   -2.055   7.6831E-02   -1.688
  1.8194E-02   1.399   7.1431E-02   0.244   1.0491E-01   0.851   9.5769E-02   -0.439   5.2839E-02   0.127
  2.7991E-01   2.836   2.0206E-01   1.333   8.5732E-02   1.480   4.6230E-02   0.064   1.3636E-02   2.659
  2.9300E-01   2.747   2.1455E-01   1.229   9.7231E-02   1.332   5.6188E-02   -0.140   4.9338E-03   2.312
  1.5712E-01   -2.025   1.2915E-01   2.845   7.6859E-02   3.070   5.8663E-02   1.449   5.4139E-02   1.488
  2.1416E-01   -2.392   1.3060E-01   2.455   2.4059E-02   1.943   2.9543E-02   -0.835   6.3752E-02   -0.877
  1.2600E-01   2.573   1.3773E-01   1.174   9.7710E 02   1.669   6.1846E-02   0.315   1.8937E-02   -0.801
  2.0519E-01   -2.946   1.2730E-01   1.914   1.9087E-02   2.582   1.9965E-02   -3.122   5.6539E-02   -2.484
  1.5361E-01   -2.149   1.3383E-01   2.741   7.6640E-02   3.068   5.3322E-02   1.457   5.0212E-02   1.267
  3.2402E-01   2.604   2.4587E-01   1.064   1.2790E-01   1.117   8.5237E-02   -0.407   2.0560E-02   -0.374
  1.5432E-01   2.120   1.5223E-01   0.751   8.5536E-02   1.198   4.3646E-02   -0.260   3.1482E-02   -1.998
  1.6134E-01   -1.647   7.8624E-02   -2.996   2.9423E-02   1.237   5.5035E-02   -0.295   7.5703E-02   0.194
  1.2311E-01   -2.679   1.2959E-01   2.148   9.7169E-02   2.649   6.6319E-02   1.300   2.8356E-02   0.847
  2.6147E-02   -1.745   3.7334E-02   0.282   9.3878E-02   0.806   9.4775E-02   -0.502   6.7787E-02   0.030
  1.5300E-01   -2.461   1.3637E-01   2.442   7.5462E-02   2.797   4.8563E-02   1.153   4.8353E-02   0.791
  4.4015E-02   -2.466   6.5005E-02   1.462   1.0190E-01   1.791   9.5264E-02   0.497   5.7900E-02   1.031
  8.1850E-02   -2.722   1.0388E-01   1.911   1.0629E-01   2.392   8.5038E-02   1.110   3.4658E-02   1.422
附表三
  序号i   MSE(i)1   MSE(i)2   R(i)
  1   0.225785   0.174436   22.74%
  2   0.222771   0.140624   36.88%
  3   0.219369   0.151782   30.81%
  4   0.285674   0.196681   31.15%
  5   0.200859   0.117462   41.52%
  6   0.241725   0.131118   45.76%
  7   0.173904   0.12585   27.63%
  8   0.265538   0.170926   35.63%
  9   0.252554   0.188398   25.40%
  10   0.200698   0.158414   21.07%
  11   0.244568   0.16654   31.90%
  12   0.266518   0.175618   34.11%
  13   0.262964   0.193943   26.25%
  14   0.205767   0.136285   33.77%
  15   0.178131   0.15117   15.14%
  16   0.242517   0.144842   40.28%
  17   0.190253   0.148196   22.11%
  18   0.178386   0.128682   27.86%
  19   0.195365   0.12737   34.80%
  20   0.166367   0.120502   27.57%
  21   0.17713   0.122623   30.77%
  22   0.238738   0.158431   33.64%
  23   0.228193   0.127405   44.17%
  24   0.200561   0.142065   29.17%
  25   0.191897   0.118061   38.48%
  26   0.221612   0.121789   45.04%
  27   0.179281   0.094282   47.41%
  28   0.180134   0.116131   35.53%
  29   0.18556   0.155508   16.19%
  30   0.238921   0.149834   37.29%
  31   0.264738   0.156879   40.74%
  32   0.185014   0.173895   6.01%
  33   0.204594   0.120216   41.24%
  34   0.231572   0.194926   15.82%
  35   0.189103   0.120844   36.10%
  36   0.266509   0.140171   47.40%
  37   0.204778   0.1762   13.96%
  38   0.196356   0.111997   42.96%
  39   0.205158   0.160156   21.94%
  序号i   MSE(i)1   MSE(i)2   R(i)
  40   0.219708   0.115601   47.38%
  41   0.285465   0.152447   46.60%
  42   0.159576   0.107874   32.40%
  43   0.246349   0.141367   42.62%
  44   0.258427   0.166202   35.69%
  45   0.180346   0.123356   31.60%
  46   0.303799   0.171047   43.70%
  47   0.208293   0.177264   14.90%
  48   0.165682   0.121439   26.70%
  49   0.189464   0.153487   18.99%
  50   0.237498   0.149028   37.25%
  51   0.199378   0.146182   26.68%
  52   0.237939   0.172508   27.50%
  53   0.186093   0.151056   18.83%
  54   0.264307   0.189193   28.42%
  55   0.225921   0.178834   20.84%
  56   0.246504   0.169229   31.35%
  57   0.190387   0.137536   27.76%
  58   0.227416   0.148789   34.57%
  59   0.234729   0.136373   41.90%
  60   0.186583   0.126433   32.24%
  61   0.209834   0.172004   18.03%
  62   0.227755   0.141199   38.00%
  63   0.189872   0.132602   30.16%
  64   0.184515   0.134046   27.35%
  65   0.20329   0.171196   15.79%
  66   0.214476   0.17969   16.22%
  67   0.227483   0.18916   16.85%
  68   0.231371   0.15906   31.25%
  69   0.225351   0.16981   24.65%
  70   0.21576   0.121369   43.75%
  71   0.269272   0.182093   32.38%
  72   0.211255   0.137489   34.92%
  73   0.15669   0.106713   31.90%
  74   0.228594   0.153705   32.76%
  75   0.300432   0.202689   32.53%
  76   0.274782   0.202039   26.47%
  77   0.355836   0.198848   44.12%
  78   0.197504   0.140051   29.09%
  79   0.193739   0.142915   26.23%
  序号i   MSE(i)1   MSE(i)2   R(i)
  80   0.261443   0.166172   36.44%
  81   0.282959   0.217548   23.12%
  82   0.203374   0.110308   45.76%
  83   0.282999   0.192667   31.92%
  84   0.221067   0.207589   6.10%
  85   0.187011   0.165053   11.74%
  86   0.180443   0.109969   39.06%
  87   0.21067   0.178795   15.13%
  88   0.239355   0.185835   22.36%
  89   0.172395   0.143033   17.03%
  90   0.259115   0.163795   36.79%
  91   0.229833   0.157998   31.26%
  92   0.207775   0.149088   28.25%
  93   0.230449   0.144534   37.28%
  94   0.200673   0.121596   39.41%
  95   0.184476   0.151907   17.65%
  96   0.256026   0.185366   27.60%
  97   0.240359   0.198669   17.34%
  98   0.289224   0.208318   27.97%
  99   0.210961   0.171762   18.58%
  100   0.312019   0.203808   34.68%
  最小   6.01102%
  最大   47.411%
  平均   30.221%
附表四
  序号i   MSE(i)1   MSE(i)2   R(i)
  1   0.946877   0.946877   6.0392E-12
  2   0.946877   0.946877   2.6518E-10
  3   1.050438   0.946877   9.8588E-02
  4   0.946877   0.946877   6.3864E-12
  5   1.124862   0.963516   1.4344E-01
  6   1.153314   0.946877   1.7899E-01
  7   1.001377   0.946877   5.4425E-02
  8   0.970591   0.946877   2.4433E-02
  9   1.002329   0.99301   9.2983E-03
  10   1.07903   1.0333   4.2380E-02
  11   1.262371   1.150208   8.8851E-02
  12   0.946877   0.946877   1.5735E-10
  13   0.946877   0.946877   3.5137E-10
  14   0.946877   0.946877   5.5187E-11
  序号i   MSE(i)1   MSE(i)2   R(i)
  15   1.041388   0.946877   9.0755E-02
  16   1.165067   0.967223   1.6981E-01
  17   1.010929   0.970339   4.0151E-02
  18   1.104709   0.951307   1.3886E-01
  19   0.946877   0.946877   1.1357E-10
  20   1.099756   0.961092   1.2609E-01
  21   0.969597   0.946877   2.3432E-02
  22   1.007818   0.946877   6.0468E-02
  23   0.946877   0.946877   1.0301E-10
  24   0.946877   0.946877   2.7839E-11
  25   0.946877   0.946877   1.1132E-10
  26   1.10711   0.986092   1.0931E-01
  27   0.946877   0.946877   8.1351E-11
  28   1.091142   0.946877   1.3221E-01
  29   1.025362   1.024721   6.2581E-04
  30   1.159731   0.946877   1.8354E-01
  31   0.991997   0.940276   5.2138E-02
  32   0.946877   0.946877   6.1983E-11
  33   1.1141   1.028148   7.7149E-02
  34   0.946877   0.946877   3.7909E-11
  35   0.946877   0.946877   1.1297E-10
  36   1.078687   0.946877   1.2219E-01
  37   0.946877   0.946877   1.0031E-10
  38   0.946877   0.946877   5.6634E-11
  39   1.13359   0.946877   1.6471E-01
  40   0.946877   0.946877   9.9736E-11
  41   0.946877   0.946877   1.3337E-10
  42   0.946877   0.946877   1.5685E-10
  43   0.946877   0.946877   2.6871E-11
  44   1.136901   0.946877   1.6714E-01
  45   0.946877   0.946877   6.8552E-11
  46   0.990461   0.941239   4.9696E-02
  47   1.172618   1.151192   1.8271E-02
  48   0.946877   0.946877   8.5137E-11
  49   0.946877   0.946877   1.0347E-10
  50   0.946877   0.946877   7.2397E-11
  51   0.946877   0.946877   5.8018E-11
  52   0.946877   0.946877   1.8200E-10
  53   0.946877   0.946877   1.0222E-10
  54   1.054654   0.946877   1.0219E-01
  序号i   MSE(i)1   MSE(i)2   R(i)
  55   0.946877   0.946877   2.6545E-11
  56   0.968318   0.946877   2.2143E-02
  57   0.946877   0.946877   1.7539E-10
  58   0.946877   0.946877   6.6099E-11
  59   0.946877   0.946877   1.1318E-10
  60   1.091091   1.088644   2.2432E-03
  61   0.946877   0.946877   9.8471E-12
  62   1.138308   0.996108   1.2492E-01
  63   0.946877   0.946877   1.0855E-10
  64   0.946877   0.946877   3.9926E-11
  65   0.946877   0.946877   7.5123E-11
  66   0.946877   0.946877   1.8470E-10
  67   0.946877   0.946877   2.4358E-10
  68   0.946877   0.946877   2.1954E-10
  69   1.02679   0.985965   3.9760E-02
  70   0.946877   0.946877   8.8524E-11
  71   0.946877   0.946877   6.5527E-12
  72   1.065111   0.952657   1.0558E-01
  73   0.946877   0.946877   9.6193E-11
  74   0.946877   0.946877   1.2302E-10
  75   1.152983   0.946877   1.7876E-01
  76   0.946877   0.946877   1.0588E-10
  77   0.946877   0.946877   1.4119E-10
  78   1.132398   0.996344   1.2015E-01
  79   1.013142   0.961945   5.0533E-02
  80   0.946877   0.946877   5.8076E-11
  81   0.946877   0.946877   2.5037E-11
  82   1.05879   1.009608   4.6451E-02
  83   0.955985   0.946877   9.5267E-03
  84   0.946877   0.946877   8.9234E-11
  85   1.021301   0.946877   7.2872E-02
  86   1.167912   0.959524   1.7843E-01
  87   1.243226   0.917475   2.6202E-01
  88   1.096489   0.946877   1.3645E-01
  89   0.946877   0.946877   1.8468E-10
  90   1.077496   0.986051   8.4868E-02
  91   1.02596   0.974687   4.9976E-02
  92   1.079016   0.946877   1.2246E-01
  93   1.020752   1.013912   6.7015E-03
  94   0.946877   0.946877   4.4630E-11
  序号i   MSE(i)1   MSE(i)2   R(i)
  95   0.946877   0.946877   1.7118E-10
  96   0.946877   0.946877   4.1489E-11
  97   1.113406   0.967828   1.3075E-01
  98   1.086811   1.026833   5.5187E-02
  99   0.946877   0.946877   1.4507E-10
  100   0.946877   0.946877   1.6820E-10
  最小   6.0392E-12
  最大   2.6202E-01
  平均   4.2689E-02
附表五
  序号i   MSE(i)1   MSE(i)2   R(i)
  1   0.946877   0.061957   9.3457E-01
  2   0.946877   0.061957   9.3457E-01
  3   0.946877   0.097737   8.9678E-01
  4   0.946877   0.205583   7.8288E-01
  5   0.946877   0.061957   9.3457E-01
  6   0.946877   0.061957   9.3457E-01
  7   0.946877   0.061957   9.3457E-01
  8   1.006911   0.061957   9.3847E-01
  9   0.946877   0.061957   9.3457E-01
  10   0.946877   0.061957   9.3457E-01
  11   0.946877   0.061957   9.3457E-01
  12   0.946877   0.061957   9.3457E-01
  13   0.946877   0.061957   9.3457E-01
  14   0.946877   0.061957   9.3457E-01
  15   0.946877   0.061957   9.3457E-01
  16   0.946877   0.061957   9.3457E-01
  17   0.946877   0.061957   9.3457E-01
  18   0.946877   0.08256   9.1281E-01
  19   0.946877   0.061957   9.3457E-01
  20   0.946877   0.061957   9.3457E-01
  21   0.946877   0.061957   9.3457E-01
  22   0.946877   0.061957   9.3457E-01
  23   0.946877   0.061957   9.3457E-01
  24   1.000337   0.061957   9.3806E-01
  25   0.946877   0.117911   8.7547E-01
  26   0.946877   0.258412   7.2709E-01
  27   0.946877   0.061957   9.3457E-01
  28   0.946877   0.061957   9.3457E-01
  29   0.946877   0.061957   9.3457E-01
  序号i   MSE(i)1   MSE(i)2   R(i)
  30   0.946877   0.061957   9.3457E-01
  31   0.946877   0.061957   9.3457E-01
  32   1.005136   0.061957   9.3836E-01
  33   0.946877   0.061957   9.3457E-01
  34   0.946877   0.061957   9.3457E-01
  35   0.946877   0.208362   7.7995E-01
  36   0.946877   0.12353   8.6954E-01
  37   0.946877   0.061957   9.3457E-01
  38   0.946877   0.17125   8.1914E-01
  39   0.946877   0.061957   9.3457E-01
  40   0.946877   0.061957   9.3457E-01
  41   0.957566   0.061957   9.3530E-01
  42   1.166407   0.061957   9.4688E-01
  43   1.177388   0.061957   9.4738E-01
  44   0.946877   0.061957   9.3457E-01
  45   0.946877   0.061957   9.3457E-01
  46   0.946877   0.061957   9.3457E-01
  47   0.987934   0.061957   9.3729E-01
  48   0.946877   0.061957   9.3457E-01
  49   0.946877   0.061957   9.3457E-01
  50   0.946877   0.061957   9.3457E-01
  51   0.955368   0.061957   9.3515E-01
  52   0.946877   0.061957   9.3457E-01
  53   0.946877   0.061957   9.3457E-01
  54   0.946877   0.194769   7.9430E-01
  55   0.946877   0.061957   9.3457E-01
  56   0.946877   0.061957   9.3457E-01
  57   1.11138   0.061957   9.4425E-01
  58   0.946877   0.061957   9.3457E-01
  59   0.946877   0.061957   9.3457E-01
  60   0.946877   0.061957   9.3457E-01
  61   0.946877   0.061957   9.3457E-01
  62   0.946877   0.061957   9.3457E-01
  63   0.946877   0.061957   9.3457E-01
  64   0.946877   0.145697   8.4613E-01
  65   0.946877   0.061957   9.3457E-01
  66   0.946877   0.061957   9.3457E-01
  67   0.977634   0.061957   9.3663E-01
  68   0.946877   0.061957   9.3457E-01
  69   1.093068   0.072028   9.3411E-01
  序号i   MSE(i)1   MSE(i)2   R(i)
  70   1.002396   0.061957   9.3819E-01
  71   0.946877   0.061957   9.3457E-01
  72   0.946877   0.061957   9.3457E-01
  73   0.946877   0.061957   9.3457E-01
  74   0.946877   0.183999   8.0568E-01
  75   0.946877   0.061957   9.3457E-01
  76   0.946877   0.061957   9.3457E-01
  77   0.946877   0.273566   7.1109E-01
  78   0.946877   0.179146   8.1080E-01
  79   0.946877   0.061957   9.3457E-01
  80   0.946877   0.061957   9.3457E-01
  81   0.946877   0.061957   9.3457E-01
  82   0.946877   0.061957   9.3457E-01
  83   0.946877   0.061957   9.3457E-01
  84   0.946877   0.129456   8.6328E-01
  85   0.946877   0.061957   9.3457E-01
  86   0.967612   0.061957   9.3597E-01
  87   0.946877   0.061957   9.3457E-01
  88   0.946877   0.061957   9.3457E-01
  89   0.946877   0.061957   9.3457E-01
  90   0.946877   0.061957   9.3457E-01
  91   0.946877   0.214195   7.7379E-01
  92   0.946877   0.061957   9.3457E-01
  93   0.946877   0.061957   9.3457E-01
  94   0.969536   0.061957   9.3610E-01
  95   0.951677   0.061957   9.3490E-01
  96   0.946877   0.061957   9.3457E-01
  97   0.946877   0.061957   9.3457E-01
  98   0.946877   0.061957   9.3457E-01
  99   1.120744   0.061957   9.4472E-01
  100   0.946877   0.181525   8.0829E-01
  最小   7.1109E-01
  最大   9.4738E-01
  平均   9.1649E-01

Claims (6)

1. 一种电力系统负荷谐波稳态模型的建立方法,建立如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
其特征是,包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据。
2. 一种电力系统负荷谐波稳态模型的建立方法,建立如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
其特征是,包括预处理输出数据的步骤,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,…,INcosθNi]作为训练及建立模型的输出数据。
3. 一种电力系统负荷谐波稳态模型的建立方法,建立如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
其特征是,包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据;
还包括预处理输出数据的步骤,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,…,INcosθNi]作为训练及建立模型的输出数据。
4. 一种电力系统负荷谐波稳态模型的仿真方法,所仿真的模型为如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
模型建立过程包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据;
其特征是,仿真过程包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为在上述建模过程中新产生的黑箱模型的输入数据。
5. 一种电力系统负荷谐波稳态模型的仿真方法,所仿真的模型为如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN ∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
模型建立过程包括预处理输出数据的步骤,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,…,INcosθNi]作为训练及建立模型的输出数据;
其特征是,仿真过程包括输出数据的反变换的步骤,在按上述建模过程新产生的黑箱模型的输出数据输出前对其进行处理,将表达形式为[ReI,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]的数据转化为[I1,I2,…IN;θ1i,θ2i,…,θNi]形式的数据作为输出数据输出,其中,ReIh=Ihcosθhi,ImIh=Ihsinθhi  h=1,2,…,N。
6. 一种电力系统负荷谐波稳态模型的仿真方法,所仿真的模型为如下式所表示的电力系统负荷谐波稳态黑箱模型:
Ih∠θhi=Fh(V1∠θ1u,V2∠θ2u,…,VN∠θNu,C)  h=1,2,…,N
上式中,Ih∠θhi为负荷吸收的第h次谐波电流向量,Ih和θhi为负荷吸收的第h次谐波电流向量的幅值和相角;V1∠θ1u,V2∠θ2u,…,VN∠θNu为供电电压中的基波和各次谐波的以幅值和相角方式表示的电压向量;C为负荷的特征参数集合向量;h从1到N,代表从基波到N次谐波;
模型建立过程包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为训练及建立模型的输入数据;
模型建立过程还包括预处理输出数据的步骤,将所述负荷吸收的第h次谐波电流向量Ih∠θhi(h=1,2,…N)转化为复数直角坐标Ih·cosθhi,Ih·sinθhi(h=1,2,…,N),以[I1cosθ1i,I2cosθ2i,…,INcosθNi]作为训练及建立模型的输出数据;
其特征是,仿真过程包括预处理输入数据的步骤,将所述供电电压中的基波和各次谐波的电压向量V1∠θ1u,V2∠θ2u,…,VN∠θNu转化为复数直角坐标Vh·cosθhu,Vh·sinθhu(h=1,2,…,N),并以[V1·cosθ1u,V1·sinθ1u,V2·cosθ2u,V2·sinθ2u,…,VN·cosθNu,VN·sinθNu,C]作为在上述建模过程中新产生的黑箱模型的输入数据;
仿真过程还包括输出数据的反变换的步骤,在新产生的黑箱模型的输出数据输出前对其预处理,将表达形式为[ReI1,ReI2,…,ReIN;ImI1,ImI2,…,ImIN]的数据转化为[I1,I2,…IN;θ1i,θ2i,…,θNi]形式的数据作为输出数据输出,其中,ReIh=Ihcosθhi,ImIh=Ihsinθhi  h=1,2,…,N。
CN2007100373596A 2007-02-09 2007-02-09 电力系统负荷谐波稳态模型的建立方法及仿真方法 Expired - Fee Related CN101242101B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007100373596A CN101242101B (zh) 2007-02-09 2007-02-09 电力系统负荷谐波稳态模型的建立方法及仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007100373596A CN101242101B (zh) 2007-02-09 2007-02-09 电力系统负荷谐波稳态模型的建立方法及仿真方法

Publications (2)

Publication Number Publication Date
CN101242101A true CN101242101A (zh) 2008-08-13
CN101242101B CN101242101B (zh) 2011-07-20

Family

ID=39933387

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007100373596A Expired - Fee Related CN101242101B (zh) 2007-02-09 2007-02-09 电力系统负荷谐波稳态模型的建立方法及仿真方法

Country Status (1)

Country Link
CN (1) CN101242101B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101860027A (zh) * 2010-05-28 2010-10-13 北京中恒博瑞数字电力科技有限公司 供电分析与运行方式相结合的电网运行时模型建立方法
CN101894190A (zh) * 2010-07-14 2010-11-24 北京交通大学 牵引变电站谐波发射水平估计方法
CN102096747A (zh) * 2011-03-18 2011-06-15 长沙高新开发区德研电气技术有限公司 电力电子系统仿真方法和装置
CN102298655A (zh) * 2010-06-23 2011-12-28 中国科学院微电子研究所 一种soi体电阻建模方法
CN104362622A (zh) * 2014-11-11 2015-02-18 安顺供电局 一种带小水电机组地区电网的负荷模型和建模方法
CN106503325A (zh) * 2016-10-18 2017-03-15 哈尔滨工业大学 一种针对正多边形拓展互承结构几何构形的方法
CN108038546A (zh) * 2017-12-29 2018-05-15 百度在线网络技术(北京)有限公司 用于压缩神经网络的方法和装置
CN110210152A (zh) * 2019-06-06 2019-09-06 福州大学 一种超高次谐波源建模方法
CN110766138A (zh) * 2019-10-21 2020-02-07 中国科学院自动化研究所 基于脑发育机制的自适应神经网络模型的构建方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2861869A1 (fr) * 2003-11-05 2005-05-06 Michelin Soc Tech Procede d'extraction d'au moins une grandeur caracterisant un signal periodique, et application.
CN100361362C (zh) * 2005-08-23 2008-01-09 湖南大学 单独注入式有源电力滤波器的分频自适应控制方法
CN100370671C (zh) * 2006-03-24 2008-02-20 武汉大学 基于α-β电流分量直接注入的有源电力滤波方法及其装置
CN100444505C (zh) * 2006-07-06 2008-12-17 北京思源清能电气电子有限公司 一种三电平单相桥逆变器的特定消谐优化方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101860027A (zh) * 2010-05-28 2010-10-13 北京中恒博瑞数字电力科技有限公司 供电分析与运行方式相结合的电网运行时模型建立方法
CN101860027B (zh) * 2010-05-28 2012-09-26 北京中恒博瑞数字电力科技有限公司 供电分析与运行方式相结合的电网运行时模型建立方法
CN102298655A (zh) * 2010-06-23 2011-12-28 中国科学院微电子研究所 一种soi体电阻建模方法
CN101894190B (zh) * 2010-07-14 2013-04-24 北京交通大学 牵引变电站谐波发射水平估计方法
CN101894190A (zh) * 2010-07-14 2010-11-24 北京交通大学 牵引变电站谐波发射水平估计方法
CN102096747A (zh) * 2011-03-18 2011-06-15 长沙高新开发区德研电气技术有限公司 电力电子系统仿真方法和装置
CN102096747B (zh) * 2011-03-18 2013-02-27 长沙高新开发区德研电气技术有限公司 电力电子系统仿真方法和装置
CN104362622A (zh) * 2014-11-11 2015-02-18 安顺供电局 一种带小水电机组地区电网的负荷模型和建模方法
CN106503325A (zh) * 2016-10-18 2017-03-15 哈尔滨工业大学 一种针对正多边形拓展互承结构几何构形的方法
CN106503325B (zh) * 2016-10-18 2019-11-29 哈尔滨工业大学 一种针对正多边形拓展互承结构几何构形的方法
CN108038546A (zh) * 2017-12-29 2018-05-15 百度在线网络技术(北京)有限公司 用于压缩神经网络的方法和装置
CN110210152A (zh) * 2019-06-06 2019-09-06 福州大学 一种超高次谐波源建模方法
CN110210152B (zh) * 2019-06-06 2021-01-29 福州大学 一种超高次谐波源建模方法
CN110766138A (zh) * 2019-10-21 2020-02-07 中国科学院自动化研究所 基于脑发育机制的自适应神经网络模型的构建方法及系统

Also Published As

Publication number Publication date
CN101242101B (zh) 2011-07-20

Similar Documents

Publication Publication Date Title
CN101242101B (zh) 电力系统负荷谐波稳态模型的建立方法及仿真方法
Merabet et al. A comparative study of harmonic currents extraction by simulation and implementation
CN106874548A (zh) 一种基于双重傅里叶变换分析逆变器的方法
Hou et al. Adaptive fuzzy backstepping control of three-phase active power filter
CN103020385B (zh) 基于RTDS电网500kV主网建模仿真系统
CN105140907A (zh) 直流微网多智能体自适应下垂一致性协调控制方法及装置
Nascimento et al. Harmonic distortion monitoring for nonlinear loads using neural-network-method
CN105006839B (zh) 并网发电系统的弱电网网源荷阻抗模型分析方法
CN102870309A (zh) 便于设计高压(hvdc)控制系统的方法和系统、hvdc系统和优化hvdc系统的方法
Shukla et al. Implementation of empirical mode decomposition based algorithm for shunt active filter
Acikgoz et al. Performance analysis of electronic power transformer based on neuro-fuzzy controller
CN104978304A (zh) 电流连续模式分数阶开关变换器的符号分析方法及装置
Rizwan et al. Implementation of ANN-based embedded hybrid power filter using HIL-topology with real-time data visualization through node-RED
Abaza et al. Implementation of an electronically based active power filter associated with a digital controller for harmonics elimination and power factor correction
Li et al. Dynamic modeling for VSG cluster by using data-physical driven method
Tao et al. Optimal power flow via quadratic power flow
Yao et al. Modeling and simulation of VSC-HVDC with dynamic phasors
Osaloni et al. Integrated AC/DC systems power flow solution using Newton-Raphson and Broyden approaches
Tang et al. Application of hardware‐in‐loop in teaching power electronic course based on a low‐cost platform
CN107482624A (zh) 一种基于upfc的电网损耗控制策略
Tozak et al. Modeling, simulation, and experimental verification of a smart grid system
CN106952547A (zh) 并网光伏发电教学实验装置
CN204012706U (zh) 一种有源电力滤波器
Singh et al. A novel reduced-rule fuzzy logic based self-supported dynamic voltage restorer for mitigating diverse power quality problems
Mahdavi et al. A droop-based diesel generator emulator for microgrid frequency control studies

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: TONGJI UNIVERSITY

Free format text: FORMER OWNER: SHANGHAI MAGNETIC SUSPENSION COMMUNICATION ENGINEERING TECHNOLOGY RESEARCH CENTE

Effective date: 20131203

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 201203 PUDONG NEW AREA, SHANGHAI TO: 200092 YANGPU, SHANGHAI

TR01 Transfer of patent right

Effective date of registration: 20131203

Address after: 200092 Shanghai City, Yangpu District Siping Road No. 1239

Patentee after: Tongji University

Address before: 201203 No. 2520, Pudong New Area, Shanghai, Longyang Road

Patentee before: Shanghai Magnetic Suspension Communication Engineering Technology Research Cente

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110720

Termination date: 20160209

CF01 Termination of patent right due to non-payment of annual fee