facat
/
pes2014-wronglowervoltagebound
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity
93 Commits 1 Branch 0 Tags 1 MiB
Commit Graph

93 Commits

Author SHA1 Message Date
dugg@lab-desk 02ddc534c7 修复了统计最大相角偏差的错误。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-03-13 08:35:32 +08:00
dugg@lab-desk 4afbbaf7fb a 2015-03-12 10:28:56 +08:00
dugg@lab-desk 94387f0bb8 a 2015-03-12 10:26:50 +08:00
dugg@lab-desk 7b43edabe7 a 2015-03-12 10:24:58 +08:00
dugg@lab-desk 9026198d61 fsdfs 2015-03-12 08:40:52 +08:00
dugg@lab-desk 32a9fba3ba fsdf 2015-03-12 08:39:07 +08:00
dugg@lab-desk d5585a6e6d dfs 2015-03-12 08:37:18 +08:00
dugg@lab-desk 27c7055747 f 2015-03-12 08:27:32 +08:00
dugg@lab-desk 5ab043255f s 2015-03-11 21:50:50 +08:00
dugg@lab-desk 584ff70634 f 2015-03-11 21:42:29 +08:00
dugg@lab-desk 1134b089ad f 2015-03-11 21:36:26 +08:00
dugg@lab-desk aa0e19f861 乱提交 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-03-11 21:02:17 +08:00
dugg@lab-desk 66fd308278 电压下界错误仿真。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-03-11 20:53:51 +08:00
dugg@lab-desk 1f9f3ff999 fasdf 2015-02-14 14:35:24 +08:00
dugg@lab-desk 5a752ca4fb 测试电压幅值下降错误的情况。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-10 15:55:19 +08:00
dugg@lab-desk 1c1e8f4b3c 可以做到从19节点开始没数据。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-07 15:40:19 +08:00
dugg@lab-desk 1bce6a7520 上一个提交提前停止迭代了，都没有收敛，互补约束不为0，现在这个可以了，加了互补变量的初始化。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-07 15:38:00 +08:00
dugg@lab-desk a964aa453d 做电压上下界错误的测试，这个参数是最好的。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-07 14:47:20 +08:00
dugg@lab-desk 1fd274a00d 修复一个电压没有测量值时上下界的错误。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-06 21:40:14 +08:00
dugg@lab-desk 2237f178a9 按时保存，好习惯。还在做测试。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-04 22:16:21 +08:00
dugg@lab-desk a091d9f502 1.加了一个画错误数据的图 
2.开始一点点做测试

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-03 21:52:51 +08:00
dugg@lab-desk beaa6855d4 又调了一下，更理想了。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-02 21:28:56 +08:00
dugg@lab-desk a2b2ba134f 把0-1变量变成0-0.1变量 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-02-02 21:23:04 +08:00
dugg@lab-desk 9a3fd4edd0 尝试了半天，这一组参数是最好的 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-01-29 21:21:03 +08:00
dugg@lab-desk e77f15c82e 互补理论调得不好。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-01-29 14:23:52 +08:00
dugg@lab-desk 72a068c4ff 修改了一下统计指标，按于尔铿的书上公式。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-01-24 22:27:51 +08:00
dugg@lab-desk 779c2d0e16 开始认真写中文期刊。 
1.把图表改成绝对值。
2.找到一个比较好的随机数值。

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2015-01-23 21:27:14 +08:00
dugg@lab-desk 466c90d8a7 1 加了算例A的收敛曲线展示。 
2 算离散结果时，重新初始化松弛变量。
3 结果对松弛变量的初衷也比较敏感。

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-09-09 16:19:29 +08:00
dugg@lab-desk 382763f559 1.为了多次计算，把原来的内点法循环单独放到一个文件中。 
2.没有量测量的时候，还是要约束一下的。

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-07-21 21:49:50 +08:00
dugg@lab-desk 0e24eed4bf 1.修复了这个地方的bug 
lPD=abs(RealPD*3*sigma);
uPD=abs(RealPD*3*sigma);
lQD=abs(RealQD*3*sigma);
uQD=abs(RealQD*3*sigma);
lVolt=abs(rVolt'*3*sigma);
uVolt=abs(rVolt'*3*sigma);
2.加入了考虑发电机的情况，其实也就是加入负数负荷。

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-07-20 21:08:16 +08:00
dugg@lab-desk a0a0d45e7d 加了真实的误差，收敛了。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-07-19 16:36:55 +08:00
dugg@lab-desk 16db56fb26 收敛了的，而且是做好了的。准备加真实的tolerance 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-07-19 16:30:16 +08:00
dugg@lab-desk 1599e58150 先算连续量，再算离散量，收敛了，而且算的结果是对的。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-25 12:13:10 +08:00
dugg@lab-desk c1727c877f 加了错误数据后，取一个比较大的互补变量的eps是可以收敛的，不过迭代次数多，达到29次。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-25 12:06:03 +08:00
dugg@lab-desk 570808f747 改变了一下Vbi,PDbi,QDbi的初值，收敛次数变为14~15次。 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-25 11:38:22 +08:00
dugg@lab-desk 240f989115 有把互补约束全部加上，收敛次数18~19次。 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
   (1-PDbi).*PDbi;
    (1-QDbi).*QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-25 11:35:26 +08:00
dugg@lab-desk 721035b580 去掉PDbi的互补约束，14次收敛 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
   % (1-PDbi).*PDbi;
    (1-QDbi).*QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-24 14:34:30 +08:00
dugg@lab-desk 00f34d732e 18~19次收敛 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
    (1-PDbi).*PDbi;
    (1-QDbi).*QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-24 14:18:29 +08:00
dugg@lab-desk f4d463540d 14次收敛 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
    (1-PDbi).*PDbi;
%     (1-QDbi).*QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-24 14:05:28 +08:00
dugg@lab-desk 9ac8691dcb 14次收敛 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
%     (1-PDbi).*PDbi;
%     (1-QDbi).*QDbi;
    ];
 2014-05-24 13:52:48 +08:00
dugg@lab-desk 431b0c8f78 1.删掉互补约束，看看是什么影响了收敛性。当前收敛次数为13 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
%     (1-Vbi).*Vbi;
%     (1-PDbi).*PDbi;
%     (1-QDbi).*QDbi;
    ];
2.更新了.gitignore

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-24 12:04:03 +08:00
dugg@lab-desk e27131d828 1.修复几个重大bug 
2.
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
    (1-PDbi).*PDbi;
    (1-QDbi).*QDbi;
    ];
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-22 11:56:06 +08:00
dugg@lab-desk 40c0a6706f Mat_G=[ 
sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
    ];
 2014-05-22 10:24:34 +08:00
dugg@lab-desk 305f059cf6 不等式约束 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-22 09:46:10 +08:00
dugg@lab-desk 6fcf53a485 等式约束 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi));
    sparse(QD(Loadi))-0.001;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-22 09:43:30 +08:00
dugg@lab-desk b7bb29928b 不等式约束 
Mat_G=[
    sparse(PD(Loadi));
    sparse(PD(Loadi))-0.001;
    sparse(QD(Loadi));
    sparse(QD(Loadi))-0.001;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-22 09:36:36 +08:00
dugg@lab-desk 57d87a95ae 加入.gitignore文件 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-21 15:47:55 +08:00
dugg@lab-desk f1a1f9be85 删掉所有asv文件 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-21 15:47:21 +08:00
dugg@lab-desk a14a61ce76 不等式是这样的 
Mat_G=[
    sparse(PD(Loadi));
    sparse(PD(Loadi))-0.001;
    sparse(QD(Loadi));
    sparse(QD(Loadi))-0.001;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-21 15:46:16 +08:00
dugg@lab-desk 28dfbb7664 修复func_deltG中的求导错误 
Signed-off-by: dugg@lab-desk <dugg@lab-desk>
 2014-05-21 15:42:16 +08:00