BÁO CÁO THỰC HÀNH MATLAB

MATLAB là một trong những hệ thống được soạn thảo hết sức chu đáo, mạch lạc với việc áp dụng nhiều thuật toán phong phú, đặc biệt như chính tên của nó

¬
6 



MATLAB là m͡t trong nhͷng h͏ th͙ng đưͫc so̩n th̫o h͇t sͱc chu đáo, m̩ch l̩c vͣi vi͏c áp dͭng nhi͉u thu̵t toán phong phú, đ̿c bi͏t như chính tên cͯa nó

Cú pháp ngôn ngͷ trong chương trình MATLAB đưͫc thi͇t l̵p r̭t tinh t͇ r̭t thu̵n lͫi cho ngưͥi s͵ dͭng Ma tr̵n đưͫc áp dͭng r͡ng rãi trong các bài toán phͱc t̩p như h͏ th͙ng đi͏n l͹c, t͹ đ͡ng hóa, các h͏ th͙ng cơ h͕c đ͡ng, cơ h͕c tĩnh v.v«

MATLAB th͹c s͹ là m͡t phương ti͏n hͷu hi͏u, cho phép nhân kh̫ năng cͯa con ngưͥi trong lĩnh v͹c h͕c t̵p và nghiên cͱ và ͱng dͭng khoa h͕c kͿ thu̵t lên nhi͉u l̯n

Chương trình MATLAB nó cho phép vӁ tҩt cҧ các đӗ thӏ vӟi các dҥng khác nhau trong không gian 2D, 3D

Lũy thӯa

Khia căn

Ký

b) Các phép toán v͉ s͙ phͱc như sau:

Khi mӝt sӕ S đưӧc biӇu thӏ dưӟi dҥng phӭc như cho trong bҧng sau:

Liên hӧp phӭc Slh = conj(S)

hase(x) Góc pha cӫa sӕ phӭc sum(v) Tәng các phҫn tӱ vector log(x) Logarit tӵ nhiên prod(v) Tích các phҫn tӱ vector log10(x) Logarit cơ sӕ 10 min(v) hҫn tӱ vector bé nhҩt sqtr(x) Căn bұc hai max)v) hҫn tӱ vector lӟn nhҩt conv(x,x) Tích chұp x cӫa chính nó mean(v) Giá trӏ trung bình cӝng Các hàm lưӧng giác

Sinh(x) Hàm sin hyperbol asinh(x) Hàm arsin hyperbol

Cosh(x) Hàm cosin hyperbol acosh(x) Hàm arccos yperbol

?ưӧc thӵc hiӋn bӣi lӋnh :

xac dinh gia tri hieu dung cua dien ap URL Cac thanh phan Ur va UL

UR,V UL,V URL,V

Gia tri hieu dung dien ap va dong dien chay trong mach

I,A Ur,V Uc,V

Mӝt khu chung cư gӗm n = 158 hӝ gia đình, công suҩt trung bình cӫa mӛi hӝ là p0

= 1.62 kW(không dùng bӃp điӋn), cosphi = 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à:

¬
6 

Chương trình Matlab như sau:

11 disp('a) Giai bai toan xac dinh phu tai tinh toan theo phuong phap so gia'); n=158;

w Phương pháp h͏ s͙ nhu c̯u:

Chương trình Matlab như sau:

11disp('b) Giai bai toan xac dinh phu tai tinh toan theo phuong phap he so nhu cau');

¬
6 

tt=[45,30,25,22,20,20,16,16,10,7.5,5];

ksd=[0.63,0.49,0.55,0.52,0.46,0.60,0.63,0.72,0.70,0.68,0.66];

cosphi=[0.75,0.71,0.72,0.78,0.68,0.8,0.8,0.82,0.8,0.75,0.74];

Chương trình Matlab như sau:

11 disp('Xac dinh phu tai tinh toan cua mot phan xuong gom 11 phu tai');

dl=[45,30,25,22,20,20,16,16,10,7.5,5;0.63,0.49,0.55,0.52,0.46,0.60,0.63,0.72,0.70,0.68,0.66;0.75,0.71,0.72,0.78,0.68,0.8,0.8,0.82,0.8,0.75,0.74];

| h̵n xét: Như v̵y tính toán phͭ t̫i đi͏n vͣi s͹ trͫ giúp cͯa chương trình

MATLAB r̭t thu̵n lͫi và cho k͇t qu̫ nhanh chóng

ñ

Hãy xác đӏnh phө tҧi tính toán cho mӝt xí nghiӋp gӗm 4 phân xưӣng sҧn xuҩt vӟi các thiӃt bӏ đӝng lӵc có các tham sӕ cho trong bҧng bt 10.3 như sau:

dl1=[45,30,25,22,20,20,16,16,10;0.63,0.49,0.55,0.52,0.46,0.60,0.63,0.72,0.70;0.75,0.71,0.72,0.78,0.68,0.80,0.80,0.82,0.80];

dl2=[25,22,20,16,10,10,7.5,6.3,5;0.55,0.52,0.46,0.72,0.70,0.62,0.61,0.55,0.52;0.72,0.78,0.68,0.82,0.80,0.72,0.71,0.72,0.78];

dl3=[2.8,8,10,16,4.5,1.5,5.5,5.5,5.5,8;0.35,0.32,0.32,0.23,0.26,0.42,0.35,0.26,0.45,0.45;0.67,0.58,0.58,0.65,0.66,0.62,0.67,0.66,0.67,0.67];

dl4=[8.5,3,7.5,6.3,10,7.5,40,35;0.37,0.3,0.41,0.45,0.47,0.45,0.53,0.45;0.66,0.58,0.63,0.67,0.7,0.63,0.9,0.58];

Chương trình Matlab như sau:

11 disp('Xac dinh phu tai tinh toan gom 4 phan xuong gom');

disp('Xac dinh phu tai tinh toan phan xuong 1 ');

dl1=[45,30,25,22,20,20,16,16,10;0.63,0.49,0.55,0.52,0.46,0.60,0.63,0.72,0.70;0.75,0.71,0.72,0.78,0.68,0.80,0.80,0.82,0.80];

ksdth=[0.2,0.3,0.4,0.5,0.6,0.7,0.8];

kb=[3,3.5,4,5,6.5,8,10];

disp('ket qua la')

disp(' S1,kVA 1,kW Q1,kVAr cosphitb1')

fprintf('%g'),disp([S1 1 Q1 cosphitb1])

disp('xac dinh phu tai tinh toan phan xuong 2');

dl2=[25,22,20,16,10,10,7.5,6.3,5;0.55,0.52,0.46,0.72,0.70,0.62,0.61,0.55,0.52;0.72,0.78,0.68,0.82,0.80,0.72,0.71,0.72,0.78];

disp('ket qua la')

disp(' S2,kVA 2,kW Q2,kVAr cosphitb2')

fprintf('%g'),disp([S2 2 Q2 cosphitb2])

disp('xac dinh phu tai tinh toan phan xuong 3');

dl3=[2.8,8,10,16,4.5,1.5,5.5,5.5,5.5,8;0.35,0.32,0.32,0.23,0.26,0.42,0.35,0.26,0.45,0.45;0.67,0.58,0.58,0.65,0.66,0.62,0.67,0.66,0.67,0.67];

¬
6 

cosphitb3=sum(dl3(1,:).*dl3(3,:))/sum(dl3(1,:));

S3=3/cosphitb3;

Q3=sqrt(S3^2-3^2);

disp('ket qua la')

disp(' S3,kVA 3,kW Q3,kVAr cosphitb3')

fprintf('%g'),disp([S3 3 Q3 cosphitb3])

disp('Xac dinh phu tai tinh toan phan xuong 4');

dl4=[8.5,3,7.5,6.3,10,7.5,40,35;0.37,0.3,0.41,0.45,0.47,0.45,0.53,0.45;0.66,0.58,0.63,0.67,0.7,0.63,0.9,0.58];

disp('ket qua la')

disp(' S4,kVA 4,kW Q4,kVAr cosphitb4')

disp('Tong phu tai tinh toan')

disp('Ket qua la')

disp(' Stt,kVA tt,kW Qtt,kVAr cosphitong')

fprintf('%g'),disp([Stt tt Qtt cosphitong])

KӃt quҧ là:

W ac dinh phu tai tinh toan gom 4 phan xuong gom

W ac dinh phu tai tinh toan phan xuong 1

Phu tai tinh toan tong cua 4 phan xuong

Tong phu tai tinh toan

2 et qua la

Stt,kVA Ptt,kW Qtt,kVAr cosphitong

391.1310 285.6631 267.1705 0.7304

|h̵n xét: Như v̵y tính toán phͭ t̫i đi͏n vͣi s͹ trͫ giúp cͯa chương trình

MATLAB r̭t thu̵n lͫi khi gi̫i các bài toán phͱc t̩p và cho k͇t qu̫ nhanh

chóng.N͇u tính b̹ng tay có th͋ m̭t r̭t nhi͉u thͥi gian

ñ

Sӕ liӋu đo đӃm trên thanh cái trҥm biӃn áp trong mӝt ngày đêm sau mӛi giӡ đưӧc

cho trong bҧng.

t=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23];

¬
6 

=[121,132,124,127,132,168,208,216,254,281,305,312,307,296,315,302,278,256,218,284,323,337,262,158];

Hãy xây dӵng đӗ thӏ phө tҧi và xác đӏnh các tham sӕ cӫa đӗ thӏ theo phương pháp

đӗ thӏ phө tҧi ngày:

Chương trình Matlab như sau:

11 disp('Xay dung do thi phu tai va xac dinh cac tham so theo phuong phap do thi phu tai ngay');

t=0:1:23;

=[121,132,124,127,132,168,208,216,254,281,305,312,307,296,315,302,278,256,218,284,323,337,262,158];

¬
6 

Xác đӏnh phө tҧi dӵ báo trong giai đoҥn đӃn năm thӭ 15:

Chương trình Matlab như sau:

11 t=8:15;

=polyval(c,t);

disp('ket qua la')

disp(' Nam 8 9 10 11 12 13 14 15') fprintf('%g'),disp([])

ÐÐ ÐÐ

Chương trình Matlab như sau:

11 disp('Tong dien nang tieu thu trong linh vuc cong nghiep Viet Nam trong 10 nam');

¬
6 

xlabel('t,Nam');

ylabel('A,GWh');

grid;

legend('Du lieu','Ham hoi quy bac hai','Sai so',2)

ssb=(A-y2).^2; % Xac dinh sai so

| h̵n xét: Như v̵y vͣi s͹ trͫ giúp cͯa chương trình MATLAB r̭t thu̵n lͫi và

cho k͇t qu̫ nhanh chóng Sai s͙ tương đ͙i ͧ mô hình Parabol là 3.81%

¬
6 

 $%&'()'

*TiӃt diӋn dây dүn cӫa mҥng điӋn phân phӕi đưӧc chӑn theo phương pháp hao

tәn điӋn áp cho phép theo biӇu thӭc:

2 1

3

,

.

.

10

mm U

U

l P F

R

m i i





ǻUR thành phҫn tác dөng cӫa hao tәn điӋn áp cho phép,V;

i công suҩt tác dөng trên đoҥn dây thӭ I;

Li chiӅu dài đoҥn dây thӭ I;

¬
6 

U điӋn áp cӫa mҥng điӋn;

Ȗ điӋn dүn suҩt(m/ȍmm2), vӟi Ȗcu = 54, ȖAL = 32, dây nhôm lõi thép Ȗ = 31.5 b) Theo phương pháp chi phí kim loҥi màu cӵc tiӇu, tiӃt diӋn đoҥn dây cuӕi cùng, đưӧc xác đӏnh theo phương pháp:

2,

m

m công suҩt tác dөng trên đoҥn dây cuӕi cùng, kW

THÔ| S͘, S͘ LI͎U CͮA M͠T S͘ LǪI DÂY D̲| CͮA CÁC C̬P ĐI͎|

ÁP

Ta lưu các thông sӕ, sӕ liӋu cӫa các loҥi dây dүn cӫa cҩp điӋn áp vào m.file khi giҧi bài toán cҫn chӑn dây dүn ta chӍ cҫn gӑi ra là đưӧc Cho nên viӋc giҧi toán thuұn tiӋn hơn khi ta không lưu vào trong m.file lúc đó ta phҧi tra bҧng nên rҩt lâu và khó khăn hơn

V tәng vӕn đҫu tư thiӃt bӏ, đ;

C tәng chi phí thưӡng liên, đ/năm;

p hӋ sӕ hiӋu quҧ sӱ dөng và khҩu hao thiӃt bӏ;

p = atc + kkh

kkh hӋ sӕ khҩu hao thiӃt bӏ

atc hӋ sӕ tiêu chuҭn sӱ dөng hiӋu quҧ vӕn, đưӧc xác đӏnh như sau:

¬
6 

1 ) 1 (

) 1 (

+

h h

Th tuәi thӑ cӫa công trình;

i hӋ sӕ chiӃt khҩu, thưӡng lҩy(0.1÷0.2) Vӕn đҫu tư:

?ӕi vӟi trҥm biӃn áp: Vd = (a+bF).l

?ӕi vӟi trҥm biӃn áp: Vb = m + nSb

Vӏ trí tӕi ưu cӫa trҥm biӃn áp

y

Trong đó xi, yi tӑa đӝ điӇm tҧi thӭ i;

i công suҩt cӫa điӇm tҧi thӭ i;

X,Y tӑa đӝ cӫa trҥm biӃn áp phân phӕi

Xác đӏnh như sau:

+ Xác đӏnh điӋn trӣ yêu cҫu

Ryc = UL/Id Trong đó : Id dòng điӋn ngҳn mҥch chҥy trong đҩt, A;

¬
6

UL điӋn áp tính toán có giá trӏ UL<Ucp;

Ucp giá trӏ điӋn áp tiӃp xúc

+ Xác đӏnh điӋn trӣ nhân tҥo:

yc tn

yc tn tao

7 4 ln 2

1 2 (ln

h l d

l l

dcR

R n

.

KXác đӏnh điӋn trӣ cӫa hӋ thӕng nӕi đҩt nhân tҥo:

tao n ng

tao n ng tao

R n

tc Q P Q

pv

b b b

c t P pv

) (

R Q Q

Q

Tҥo m.file voi noi dung nhu sau :

%day dan U=10kV

¬
6

ket qua la

Fc dUt dUcp

70.0000 216.9055 225.0000

| h̵n xét: Như v̵y s͹ trͫ giúp cͯa chương trình MATLAB r̭t thu̵n lͫi và cho k͇t

qu̫ nhanh chóng T͝n th̭t đi͏n áptính toán nh͗ hơn so vͣi gái tr͓ cho phép nên đ̩t yêu c̯u

ñ

Hãy chӑn tiӃt diӋn cho mҥng điӋn áp U=0.38 kV, theo phương pháp chi phí kim loҥi màu cӵc tiӇu, biӃt hao tәn điӋn áp cho phép là dUpt = 7.5% các phө tҧi và chiӅu dài đưӡng dây cho trong bҧng sau:

?oҥn = [ab,bc,bd];

L = [0.33,0.22,0.30];

?iӇm tҧi=[B,C,D];

S=[35+j*21,32+j*18,37+j*22]; gama=32;

Tҥo m.file vӟi nӝi dung như sau :

%day dan U=0.4 kV

[A25]=[1.28+j*0.35];

[A35]=[0.82+j*0.33];

lab

lbc

lbd

C

D

S2

S3S

1

Lưu vӟi tên là DdU04.m

Chương trình Matlab như sau:

11 disp('chon tiet dien day dan cho mang dien ha ap');

¬
6

Fab=real(Sab)*l(1)*10^3/(gama*U*dUrab)

F1=input('Chon ma hieu day A = ');

F1d=input('Go tiet dien day dan Fab = ');

Fab=F1d;

dUabr=real(Sab)*real(F1)*l(1)/U;

dUbcr=dUr-dUabr;

Fbc=real(S(2))*l(2)*10^3/(gama*U*dUbcr)

F2=input('Chon ma hieu day dan A = ');

F2d=input('Go tiet dien day Fbc =');

Fbc=F2d;

Fbd=real(S(3))*l(3)*10^3/(gama*U*dUbcr)

F3=input('Chon ma hieu day dan A = ');

F3d=input('Go tiet dien day Fbd =');

disp('ket qua la')

disp(' Fab Fbc Fbd dU dUcp ')

fprintf('%g'),disp([Fab Fbc Fbd dU dUcp])

Chon ma hieu day A = A240

Go tiet dien day dan Fab = 240

Fbc =

44.10

Chon ma hieu day dan A = A70

Go tiet dien day Fbc =70

Fbd =

69.53

Chon ma hieu day dan A = A70

Go tiet dien day Fbd =70

h, giá thành tәn thҩt cǻ = 1000đ/kWh Tuәi thӑ công trình Th = 25 năm, hӋ sӕ chiӃt khҩu I = 0.1

¬
6

Chương trình Matlab như sau:

11 disp('Xac dinh chi phi quy dan cua duong day 10kV');

DdU10;

U=10; a=158.01;b=0.89;l=8.5;S=650;F=95;Th=25;Tm=4760;kkhd=0.04;ii=0.1; cd=0.001;

Hãy xác đӏnh chi phí quy dүn cӫa trҥm biӃn áp 10/0.4 kV có công suҩt đӏnh mӭc Sn

= 400 kVA, công suҩt cӫa phө tҧi là S = 332 kVA; thӡi gian sӱ dөng công suҩt cӵc

¬
6

đҥi Tm = 4760 h, giá thành tәn thҩt cǻ = 100đ/kWh; tuәi thӑ trҥm biӃn áp Th = 25 năm, hӋ sӕ chiӃt khҩu i =0.1 Các tham sӕ kinh tӃ: m = 19.04*10^6 VND, n = 0.18*10^6 VND/kVA, hӋ sӕ khҩu hao máy biӃn áp kkhb = 6.5%

Chương trình Matlab như sau:

11 disp('Xac dinh chi phi quy dan cua tram bien áp 10/0.4kV');

DdU10;

U=10;Sn=400;m=19.04;n=0.18;cd=0.001;Th=25;ii=0.1;Tm=4760;Spt=332;kkhb=0.065;

disp('ket qua cap dien ap 10kV la')

disp(' Fc10 dUt10 dUcp10 ')

fprintf('%g'),disp([Fc10 dUt10 dUcp10]);

¬
6

disp('ket qua cap dien ap 22kV la')

disp(' Fc22 dUt22 dUcp22 ')

fprintf('%g'),disp([Fc22 dUt22 dUcp22])

disp('ket qua cap dien ap 35kV la')

disp(' Fc35 dUt35 dUcp35 ')

fprintf('%g'),disp([Fc35 dUt35 dUcp35])

KӃt quҧ màn hình như sau:

| h̵n xét: ãng vͣi c͡t thͱ nh̭t, thͱ hai và thͱ ba là chi phí quy d̳n cͯa các c̭p

đi͏n áp 10; 22 và 35 kV Các s͙ li͏u tính toán cho th̭y chi phí quy d̳n cͯa tr̩m bi͇n áp 10/0.4kV có giá tr͓ th̭p nh̭t T͝ng chi phí quy d̳n cͯa phương án10kV có giá tr͓ nh͗ nh̭t(Z 10 = 150.47.10 6 VND)=> nên phương án 10kV là t͙i ưu

Chương trình Matlab như sau:

11 disp('Xac dinh vi tri toi uu cua tram bien ap');

¬
6 

Mӝt xí nghiӋp sҧn xuҩt có công suҩt tính toán năm hiӋn tҥi là S0 = 655kVA, suҩt tăng phө tҧi trung bình hàng năm là a = 5%, hàm dӵ báo dҥng tuyӃn tính theo thӡi gian:St = S0*(1+a*(t-t0)); hӋ sӕ công suҩt coi là không đәi và bҵng cosij = 0.83 trong đó phө tҧi loҥi I và loҥi II chiӃm 72% (m1+2 = 0.75), thӡi gian sӱ dөng công suҩt cӵc đҥi TM = 5750 h, hӋ sӕ chiӃt khҩu i = 0.1; giá thành tәn thҩt điӋn năng cǻ = 1000đ/kWh; suҩt thiӋt hҥi do mҩt điӋn gth = 5000đ/kWh, chu kǤ tính toán T = 7 năm Hãy chӑn sӕ lưӧng và công suҩt máy biӃn áp cӫa trҥm biӃn áp 22/0.4kV cung cҩp cho xí nghiӋp

Chương trình Matlab như sau:

11 disp('Chon so luong va cong suat may bien ap cua tram bien ap 22/0.4 kV'); disp('Xac dinh phu tai du bao nam thu 7');

¬
6 

Chon so luong va cong suat may bien ap cua tram bien ap 22/0.4 kV

W ac dinh phu tai du bao nam thu 7

S7,kVA

525.4200

W ac dinh chi phi quy dan cua cac phuong an

2et qua la:

Chương trình Matlab như sau:

11 disp('Tinh toan noi dat cho tram bien ap 22/0.4 kV');

nuyc=input('he so su dung nuyc = ');

nuyng=input(' he so su dung nuyng = ');

Chương trình Matlab như sau:

11 disp('Tinh toan noi dat cho tram bien ap 110/22 kV');

¬
6 

la=L/ndc;

ka=la/lc;

nuyc=input('he so su dung nuyc(he so su dung dien cuc) = ');

nuyc=input('he so su dung nuyng(he so su dung thanh noi) = ');

Chon gia tri(so luong dien cuc) n2 = 22

he so su dung nuyc(he so su dung dien cuc) = 0.3

he so su dung nuyng(he so su dung thanh noi) = 0.2

n Rd,ohm Fmin,mm2 F,mm2

22.00 0.37 11.23 400.00

ñ

tӍ lӋ chi phí khҩu hao và thu hӗi vӕn p = 0.12; thӡi gian vұn hành trong năm t = 8760; thӡi gian hao tәn cӵc đҥi IJ = 3580 h

Chương trình Matlab như sau:

11 disp('Xac dinh dung luong bu toi uu cho mgna dien ha ap 0.38 kV');

Matlab tính toán phân phӕi công suҩt phҧn kháng bù cho các phân xưӣng, biӃt tәng tәn thҩt công suҩt cҫn bù là Qb = 155 kVAr, coi hӋ sӕ đӗng thӡi bҵng 1

Chương trình Matlab như sau:

11 disp('Tinh toan cong suat phan khang bu cho cac phan xuong mang dien ap 0.38 kV');

2

Q b2 3

Q b3

1

P + jQ

2 2

P + jQ

3 3

Q b¬

A

l

l

l 1

2

3

¬
6 

Chương trình Matlab như sau:

11 disp('Xac dinh cac chi tieu danh gia du an');

Vtl=(Vv-Vtv).*Is; % tra von = C hang nam

Ckh=V/(length(t)-1); % chi phi khau hao

sӕ vӕn cӫa dӵ án V™ = 2058*10^6 VND, trong đó vӕn tӵ có là Vtc = 1058*10^6,

HӋ sӕkO&M = 2% Thӡi gian tính toán công trình N = 10 năm

Chương trình Matlab như sau:

11 disp('Xac dinh cac chi tieu kinh te tai chinh cua du an xay dung mang dien');

m=468.83;Tm=4750;a=0.05;dA=0.076;V=2058;st=0.25;ii=0.1;cm=0.00058; cb=0.0008;;kom=0.02;Cvh=105;N=10;Vtc=1058;ls=0.08;

¬
6 

ñ

Cũng vӟi sӕ liӋu bài tұp 15.13, hãy xác đӏnh các chӍ tiêu kinh tӃ tài chính và cho nhұn xét trong các trưӡng hӧp:

Bài toán gi̫i như bài toán 15.13 2hi gi̫i vͣi trưͥng hͫp b thì thay Vtc =0

a) Toàn b͡ v͙n đ̯u tư là v͙n t͹ có

W ac dinh cac chi tieu kinh te tai chinh cua du an xay dung mang dien

NPV,10^6 VND R

1122.45 1.08

b) Toàn b͡ v͙n đ̯u tư là v͙n vay

W ac dinh cac chi tieu kinh te tai chinh cua du an xay dung mang dien

Chương trình Matlab sau:

11 disp(''Xac dinh dung luong bu... data-page="71">

Matlab tính tốn phân phӕi cơng suҩt phҧn kháng bù cho phân xưӣng, biӃt tәng tәn thҩt công suҩt cҫn bù Qb = 155 kVAr, coi hӋ sӕ đӗng thӡi bҵng

Chương trình Matlab sau: ...

Chon ma hieu day dan A = A70

Go tiet dien day Fbd =70

h, giá thành tәn thҩt cǻ = 1000đ/kWh Tuәi thӑ cơng trình Th = 25 năm, hӋ sӕ chiӃt khҩu I = 0.1