Bài thí nghiệm 1 | Lý thuyết mạch 1 | Đại học Bách Khoa Hà Nội
Bài thí nghiệm 1 | Lý thuyết mạch 1 | Đại học Bách Khoa Hà Nội. Tài liệu được biên soạn giúp các bạn tham khảo, củng cố kiến thức, ôn tập và đạt kết quả cao kết thúc học phần. Mời các bạn đọc đón xem!
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