潮流计算Word下载.docx
- 文档编号:20773876
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:15
- 大小:16.66KB
潮流计算Word下载.docx
《潮流计算Word下载.docx》由会员分享,可在线阅读,更多相关《潮流计算Word下载.docx(15页珍藏版)》请在冰豆网上搜索。
Y(i,j)=-y(i,j);
%求自导纳
Y(i,i)=sum(y(i,:
));
Y%Y为导纳矩阵
G=real(Y);
B=imag(Y);
y(1,2)=-1/(0.04+0.25i);
y(1,3)=-1/(0.1+0.35i);
y(1,4)=0;
y(1,5)=0;
y(2,3)=-1/(0.08+0.03i);
y(2,4)=-1/0.05i;
y(2,5)=0;
y(3,4)=0;
y(3,5)=-1/0.03i;
%原始节点功率
S
(1)=1.6+0.8i;
S
(2)=2+0.1i;
S(3)=3.7+1.3i;
S(4)=0;
S(5)=0;
P=real(S);
Q=imag(S);
%赋初值
U=ones(1,5);
U(5)=1.05;
e=zeros(1,5);
ox=ones(8,1);
fx=ones(8,1);
count=0
%计算迭代次数
whilemax(fx)>
1e-5
4
H(i,j)=0;
N(i,j)=0;
M(i,j)=0;
L(i,j)=0;
oP(i)=0;
oQ(i)=0;
5oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
oP(i)=oP(i)+P(i);
oQ(i)=oQ(i)+Q(i);
fx=[oP,oQ]'
;
%求雅克比矩阵
%当i~=j时候求H,N,M,L如下:
ifi~=jH(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
L(i,j)=H(i,j);
M(i,j)=-N(i,j);
H,N,M,L
%当i=j时H,N,M,L如下:
H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
N(i,i)=N(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
N(i,i)=N(i,i)-2*(U(i))^2*G(i,i);
L(i,i)=L(i,i)+2*(U(i))^2*B(i,i);
J=[H,N;
M,L]%J为雅克比矩阵
ox=-((inv(J))*fx);
oe(i)=ox(i);
oU(i)=ox(i+4)*U(i);
e(i)=e(i)+oe(i);
U(i)=U(i)+oU(i);
count=count+1;
ox,U,e,count
%求节点注入的净功率
i=5;
P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)))+P(i);
Q(i)=U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)))+Q(i);
S(5)=P(5)+Q(5)*sqrt(-1);
S
%求节点注入电流
I=Y*U'
U=ones(4,5);
e=zeros(4,5);
count=0%计算迭代次数
Y=
9.7907-31.4002i-0.6240+3.9002i-1.6667+5.0000i-2.5000+7.5000i-5.0000+15.0000i
-0.6240+3.9002i11.8740-37.6502i-10.0000+30.0000i0-1.2500+3.7500i
-1.6667+5.0000i-10.0000+30.0000i12.9167-38.7500i-1.2500+3.7500i0
-2.5000+7.5000i0-1.2500+3.7500i3.7500-11.2500i0
-5.0000+15.0000i-1.2500+3.7500i006.2500-18.7500i
count=
0
H=
03.90025.00007.5000
3.9002030.00000
5.000030.000003.7500
7.500003.75000
N=
00.62401.66672.5000
0.6240010.00000
1.666710.000001.2500
2.500001.25000
M=
0-0.6240-1.6667-2.5000
-0.62400-10.00000
-1.6667-10.00000-1.2500
-2.50000-1.25000
L=
J=
-32.15023.90025.00007.5000-9.54070.62401.66672.5000
3.9002-37.837730.000000.6240-11.811510.00000
5.000030.0000-38.75003.75001.666710.0000-12.91671.2500
7.500003.7500-11.25002.500001.2500-3.7500
10.0407-0.6240-1.6667-2.5000-30.65023.90025.00007.5000
-0.624011.9365-10.000003.9002-37.462730.00000
-1.6667-10.000012.9167-1.25005.000030.0000-38.75003.7500
-2.50000-1.25003.75007.500003.7500-11.2500
06.37667.866612.0338
6.8838058.38530
9.439059.623707.3107
12.769106.48620
02.64585.24285.2368
-0.5241021.52570
0.525717.810601.0626
3.030903.53630
0-2.6458-5.2428-5.2368
0.52410-21.52570
-0.5257-17.81060-1.0626
-3.03090-3.53630
-47.18346.37667.866612.0338-16.63832.64585.24285.2368
6.8838-70.937358.38530-0.5241-25.486521.52570
9.439059.6237-76.37347.31070.525717.8106-32.76821.0626
12.769106.4862-19.25533.030903.5363-6.4054
14.5595-2.6458-5.2428-5.2368-52.87226.37667.866612.0338
0.524120.0511-21.525706.8838-73.453258.38530
-0.5257-17.810619.3989-1.06269.439059.6237-80.12797.3107
-3.03090-3.53636.567212.769106.4862-19.6624
05.45976.856910.2936
5.7778047.68890
7.832848.382605.9781
10.772605.48390
01.89283.91214.2294
-0.0948017.05250
0.984514.971301.1691
2.792602.65160
0-1.8928-3.9121-4.2294
0.09480-17.05250
-0.9845-14.97130-1.1691
-2.79260-2.65160
-41.93725.45976.856910.2936-14.86151.89283.91214.2294
5.7778-58.642347.68890-0.0948-20.592117.05250
7.832848.3826-62.19345.97810.984514.9713-25.23441.1691
10.772605.4839-16.25652.792602.6516-5.4175
11.9311-1.8928-3.9121-4.2294-43.99055.45976.856910.2936
0.094816.5571-17.052505.7778-59.150447.68890
-0.9845-14.971317.1249-1.16917.832848.3826-64.88455.9781
-2.79260-2.65165.444210.772605.4839-16.3284
05.31866.697010.0289
5.6161046.17140
7.603746.775905.7849
10.482005.33130
01.80463.74354.0981
-0.0550016.39800
1.023414.584501.1724
2.738902.53300
0-1.8046-3.7435-4.0981
0.05500-16.39800
-1.0234-14.58450-1.1724
-2.73890-2.53300
-41.14555.31866.697010.0289-14.67651.80463.74354.0981
5.6161-56.880346.17140-0.0550-19.969916.39800
7.603746.7759-60.16455.78491.023414.5845-24.19611.1724
10.482005.3313-15.81332.738902.5330-5.2709
11.4834-1.8046-3.7435-4.0981-42.75325.31866.697010.0289
0.055015.9728-16.398005.6161-57.086946.17140
-1.0234-14.584516.7803-1.17247.603746.7759-62.76485.7849
-2.73890-2.53305.271910.482005.3313-15.8151
05.31606.694210.0240
5.6132046.14430
7.599546.746205.7813
10.476705.32860
01.80313.74034.0958
-0.0544016.38470
1.024214.578801.1726
2.737802.53070
0-1.8031-3.7403-4.0958
0.05440-16.38470
-1.0242-14.57880-1.1726
-2.73780-2.53070
-41.13125.31606.694210.0240-14.67431.80313.74034.0958
5.6132-56.848646.14430-0.0544-19.960316.38470
7.599546.7462-60.12715.78131.024214.5788-24.17571.1726
10.476705.3286-15.80532.737802.5307-5.2684
11.4743-1.8031-3.7403-4.0958-42.73125.31606.694210.0240
0.054415.9603-16.384705.6132-57.048646.14430
-1.0242-14.578816.7757-1.17267.599546.7462-62.72715.7813
-2.73780-2.53075.268410.476705.3286-15.8053
ox=
1.0e-06*
0.0465
0.0225
-0.0186
0.0482
-0.0799
-0.1160
-0.1154
-0.1080
U=
1.15561.22991.25911.18531.0500
e=
0.22590.39430.41370.29210
5
S=
1.6000+0.8000i2.0000+0.1000i3.7000+1.3000i0-6.4501+0.4422i
I=
0.2349-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 潮流 计算