Phương pháp giải mạch với Matlab

2;R2=4;R3=2;R4=5;R5=2;Va=24;Vb=36; 
R=zeros(2,2); 
R=A_array(1,0,R1+R2,R); 
R=A_array(2,0,R4+R5,R); 
R=A_array(1,2,R3,R); 
V=zeros(2,1); 
V=b_array(1,0,Va,V); 
V=b_array(2,0,-Vb,V); 
I=R\V; 
disp(['dong dien vong 1=',num2str(I(1)),' A']); 
disp(['dong dien vong 2=',num2str(I(2)),' A']); 
disp(['dong dien I1=',num2str(I(1)),' A']); 
disp(['dong dien I2=',num2str(I(2)),' A']); 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 9 
disp(['dong dien I3=',num2str(I(1)-I(2)),' A']); 
kết quả 
dong dien vong 1=-2.1176 A 
dong dien vong 2=3.5294 A 
dong dien I1=-2.1176 A 
dong dien I2=3.5294 A 
dong dien I3=-5.6471 A 
Bài 11 
I3
Vd=24V
Vb=6V
R2=1k
Va=12V
R1=2k
R3=3k
I1
Vc=18V
R7=1k
R6=2k
I2R5=4k
R4=1k
 function bai11 
R1=2;R2=1;R3=3;R4=1;R5=4;R6=2;R7=1; 
Va=12;Vb=6;Vc=18;Vd=24; 
R=zeros(3,3); 
R=A_array(1,0,R1+R2,R); 
R=A_array(2,0,R5,R); 
R=A_array(3,0,R6+R7,R); 
R=A_array(1,2,R3,R); 
R=A_array(1,3,R4,R); 
V=zeros(3,1); 
V=b_array(1,0,-Va-Vd,V); 
V=b_array(2,0,-Vb-Vc,V); 
V=b_array(3,0,Vd+Vc,V); 
I=R\V; 
disp(['dong dien vong 1= ',num2str(I(1)),' A']); 
disp(['dong dien vong 2= ',num2str(I(2)),' A']); 
disp(['dong dien vong 3= ',num2str(I(3)),' A']); 
kết quả 
dong dien vong 1= -6.549 A 
dong dien vong 2= -6.2353 A 
dong dien vong 3= 8.8627 A 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 10 
Bài 12 
R5
2k
I3
R6
4k
I1
Va=6V R3
3k
Ve=24V
I2
Vd=18V
I=?
Vc=12V
R4
2k
R8
1k
R2
3k
Vb=10V
R7
6k
R1
5k
I4
R9
7k
function bai12 
R1=5;R2=3;R3=1;R4=2;R5=2;R6=4;R7=6;R8=1;R9=7; 
Va=6;Vb=10;Vc=12;Vd=18;Ve=24; 
R=zeros(4,4); 
R=A_array(1,0,R1,R); 
R=A_array(2,0,R3+R4,R); 
R=A_array(3,0,R9,R); 
R=A_array(4,0,R6+R7,R); 
R=A_array(1,2,R2,R); 
R=A_array(2,3,R5,R); 
R=A_array(3,4,R8,R); 
V=zeros(4,1); 
V=b_array(1,0,-Va-Vc,V); 
V=b_array(2,0,Vb,V); 
V=b_array(3,0,Ve-Vd,V); 
V=b_array(4,0,Vc+Vd,V); 
I=R\V; 
disp(['dong dien vong 1= ',num2str(I(1)),' A']); 
disp(['dong dien vong 2= ',num2str(I(2)),' A']); 
disp(['dong dien vong 3= ',num2str(I(3)),' A']); 
disp(['dong dien vong 4= ',num2str(I(4)),' A']); 
disp(['dong dien I= ',num2str(I(1)-I(2)),' A']); 
kết quả 
dong dien vong 1= 1.9596 A 
dong dien vong 2= -0.77441 A 
dong dien vong 3= -1.037 A 
dong dien vong 4= -2.8215 A 
dong dien I= 2.734 A 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 11 
Bài 13 
I3
R1=6k
R3
1k
R2=3k
R5
2k
Vo=?
R4=4k
I1
R6
3k
R7
6k
Va=72V
I2
function bai13 
R1=6;R2=3;R3=1;R4=4;R5=2;R6=3;R7=6; 
Va=72; 
R=zeros(3,3); 
R=A_array(1,0,R1+R4,R); 
R=A_array(2,0,0,R); 
R=A_array(3,0,R2+R3,R); 
R=A_array(1,2,R6,R); 
R=A_array(2,3,R7,R); 
R=A_array(1,3,R5,R); 
V=zeros(3,1); 
V=b_array(1,0,-Va,V); 
V=b_array(2,0,0,V); 
V=b_array(3,0,0,V); 
I=R\V; 
disp(['dong dien vong 1= ',num2str(I(1)),' A']); 
disp(['dong dien vong 2= ',num2str(I(2)),' A']); 
disp(['dong dien vong 3= ',num2str(I(3)),' A']); 
disp(['dong dien I= ',num2str(I(3)),' A']); 
disp(['dien ap Vo= ',num2str((I(3)*R3)),' V']); 
kết quả 
dong dien vong 1= 6 A 
dong dien vong 2= 4 A 
dong dien vong 3= 3 A 
dong dien I= 3 A 
dien ap Vo= 3 V 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 12 
Bài 14 
R4
3k
I2Va=24V
Vo=?
R3=1k
I1
Vb=12V
R2
6k
I=?
R1=3k
 function bai14 
R1=3;R2=6;R3=1;R4=3; 
Va=24;Vb=12; 
R=zeros(2,2); 
R=A_array(1,0,R1,R); 
R=A_array(2,0,R3+R4,R); 
R=A_array(1,2,R2,R); 
V=zeros(2,1); 
V=b_array(1,0,-Va-Vb,V); 
V=b_array(2,0,Vb,V); 
I=R\V; 
disp(['dong dien vong 1= ',num2str(I(1)),' A']); 
disp(['dong dien vong 2= ',num2str(I(2)),' A']); 
disp(['dong dien I= ',num2str(I(1)-I(2)),' A']); 
disp(['dien ap Vo= ',num2str(I(2)*R4),' A']); 
kết quả 
dong dien vong 1= 5.3333 A 
dong dien vong 2= 2 A 
dong dien I= 3.3333 A 
dien ap Vo= 6 A 
Bài 15 
R3
2k
Vb=12V
R4
8k
Va=24V
Vo=?
R2=4k
I1
R1=1k
I2
function bai15 
R1=1;R2=4;R3=2;R4=8; 
Va=24;Vb=12; 
R=zeros(2,2); 
R=A_array(1,0,R1+R4,R); 
R=A_array(2,0,R2,R); 
R=A_array(1,2,R3,R); 
V=zeros(2,1); 
V=b_array(1,0,Va,V); 
V=b_array(2,0,Vb-Va,V); 
I=R\V; 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 13 
disp(['dong dien vong 1= ',num2str(I(1)),' A']); 
disp(['dong dien vong 2= ',num2str(I(2)),' A']); 
disp(['dong dien I= ',num2str(I(1)-I(2)),' A']); 
disp(['dien ap Vo= ',num2str((I(1)-I(2))*R3),' V']); 
kết quả 
dong dien vong 1= -1.9355 A 
dong dien vong 2= 1.3548 A 
dong dien I= -3.2903 A 
dien ap Vo= -6.5806 V 
Bài 16 
R6=0.5k
Vo=?
I2
Va=6V
R4
0.5k
R3
1k
R2=0.5k
I3I1
I4
R5
0.5k
R1=1k Vb=12V
I=?
 function bai16 
R1=1;R2=0.5;R3=1;R4=0.5;R5=0.5;R6=0.5; 
Va=6;Vb=12; 
R=zeros(4,4); 
R=A_array(1,0,R1,R); 
R=A_array(2,0,0,R); 
R=A_array(3,0,R5,R); 
R=A_array(4,0,R6,R); 
R=A_array(1,2,R3,R); 
R=A_array(2,3,R4,R); 
R=A_array(3,4,R2,R); 
V=zeros(4,1); 
V=b_array(1,0,-Va,V); 
V=b_array(2,0,-Vb,V); 
V=b_array(3,0,0,V); 
V=b_array(4,0,Vb,V); 
I=R\V; 
disp(['dong dien vong 1= ',num2str(I(1)),' A']); 
disp(['dong dien vong 2= ',num2str(I(2)),' A']); 
disp(['dong dien vong 3= ',num2str(I(3)),' A']); 
disp(['dong dien vong 4= ',num2str(I(4)),' A']); 
disp(['dong dien I= ',num2str(I(3)),' A']); 
disp(['dien ap Vo= ',num2str((I(3)*R5)),' V']); 
kết quả 
dong dien vong 1= 10.875 A 
dong dien vong 2= 15.75 A 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 14 
dong dien vong 3= 1.5 A 
dong dien vong 4= -11.25 A 
dong dien I= 1.5 A 
dien ap Vo= 0.75 V 
Bài 17 
Va=60V
V2I3
R2
12k
R1
2k
R3
2k
I5
I4
R6
1k
R4
6k
V1 V3R5
2k
I2
I1
Vo=?
 function bai17 
R1=2;R2=12;R3=2;R4=6;R5=2;R6=1; 
G1=1/R1;G2=1/R2;G3=1/R3;G4=1/R4;G5=1/R5;G6=1/R6; 
Va=60; 
G=zeros(3,3); 
G=A_array(1,0,G1+G2,G); 
G=A_array(1,2,G3,G); 
G=A_array(2,0,G4,G); 
G=A_array(2,3,G5,G); 
G=A_array(3,0,G6,G); 
I=zeros(3,1); 
I=b_array(0,1,Va/R1,I); 
V=G\I; 
disp(['Dien ap V1= ',num2str(V(1)),' V']); 
disp(['Dien ap V2= ',num2str(V(2)),' V']); 
disp(['Dien ap V3= ',num2str(V(3)),' V']); 
disp(['Dong dien I1= ',num2str((Va-V(1))/R1),' A']); 
disp(['Dong dien I2= ',num2str(V(1)/R2),' A']); 
disp(['Dong dien I3= ',num2str((V(1)-V(2))/R3),' A']); 
disp(['Dong dien I4= ',num2str(V(2)/R4),' A']); 
disp(['Dong dien I5= ',num2str((V(2)-V(3))/R5),' A']); 
disp(['Dien ap Vo= ',num2str(((V(2)-V(3))/R5)*R6),' A']); 
Dien ap V1= 36 V 
Dien ap V2= 18 V 
Dien ap V3= 6 V 
Dong dien I1= 12 A 
Dong dien I2= 3 A 
Dong dien I3= 9 A 
Dong dien I4= 3 A 
Dong dien I5= 6 A 
Dien ap Vo= 6 A 
Bài 18 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 15 
R1
3k
Va=24V
V1
I2
Ia 5A
I1
I3
R2
6k
R3
12k
 function bai18 
R1=3;R2=6;R3=12; 
G1=1/R1;G2=1/R2;G3=1/R3; 
Va=24;Ia=5; 
G=zeros(1,1); 
G=A_array(1,0,G1+G2+G3,G); 
I=zeros(1,1); 
I=b_array(0,1,(Va/R3)+Ia,I); 
V=G\I; 
disp(['Dien ap V1= ',num2str(V(1)),' V']); 
disp(['Dong dien I1= ',num2str(V(1)/R1),' A']); 
disp(['Dong dien I2= ',num2str(V(1)/R2),' A']); 
disp(['Dong dien I3= ',num2str((Va-V(1))/R3),' A']); 
Dien ap V1= 12 V 
Dong dien I1= 4 A 
Dong dien I2= 2 A 
Dong dien I3= 1 A 
Bài 19 
I4
R2 2k
Ib 5A
R1
4k
R3
4k
Vb=?
Va=?
I1I2 V2
I3
V1
Ia
3A
Ic 6A
 function bai19 
R1=4;R2=2;R3=4; 
G1=1/R1;G2=1/R2;G3=1/R3; 
Ia=3;Ib=5;Ic=6; 
G=zeros(2,2); 
G=A_array(1,0,G2,G); 
G=A_array(2,0,G3,G); 
G=A_array(1,0,G1,G); 
I=zeros(2,1); 
I=b_array(0,1,Ia,I); 
I=b_array(1,2,Ib,I); 
I=b_array(0,2,Ic,I); 
V=G\I; 
Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 16 
disp(['Dien ap V1= ',num2str(V(1)),' V']); 
disp(['Dien ap V2= ',num2str(V(2)),' V']); 
disp(['Dong dien I2= ',num2str((V(2)-V(1))/R1),' A']); 
disp(['Dong dien I4= ',num2str(-V(1)/R2),' A']); 
disp(['Dien ap Va= ',num2str(((V(2)-V(1))/R1)*R1),' V']); 
disp(['Dien ap Vb= ',num2str(V(1)),' V']); 
Dien ap V1= -2.6667 V 
Dien ap V2= 44 V 
Dong dien I2= 11.6667 A 
Dong dien I4= 1.3333 A 
Dien ap Va= 46.6667 V 
Dien ap Vb= -2.6667 V 
Bài 20 
function bai20 
ZC=-2j; R2=1; ZL=j; R4=0.5; Ia=1+j; 
G1=1/ZC; G2=1/R2; G3=1/ZL; G4=1/R4; 
G=zeros(2,2); 
G=A_array(1,0,G1+G2,G); 
G=A_array(1,2,G3,G); 
G=A_array(2,0,G4,G); 
I=zeros(2,1); 
I=b_array(0,2,Ia,I); 
U=G\I; 
disp(['Dien the U1=',num2str(U(1)), '(V)']); 
disp(['Do lon : ',num2str(abs(U(1))), '(V)']); 
disp(['Goc pha : ',num2str(rad2deg(angle(U(1)))), '(do)']); 
disp(['Dien the U2=',num2str(U(2)), '(V)']); 
disp(['Do lon :',num2str(abs(U(2))), '(V)']); 
disp(['Goc pha : ',num2str(rad2deg(angle(U(2)))), '(do)']); 
disp(['Dong dien I1=',num2str(U(1)/ZC), '(A)']); 
disp(['Do lon :',num2str(abs(U(1)/ZC)), '(A)']); 
disp(['Goc pha :',num2str(rad2deg(angle(U(1)/ZC))), '(do)']); 
disp(['Dong dien I2=',num2str(U(1)/R2), '(A)']); 
disp(['Do lon :',num2str(abs(U(1)/R2)), '(A)']); 
disp(['Goc pha :',num2str(rad2deg(angle(U(1)/R2))), '(do)']); 
disp(['Dong dien I3=',num2str((U(1)-U(2))/ZL), '(A)']); 
disp(['Do lon :',num2str(abs((U(1)-U(2))/ZL)), '(A)']); 
disp(['Goc pha :',num2str(rad2deg(angle((U(1)-U(2))/ZL))), '(do)']); 
disp(['Dong dien I4=',num2str(U(2)/R4), '(A)']); 
disp(['Do lon :',num2str(abs(U(2)/R4)), '(A)']); 
disp(['Goc pha :',num2str(rad2deg(angle(U(2)/R4))), '(do)']); 
3I

1I

-2j  1aI j A

  
4I

R2 
1 
j 
1U

R4 
0.5
0 
2I

2U

Phương pháp giải Mạch 1 với Matlab Biên soạn: Ths Phan Như Quân 
 Trang 17 
kết quả 
Dien the U1=0.43902-0.04878i(V) 
Do lon : 0.44173(V) 
Goc pha : -6.3402(do) 
Dien the U2=0.26829+0.41463i(V) 
Do lon :0.49386(V) 
Goc pha : 57.0948(do) 
Dong dien I1=0.02439+0.21951i(A) 
Do lon :0.22086(A) 
Goc pha :83.6598(do) 
Dong dien I2=0.43902-0.04878i(A) 
Do lon :0.44173(A) 
Goc pha :-6.3402(do) 
Dong dien I3=-0.46341-0.17073i(A) 
Do lon :0.49386(A) 
Goc pha :-159.7751(do) 
Dong dien I4=0.53659+0.82927i(A) 
Do lon :0.98773(A) 
Goc pha :57.0948(do)