


Preview text:
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 Bài 1: Code Mathlab:
j = sqrt(-1); E1 = 100; E2 = 220*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.4; Z6 = 10+j*20;
Z35 = j*2*pi*60*0.6*sqrt(0.2*0.4); 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]; Phương pháp dòng vòng Phương pháp thế nút
A = [1 0 1 1 0 0;0 1 1 0 1 0;0 0 0 -1 1 1];
B = [-1 0 0 1 0 1;0 0 1 -1 -1 0;1 1 -1 0 0 0]; Zv = A*Znh*A'; Ynh = inv(Znh); Ev = A*(Enh - Znh*Jnh); Ynut = B*Ynh*B'; Iv = Zv\Ev; Jnut = B*(Jnh - Ynh*Enh); Inh = A'*Iv Vnut = Ynut\Jnut; Unh = Znh*(Inh + Jnh) - Enh Unhn = B'*Vnut
Sng = (Inh + Jnh)'*Enh + Jnh'*Unh
Inhn = Ynh*(Unhn + Enh) – Jnh
Sng = (Inhn + Jnh)'*Enh + Jnh'*Unhn Kết quả thu được : Kết quả thu được : Inh = Inhn = -2.8620 - 3.0433i -2.8620 - 3.0433i 3.9151 + 2.8309i 3.9151 + 2.8309i 1.0531 - 0.2124i 1.0531 - 0.2124i 1.2749 - 0.0992i 1.2749 - 0.0992i -0.2218 - 0.1132i -0.2218 - 0.1132i -4.1369 - 2.9441i -4.1369 - 2.9441i Unh = Unhn = 1.0e+02 * 1.0e+02 * -0.6413 - 2.0578i -0.6413 - 2.0578i -0.6001 - 0.9475i -0.6001 - 0.9475i 0.3379 + 0.6309i 0.3379 + 0.6309i 0.3034 + 1.4269i 0.3034 + 1.4269i 0.2622 + 0.3167i 0.2622 + 0.3167i 0.0412 + 1.1103i 0.0412 + 1.1103i Sng = Sng = 1.2746e+03 + 1.6798e+03i 1.2746e+03 + 1.6798e+03i Bài 2: Code Mathlab:
j = sqrt(-1); E1 = 220*sqrt(2)*exp(j*0); E5=200; Enh = [E1;0;0;0;E5]; Jnh = [0;0;0;0;0]; Z1 = 200; Z2 = 200; Z3 = 10; Z4 = -j*100; Z5 = j*100;
Znh = [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]; Phương pháp dòng vòng Phương pháp thế nút
A = [1 1 0 0 0;0 1 -1 1 0;0 0 0 -1 1]; B = [-1 1 1 0 0;1 -1 0 1 1]; Zv = A*Znh*A'; Ynh = inv(Znh); Ev = A*(Enh - Znh*Jnh); Ynut = B*Ynh*B'; Iv = Zv\Ev; Jnut = B*(Jnh-Ynh*Enh); Inh = A'*Iv; Vnut = Ynut\Jnut; Unh = Znh*(Inh + Jnh) - Enh; Unhn = B'*Vnut;
Inhn = Ynh*(Unhn + Enh) - Jnh; I3 = Inh(3) I3 = Inhn(3) UBC = Unh(4) UBC = Vnut(2) Pe1 = E1*Inh(1) Pe1 = E1*Inhn(1) Pe5 = E5*Inh(5) Pe5 = E5*Inhn(5) Kết quả thu được: Kết quả thu được: I3 = I3 = 0.0000+ 2.0000i 0.0000+ 2.0000i UBC = UBC = -1.5556e+02 + 2.2000e+02i -1.5556e+02 + 2.2000e+02i Pe1 = Pe1 = 2.4200e+02 + 3.1113e+02i 2.4200e+02 + 3.1113e+02i Pe5 = Pe5 = 4.4000e+02 - 8.8873e+01i 4.4000e+02 - 8.8873e+01i