Bài tập lớn giải tích mô phỏng hệ thống điện
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (439.47 KB, 41 trang )
TRƯỜNG ĐẠI HỌC QUY NHƠN
KHOA KỸ THUẬT VÀ CÔNG NGHỆ
BÀI TẬP LỚN
MÔN HỌC :
GIẢI TÍCH VÀ MÔ PHỎNG TRONG HỆ THỐNG ĐIỆN
GVHD
: TS. TRẦN MẠNH HÀ
Sinh viên thực hiện
: TRẦN NGỌC SƠN
Lớp
: Kỹ Thuật Điện- Điện tử ĐH54
Khoa
: ĐIỆN.
Bình Định / Ngày.......Tháng.........Năm.......
1
Đề bài:
Cho sơ đồ hệ thống điện mẫu IEEE 30 nút như hình vẽ sau. Các dữ liệu nút và
các dữ liệu nhánh của hệ thống điện được cho trong bảng 1 và 2 tương ứng. Các thông
số được cho trong hệ đơn vị tương đối ứng với lượng công suất cơ bản bằng 100MVA.
Yê u cầu sinh viên thực hiện các nội dung lớn sau đây:
Sinh viên: Nguyễn Văn Lâm
Lớp
: Lớp KTĐ-ĐT K38A
+ Hệ số tải Kt (%) =108
+ Mất đường dây: 10 - 20
Hình 1. Sơ đồ môt sợi hệ thống điện IEEE 30 nút.
I. ỨNG DỤNG POWERWORLD:
1.Xây dựng sơ đồ mô phỏng
2.Xuất kết quả ma trận Ybus
3.Xuất kết quả điện áp tại các nút
4.Xuất kết quả phân bổ công suất trong hệ thống điện
II. ỨNG DỤNG MATLAB TÍNH TOÁN MA TRẬN Ybus
1.Viết hàm tính toán
2.Xuất kết quả ma trậnYbus
3.So sánh với kết quả từ phần mềm Powerworld
III. ỨNG DỤNG MATLAB ( PHƯƠNG PHÁP GAUSS-SEIDEL)
1.Viết đoạn lệnh tính toán chế độ xác lập
2.Xuất kết quả điện áp tại các nút
3.Xuất kết quả phân bố công suất trong hệ thống điện
IV. ỨNG DỤNG MATLAB ( PHƯƠNG PHÁP NEWTON-RAPHSON)
1.Viết đoạn lệnh tính toán chế độ xác lập
2.Xuất kết quả điện áp tại các nút
3.Xuất kết quả phân bố công suất trong hệ thống điện
Bảng 1. Dữ liệu nút của hệ thống điện IEEE 30 nút.
Nút
Mã
nút
Điện
áp
(pu)
Góc
pha
(độ)
Pt
Qt
Pf
Qf
Qmin
Qmax
Qshunt
(MW)
(Mvar)
(MW)
(Mvar)
(Mvar)
(Mvar)
(Mvar)
1
1
1.06
0
0
0
0
0
0
0
0
2
2
1
0
23,436
13,716
40
48.822
0
0
0
3
0
1
0
2,592
1,296
0
0
0
0
0
4
0
1.06
0
8,208
1,728
0
0
0
0
0
5
2
1
0
101,73
6
20,52
0
35.975
0
0
0
6
0
1
0
0
0
0
0
0
0
0
7
0
1
0
24,624
11,772
0
0
0
0
0
8
2
1
0
32,4
32,4
0
30.826
0
0
0
9
0
1
0
0
0
0
0
0
0
0
10
0
1
0
6,264
2,16
0
0
0
0
19
11
2
1
0
0
0
0
16.119
0
0
0
12
0
1
0
12,096
8,1
0
0
0
0
0
13
2
1
0
0
0
0
10.423
0
0
0
14
0
1
0
6,696
1,728
0
0
0
0
0
15
0
1
0
8,856
2,7
0
0
0
0
0
16
0
1
0
3,78
1,944
0
0
0
0
0
17
0
1
0
9,72
6,264
0
0
0
0
0
18
0
1
0
3,456
0,972
0
0
0
0
0
19
0
1
0
10,26
3,672
0
0
0
0
0
20
0
1
0
2,376
0,756
0
0
0
0
0
21
0
1
0
18,9
12,096
0
0
0
0
0
22
0
1
0
0
0
0
0
0
0
0
23
0
1
0
3,456
1,728
0
0
0
0
0
24
0
1
0
9,396
7,236
0
0
0
0
4.3
25
0
1
0
0
0
0
0
0
0
0
26
0
1
0
3,78
2,484
0
0
0
0
0
27
0
1
0
0
0
0
0
0
0
0
28
0
1
0
0
0
0
0
0
0
0
29
0
1
0
2,592
0,972
0
0
0
0
0
30
0
1
0
11,448
2,052
0
0
0
0
0
Bảng 2. Dữ liệu nhánh của hệ thống điện IEEE 30 nút
Nút đi
Nút đến
R (pu)
X (pu)
B/2 (pu)
k
1
2
0.0192
0.0575
0.0264
1
1
3
0.0452
0.1852
0.0204
1
2
4
0.057
0.1737
0.0184
1
3
4
0.0132
0.0379
0.0042
1
2
5
0.0472
0.1983
0.0209
1
2
6
0.0581
0.1763
0.0187
1
4
6
0.0119
0.0414
0.0045
1
5
7
0.046
0.116
0.0102
1
6
7
0.0267
0.082
0.0085
1
6
8
0.012
0.042
0.0045
1
6
9
0
0.208
0
0.978
6
10
0
0.556
0
0.969
9
11
0
0.208
0
1
9
10
0
0.11
0
1
4
12
0
0.256
0
0.932
12
13
0
0.14
0
1
12
14
0.1231
0.2559
0
1
12
15
0.0662
0.1304
0
1
12
16
0.0945
0.1987
0
1
14
15
0.221
0.1997
0
1
16
17
0.0824
0.1923
0
1
15
18
0.1073
0.2185
0
1
18
19
0.0639
0.1292
0
1
19
20
0.034
0.068
0
1
10
20
0.0936
0.209
0
1
10
17
0.0324
0.0845
0
1
10
21
0.0348
0.0749
0
1
10
22
0.0727
0.1499
0
1
21
22
0.0116
0.0236
0
1
15
23
0.1
0.202
0
1
22
24
0.115
0.179
0
1
23
24
0.132
0.27
0
1
24
25
0.1885
0.3292
0
1
25
26
0.2544
0.38
0
1
25
27
0.1093
0.2087
0
1
28
27
0
0.396
0
0.968
27
29
0.2198
0.4153
0
1
27
30
0.3202
0.6027
0
1
29
30
0.2399
0.4533
0
1
8
28
0.0636
0.2
0.0214
1
6
28
0.0169
0.0599
0.065
1
I. ỨNG DỤNG POWERWORLD:
1.
Xây dựng sơ đồ mô phỏng
Hình 2: Xây dựng sơ đồ mô phỏng trên phần mềm Power world 20
Bảng 3: Dữ liệu nút của sơ đồ xây dụng trên Case Information
Bảng 4: Dữ liệu nhánh của sơ đồ xây dụng trên Case Information
2.
Xuất kết quả ma trận Ybus
j = sqrt(-1);
Ybus = sparse(30);
Ybus( 1, 1) = 6.4684+ j*( -20.6959);
Ybus( 1, 2) = -5.2246+ j*( 15.6467);
Ybus( 1, 3) = -1.2437+ j*( 5.0960);
Ybus( 2, 1) = -5.2246+ j*( 15.6467);
Ybus( 2, 2) = 9.7523+ j*( -30.6487);
Ybus( 2, 4) = -1.7055+ j*( 5.1974);
Ybus( 2, 5) = -1.1360+ j*( 4.7725);
Ybus( 2, 6) = -1.6861+ j*( 5.1165);
Ybus( 3, 1) = -1.2437+ j*( 5.0960);
Ybus( 3, 3) = 9.4392+ j*( -28.6023);
Ybus( 3, 4) = -8.1954+ j*( 23.5309);
Ybus( 4, 2) = -1.7055+ j*( 5.1974);
Ybus( 4, 3) = -8.1954+ j*( 23.5309);
Ybus( 4, 4) = 16.3141+ j*( -54.9186);
Ybus( 4, 6) = -6.4131+ j*( 22.3112);
Ybus( 4, 12) = -0.0000+ j*( 3.9062);
Ybus( 5, 2) = -1.1360+ j*( 4.7725);
Ybus( 5, 5) = 4.0900+ j*( -12.1906);
Ybus( 5, 7) = -2.9540+ j*( 7.4493);
Ybus( 6, 2) = -1.6861+ j*( 5.1165);
Ybus( 6, 4) = -6.4131+ j*( 22.3112);
Ybus( 6, 6) = 22.3416+ j*( -82.4350);
Ybus( 6, 7) = -3.5902+ j*( 11.0261);
Ybus( 6, 8) = -6.2893+ j*( 22.0126);
Ybus( 6, 9) = -0.0000+ j*( 4.8077);
Ybus( 6, 10) = -0.0000+ j*( 1.7986);
Ybus( 6, 28) = -4.3628+ j*( 15.4636);
Ybus( 7, 5) = -2.9540+ j*( 7.4493);
Ybus( 7, 6) = -3.5902+ j*( 11.0261);
Ybus( 7, 7) = 6.5442+ j*( -18.4567);
Ybus( 8, 6) = -6.2893+ j*( 22.0126);
Ybus( 8, 8) = 7.7333+ j*( -26.5275);
Ybus( 8, 28) = -1.4440+ j*( 4.5408);
Ybus( 9, 6) = -0.0000+ j*( 4.8077);
Ybus( 9, 9) = 0.0000+ j*( -18.7063);
Ybus( 9, 10) = -0.0000+ j*( 9.0909);
Ybus( 9, 11) = -0.0000+ j*( 4.8077);
Ybus( 10, 6) = -0.0000+ j*( 1.7986);
Ybus( 10, 9) = -0.0000+ j*( 9.0909);
Ybus( 10, 10) = 16.7791+ j*( -48.3791);
Ybus( 10, 17) = -3.9560+ j*( 10.3174);
Ybus( 10, 21) = -10.2037+ j*( 21.9614);
Ybus( 10, 22) = -2.6193+ j*( 5.4008);
Ybus( 11, 9) = -0.0000+ j*( 4.8077);
Ybus( 11, 11) = 0.0000+ j*( -4.8077);
Ybus( 12, 4) = -0.0000+ j*( 3.9062);
Ybus( 12, 12) = 6.5740+ j*( -24.4242);
Ybus( 12, 13) = -0.0000+ j*( 7.1429);
Ybus( 12, 14) = -1.5266+ j*( 3.1734);
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
Ybus(
12,
12,
13,
13,
14,
14,
14,
15,
15,
15,
15,
15,
16,
16,
16,
17,
17,
17,
18,
18,
18,
19,
19,
19,
20,
20,
21,
21,
21,
22,
22,
22,
22,
23,
23,
23,
24,
24,
24,
24,
25,
25,
25,
25,
26,
26,
27,
27,
27,
27,
27,
28,
28,
15) =
16) =
12) =
13) =
12) =
14) =
15) =
12) =
14) =
15) =
18) =
23) =
12) =
16) =
17) =
10) =
16) =
17) =
15) =
18) =
19) =
18) =
19) =
20) =
19) =
20) =
10) =
21) =
22) =
10) =
21) =
22) =
24) =
15) =
23) =
24) =
22) =
23) =
24) =
25) =
24) =
25) =
26) =
27) =
25) =
26) =
25) =
27) =
28) =
29) =
30) =
6) =
8) =
-3.0954+ j*( 6.0973);
-1.9520+ j*( 4.1044);
-0.0000+ j*( 7.1429);
0.0000+ j*( -7.1429);
-1.5266+ j*( 3.1734);
4.0175+ j*( -5.4243);
-2.4910+ j*( 2.2509);
-3.0954+ j*( 6.0973);
-2.4910+ j*( 2.2509);
9.3655+ j*( -16.0116);
-1.8108+ j*( 3.6874);
-1.9683+ j*( 3.9761);
-1.9520+ j*( 4.1044);
3.8346+ j*( -8.4979);
-1.8826+ j*( 4.3935);
-3.9560+ j*( 10.3174);
-1.8826+ j*( 4.3935);
5.8386+ j*( -14.7110);
-1.8108+ j*( 3.6874);
4.8865+ j*( -9.9062);
-3.0757+ j*( 6.2188);
-3.0757+ j*( 6.2188);
8.9580+ j*( -17.9835);
-5.8824+ j*( 11.7647);
-5.8824+ j*( 11.7647);
5.8824+ j*( -11.7647);
-10.2037+ j*( 21.9614);
26.9783+ j*( -56.0891);
-16.7746+ j*( 34.1277);
-2.6193+ j*( 5.4008);
-16.7746+ j*( 34.1277);
21.9345+ j*( -43.4829);
-2.5405+ j*( 3.9544);
-1.9683+ j*( 3.9761);
3.4298+ j*( -6.9653);
-1.4614+ j*( 2.9892);
-2.5405+ j*( 3.9544);
-1.4614+ j*( 2.9892);
5.3118+ j*( -9.2313);
-1.3099+ j*( 2.2876);
-1.3099+ j*( 2.2876);
4.4957+ j*( -7.8650);
-1.2165+ j*( 1.8171);
-1.9693+ j*( 3.7602);
-1.2165+ j*( 1.8171);
1.2165+ j*( -1.8171);
-1.9693+ j*( 3.7602);
3.6523+ j*( -32.1877);
-0.0001+ j*( 25.2525);
-0.9955+ j*( 1.8810);
-0.6875+ j*( 1.2940);
-4.3628+ j*( 15.4636);
-1.4440+ j*( 4.5408);
Ybus( 28, 27) = -0.0001+ j*( 25.2525);
Ybus( 28, 28) = 5.8069+ j*( -45.1705);
Ybus( 29, 27) = -0.9955+ j*( 1.8810);
Ybus( 29, 29) = 1.9076+ j*( -3.6044);
Ybus( 29, 30) = -0.9121+ j*( 1.7234);
Ybus( 30, 27) = -0.6875+ j*( 1.2940);
Ybus( 30, 29) = -0.9121+ j*( 1.7234);
Ybus( 30, 30) = 1.5995+ j*( -3.0173);
V( 1) = 1.060000 + j*( 0.000000);
V( 2) = 0.995182 + j*(-0.098041);
V( 3) = 0.969361 + j*(-0.148009);
V( 4) = 0.947499 + j*(-0.176061);
V( 5) = 0.914801 + j*(-0.265331);
V( 6) = 0.932580 + j*(-0.208049);
V( 7) = 0.914207 + j*(-0.239650);
V( 8) = 0.925352 + j*(-0.221626);
V( 9) = 0.919149 + j*(-0.264609);
V( 10) = 0.894842 + j*(-0.289569);
V( 11) = 0.951679 + j*(-0.273974);
V( 12) = 0.879909 + j*(-0.286524);
V( 13) = 0.894654 + j*(-0.291325);
V( 14) = 0.853239 + j*(-0.300140);
V( 15) = 0.844490 + j*(-0.299521);
V( 16) = 0.876134 + j*(-0.291349);
V( 17) = 0.881698 + j*(-0.292585);
V( 18) = 0.799505 + j*(-0.318362);
V( 19) = 0.778231 + j*(-0.326775);
V( 20) = 0.776154 + j*(-0.327650);
V( 21) = 0.883733 + j*(-0.290959);
V( 22) = 0.882822 + j*(-0.290266);
V( 23) = 0.846208 + j*(-0.294144);
V( 24) = 0.860788 + j*(-0.282906);
V( 25) = 0.889692 + j*(-0.255032);
V( 26) = 0.866980 + j*(-0.257574);
V( 27) = 0.918473 + j*(-0.234502);
V( 28) = 0.924169 + j*(-0.223749);
V( 29) = 0.889217 + j*(-0.252942);
V( 30) = 0.871108 + j*(-0.266416);
3.
Xuất kết quả điện áp tại các nút
Nom kV
132.00
132.00
132.00
132.00
132.00
132.00
132.00
132.00
1.00
Volt
PU Volt (kV)
106.000 139.920
100.000 132.000
0.97301
128.437
0.95498
126.058
0.94490
124.727
0.94444
124.666
0.93490
123.406
0.93958
124.024
0.94488 0.945
33.00
11.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
33.00
132.00
33.00
33.00
0.92851
0.97913
0.91445
0.93013
0.89310
0.88437
0.91181
0.91703
0.84833
0.83157
0.82997
0.91821
0.91709
0.88365
0.89328
0.91170
0.89027
0.93368
0.93655
0.90982
0.89603
30.641
10.770
30.177
30.694
29.472
29.184
30.090
30.262
27.995
27.442
27.389
30.301
30.264
29.161
29.478
30.086
29.379
30.811
123.625
30.024
29.569
4.
Xuất kết quả phân bố công suất trong hệ thống điện
Bảng 5: Kết quả phân bổ công suất của lưới điện
II. ỨNG DỤNG MATLAB TÍNH TOÁN MA TRẬN
clc;
clear all;
disp('
___________
BÀI TẬP LỚNN
___________________
');
disp('
|
GIẢI TÍCH VÀ MÔ PHỎNG HỆ THỐNG ĐIỆN
|
');
disp('
|
GVHD: TS.NGÔ MINH KHOA
|
');
disp('
|
SVTH: NGUYỄN VĂN LÂM
|
');
disp('
|
LỚP : KTD-DT K38A MÃ SV:3851070091
|
');
disp('
|_______________________________________________________|
');
%-------------------------------------------------------------------------%
Scb = 100; accuracy = 0.001; accel = 1.8; maxiter = 100;
%
nut manut dienap
gocpha Pt
Qt
Pf
Qf Qmin Qmax
dlieunut=[1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
%
tsoday=
1
2
0
0
2
0
0
2
0
0
2
0
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1.06
1.0
1.0
1.06
1.0
1.0
1.0
1.0
1.0
1.0
1.0
1.0
1.0
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
nutdi nutden R(pu)
[ 1
2
0.0192
1
3
0.0452
2
4
0.057
3
4
0.0132
2
5
0.0472
2
6
0.0581
4
6
0.0119
5
7
0.046
6
7
0.0267
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0.0
0.0
23.44 13.72
2.59 1.30
8.21 1.73
101.74 20.52
0.0 0.0
24.62 11.77
32.4 32.4
0.0 0.0
6.26 2.16
0.0 0.0
12.1 8.10
0
0.0
6.7
1.73
8.86 2.7
3.78 1.94
9.72 6.26
3.46 0.97
10.26 3.67
2.38
0.76
18.90 12.1
0
0.0
3.46
1.73
9.40
7.24
0
0.0
3.78
2.48
0
0.0
0
0.0
2.59 0.97
11.45 2.05
X(pu)
0.0575
0.1852
0.1737
0.0379
0.1983
0.1763
0.0414
0.116
0.082
B/2(pu)
0.0264
0.0204
0.0184
0.0042
0.0209
0.0187
0.0045
0.0102
0.0085
0.0
40.0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
k
1
1
1
1
1
1
1
1
1
0.0
48.822
0.0
0.0
35.975
0.0
0.0
30.826
0.0
0.0
16.119
0
10.423
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Qshunt
100
0
0
0
0
0
0
0
0
19
0
0
0
0
0
0
0
0
0
0
0
0
0
4.3
0
0
0
0
0
0];
6
8
0.012
0.042
0.0045
1
6
9
0
0.208
0
0.978
6
10 0
0.556
0
0.969
9
11 0
0.208
0
1
9
10 0
0.11
0
1
4
12 0
0.256
0
0.932
12 13 0
0.14
0
1
12 14 0.1231 0.2559 0
1
12 15 0.0662 0.1304 0
1
12 16 0.0945 0.1987 0
1
14 15 0.221
0.1997 0
1
16 17 0.0824 0.1923 0
1
15 18 0.1073 0.2185 0
1
18 19 0.0639 0.1292 0
1
19 20 0.034
0.068
0
1
10 17 0.0324 0.0845 0
1
10 21 0.0348 0.0749 0
1
10 22 0.0727 0.1499 0
1
21 22 0.0116 0.0236 0
1
15 23 0.1
0.202
0
1
22 24 0.115
0.179
0
1
23 24 0.132
0.27
0
1
24 25 0.1885 0.3292 0
1
25 26 0.2544 0.38
0
1
25 27 0.1093 0.2087 0
1
28 27 0
0.396
0
0.968
27 29 0.2198 0.4153 0
1
27 30 0.3202 0.6027 0
1
29 30 0.2399 0.4533 0
1
8
28 0.0636 0.2
0.0214
1
6
28 0.0169 0.0599 0.065
1];
j=sqrt(-1); i = sqrt(-1);
nl = tsoday(:,1); nr = tsoday(:,2); R = tsoday(:,3);
X = tsoday(:,4); Bc = j*tsoday(:,5); a = tsoday(:, 6);
nbr=length(tsoday(:,1)); nbus = max(max(nl), max(nr));
Z = R + j*X; y= ones(nbr,1)./Z;
for n = 1:nbr
if a(n) <= 0 a(n) = 1; else end
Ynut=zeros(nbus,nbus);
for k=1:nbr;
Ynut(nl(k),nr(k))=Ynut(nl(k),nr(k))-y(k)/a(k);
Ynut(nr(k),nl(k))=Ynut(nl(k),nr(k));
end
end
for n=1:nbus
for k=1:nbr
if nl(k)==n
Ynut(n,n) = Ynut(n,n)+y(k)/(a(k)^2) + Bc(k);
elseif nr(k)==n
Ynut(n,n) = Ynut(n,n)+y(k) +Bc(k);
else, end
end
end
clear Pgg
A=Ynut
2.
Xuất kết quả ma trận Ybus
A=
Columns 1 through 7
6.4684 -20.6959i -5.2246 +15.6467i -1.2437 + 5.0960i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-5.2246 +15.6467i
9.7523 -30.6487i
0.0000 + 0.0000i -1.7055 + 5.1974i
-1.1360 + 4.7725i -1.6861 + 5.1165i 0.0000 + 0.0000i
-1.2437 + 5.0960i
0.0000 + 0.0000i
9.4392 -28.6023i -8.1954 +23.5309i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -1.7055 + 5.1974i -8.1954 +23.5309i 16.3141 -55.5094i
0.0000 + 0.0000i -6.4131 +22.3112i 0.0000 + 0.0000i
0.0000 + 0.0000i -1.1360 + 4.7725i
0.0000 + 0.0000i
0.0000 + 0.0000i
4.0900 -12.1906i 0.0000 + 0.0000i -2.9540 + 7.4493i
0.0000 + 0.0000i -1.6861 + 5.1165i
0.0000 + 0.0000i -6.4131 +22.3112i
0.0000 + 0.0000i 22.3416 -82.7706i -3.5902 +11.0261i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
-2.9540 + 7.4493i -3.5902 +11.0261i 6.5442 -18.4567i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -6.2893 +22.0126i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 4.9158i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 1.8561i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 4.1913i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -4.3628 +15.4636i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
Columns 8 through 14
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 4.1913i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-6.2893 +22.0126i
0.0000 + 4.9158i
0.0000 + 1.8561i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
7.7333 -26.5275i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 -18.7063i
0.0000 + 9.0909i
0.0000 + 4.8077i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 9.0909i 11.6772 -37.5884i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 4.8077i
0.0000 + 0.0000i
0.0000 - 4.8077i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
6.5740 -24.4242i 0.0000 + 7.1429i -1.5266 + 3.1734i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 7.1429i 0.0000 - 7.1429i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
-1.5266 + 3.1734i 0.0000 + 0.0000i 4.0175 - 5.4243i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
-3.0954 + 6.0973i 0.0000 + 0.0000i -2.4910 + 2.2509i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
-1.9520 + 4.1044i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -3.9560 +10.3174i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -5.1019 +10.9807i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -2.6193 + 5.4008i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-1.4440 + 4.5408i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
Columns 15 through 21
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -3.9560 +10.3174i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -5.1019 +10.9807i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-3.0954 + 6.0973i -1.9520 + 4.1044i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-2.4910 + 2.2509i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
9.3655 -16.0116i
0.0000 + 0.0000i
0.0000 + 0.0000i -1.8108 + 3.6874i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
3.8346 - 8.4979i -1.8826 + 4.3935i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -1.8826 + 4.3935i
5.8386 -14.7110i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-1.8108 + 3.6874i
0.0000 + 0.0000i
0.0000 + 0.0000i
4.8865 - 9.9062i
-3.0757 + 6.2188i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -3.0757 + 6.2188i
8.9580 -17.9835i -5.8824 +11.7647i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
-5.8824 +11.7647i 5.8824 -11.7647i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 21.8765 -45.1084i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -16.7746 +34.1277i
-1.9683 + 3.9761i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
Columns 22 through 28
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -4.3628 +15.4636i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -1.4440 + 4.5408i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-2.6193 + 5.4008i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -1.9683 + 3.9761i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-16.7746 +34.1277i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
21.9345 -43.4829i
0.0000 + 0.0000i -2.5405 + 3.9544i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
3.4298 - 6.9653i -1.4614 + 2.9892i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
-2.5405 + 3.9544i -1.4614 + 2.9892i
5.3118 - 9.2313i -1.3099 + 2.2876i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -1.3099 + 2.2876i
4.4957 - 7.8650i
-1.2165 + 1.8171i -1.9693 + 3.7602i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -1.2165 + 1.8171i
1.2165 - 1.8171i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -1.9693 + 3.7602i
0.0000 + 0.0000i 3.6523 - 9.4604i 0.0000 + 2.6087i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 2.6087i 5.8068 -22.6130i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -0.9955 + 1.8810i 0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i -0.6875 + 1.2940i 0.0000 + 0.0000i
0.0000 + 0.0000i
Columns 29 through 30
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i
-0.9955 + 1.8810i -0.6875 + 1.2940i
0.0000 + 0.0000i 0.0000 + 0.0000i
1.9076 - 3.6044i -0.9121 + 1.7234i
-0.9121 + 1.7234i 1.5995 - 3.0173i
>>
3.So sánh với kết quả từ phần mềm Powerworld
Từ kết quả ma trận Ybus xuất từ phần mềm Powerworld và kết quả tính toán trong
Matlab ta thấy kết quả tính toán từ hai phầm mềm tương đương nhau.
III. ỨNG DỤNG MATLAB ( PHƯƠNG PHÁP GAUSS-SEIDEL)
1.
Viết đoạn lệnh tính toán chế độ xác lập
Vm=0; delta=0; yload=0; deltad =0;
nbus = length(dlieunut(:,1));
for k=1:nbus
n=dlieunut(k,1);
kb(n)=dlieunut(k,2); Vm(n)=dlieunut(k,3); delta(n)=dlieunut(k, 4);
Pd(n)=dlieunut(k,5); Qd(n)=dlieunut(k,6); Pg(n)=dlieunut(k,7); Qg(n) =
dlieunut(k,8);
Qmin(n)=dlieunut(k, 9); Qmax(n)=dlieunut(k, 10);
Qsh(n)=dlieunut(k, 11);
if Vm(n) <= 0 Vm(n) = 1.0; V(n) = 1 + j*0;
else delta(n) = pi/180*delta(n);
V(n) = Vm(n)*(cos(delta(n)) + j*sin(delta(n)));
P(n)=(Pg(n)-Pd(n))/Scb;
Q(n)=(Qg(n)-Qd(n)+ Qsh(n))/Scb;
S(n) = P(n) + j*Q(n);
end
DV(n)=0;
end
num = 0; AcurBus = 0; converge = 1;
Vc = zeros(nbus,1)+j*zeros(nbus,1); Sc = zeros(nbus,1)+j*zeros(nbus,1);
while exist('accel')~=1
accel = 1.3;
end
while exist('accuracy')~=1
accuracy = 0.001;
end
while exist('Scb')~=1
Scb= 100;
end
while exist('maxiter')~=1
maxiter = 100;
end
iter=0;
maxerror=10;
while maxerror >= accuracy & iter <= maxiter
iter=iter+1;
for n = 1:nbus;
YV = 0+j*0;
for L = 1:nbr;
if nl(L) == n, k=nr(L);
YV = YV + Ynut(n,k)*V(k);
elseif nr(L) == n, k=nl(L);
YV = YV + Ynut(n,k)*V(k);
end
end
Sc = conj(V(n))*(Ynut(n,n)*V(n) + YV) ;
Sc = conj(Sc);
DP(n) = P(n) - real(Sc);
DQ(n) = Q(n) - imag(Sc);
if kb(n) == 1
S(n) =Sc; P(n) = real(Sc); Q(n) = imag(Sc); DP(n) =0; DQ(n)=0;
