ĐẠI HỌC BÁCH KHOA HÀ NỘI
TRƯỜNG ĐIỆN – ĐIỆN TỬ
BÁO CÁO THÍ NGHIỆM
LÍ THUYT MẠCH I
HỌC KỲ 2024.1
Họ và tên: Nguyễn Khắc Long
MSSV: 20232169
Lớp: Kỹ thuật điều khiển và tự động hóa – EE2-05
Mã học phần : EE2021
Mã lớp: 749363
Hà Nội 12/2024
BÀI THÍ NGHIỆM LTM : 01
TÍNH CHẾ ĐỘ XÁC LẬP ĐIỀU HÒA TRONG MẠCH ĐIỆN
TUYẾN TÍNH BẰNG MÁY TÍNH DÙNG PHẦN MỀM MATLAB
I. Mục đích thí nghiệm
Bài thí nghiệm này giúp cho sinh viên biết sử dụng chương trình
MATLAB để tính chế độ xác lập điều hòa trong mạch điện tuyến tính
bằng.
1) Phương pháp dòng vòng
2)Phương pháp thế nút
II. Nội dung thí nghiệm
Bài 1: Mạch hình bên:
Z1 = 30+j40(Ω);
Z2 = 20+j10(Ω)
R3=10; R4=15; R5=20 (Ω)
L3= 0.2; L4=0.3; L5=0.4 (H)
Z6=10+20j (Ω)
M35=M53=0.6
√
L3L5
˙
J6
= 10∠30° ( A)
˙
E1
= 100 (V)
˙
E2
= 200∠60°(V)
f =60(Hz)
Code MATLAB:
B = [1 0 1 1 0 0;0 1 1 0 1 0;0 0 0 -1 1 1];
j = sqrt(-1);
pi = 3.1415;
E1 = 100; E2 = 200*exp(j*pi/3);
Enh = [E1;E2;0;0;0;0];
J6 = 10*exp(j*pi/6);
Jnh = [0;0;0;0;0;-J6];
Z1 = 30+j*40;
Z2 = 20+j*10;
Z3 = 10+j*2*pi*60*0.2;
Z4 = 15+j*2*pi*60*0.3;
Z5 = 20+j*2*pi*60*0.5;
Z6 = 10+j*20;
Z35 = -j*2*pi*60*0.6*sqrt(0.2*0.5);
Z53=Z35;
Znh = [Z1 0 0 0 0 0 ; 0 Z2 0 0 0 0 ; 0 0 Z3 0 Z35 0 ; 0 0 0 Z4 0 0 ; 0
0 Z53 0 Z5 0; 0 0 0 0 0 Z6];
Zv = B*Znh*B';
Ev = B*(Enh + Znh*Jnh);
Iv = Zv\Ev;
Inh = B'*Iv
Unh = Znh*(Inh + Jnh) - Enh
Sng = (Inh + Jnh)'*Enh + Jnh'*Unh
Kết quả hiện ở Command Window
Inh =
-2.5683 - 3.0815i
4.2296 + 2.8496i
1.6613 - 0.2318i
1.0440 - 0.1117i
0.6173 - 0.1201i
-3.6123 - 2.9698i
Unh =
1.0e+02 *
-0.5379 - 1.9518i
-0.4391 - 0.7391i
0.2550 + 0.7878i
0.2829 + 1.1639i
0.1841 - 0.0487i
0.3667 - 3.2515i
Sng =
1.9679e+03 + 3.7550e+03i
Bài 2: Mạch hình bên:
R1 = R2 = 200; R3 = 10 (Ω)
ωL =
1
ωC
= 100 (Ω)
f = 50 (Hz)
E5 = 200 (V) = const
E1 = 200
√2
.sin(ωt) (V)
Code MATLAB:
j=sqrt(-1);
R1 = 200; R2 = 200; R3 =10;
ZL = 100*j; ZC = -100*j;
E1 = 200; %phuc hoa nguon e1
E5 = 200; %nguon mot chieu
disp('xet chi co E1 tac dung')
A = [1 0;-1 0;-1 1;0 -1;0 1];
Z1 = R1;Z2 = R2;Z3 = R3; Z4 = -100*j;Z5 = 100*j;
Z = [Z1 0 0 0 0;
0 Z2 0 0 0;
0 0 Z3 0 0;
0 0 0 Z4 0;
0 0 0 0 Z5];
En = [E1;0;0;0;0];
Jd = [0;0];
Jn = [0;0;0;0;0];
Ynh = A'*inv(Z)*A;
Inh = A'*inv(Z)*En+Jd;
disp('dien the tai cac nut')
Vd = Ynh\Inh
Un = -A*Vd;
In = inv(Z)*(Un+En);
disp('dong dien qua R3')
I31 = In(3)
disp('hieu dien the Ubc')
Ubc = Un(3)
disp('cong suat Pe1')
Pe1 = In'*En
% Xet chi co E5 tac dung
disp('xet chi co E5 tac dung')
Z12=(R1*R2)/(R1+R2);
I32=-E5/(R3+Z12)
Pe5=-I32*E5
Kết quả hiện ở Command Window
xet chi co E1 tac dung
dien the tai cac nut
Vd =
100.0000
100.0000
dong dien qua R3
I31 =
-1.4211e-15
hieu dien the Ubc
Ubc =
-1.4211e-14
cong suat Pe1
Pe1 =
100.0000
xet chi co E5 tac dung
I32 =
-1.8182
Pe5 =
363.6364