Bai tap va loi giai matlab power tool

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.

Bà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

Bà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

Bà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);

c) 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)

Bà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

Bà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:

Bà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

legend('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

title('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

Bà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 ])

Bài tập 6.1: Tìm nghiệm của phương trình bậc bốn:

,2

)cos(

4,2)

(

1 2

5 , 1

2

2 2

4 1

−+

= −

x x

e

x x

x x

xá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

pp=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

format 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

Bà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

a1r=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')

Z=[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

Iab = 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;

disp('(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

Uabcr=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

ts = [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

Q=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