【SVM分类】基于支持向量机实现电力节点分类附matlab代码

1 简介

2 部分代码

function mpc = case39%CASE39 Power flow data for 39 bus New England system.%   Please see CASEFORMAT for details on the case file format.%   Data taken from [1] with the following modifications/additions:%%       - renumbered gen buses consecutively (as in [2] and [4])%       - added Pmin = 0 for all gens%       - added Qmin, Qmax for gens at 31 & 39 (copied from gen at 35)%       - added Vg based on V in bus data (missing for bus 39)%       - added Vg, Pg, Pd, Qd at bus 39 from [2] (same in [4])%       - added Pmax at bus 39: Pmax = Pg + 100%       - added line flow limits and area data from [4]%       - added voltage limits, Vmax = 1.06, Vmin = 0.94%       - added identical quadratic generator costs%       - increased Pmax for gen at bus 34 from 308 to 508%         (assumed typo in [1], makes initial solved case feasible)%       - re-solved power flow %   Notes:%       - Bus 39, its generator and 2 connecting lines were added%         (by authors of [1]) to represent the interconnection with%         the rest of the eastern interconnect, and did not include%         Vg, Pg, Qg, Pd, Qd, Pmin, Pmax, Qmin or Qmax.%       - As the swing bus, bus 31 did not include and Q limits.%       - The voltages, etc in [1] appear to be quite close to the%         power flow solution of the case before adding bus 39 with%         it's generator and connecting branches, though the solution%         is not exact.%       - Explicit voltage setpoints for gen buses are not given, so%         they are taken from the bus data, however this results in two%         binding Q limits at buses 34 & 37, so the corresponding%         voltages have probably deviated from their original setpoints.%       - The generator locations and types are as follows:%           1   30      hydro%           2   31      nuke01%           3   32      nuke02%           4   33      fossil02%           5   34      fossil01%           6   35      nuke03%           7   36      fossil04%           8   37      nuke04%           9   38      nuke05%           10  39      interconnection to rest of US/Canada%%   MATPOWER​%% MATPOWER Case Format : Version 2mpc.version = '2';​%%-----  Power Flow Data  -----%%%% system MVA basempc.baseMVA = 100;​%% bus data%  bus_i  type  Pd  Qd  Gs  Bs  area  Vm  Va  baseKV  zone  Vmax  Vminmpc.bus = [  1  1  97.6  44.2  0  0  2  1.0393836  -13.536602  345  1  1.06  0.94;  2  1  0  0  0  0  2  1.0484941  -9.7852666  345  1  1.06  0.94;  3  1  322  2.4  0  0  2  1.0307077  -12.276384  345  1  1.06  0.94;  4  1  500  184  0  0  1  1.00446  -12.626734  345  1  1.06  0.94;  5  1  0  0  0  0  1  1.0060063  -11.192339  345  1  1.06  0.94;  6  1  0  0  0  0  1  1.0082256  -10.40833  345  1  1.06  0.94;  7  1  233.8  84  0  0  1  0.99839728  -12.755626  345  1  1.06  0.94;  8  1  522  176.6  0  0  1  0.99787232  -13.335844  345  1  1.06  0.94;  9  1  6.5  -66.6  0  0  1  1.038332  -14.178442  345  1  1.06  0.94;  10  1  0  0  0  0  1  1.0178431  -8.170875  345  1  1.06  0.94;  11  1  0  0  0  0  1  1.0133858  -8.9369663  345  1  1.06  0.94;  12  1  8.53  88  0  0  1  1.000815  -8.9988236  345  1  1.06  0.94;  13  1  0  0  0  0  1  1.014923  -8.9299272  345  1  1.06  0.94;  14  1  0  0  0  0  1  1.012319  -10.715295  345  1  1.06  0.94;  15  1  320  153  0  0  3  1.0161854  -11.345399  345  1  1.06  0.94;  16  1  329  32.3  0  0  3  1.0325203  -10.033348  345  1  1.06  0.94;  17  1  0  0  0  0  2  1.0342365  -11.116436  345  1  1.06  0.94;  18  1  158  30  0  0  2  1.0315726  -11.986168  345  1  1.06  0.94;  19  1  0  0  0  0  3  1.0501068  -5.4100729  345  1  1.06  0.94;  20  1  680  103  0  0  3  0.99101054  -6.8211783  345  1  1.06  0.94;  21  1  274  115  0  0  3  1.0323192  -7.6287461  345  1  1.06  0.94;  22  1  0  0  0  0  3  1.0501427  -3.1831199  345  1  1.06  0.94;  23  1  247.5  84.6  0  0  3  1.0451451  -3.3812763  345  1  1.06  0.94;  24  1  308.6  -92.2  0  0  3  1.038001  -9.9137585  345  1  1.06  0.94;  25  1  224  47.2  0  0  2  1.0576827  -8.3692354  345  1  1.06  0.94;  26  1  139  17  0  0  2  1.0525613  -9.4387696  345  1  1.06  0.94;  27  1  281  75.5  0  0  2  1.0383449  -11.362152  345  1  1.06  0.94;  28  1  206  27.6  0  0  3  1.0503737  -5.9283592  345  1  1.06  0.94;  29  1  283.5  26.9  0  0  3  1.0501149  -3.1698741  345  1  1.06  0.94;  30  2  0  0  0  0  2  1.0499  -7.3704746  345  1  1.06  0.94;  31  3  9.2  4.6  0  0  1  0.982  0  345  1  1.06  0.94;  32  2  0  0  0  0  1  0.9841  -0.1884374  345  1  1.06  0.94;  33  2  0  0  0  0  3  0.9972  -0.19317445  345  1  1.06  0.94;  34  2  0  0  0  0  3  1.0123  -1.631119  345  1  1.06  0.94;  35  2  0  0  0  0  3  1.0494  1.7765069  345  1  1.06  0.94;  36  2  0  0  0  0  3  1.0636  4.4684374  345  1  1.06  0.94;  37  2  0  0  0  0  2  1.0275  -1.5828988  345  1  1.06  0.94;  38  2  0  0  0  0  3  1.0265  3.8928177  345  1  1.06  0.94;  39  2  1104  250  0  0  1  1.03  -14.535256  345  1  1.06  0.94;];​%% generator data%  bus  Pg  Qg  Qmax  Qmin  Vg  mBase  status  Pmax  Pmin  Pc1  Pc2  Qc1min  Qc1max  Qc2min  Qc2max  ramp_agc  ramp_10  ramp_30  ramp_q  apfmpc.gen = [  30  250  161.762  400  140  1.0499  100  1  1040  0  0  0  0  0  0  0  0  0  0  0  0;  31  677.871  221.574  300  -100  0.982  100  1  646  0  0  0  0  0  0  0  0  0  0  0  0;  32  650  206.965  300  150  0.9841  100  1  725  0  0  0  0  0  0  0  0  0  0  0  0;  33  632  108.293  250  0  0.9972  100  1  652  0  0  0  0  0  0  0  0  0  0  0  0;  34  508  166.688  167  0  1.0123  100  1  508  0  0  0  0  0  0  0  0  0  0  0  0;  35  650  210.661  300  -100  1.0494  100  1  687  0  0  0  0  0  0  0  0  0  0  0  0;  36  560  100.165  240  0  1.0636  100  1  580  0  0  0  0  0  0  0  0  0  0  0  0;  37  540  -1.36945  250  0  1.0275  100  1  564  0  0  0  0  0  0  0  0  0  0  0  0;  38  830  21.7327  300  -150  1.0265  100  1  865  0  0  0  0  0  0  0  0  0  0  0  0;  39  1000  78.4674  300  -100  1.03  100  1  1100  0  0  0  0  0  0  0  0  0  0  0  0;];​%% branch data%  fbus  tbus  r  x  b  rateA  rateB  rateC  ratio  angle  status  angmin  angmaxmpc.branch = [  1  2  0.0035  0.0411  0.6987  600  600  600  0  0  1  -360  360;  1  39  0.001  0.025  0.75  1000  1000  1000  0  0  1  -360  360;  2  3  0.0013  0.0151  0.2572  500  500  500  0  0  1  -360  360;  2  25  0.007  0.0086  0.146  500  500  500  0  0  1  -360  360;  2  30  0  0.0181  0  900  900  2500  1.025  0  1  -360  360;  3  4  0.0013  0.0213  0.2214  500  500  500  0  0  1  -360  360;  3  18  0.0011  0.0133  0.2138  500  500  500  0  0  1  -360  360;  4  5  0.0008  0.0128  0.1342  600  600  600  0  0  1  -360  360;  4  14  0.0008  0.0129  0.1382  500  500  500  0  0  1  -360  360;  5  6  0.0002  0.0026  0.0434  1200  1200  1200  0  0  1  -360  360;  5  8  0.0008  0.0112  0.1476  900  900  900  0  0  1  -360  360;  6  7  0.0006  0.0092  0.113  900  900  900  0  0  1  -360  360;  6  11  0.0007  0.0082  0.1389  480  480  480  0  0  1  -360  360;  6  31  0  0.025  0  1800  1800  1800  1.07  0  1  -360  360;  7  8  0.0004  0.0046  0.078  900  900  900  0  0  1  -360  360;  8  9  0.0023  0.0363  0.3804  900  900  900  0  0  1  -360  360;  9  39  0.001  0.025  1.2  900  900  900  0  0  1  -360  360;  10  11  0.0004  0.0043  0.0729  600  600  600  0  0  1  -360  360;  10  13  0.0004  0.0043  0.0729  600  600  600  0  0  1  -360  360;  10  32  0  0.02  0  900  900  2500  1.07  0  1  -360  360;  12  11  0.0016  0.0435  0  500  500  500  1.006  0  1  -360  360;  12  13  0.0016  0.0435  0  500  500  500  1.006  0  1  -360  360;  13  14  0.0009  0.0101  0.1723  600  600  600  0  0  1  -360  360;  14  15  0.0018  0.0217  0.366  600  600  600  0  0  1  -360  360;  15  16  0.0009  0.0094  0.171  600  600  600  0  0  1  -360  360;  16  17  0.0007  0.0089  0.1342  600  600  600  0  0  1  -360  360;  16  19  0.0016  0.0195  0.304  600  600  2500  0  0  1  -360  360;  16  21  0.0008  0.0135  0.2548  600  600  600  0  0  1  -360  360;  16  24  0.0003  0.0059  0.068  600  600  600  0  0  1  -360  360;  17  18  0.0007  0.0082  0.1319  600  600  600  0  0  1  -360  360;  17  27  0.0013  0.0173  0.3216  600  600  600  0  0  1  -360  360;  19  20  0.0007  0.0138  0  900  900  2500  1.06  0  1  -360  360;  19  33  0.0007  0.0142  0  900  900  2500  1.07  0  1  -360  360;  20  34  0.0009  0.018  0  900  900  2500  1.009  0  1  -360  360;  21  22  0.0008  0.014  0.2565  900  900  900  0  0  1  -360  360;  22  23  0.0006  0.0096  0.1846  600  600  600  0  0  1  -360  360;  22  35  0  0.0143  0  900  900  2500  1.025  0  1  -360  360;  23  24  0.0022  0.035  0.361  600  600  600  0  0  1  -360  360;  23  36  0.0005  0.0272  0  900  900  2500  1  0  1  -360  360;  25  26  0.0032  0.0323  0.531  600  600  600  0  0  1  -360  360;  25  37  0.0006  0.0232  0  900  900  2500  1.025  0  1  -360  360;  26  27  0.0014  0.0147  0.2396  600  600  600  0  0  1  -360  360;  26  28  0.0043  0.0474  0.7802  600  600  600  0  0  1  -360  360;  26  29  0.0057  0.0625  1.029  600  600  600  0  0  1  -360  360;  28  29  0.0014  0.0151  0.249  600  600  600  0  0  1  -360  360;  29  38  0.0008  0.0156  0  1200  1200  2500  1.025  0  1  -360  360;];​%%-----  OPF Data  -----%%%% generator cost data%  1  startup  shutdown  n  x1  y1  ...  xn  yn%  2  startup  shutdown  n  c(n-1)  ...  c0mpc.gencost = [  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;  2  0  0  3  0.01  0.3  0.2;];​

3 仿真结果

4 参考文献

[1]屈玉涛, 邓万宇. 基于matlab的svm分类预测实现[J]. 信息通信, 2017(3):2.​

博主简介:擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真,相关matlab代码问题可私信交流。

部分理论引用网络文献,若有侵权联系博主删除。

文章出处登录后可见!

已经登录?立即刷新

共计人评分,平均

到目前为止还没有投票!成为第一位评论此文章。

(0)
上一篇 2022年5月23日 下午2:08
下一篇 2022年5月23日 下午2:14

相关推荐

本站注重文章个人版权,不会主动收集付费或者带有商业版权的文章,如果出现侵权情况只可能是作者后期更改了版权声明,如果出现这种情况请主动联系我们,我们看到会在第一时间删除!本站专注于人工智能高质量优质文章收集,方便各位学者快速找到学习资源,本站收集的文章都会附上文章出处,如果不愿意分享到本平台,我们会第一时间删除!