CN107436971A - 适用于非正定型相关性控制的改进拉丁超立方抽样方法 - Google Patents
适用于非正定型相关性控制的改进拉丁超立方抽样方法 Download PDFInfo
- Publication number
- CN107436971A CN107436971A CN201710549887.3A CN201710549887A CN107436971A CN 107436971 A CN107436971 A CN 107436971A CN 201710549887 A CN201710549887 A CN 201710549887A CN 107436971 A CN107436971 A CN 107436971A
- Authority
- CN
- China
- Prior art keywords
- matrix
- correlation
- sample
- positive definite
- latin hypercube
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明公开了一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,包括:S1,获得输入变量的累积分布函数或大量实测的离散数据及各变量间的相关系数矩阵;S2,分布函数已知时采用重要抽样法抽取样本,未知时采用改进的拉丁超立方抽样法,得到样本矩阵;S3,基于正定谱分解法,对相关系数矩阵进行修正;S4,对修正后随机顺序矩阵的相关系数矩阵,进行Cholesky分解和相关变换;S5,通过指定的相关系数矩阵计算顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵;S6,将样本矩阵带入节点,进行潮流计算,得到节点电压与支路功率,计算相对误差指标。本发明解决了输入变量分布函数未知和非正定矩阵无法分解的问题,扩大了拉丁超立方抽样方法的适用范围。
Description
技术领域
本发明涉及一种拉丁超立方抽样方法,具体涉及配电网仿真计算中适用于非正定型相关性控制的改进拉丁超立方抽样技术,属于分布式电源接入后的电力系统配电网随机潮流采样仿真技术领域。
背景技术
基于光伏、风力等可再生能源的分布式发电技术,能够有效地缓解能源紧张、环境污染等问题,得到了广泛的应用。但是配电网的潮流分布会受到分布式电源并网运行的影响。光伏发电和风力发电的输出功率由于自然环境因素的影响具有较强的随机性。当分布式电源采取多点并网时,多个电源之间的相关性问题也需要考虑。例如实际中光伏出力由于同地域光照的强相关性往往呈现正相关。同理,分布式风电中相邻风机一般也呈现正相关。风光由于夜晚风大、白天风小呈现互补性,即负相关。
模拟法是计算随机潮流的有效方法,拉丁超立方抽样法作为改进的模拟法,通过分层抽样,使得样本点均匀、全面地覆盖到变量的分布范围,有效地减小方差,拉丁超立方抽样法包括抽样和相关性控制两个部分。目前抽样方法主要包括区间随机抽样、点阵抽样和重要抽样等。但是这些方法都要求已知输入变量的累积分布函数,无法解决CDF较难准确获得或未知的情况。相关性控制的方法主要包括:遗传算法、模拟退火算法和Cholesky分解。Cholesky分解,计算简便速度快,但仅适用于相关系数矩阵为正定矩阵的情况。一般在相关性控制中将相关系数矩阵处理为正定矩阵。但随着分布式电源并网规模的增加,矩阵的维数增加,正定的可能性大大降低,这种处理方式不再合适。因此本发明提出一种基于修正相关系数矩阵的改进拉丁超立方抽样法。有效解决了输入变量分布函数未知和非正定矩阵无法分解的问题,扩大了拉丁超立方抽样法的适用范围。
发明内容
本发明正是针对现有技术中存在的技术问题,提供一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,该方法解决了输入变量的累积分布函数较难准确获得或未知的情况,以及相关性控制中非正定矩阵无法cholesky分解的问题,扩大了拉丁超立方抽样法的适用范围。
为了解决上述技术问题,本发明的技术方案如下:一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,包括:
步骤S1,首先获得输入变量的累积分布函数或输入变量大量实测的离散数据及各变量间的相关系数矩阵;
步骤S2,分布函数已知时采用重要抽样法抽取样本,未知时采用改进的拉丁超立方抽样法抽取样本,得到样本矩阵;
步骤S3,基于正定谱分解法,对相关系数矩阵进行修正;
步骤S4,对修正后的随机顺序矩阵的相关系数矩阵,进行Cholesky分解和相关变换;
步骤S5,对指定的相关系数矩阵修正分解后确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵;
步骤S6,将样本矩阵带入节点,进行潮流计算,得到节点电压与支路功率,计算相对误差指标。
进一步地,所述步骤S1中具体为:确定各节点接入的新能源类型,当输入变量的累积分布函数已知时,建立输出功率的模型,确定输入变量的累积分布函数;当输入变量的累积分布函数无法准确获得或者未知时,获取该变量大量实测的离散数据,计算各变量间的相关系数矩阵。
进一步地,所述步骤S2中具体为:分布函数已知时采用重要抽样法抽取样本,抽样公式为:分布函数未知时采用改进的拉丁超立方抽样法抽取样本,步骤为:对每个变量的M个离散样本数据按照从小到大排序,得到排序后的矩阵T'N×M;将累积分布函数的取值范围其均匀分为A等份;T”N×M即为抽样后的样本,T”k,a计算公式:T”k,a=T'k,j j=round((a·M)/A)。
进一步地,所述步骤S3中具体为:所述步骤S3具体为:基于正定谱分解法,对相关系数矩阵进行修正,修正步骤为:求得非正定矩阵C的特征值与特征向量后,重新定义对角矩阵Λ',根据变换后得到修正矩阵C',具体步骤如下:1)C为n阶方阵,λi与列向量Si为它的一个特征值与相应的特征向量,满足:
CSi=λiSi
Si=[si1,si2,...sin]T;
2)选取一个很小的正值ε(取值将影响误差大小,为保证算法的精度,一般小于10-4),将对角矩阵Λ重新定义为对角矩阵Λ':
Λ=diag(λi)
3)根据特征向量矩阵S,定义对角标度矩阵T:
4)根据变换得到修正矩阵C':
C'=BBT
C'即为修正后的正定矩阵。
进一步地,所述步骤S4中具体为:用改进的Spearman秩相关系数,用样本顺序代替样本值,进行相关性控制,对修正后的随机顺序矩阵相关系数矩阵进行Cholesky分解,得到下三角矩阵,变换后得到D,将初始样本中参照D排序,此时的相关性大大减小。
进一步地,所述步骤S5中具体为:对指定的相关系数矩阵Pset修正分解,得到下三角矩阵Qset,再对D变换:Dset=QsetD,求得Dset;将顺序矩阵的每一行元素按照Dset中的每一行元素重新排序,确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵。
进一步地,所述步骤S6中具体为:将样本矩阵代入节点,进行随机潮流计算,得到各节点电压与支路功率,得出分布情况,根据 计算相对误差指标。
相对于现有技术,本发明的有益效果如下:1)本发明提供的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,以变量的离散数据为基础改进抽样方法,解决了输入变量分布函数未知的问题,基于正定谱分解法,通过对相关系数矩阵进行简单的变换修正,克服了非正定矩阵无法分解的问题,仿真证明修正后的矩阵为正定矩阵,均可进行Cholesky分解,计算速度快且误差较小;该方法由仅适用于分布函数已知推广到分布函数未知,由仅适用于正定推广到非正定,扩大了拉丁超立方采样法的适用范围;2)本发明提供的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,该算法得到输出功率的样本的准确性和收敛性要优于蒙特卡罗模拟法。通过电压、功率等输出变量的误差分析,可以看出在随机潮流的计算中,该算法的准确性和收敛性都要优于蒙特卡罗模拟法;相同计算精度下需要的规模和时间远小于蒙特卡罗模拟法。
附图说明
图1为本发明所述的改进拉丁超立方抽样方法的流程图;
图2所示为本发明实施案例一中两种算法下69节点输出变量的相对误差指标;
其中:MCSRS表示蒙特卡罗模拟法,MCILHS表示拉丁超立方抽样法。
具体实施方式
本发明提供一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其方法流程如图1所示,由图1可知,该方法包括:
步骤S1,首先获得输入变量的累积分布函数或输入变量大量实测的离散数据及各变量间的相关系数矩阵;
步骤S2,分布函数已知时采用重要抽样法抽取样本,未知时采用改进的拉丁超立方抽样法抽取样本,得到样本矩阵;
步骤S3,基于正定谱分解法,对相关系数矩阵进行修正;
步骤S4,对修正后的随机顺序矩阵的相关系数矩阵,进行Cholesky分解和相关变换;
步骤S5,对指定的相关系数矩阵修正分解后确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵;
步骤S6,将样本矩阵带入节点,进行潮流计算,得到节点电压与支路功率,计算相对误差指标。
所述步骤S1中具体为:确定各节点接入的新能源类型,当输入变量的累积分布函数已知时,建立输出功率的模型,确定输入变量的累积分布函数;当输入变量的累积分布函数无法准确获得或者未知时,获取该变量大量实测的离散数据,计算各变量间的相关系数矩阵。
所述步骤S2中具体为:分布函数已知时采用重要抽样法抽取样本,抽样公式为:分布函数未知时采用改进的拉丁超立方抽样法抽取样本,步骤为:对每个变量的M个离散样本数据按照从小到大排序,得到排序后的矩阵T'N×M;将累积分布函数的取值范围其均匀分为A等份;T”N×M即为抽样后的样本,T”k,a计算公式:T”k,a=T'k,j j=round((a·M)/A)。
所述步骤S3中具体为:所述步骤S3具体为:基于正定谱分解法,对相关系数矩阵进行修正,修正步骤为:求得非正定矩阵C的特征值与特征向量后,重新定义对角矩阵Λ',根据变换后得到修正矩阵C',具体步骤如下:1)C为n阶方阵,λi与列向量Si为它的一个特征值与相应的特征向量,满足:
CSi=λiSi
Si=[si1,si2,...sin]T;
2)选取一个很小的正值ε(取值将影响误差大小,为保证算法的精度,一般小于10-4),将对角矩阵Λ重新定义为对角矩阵Λ':
Λ=diag(λi)
3)根据特征向量矩阵S,定义对角标度矩阵T:
4)根据变换得到修正矩阵C':
C'=BBT
C'即为修正后的正定矩阵。
所述步骤S4中具体为:用改进的Spearman秩相关系数,用样本顺序代替样本值,进行相关性控制,对修正后的随机顺序矩阵相关系数矩阵进行Cholesky分解,得到下三角矩阵,变换后得到D,将初始样本中参照D排序,此时的相关性大大减小。
所述步骤S5中具体为:对指定的相关系数矩阵Pset修正分解,得到下三角矩阵Qset,再对D变换:Dset=QsetD,求得Dset;将顺序矩阵的每一行元素按照Dset中的每一行元素重新排序,确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵。
所述步骤S6中具体为:将样本矩阵代入节点,进行随机潮流计算,得到各节点电压与支路功率,得出分布情况,根据 计算相对误差指标。
实施例1:
本发明提供一种适用于非正定型相关性控制的改进拉丁超立方抽样方法的实施案例中:
步骤S1中累积分布函数和离散数据矩阵的确定包括:
在PG&E69节点系统的5个节点接入分布式电源。以光伏为例,光伏模型采用beta分布模型,形状参数的选取为:α=0.9,β=0.85,容量为100kVA。确定输入变量的累积分布函数为5个光伏输出功率样本之间的相关系数矩阵为P5。
步骤S2中分布函数已知时采用重要抽样法抽取样本包括:
抽样原理第k个变量第a个值为:当x=xka时,Zk(xka)具有中的最大值。
步骤S3中基于正定谱分解法,对相关系数矩阵进行修正包括:
C为n阶方阵,λi与列向量Si为它的一个特征值与相应的特征向量。选取一个很小的正值ε,将对角矩阵Λ重新定义为对角矩阵Λ':Λ':
根据特征向量矩阵S,定义对角标度矩阵T:T:
根据变换得到修正矩阵C':C'=BBT,
修正误差有关于所有元素和特征值两种:
修正前后的相关系数矩阵和误差为:
步骤S4中对修正后的随机顺序矩阵的相关系数矩阵,进行Cholesky分解和相关变换包括:
对输入变量X1,X2,…,XN按照分布情况抽样,规模为1*M,得到样本矩阵FN*M。生成随机的顺序矩阵AN*M,求出AN*M的相关系数矩阵PA,对PA进行Cholesky分解,得下三角矩阵Q,按D=Q-1A对顺序矩阵进行变换。将初始样本FN*M中的各行元素参照D中对应的各行元素大小重新排序,求出FD。此时FD的相关性削弱。
步骤S5中对指定的相关系数矩阵修正分解后确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵包括:
当输入变量间需要达到必要的相关系数时,令指定的相关系数矩阵为Pset,对Pset进行Cholesky分解,得到下三角矩阵Qset。Dset=QsetD,求得Dset;将顺序矩阵A的每一行元素按照Dset中的每一行元素重新排序,得到新的顺序矩阵Aset,将初始样本F按Aset排序,得到Fset,Fset的相关系数满足设定的要求。
步骤S6中将样本矩阵带入节点,进行潮流计算,得到节点电压与支路功率,计算相对误差指标包括:
样本拟合的效果用相对误差指标来表示,计算公式为:
c表示形状参数α,β,cb表示给定值。
为表示计算精确度,用平均值与标准差的相对误差指标计算,公式如下:
U表示电压,与之对应的是功率P。为选取的电压的基准值。图2为本发明实施案例一中两种算法下69节点输出变量的相对误差指标。
以上对本申请进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐释,以上实施例的说明只是用于帮助理解本申请的方法及核心思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明不应理解为对本申请的限制。
Claims (7)
1.一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述方法包括以下步骤:
步骤S1,首先获得输入变量的累积分布函数或输入变量大量实测的离散数据及各变量间的相关系数矩阵;
步骤S2,分布函数已知时采用重要抽样法抽取样本,未知时采用改进的拉丁超立方抽样法抽取样本,得到样本矩阵;
步骤S3,基于正定谱分解法,对相关系数矩阵进行修正;
步骤S4,对修正后的随机顺序矩阵的相关系数矩阵,进行Cholesky分解和相关变换;
步骤S5,对指定的相关系数矩阵修正分解后确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵;
步骤S6,将样本矩阵带入节点,进行潮流计算,得到节点电压与支路功率,计算相对误差指标。
2.如权利要求1所述的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述步骤S1具体为:确定各节点接入的新能源类型,当输入变量的累积分布函数已知时,建立输出功率的模型,确定输入变量的累积分布函数;当输入变量的累积分布函数无法准确获得或者未知时,获取该变量大量实测的离散数据,计算各变量间的相关系数矩阵。
3.如权利要求1所述的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述步骤S2具体为:分布函数已知时采用重要抽样法抽取样本,抽样公式为:其中,Zk为输入随机变量的概率密度函数,A为抽样的分层的总层数,xk,a和yk,a是第a层左边界和右边界,表示求反函数;分布函数未知时采用改进的拉丁超立方抽样法抽取样本,步骤为:对每个变量的M个离散样本数据按照从小到大排序,得到排序后的矩阵T′N×M;将累积分布函数的取值范围其均匀分为A等份,即M为A的整数倍;T″N×M即为抽样后的样本,T″k,a计算公式:T″k,a=T′k,j j=round((a·M)/A)。
4.如权利要求1所述的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述步骤S3具体为:基于正定谱分解法,对相关系数矩阵进行修正,修正步骤为:求得非正定矩阵C的特征值与特征向量后,重新定义对角矩阵Λ',根据变换后得到修正矩阵C',具体步骤如下:1)C为n阶方阵,λi与列向量Si为它的一个特征值与相应的特征向量,满足:
CSi=λiSi
Si=[si1,si2,...sin]T;
2)选取一个很小的正值ε(取值将影响误差大小,为保证算法的精度,一般小于10-4),将对角矩阵Λ重新定义为对角矩阵Λ':
Λ=diag(λi)
Λ':
3)根据特征向量矩阵S,定义对角标度矩阵T:
T:
4)根据变换得到修正矩阵C':
C'=BBT
<mrow>
<mi>B</mi>
<mo>=</mo>
<msqrt>
<mi>T</mi>
</msqrt>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mi>S</mi>
<msqrt>
<msup>
<mi>&Lambda;</mi>
<mo>&prime;</mo>
</msup>
</msqrt>
</mrow>
C'即为修正后的正定矩阵。
5.如权利要求1所述的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述步骤S4具体为:采用改进的Spearman秩相关系数,以样本顺序代替样本值进行相关性控制,对修正后的随机顺序矩阵相关系数矩阵进行Cholesky分解,得到下三角矩阵,变换后得到D,将初始样本参照D排序,此时的相关性大大减小。
6.如权利要求1所述的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述步骤S5具体为:对指定的相关系数矩阵Pset修正分解,得到下三角矩阵Qset,再对D变换:Dset=QsetD,求得Dset;将顺序矩阵的每一行元素按照Dset中的每一行元素重新排序,确定最终的顺序矩阵,按顺序矩阵排序,确定最终的样本矩阵。
7.如权利要求1所述的一种适用于非正定型相关性控制的改进拉丁超立方抽样方法,其特征在于,所述步骤S6具体为:将样本矩阵代入节点,进行随机潮流计算,得到各节点电压与支路功率,得出分布情况,根据 计算相对误差指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710549887.3A CN107436971B (zh) | 2017-07-07 | 2017-07-07 | 适用于非正定型相关性控制的改进拉丁超立方抽样方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710549887.3A CN107436971B (zh) | 2017-07-07 | 2017-07-07 | 适用于非正定型相关性控制的改进拉丁超立方抽样方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107436971A true CN107436971A (zh) | 2017-12-05 |
CN107436971B CN107436971B (zh) | 2020-10-02 |
Family
ID=60459795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710549887.3A Active CN107436971B (zh) | 2017-07-07 | 2017-07-07 | 适用于非正定型相关性控制的改进拉丁超立方抽样方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107436971B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108123437A (zh) * | 2017-12-08 | 2018-06-05 | 云南电网有限责任公司 | 一种用交替投影法修正配电网相关性模型的随机潮流算法 |
CN109494747A (zh) * | 2018-11-09 | 2019-03-19 | 杭州电子科技大学 | 一种基于交替梯度算法的电网概率潮流计算方法 |
CN109598454A (zh) * | 2019-01-22 | 2019-04-09 | 广东电网有限责任公司 | 一种温控负荷可调度容量分析方法、装置和设备 |
CN109657385A (zh) * | 2018-12-26 | 2019-04-19 | 国网江苏省电力有限公司 | 一种适用于输电线路随机状态生成的仿真计算方法 |
CN110334366A (zh) * | 2019-03-14 | 2019-10-15 | 华电电力科学研究院有限公司 | 一种基于运用拉丁超立方抽样的蒙特卡罗法的建筑瞬时冷负荷预测方法 |
CN111881591A (zh) * | 2020-07-31 | 2020-11-03 | 贵州大学 | 基于svd的大规模相关性随机变量两次排序抽样方法 |
CN111985092A (zh) * | 2020-07-30 | 2020-11-24 | 哈尔滨工业大学 | 一种智能汽车仿真测试矩阵生成方法 |
CN113408061A (zh) * | 2021-07-08 | 2021-09-17 | 中汽院智能网联科技有限公司 | 基于改进拉丁超立方抽样的虚拟驾驶场景要素重组方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104682387A (zh) * | 2015-03-10 | 2015-06-03 | 东南大学 | 一种基于多区域交互迭代的概率潮流计算方法 |
CN106229986A (zh) * | 2016-08-30 | 2016-12-14 | 上海交通大学 | 一种电力系统的概率潮流计算方法 |
-
2017
- 2017-07-07 CN CN201710549887.3A patent/CN107436971B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104682387A (zh) * | 2015-03-10 | 2015-06-03 | 东南大学 | 一种基于多区域交互迭代的概率潮流计算方法 |
CN106229986A (zh) * | 2016-08-30 | 2016-12-14 | 上海交通大学 | 一种电力系统的概率潮流计算方法 |
Non-Patent Citations (4)
Title |
---|
M.N.KABIR: ""Probabilistic load flow for distribution systems with uncertain PV generation"", 《APPLIED ENERGY》 * |
云永琥 等: ""改进Kriging的热结构耦合梁共振非概率可靠性分析"", 《哈尔滨工业大学学报》 * |
李俊芳: ""基于进化算法改进拉丁超立方抽样的概率潮流计算"", 《中国电机工程学报》 * |
黄煜 等: ""基于拉丁超立方采样技术的半不变量法随机潮流计算"", 《电力自动化设备》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108123437A (zh) * | 2017-12-08 | 2018-06-05 | 云南电网有限责任公司 | 一种用交替投影法修正配电网相关性模型的随机潮流算法 |
CN109494747A (zh) * | 2018-11-09 | 2019-03-19 | 杭州电子科技大学 | 一种基于交替梯度算法的电网概率潮流计算方法 |
CN109494747B (zh) * | 2018-11-09 | 2021-09-28 | 杭州电子科技大学 | 一种基于交替梯度算法的电网概率潮流计算方法 |
CN109657385A (zh) * | 2018-12-26 | 2019-04-19 | 国网江苏省电力有限公司 | 一种适用于输电线路随机状态生成的仿真计算方法 |
CN109598454A (zh) * | 2019-01-22 | 2019-04-09 | 广东电网有限责任公司 | 一种温控负荷可调度容量分析方法、装置和设备 |
CN110334366A (zh) * | 2019-03-14 | 2019-10-15 | 华电电力科学研究院有限公司 | 一种基于运用拉丁超立方抽样的蒙特卡罗法的建筑瞬时冷负荷预测方法 |
CN111985092A (zh) * | 2020-07-30 | 2020-11-24 | 哈尔滨工业大学 | 一种智能汽车仿真测试矩阵生成方法 |
CN111881591A (zh) * | 2020-07-31 | 2020-11-03 | 贵州大学 | 基于svd的大规模相关性随机变量两次排序抽样方法 |
CN113408061A (zh) * | 2021-07-08 | 2021-09-17 | 中汽院智能网联科技有限公司 | 基于改进拉丁超立方抽样的虚拟驾驶场景要素重组方法 |
CN113408061B (zh) * | 2021-07-08 | 2023-05-05 | 中汽院智能网联科技有限公司 | 基于改进拉丁超立方抽样的虚拟驾驶场景要素重组方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107436971B (zh) | 2020-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107436971A (zh) | 适用于非正定型相关性控制的改进拉丁超立方抽样方法 | |
CN108306303B (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
Xu et al. | Robust dispatch of high wind power-penetrated power systems against transient instability | |
CN110738344B (zh) | 电力系统负荷预测的分布式无功优化方法及装置 | |
CN111008728A (zh) | 一种用于分布式光伏发电系统短期出力的预测方法 | |
CN103106544B (zh) | 一种基于t-s型模糊神经网络的光伏发电预测系统 | |
CN110556847B (zh) | 含光伏配电网中储能系统规划运行联合优化方法及系统 | |
CN105207573A (zh) | 基于离散概率模型的风光互补供电系统定量优化配置方法 | |
CN105512808A (zh) | 一种面向大数据的电力系统暂态稳定评估方法 | |
CN107834593B (zh) | 一种下垂控制型孤岛微电网静态电压稳定概率评估方法 | |
CN103986193B (zh) | 一种最大风电并网容量获取的方法 | |
CN108599239A (zh) | 一种下垂控制型孤岛微电网电压质量概率评估方法 | |
CN106786606B (zh) | 一种基于多种随机变量的电力系统概率潮流的计算方法 | |
CN106532772A (zh) | 一种基于改进正交优化群智能算法的分布式电源规划方法 | |
CN106169910A (zh) | 基于群搜索算法的光伏电池参数辨识方法 | |
CN107069708B (zh) | 一种基于极限学习机的输电网线路有功安全校正方法 | |
CN113328467B (zh) | 一种概率电压稳定性评估方法、系统、终端设备及介质 | |
CN102709955A (zh) | 一种基于多断面潮流控制的方法 | |
CN113536694A (zh) | 综合能源系统鲁棒优化运行方法、系统、装置及存储介质 | |
CN107276093A (zh) | 基于场景削减的电力系统概率潮流计算方法 | |
CN110336322B (zh) | 基于日最小负荷置信区间的光伏发电准入容量确定方法 | |
CN108460177A (zh) | 一种大规模多状态串并联系统的可靠性近似计算方法 | |
Cheng et al. | Short-term power load forecasting based on improved online ELM-K | |
CN115173421A (zh) | 一种基于渐进拉丁超立方采样的概率最优潮流计算方法 | |
CN107086603A (zh) | 一种含dfig电力系统随机模糊连续潮流方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |