MÔI TRƯỜNG MATLAB
4.2.2. Nội dung chương trình Giải :
Giải:
Sau khi mở chương trình Matlab. Trong cửa sổ lệnh ta bắt ựầu việc tắnh toán.
Hình 4.3. Màn hình chương trình trước khi thao tác
Trước hết ta thiết lập các ma trận dữ liệụ Ở ựây không xét ựến máy phát và máy biến áp, nên các ma trận Mp và MBA là các ma trận zerọ Vì các ựường dây kép, nên các giá trị ựiện trở ựược chia ựôi, còn giá trị dung dẫn ựược nhân ựôị
>> clear Ucb=220; Mp=[0 0 0 0 0]; MBA=[0 0 0 0 0]; Dd=[150 220 0.10/2 0.424/2 2*2.68*10^-6 130 220 0.10/2 0.424/2 2*2.68*10^-6 97 220 0.10/2 0.424/2 2*2.68*10^-6 145 220 0.12/2 0.430/2 2*2.64*10^-6 108 220 0.10/2 0.424/2 2*2.68*10^-6 75 220 0.12/2 0.430/2 2*2.64*10^-6];
Trường đại học Nông nghiệp Hà Nội Ờ Luận văn thạc sĩ khoa học kỹ thuật ẦẦẦẦẦẦẦẦẦẦ.. 57
Sau khi khai báo các tham số ban ựầu của mạng ựiện ta gọi chương trình chuyển ựổi ựơn vị sang hệ ựơn vị tương ựối tsoptuHTD :
tsoptuHTD(Ucb,Mp,MBA,Dd);
Hình 4.5 : Gọi chương trình tsoptuHTD
Kết quả hiển thị trên màn hình là:
Tham so duong day:
Rđ(pu) Xđ Bđ 0.0155 0.0657 0.1946 0.0134 0.0569 0.1686 0.0100 0.0425 0.1258 0.0180 0.0644 0.1853 0.0112 0.0473 0.1401 0.0093 0.0333 0.0958
Trường đại học Nông nghiệp Hà Nội Ờ Luận văn thạc sĩ khoa học kỹ thuật ẦẦẦẦẦẦẦẦẦẦ.. 58
Hình 4.6 :Kết quả tắnh toán tham số máy phát, MBA, ựường dây
Bước tiếp theo ta thiết lập chương trình tắnh toán trên MATLAB với việc khai báo các tham số Scb ,ựộ chắnh xác accuracy, bước lặp tối ựa maxit và các mảng dữ liệu về nút và ựường dây ựã ựược xử lý ở trên :
>> clear
Scb = 100; accuracy = 0.0001; maxit = 10;
% No code Mag. Degree MW Mvar MW Mvar Qmin Qmax Mvar dln=[1 1 1.05 0.0 0 0 0 0 10 50 0
2 2 1.035 0.0 15 8 60 35 10 50 0 3 0 1.000 0.0 40 25 0 0 0 0 0 3 0 1.000 0.0 40 25 0 0 0 0 0 4 2 1.04 0.0 20 12 40 20 0 0 0 5 0 1.00 0.0 25 16 0 0 0 0 0]; % Du lieu duong day
Trường đại học Nông nghiệp Hà Nội Ờ Luận văn thạc sĩ khoa học kỹ thuật ẦẦẦẦẦẦẦẦẦẦ.. 59 % Nut Nut R X 1/2 B Dd % dau cuoi p.ụ p.ụ p.ụ dld=[1 2 0.0155 0.0657 0.1946 1 1 5 0.0134 0.0569 0.1686 1 2 3 0.0100 0.0425 0.1258 1 2 4 0.0180 0.0644 0.1853 1 3 4 0.0112 0.0473 0.1401 1 4 5 0.0093 0.0333 0.0958 1];
Hình 4.7 : Khai báo tham số sau khi chuyển về hệ ựơn vị tương ựối
Sau ựó ta tiến hành gọi chương trình ynut, gtppNewton, hienthikq
ynut % Ma tran dien dan
gtppNewton % Giai tich bang phuong phap Newton-Raphson hienthikq % Hien thi ket qua tren man hinh
Trường đại học Nông nghiệp Hà Nội Ờ Luận văn thạc sĩ khoa học kỹ thuật ẦẦẦẦẦẦẦẦẦẦ.. 60
Giai tich bang pp Newton-Raphson
Sai so cuc dai = 9.45948e-005 Nọ buoc tinh = 7
Nọ Dien ap Gocfa ---CS tai--- ---CS phat---- CS bu--- Nut Modun Do MW Mvar MW Mvar Mvar 1 1.050 0.000 0.000 0.000 0.624 -54.805 0.000 2 1.065 0.111 15.000 8.000 60.000 37.781 0.000 3 1.052 -0.343 40.000 25.000 0.000 0.000 0.000 4 1.040 0.157 20.000 12.000 40.000 -120.757 0.000 5 1.045 -0.220 25.000 16.000 0.000 0.000 0.000 Total 100.000 61.000 100.624 -137.781 0.000
Trường đại học Nông nghiệp Hà Nội Ờ Luận văn thạc sĩ khoa học kỹ thuật ẦẦẦẦẦẦẦẦẦẦ.. 61
Nhận xét kết quả thu ựược : ta có thể thấy bài toán hội tụ sau 7 bước lặp với sai số cho phép không quá 9.459.10-5 , ựồng thời các số liệu thu ựược tương ựối chắnh xác với các ựiều kiện thực tế.