Các bài tập về Matlab power tool dành cho các sinh viên ngành hệ thống điện, từ bài tập đơn giản đến các bài tập phức tapj, có code sẵn.Việc mô phỏng hệ thống điện liên quan đến một loạt các khung thời gian, các yêu cầu về thời gian thực có thể khác nhau, bắt đầu từ những trường hợp có thời gian cực nhỏ như micro giây khi mô phỏng nhanh điện transients và kéo dài đến một vài năm trong các nghiên cứu quy hoạch hệ thống. Hệ thống có thể có được mô hình và giải quyết trong nhiều cách khác nhau, phụ thuộc vào lĩnh vực nghiên cứu. Hiện nay việc xây dựng phần mềm ứng dụng cho từng công việc trong các khâu của sản xuất điện năng đã trở thành phổ biến. Tuy nhiên trong lĩnh vực nghiên cứu và áp dụng cho giáo dục mong muốn có một phần mềm duy nhất làm nền tảng cung cấp được các chức năng mà từ đó có thể áp dụng phân tích một số bài toán hệ thống điện. Trong những thập kỷ qua, một số ngôn ngữ lập trình cao cấp ví dụ như Matlab, Mathcad, Mathematic đã trở nên phổ biến hơn cho cả hai lĩnh vực gồm mục đích nghiên cứu và mục đích giáo dục. Trong các ngôn ngôn ngữ thì Matlab được người dùng lựa chọn tốt nhất. Các tính năng quan trọng của Matlab là việc xử lý các ma trận, biểu diễn đồ thị và cung cấp cho người dùng môi trường đồ họa (Simulink), với các mô hình điều khiển rất đơn giản lược đồ được thiết kế trong thư viện matlab hoặc được các chuyên gia cập nhật thường xuyên. Các vấn đề được trình bày trong bài báo này bao gồm, phần đầu trình bày các tính năng của một số Matlab toolboxes được sử dụng trong phân tích hệ thống điện, như MatPower toolbox (MPT), Power system analysic toolbox (PSAT) và Voltage Stability Toolbox (VST). Phần thứ hai trình bày một trường hợp nghiên cứu dựa trên sơ đồ lưới điện, sử dụng cho MatPower toolbox (MPT) để minh họa các khả năng của công cụ này cho mục đích nghiên cứu và mục đích giáo dục.
Trang 1Bài tập và lời giải
Hãy sử dụng các lệnh của MATLAB để thực hiện các phép tính sau:
6 Xác định mô dun và góc pha của đại lượng S (bài 5)
7 Hãy biểu thị kết quả xác định diện tích mặt cầu bán kính r = 12,4 m với:
a) Số nguyên;
b) Hai chữ số sau dấu phảy tĩnh;
c) Ba chữ số sau dấu phảy động
6 abs(S) và phase(S) S= 1.8385; gam= 1.1802
7.a fprintf('S=%7.0f m2 \n', 4*pi*12.8^2) S= 2059 m2
7.b fprintf('S=%7.2f m2 \n', 4*pi*12.8^2) S=2058.87 m2
7.c fprintf('S=%7.3e m2 \n', 4*pi*12.8^2) S=2.059e+003 m2
Bài tập 2.1: Hãy xây dựng hàm xác định điện trở của mạch gồm 4 nhánh song với các
Trang 2Bài tập 2.4: Hãy vẽ đồ thị hàm số y = 1,82x3 + 2 và lấy 10 dữ liệu từ đồ thị ;
>> x=1:10;
y=1.82.*x.^3+2;
plot(x,y), grid;
>> ginput(10)
Bài tập 2.5 : Hãy sử dụng vòng lặp if để giải bài toán xác định phụ tải tính toán của hai
hộ dùng điện với công suất tính toán là: P1=1020; P2=2115 kW
Bài tập 2.6: Hãy dùng vòng lặp if – elseif để viết chương trình giải bài toán xác định
một số dương có chia hết cho 4 hoặc cho 6 không Ví dụ với trường hợp n=22
Trang 3Bài tập 2.7: Hãy sử dụng vòng lặp while cho bài toán tìm giá trị lớn nhất của số n mà có
Bài tập 2.8: Hãy sử dụng vòng lặp for cho chương trình giải bài toán xác định vectơ hệ
số nhu cầu phụ thuộc vào số hộ dùng điện:
Bài tập 2.9: Hãy xây dựng hàm xác định tổn thất điện năng trên đường dây và áp dụng
hàm m.file để tính tổn thất trên đường dây 10 kV dài 13 km làm bằng dây dẫn AC-50 (r0=0,65 và x0=0,392 Ω/km), phụ tải của mạng điện là S= 185 kVA, hệ số cosφ = 0,84; thời gian tổn thất cực đại τ=2150 h
7 1 3
4 2 6
=
A
Hãy xác định:
a) các giá trị cực đại của mỗi cột (C1);
b) các giá trị cực tiểu của mỗi cột (C2);
Trang 4c) các giá trị cực đại của mỗi dòng (C3);
d) các giá trị cực tiểu của mỗi dòng (C4)
e) các giá trị cực đại của cột và chỉ số của chúng (C1, I)
>> A=[6 2 4; 3 1 7; 5 8 9];
C1= max(A) % Gia tri cua phan tu cuc dai cua moi cot
C2= min(A) % Gia tri cuc tieu cua moi cot
C3= max(A,[ ],2) % Gia tri cua phan tu cuc dai cua moi dong
C4= min(A,[ ],2) % Gia tri cua phan tu cuc dai cua moi dong
[C4,I]=max(A) % Gia tri cua phan tu cuc dai cua moi cot va chi so cua chung
Bài tập 3.2: Số liệu thống kê về các đại lượng x và y cho trong bảng:
x 24,8 23,5 21,6 27,2 33,5 22,8 19,5 26,43 30,6
Hãy sử dụng các hàm trong Matlab để xác định các đại lượng:
a) Giá trị trung bình của đại lượng x và y: X_tb, Y_tb;
b) Phương sai của các đại lượng x và y: Cx, Cy;
c) Độ lệch chuẩn của đại lượng x và y; xigma1x, xigma1y
d) Độ lệch trung bình bình phương của đại lượng x và y; xigma2x, xigma2y;e) Hệ số tương quan của các giá trị quan sát giữa các đại lượng Rx_y
>> x=[24.8 23.5 21.6 27.2 33.5 22.8 19.5 26.43 30.6];
y=[40.6 40.0 38.7 44.5 52.4 39.4 36.7 40.3 43.6];
X_tb = mean(x) % Gia tri trung binh cua dai luong x
Y_tb = mean(y) % Gia tri trung binh cua dai luong y
Cx=cov(x);
Cy=cov(y)
xigmax1 = std(x,0) % Do lech chuan cua dai luong x
xigmax2 = std(x,1) % Do lech trung binh binh phuong
xigmay2 = std(y,1) % Do lech trung binh binh phuong
Bài tập 3.3: Cho dãy số liệu thống kê như bài 3.2:
x 24,8 23,5 21,6 27,2 33,5 22,8 19,5 26,43 30,6
a) Hãy sắp xếp các phần tử của các đại lượng x và y theo thứ tự tăng dần:
b) Hãy sắp xếp các phần tử của đại lượng y theo thứ tự tăng dần kèm thoe chỉ số;
>> x=[24.8 23.5 21.6 27.2 33.5 22.8 19.5 26.43 30.6];
y=[40.6 40.0 38.7 44.5 52.4 39.4 36.7 40.3 43.6];
sort(x)
[X, INDEX] = sort(x)
Trang 5Bài tập 3.4 : số liệu đo đếm của x và y cho trong bảng sau:
x 24,8 23,5 21,6 27,2 33,5 22,8 19,5 26,43 30,6
a) Hãy xây dựng hàm hồi quy thực nghiệm bậc nhất
b) Hãy xây dựng hàm hồi quy thực nghiệm bậc hai
b) Hãy xác định hàm tương quan hồi quy dạng mũ: y = a0 + a1ex + a2xe-x
và đánh giá sai số tương đối;
XX = [ones(size(x)) exp(-x) x.*exp(-x)];
a = XX\y % He so hoi quy
Trang 6Bài tập 3.8: Hãy xác định hệ số nhu cầu của phụ tải động lực ứng với số lượng thiết bị
là n=12, biết sự phụ thuộc của hệ số này và số lượng thiết bị dùng điện như sau:
Bài tập 3.10: Hãy xác định các tham số của máy biến áp 22/0,4 kV, công suất 250 kVA.
Biết số liệu về các tham số của các máy biến áp như sau:
Trang 7Bài tập 3.11: Hãy xây dựng hàm y = 1,85sin(x) trong khoảng x = 0÷10 và các hàm nội
suy: bậc thang, tuyến tính hóa, lập phương mô phỏng và hàm lập phương và đánh giá sai số so với hàm gốc
Trang 8legend('data',
'function',
'bac thang',
'tuyen tinh hoa',
'lap phuong mo phong',
'lap phuong', 2);
sbt=max (abs(yn - yy));
stt=max (abs(yl - yy));
slpm=max (abs(yc - yy));
slp=max (abs(ys - yy));
disp(' sbt stt slpm slp')
fprintf('%g'), disp([sbt,stt,slpm,slp])
Bài tập 4.1 : Hãy biểu thị ma trận A gồm 2 dòng và 4 cột trên màn hình :
32 14 2
8 6 4 2 π π
=
A
>> A=[2,4,6,8;pi,2*pi,14,32]
Bài tập 4.4: Cho một vectơ-cột v1=[3;5;8], hãy thiết lập một ma trận kích thước nhân
đôi từ các phần tử của vectơ này
>> v1=[3;5;8],
U1=v1(:,[1 1])
Bài tập 4.5: Cho một vectơ-dòng v2=[3,5,8], hãy thiết lập một ma trận kích thước nhân
đôi từ các phần tử của vectơ này
data function bac thang tuyen tinh hoa lap phuong mo phong lap phuong
Trang 9title('Do thi y=f(x) theo lenh plot(x,y)')
Bài tập 5.3 Vẽ đồ thị hàm số y=2x2+1,6x-0,83 bằng lệnh fplot trong khoảng x=1÷2
Trang 10Bài tập 5.4: Vẽ đồ thị hàm số y= 1,3sinx2 trong khoảng -3 ≤ x ≤ 3 bằng lệnh ezplot;
>> ezplot('1.3*(sin(x))^2', -3, 3), grid;
Bài tập 5.5 Vẽ đồ thị hàm số y = 0,23x3 trong khoảng x = 1 ÷ 1000:
a) Với trục các trục tọa độ chia theo tỷ lệ log cơ số 10:
b) Với trục y chia theo tỷ lệ log cơ số 10:
% a)
>> x=2:1:10^3;
y=0.23.*x;
loglog (x,y),grid
xlabel('x,log10'); ylabel('Y, log10')
title('Do thi y=f(x) lenh loglog(x,y) ')
% b)
>> x=2:1:10^3;
y=0.23.*x;
semilogy(x,y),grid
xlabel('x '); ylabel('Y, log10')
title('Do thi y=f(x) lenh semilogy')
Bài tập 5.6 Xây dựng đồ thị hàm d = sin(5ϕ) trong hệ tọa độ cực trong khoảng biến thiên ϕ = 0 ÷ 2π
>> phi= 0 : 0.01 : 2*pi; d = sin( 5* phi );
polar( phi , d )
title('Do thi ham d=sin(5phi) trong toa do cuc')
Bài tập 5.7 Hãy xây dựng đồ thị phụ tải dạng bậc thang ứng với số liệu cho trước
title('Bieu do dang bac thang')
Bài tập 5.17: Vẽ đồ thị 3D bằng lệnh ezplot3 của hàm: x=cost+1; y=tsint và z=1,5t trong khoảng 0≤t≤10π
>> ezplot3 ('3*t*sin(t)', '1.5*t', 't ', [0, 10*pi ])
Trang 11Bài tập 6.1: Tìm nghiệm của phương trình bậc bốn:
Trang 12,2
)cos(
4,2)
(
1 2
5 , 1
2
2 2
4 1
−+
= −
x x
e
x x
x x
Trang 13xác định các giá trị dòng điện chạy trong mạch
>> clear
E=[100; 0; 90];
R=[12; 21; 25];
Z=[1 1 1; R(1) -R(2) 0; 0 R(2) -R(3)];
U=[0; E(1)-E(2); E(2)-E(3)];
I=Z\U; % Dong dien tren cac nhanh
disp('Ket qua la:')
fprintf('%g'), disp([I(1), I(2), I(3)]);
Ket qua la:
Bài tập 9.3 Hãy xác định dòng điện và công suất của mạch điện (hình 9.27), biết điện
áp nguồn u(t)= 120cosωt và phụ tải Z=3,75∠520Ω Vẽ đồ thị của các đại lượng
Um=120; tetau=0; Z=3.75; gam=52;
tetai=tetau-gam; teta =(tetau-tetai)*pi/180;
-100 0 100 200 Dien ap va dong dien
wt, do
0 100 200 300 400 -2000
0 2000 4000 Cong suat toan phan
wt, do
0 100 200 300 400 0
1000 2000 3000 Cong suat tac dung
0 100 200 300 400 -2000
-1000 0 1000 2000 Cong suat phan khang
Trang 14pp=P*ones(1,length(wt));
xline=zeros(1, length(wt));
wt=180/pi*wt;
subplot(2,2,1), plot(wt,u,wt,i,wt,xline), grid
title('Dien ap va dong dien')
xlabel('wt, do'), ylabel('U, I, ')
subplot(2,2,2), plot(wt,p,wt,xline),grid
title('Cong suat toan phan')
xlabel('wt, do'), ylabel('S, VA')
subplot(2,2,3), plot(wt,pr), grid
title('Cong suat tac dung')
xlabel('wt, do'), ylabel('p, W ')
subplot(2,2,4), plot(wt,px),grid
title('Cong suat phan khang')
xlabel('wt, do'), ylabel('q, VAr')
Bài tập 9.4: Nguồn điện xoay chiều u(t) = Umsin(ωt+Ψu), tần số 50 Hz được nối với điện trở R = 5,6 Ω và cuộn dây L = 22,3 mH, điện trở shun Rsh=0,15 Ω, Dòng điện chạy trong mạng Ish=1,5A (hình 4.28) Hãy xác định các giá trị hiệu dụng của điện áp URL, Các thành phần UR và UL Vẽ
đồ thị của các đại lượng này
R = 5.6; Rsh = 0.15; Ish=1.5; f=50; L=22.3*10^-3;
t=0:0.0001:0.04;
omeg=2*pi*f; % Bieu thi gia tri cua omega
X=omeg*L; % Bieu thi gia tri cua dien tro khang
Ushr=Ish*R; % Do roi dien ap tren dien tro shun
Ushx = Ish*X; % Do roi dien ap tren dien tro X
UR=Ushr/sqrt(2); % Gia tri hieu dung cua dien ap tren phan tu dien tro R
UL=Ushx/sqrt(2); % Gia tri hieu dung cua dien ap tren phan tu cam khang
URL=UR*sqrt(2); % Mo dun dien ap tong
UshRL=Ushr*sqrt(2); % Modun dien ap tren dien tro shun
uR=Ushr*sin(omeg*t); % Ham bien thien cua dien ap tren phan tu R
uL= Ushx*sin(omeg*t+pi/2); % Ham bien thien cua dien ap tren phan tu L
Trang 15format bank
disp('Ket qua la:')
fprintf('%g'), disp([UR, UL, URL])
plot(t,uR, t,uL, t,uRL), grid
xlabel('t'), ylabel('U, V')
title('Do thi dien ap')
legend('uR', 'uL', 'uRL')
Bài tập 9.5: Cho mạch điện hình 9.29 vớí nguồn điện áp hình sin tần số 50 Hz:
u=Umsin(ωt) V Hãy xác định các giá trị hiệu dụng của điện áp và dòng điện chạy trong mạch, biết Um=250; R=15,6 Ω và C = 116,42 µF Vẽ đồ thị biến thiên của điện áp
i=Im*sin(omega*t+phi); % Gia tri tuc thoi cua dong dien
ur=UmR*sin(omega*t+phi); % Gia tri tuc thoi cua dien ap Ur
uc=UmC*sin(omega*t+phi-pi/2); % Gia tri tuc thoi cua dien ap Uc
I = Im/sqrt(2); % Gia tri hieu dung cua dong dien
Ur = UmR/sqrt(2); % Gia tri hieu dung cua dien ap
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 -15
-10 -5 0 5 10 15
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 -250
-200 -150 -100 -50 0 50 100 150 200 250
uR uC
Trang 16Bài tập 9.6: Hãy xác định giá trị dòng điện và công suất tác dụng trên đầu vào mạng hai
cực hình 9.30, biết điện áp trên đầu vào là U = 24 V, các giá trị điện dẫn tác dụng và phản kháng tương ứng là : G = 0,018 1/Ω và B=0,022 1/Ω
Ví dụ 9.7: Giá trị hiệu dụng của điện áp trên đầu vào mạch điện hình 9.31 là U=220∠32
V Hãy xác định các giá trị hiệu dụng của dòng điện chạy trong mạch, biết C = 115 µF, R=3,5Ω, L = 210 mH, tần số f = 50 Hz
e1=240; a1 = -12; e2 = 220; a2 = 3.5; R=0.82; X=3.53;
Z=R+j*X; % Dien tro duong day
a1=(-20+a1:5:20+a1)'; % Cho goc a1 bien thien tu -20 den 20
Hình 9.30 Sơ đồ mạch
điện bài tập 9.6
I U
Trang 17a1r=a1*pi/180; % Doi goc a1 sang radian
k=length(a1); % Kich thuoc cua day a1 (so phan tu)
a2=ones(k,1)*a2; % Tao mang voi do dai cua a2
a2r=ones*pi/180; % Chuyen doi do ra radian
U1=e1*cos(a1r)+j*e1.*sin(a1r);
U2=e2*cos(a2r)+j*e2.*sin(a2r);
I12=(U1-U2)./Z; I21=-I12;
S1=U1.*conj(I12); P1=real(S1); Q1=imag(S1);
S2=U2.*conj(I21); P2=real(S2); Q2=imag(S2);
plot(a1,P1, a1,P2, a1,PL), grid;
xlabel('Goc pha a1'), ylabel('P, W')
Trang 18Z=[Z1, Z2, Z3];
[I] = giaimach(E,Z);
Im = abs(I);
theta = angle(I);
Ví dụ 9.10: Mạch điện ba pha hình 9.34 với nguồn điện đối xứng điện áp U=380V
cung cấp cho hai phụ tải: phụ tải thứ nhất mắc theo hình sao điện trở của mỗi pha là
Zs=15+j34 Ω; phụ tải thứ hai mắc theo hình tam giác, điện trở mỗi pha là Ztg=30-j22 Ω Điện trở của đường dây là Zd= 2,3+j5,8 Ω Hãy xác định:
a) Dòng điện, công suất tác dụng và công suất phản kháng cấp đến từ nguồn;
b) Điện áp dây tại đầu vào của phụ tải hỗn hợp;
c) Dòng điện pha của mỗi phụ tải;
d) Công suất tác dụng và phản kháng của tải và đường dây
>> clear all
U=380; Zd=2.3+j*5.8;
Ztg=30-j*22; % Dien tro cua phu tai tam giac
Zs=15+j*34; % Dien tro cua phu tai sao
Z2 = Ztg/3; % Dien tro mot pha cua phu tai tam giac
U1=U/sqrt(3); % Dien áp pha
Z = Zd+(Z2*Zs)/(Z2+Zs); % Tong tro
I=U1/Z; % Dong dien chay tren day dan
S=3*U1*conj(I); % Cong suat ba pha
dU=Zd*I; % Ton that dien ap tu nguon den cac diem tai
U2=U1-dU; % Dien ap tren dau vao phu tai
a30=30*pi/180; % Chuyen doi goc 30 do ra radian
U2d = sqrt(3)*exp(j*a30)*U2; % Dien ap day o cuoi phu tai
Is = U2/Zs; % Dong dien pha cua phu tai sao
Itg = U2/Z2;% Dong dien pha cua phu tai tam giac
Trang 19Iab = Itg/(sqrt(3)* exp(j*a30)); % Dong dien chay trong cuon day tam giac
S1=3*U2*conj(Is); % Cong suat phu tai sao
S2=3*U2*conj(Itg); % Cong suat phu tai tam giac
dS = 3*Zd*I^2; % Ton that cong suat tren duong day
Ssum = S1+S2+dS; % Tong cong suat tu nguon
disp(' I, A S, VA U2, V Is Itg Iab S1 S2 dS Ssum')
Kq = [I,S,U2,Is,Itg,Iab,S1,S2,dS,Ssum]
Bài tập 9.11 Mạng điện ba pha đối xứng điện áp pha 220 V (hình 9.35), điện trở mỗi
pha Zs=1,6+j4,3 Ω, điện trở hỗ cảm giữa các pha Zm= j4,15 Ω Điện trở dây trung tính
Zn=2,5+j*3,4 Hãy xác định dòng điện chạy trong các pha:
fprintf('%g'), disp([abs(Iabc), angle(Iabc)*180/pi])
Bài tập 9.12: Nguồn điện ba pha đối xứng 127 V (hình 9.36) cung cấp cho phụ tải đấu
theo hình sao, trung tính cách ly Mỗi pha có điện trở Zs=3,6+j1,3 Ω, điện trở hỗ cảm giữa các pha Zm= j2,4 Ω Hãy xác định dòng điện chạy trên các pha theo phương pháp:a) Định luật Kirchhoff;
Trang 20disp('(a) Phuong phap DL Kirchhoff:');
a=cos(2*pi/3)+ j*sin(2*pi/3); % Toan tu quay
A=[1 1 1; 1 a^2 a; 1 a a^2];
125
97
140 130
10 210
(a) Ma trận điện trở của phụ tải Z012=A-1ZabcA
(b) Các thành phần đối xứng của điện áp
(c) Các thành phần đối xứng của dòng điện
(d) Dòng điện pha của phụ tải
(e) Công suất toàn phần cung cấp cho phụ tải dưới dạng các thành phần đối xứng(f) Công suất toàn phần cung cấp cho phụ tải dưới dạng công suất theo từng pha.clear
% Ma tran dien ap
Uabc= [210, 10
130, -140
Trang 21Uabcr=Uabc(:, 1).*(cos(pi/180*Uabc(:, 2))+j*sin(pi/180*Uabc(:,2)));
% Cong suat theo tung pha
S3ph=(Uabcr.')*conj(Iabc)
Bài 9.14 Cho mạch điện RLC (hình 9.37) với các tham số L=1.25 Henri; R=0.86 Ω;
C=1.5 Fara Hãy xác định và vẽ đồ thị dòng điện trong mạch và điện áp trên tụ, biểu thị mối quan hệ giữa điện áp tụ và dòng điện trong khoảng thời gian 15s kể từ khi đóng công tắc Biết ở thời điểm ban đầu dòng điện IL(0)=0 và điện áp tại tụ Uc(0)=0,5V, điện
áp nguồn US=2,4V
function xdot = mqdo9_14(t,x);
E = 2.4; R =0.86; L = 1.25; C = 1.5;
xdot = [x(2)/C ; 1/L*( E - x(1) - R*x(2) ) ];
Sau khi scrip này được cất giữ trong m.file ta thực hiện các lệnh:
>> t0= 0; tc =15; % khoang thoi gian
x0 = [0.3, 0]; % Dieu kien ban dau
R
i(t)
L
+_ US
UC
Hình 9.37 Sơ đồ mạch điện
Trang 22ts = [t0, tc];
[t,x] = ode45('mqdo9_14', ts, x0);
plot(t,x( :,2)), grid
title('Bieu do bien thien cua dong dien trong qua trinh qua do')
xlabel('t, sec'), ylabel('I,A')
Bài tập 10.1 Một khu chung cư gồm n=158 hộ gia đình, công suất tiêu thụ trung bình
của mỗi hộ là p0=1,62 kW (không dùng bếp điện), cosϕ = 0,92; Phụ tải động lực gồm các động cơ vệ sinh kỹ thuật với công suất định mức tương ứng là:
Hãy xác định phụ tải tính toán theo hai phương pháp:
a) Phương pháp số gia;
b) Phương pháp hệ số nhu cầu ;
% Xac dinh phu tai tinh toan
>> SH=[158 1.62 0.92]; % Ma tran phu tai sinh hoat
DL=[25 16 10 7.5 5.6 3; % Ma tran phu tai dong luc
0.74 0.78 0.72 0.67 0.80 0.78]; % He so cong suat cua phu tai dong luc
kncvs=interp1(nvs,knc,length(nvs)); % He so nhu cau cua phu tai d.luc
Psh= kdtsh*SH(1)* SH(2); % Cong suat tinh toan phu tai sinh hoat
Pdl=kncvs*sum(DL(1,:)); % Cong suat tinh toan phu tai dong luc
% Ap dung phuong phap so gia
Trang 23Q=sqrt(S^2-P^2);
format bank
disp('a) Ket qua tinh toan theo p.p so gia:')
disp(' Psh, kW Pdl P S, kVA Q,kVAr')
Pdli=[25 16 10 7.5 5.6 3]; % Ma tran phu tai dong luc
cofi =[0.74 0.78 0.72 0.67 0.80 0.78]; % He so cong suat cua phu tai dong lucnvs=[2 3 5 8 10 15 20 30 50];
knc=[1 0.9 0.8 0.75 0.7 0.65 0.63 0.6 0.55];
kncvs=interp1(nvs,knc,length(nvs)); % He so nhu cau cua phu tai d.luc
Psh= kdtsh*SH(:,1)*SH(:,2); % Cong suat tinh toan phu tai sinh hoat
Pdl=kncvs*sum(DL(1,:)); % Cong suat tinh toan phu tai dong luc
disp('b) Ket qua tinh toan theo p.p he so nhu cau:')
disp(' Psh, kW Pdl P S, kVA Q,kVAr')
fprintf('%g'), disp([Psh,Pdl,P,S,Q])
b) Ket qua tinh toan theo p.p he so nhu cau:
Psh, kW Pdl P S, kVA Q,kVAr
81.50 48.65 117.13 137.03 71.12