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

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 1


BÁO CÁO THӴC HANH MATLAP
Sinh viên:Bùi Bҧo Linh.
STT:33.


















TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 2

I. Mӣ đҫu.

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.













TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 3

Phҫn I: Thӵc hành MATLAB cơ bҧn
BÀI THӴC HÀNH SӔ 1

Chương 1: Cơ sӣ MATLAB
I. Cơ sӣ lý thuyӃt
Ngôn ngӳ cӫa hӋ thӕng MATLAB rҩt phong phú, nó gҫn như áp dөng tҩt cҧ các
phương tiӋn lұp trình đã biӃt, kӇ cҧ mô phӓng, hoҥt hình, đӗ hӑa«
MATLAB chӭa khӕi lưӧng khәng lӗ các phép tính, thuұt giҧi và các hàm cho
phép giҧi rҩt nhiӅu các bài toán phӭc tҥp như nghӏch đҧo ma trұn, đҥo hàm, vi
phân, tích phâ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.
1.1) Sӕ và phép toán sơ cҩp
a) Các phép toán thông dͭng thưͥng dùng căn b̫n như sau:
Phép
toán
Cӝng
Trӯ
Nhân
Chia
Chia
trái
Lũy
thӯa
Khia
căn
Ký
hiӋu
+
-
*
/
\

^
sqrt

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:
Tham sӕ
Cú pháp
Tham sӕ
Cú pháp
Phҫn thӵc
P = real(S)
Môdul
Sm = abs(S)
Phҫn ҧo
Q = imag(S)
Góc pha
Theta = angle(S)
Liên hӧp phӭc
Slh = conj(S)




TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 4

c) Các bi͇n và hàm trong MATLAB
Các hàm toán hӑc
LӋnh

Hàm
LӋnh
Hàm
abs(x)
Giá trӏ tuyӋt đӕi
rem(x,y)
Sӕ dư cӫa phép chia x/y
exp(x)
Hàm mũ cơ sӕ e
round(x)
Làm tròn sӕ
imag(x)
Phҫn ҧo
Ceil(x)
Làm tròn lên
real(x)
Phҫn thӵc
floor(x)
Làm tròn xuӕng
Phase(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)

Phҫn tӱ vector bé nhҩt
sqtr(x)
Căn bұc hai
max)v)
Phҫ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
Sin(x)
Hàm sin
atan(x)
Hàm artg ±90
0

Cos(x)
Hàm cosin
atan2(x,y)
Hàm artg ±180
0

Tan(x)
Hàm tg
sinc(x)
Hàm sin(ɉx)/(ɉx)
Asin(x)
Hàm arsin
acos(x)
Hàm arcos

Sinh(x)
Hàm sin hyperbol
asinh(x)
Hàm arsin hyperbol
Cosh(x)
Hàm cosin hyperbol
acosh(x)
Hàm arccos yperbol

1.2 Các dҥng xuҩt dӳ liӋu
Đưӧc thӵc hiӋn bӣi lӋnh :
fprintf('S=%v.uf d\n',F)
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 5

trong đó S là tham sӕ ra, v là chӳ sӕ cӵc đҥi cӕ thӇ, u là chӳ sӕ sau dҩu phҭy, f
hiӋn thӏ dҩu phҭy tĩnh, (e hiӋn thӏ dҩu phҭy đӝng), d hiӋn thӏ đӧn vӏ tính và F hàm
tính toán.
1.3 Văn bҧn trong MATLAB
ĐӇ xuҩt chuӛi ký tӵ nào đó ra màn hình thì ta sӱ dөng lӋnh như sau:
disp('string')
Bài tұp

Bai 1.1:Hãy sӱ dөng các lӋnh trong matlap dӇ thӵc hiӋ các phép tính.
a)25,4+17(34/4,2)-2,5;
b)cos(x)+sin(y);vӟi 1,5ʌ,y=6,4;
c)sin
2
(x)-cos

1/2
(y);vӟi x=1,5ʌ,y=6,4.
d)e
2,5
+tg(x);vӟi x=1,2ʌ.
e)S=(2,7+j3,2)-(2+j1,5).
a)25.4+17*(34/4.2)-2.5
b)x=5.3;
y=1.2;
cos(x)+sin(y)
c) x=1.5*pi;
y=6.4;
(sin(x))^2-(cos(x))^1/2
d) x=1.2*pi;
exp(2.5)+tan(x)
e) (2.7+3.2j)-(2+1.5j)
f) a=23.76;b=16.33;c=5.67;d=6.89;
x=1.35;
y=((cos(a^3+x^2))^2-log(b)+d)/sqrt(c^2+2*x^2);
fprintf(
'y=%7.8f\n',y)
Bài 1.2:Mӝt mҥch đҥi lưӧng xoay chiӅu có các đҥi lưӧng điӋn áp U=U
R
+jU
X
V
và giá trӏ tәng trӣ Z=R+jX ;
a)Thao tác hiӇn thӏ tên ngưӡi giҧ trên màn hình;
b)Xác đӏnh giá trӏ dòng điӋn chҥy trong mҥch;
c)BiӇu thӏ giá trӏ tính toán dưӟi dҥng I

m
jĮ
,vӟi I
m
là modun và Į là góc pha;

ten=
'Bui Bao Linh'
UR=122.5;UX=117.5;
R=4.25;X=5.34;
U=UR+i*UX;
Z=R+i*X;
I=U/Z;
Im=abs(I);
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 6

anpha=angle(I);
fprintf(
'c=%s,Im=%f,anpha=%f\n',c,Im,anpha)

Bài 1.3
P=66.78;Q=35.67;
U=35;R=4.75;
to=4620;n=2;
DeltaP=((P^2+Q^2)/U^2)*R;
DeltaA=DeltaP*to;
fprintf(
'DeltaA=%7.2f kWh \n',DeltaA)

KӃt quҧ:
DeltaA=102683.20 kWh
Bài 1.4
a)Z1=3.21e
15i
chuyӇn vӅ dҥng đҥi sӕ và kiӇm tra kӃt quҧ:
3.21*exp(15i)
KӁt quҥ
ans =

-2.4386 + 2.0874i
b)Z2=4+3i chuyӇn vӅ dҥng sӕ mũ và kiӇm tra kӃt quҧ:

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 7

BAI THӴC HÀNH SӔ 2.
LҰP TRÌNH TRÊN MATLAB

Bai 2.1
a)Lұp trình bài giҧi bài toán 1.3.
Xây dӵng hàm lҩy tên là tonthatcongsuat voi cu pháp như sau:
function[dentaA]=tonthatcongsuat(U,P,Q,R,to)
U=input(
'nhap gia tri U=');
P=input(
'nhap gia tri P=');
Q=input(
'nhap gia tri Q=');
R=input(
'nhap gia tri R=');
to=input(
'nhap gia tri to=');
DeltaP=((P^2+Q^2)/U^2)*R;

DeltaA=DeltaP*to;
fprintf(
'DeltaA=%7.2f kWh \n',DeltaA)
Lưu dưӟi dҥng m.file.Khi đánh lӋnh:
tonthatcongsuat(U,P,Q,R,to)
sӁ hiӋn ra cac tӯ khóa đӇ nhұp U,P,Q,R,to và hiӇn thӏ kӃt quҧ:
nhap gia tri U=35
nhap gia tri P=66.78
nhap gia tri Q=35.67
nhap gia tri R=4.75
nhap gia tri to=4620
DeltaA=102683.20 kWh
Bài 2.2 Xây dӵng hàm xác đӏnh tәng trӣ gӗm n nhánh song song và áp dөng
hàm vӯa xây dӵng đӇ tính điӋn trӣ tương đương gӗm 5 nhánh vӟi:
R1=4.25 R2=0 R3=4.25 R4=3.78.
A=input(
'nhap ma tran tong tro [R1 R2 Rn]');
tongtro=sum(A.^-1)
KӃt quҧ.
nhap ma tran tong tro [R1 R2 Rn][4.25 6.45 4.25 5.34 3.78]

tongtro =

1.0774
Bài 2.3 Hãy biӇu thӏ hàm f1(x)=x
4
+3x
3
-6x
2

+1,7x
1/2
dưӟi dҥng inline objects và
xác đӏnh giá trӏ cӫa hàm này ӭng vӟi x=0,3.
f1=inline('x^4+3*x^3-6*x^2+1.7*x^(1/2)','x')
y1=f1(0.3)
kӃt quҧ:
f1 =

Inline function:
f1(x) = x^4+3*x^3-6*x^2+1.7*x^(1/2)

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 8

y1 =

0.4802

Bài 2.4 BiӇu thӏ hàm vecto F=[f1 f3] vӟi f2=x
1
4
+x
2
3
-7;x
1
2
.x

2
-x
1
x
2
,
F3=x
2
3
-1;4x
1
2
x
2
-x
1
x
2
+lnx
1
.Vӟi x
1
=1,5;x
2
=2,4;

F=inline(
'[x1^4+x2^3-7 (x1^2)*x2-x1*x2;x2^3-1 4*x1^2*x2-
x1*x2+log(x1)]'
,'x1','x2')

f=F(1.5,2.4)
KӃt quҧ:
F =

Inline function:
F(x1,x2) = [x1^4+x2^3-7 (x1^2)*x2-x1*x2;x2^3-1 4*x1^2*x2-
x1*x2+log(x1)]


f =

11.8865 1.8000
12.8240 18.4055

Bài 2.5 Hãy sӱ dөng vòng lһp while cho bài toán tìm giá trӏ lӟn nhҩt cӫa sӕ n
mà có tәng S=1
3
+2
3
+«+n
3
nhӓ hơn 250.
S=1;n=1;
while S+(n+1)^3<250
n=n+1;S=S+n^3;
end
[n,S]
KӃt quҧ.
ans =

5 225







TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 9

BAI THӴC HÀNH SӔ 3.
MA TRҰN

Bài 3.1 ThiӃt lұp ma trұn A vӟi các phҫn tӱ bҩt kì.Sӕ dòng bҵng sӕ kí tӵ cӫa
hӑ=3.Sӕ cӝt bҵng sӕ kí tӵ cӫa tên=4.
-B có sӕ cӝt gҩp đôi A.
-C gӗm các phҫn tӱ dòng thӭ nhҩt đӃn dòng thӭ 2,cӝt thӭ nhҩt tӟi cӝt thӭ 3 A.
-D ngүu nhiên cùng kích thưӟc vӟi A.
-gӑi các phҫn tӱ cӫa ma trұn gӕc
Dòng 1 cӝt 2,cӝt 3,dòng 1.
A=[1 2 3 4;2 7 1 0;0 3 0 9;]
B=repmat(A,1,2)
C=[A(1:2,1) A(1:2,2) A(1:2,3)]
D=[randperm(5);randperm(5);randperm(5);randperm(5)]
S=A(1,2)
Y=A(:,3)
X=A(1,:)
KӃt quҧ

A =
1 2 3 4
2 7 1 0
0 3 0 9
B =

1 2 3 4 1 2 3 4
2 7 1 0 2 7 1 0
0 3 0 9 0 3 0 9
C =

1 2 3
2 7 1
D =

2 4 3 5 1
3 5 2 1 4
5 1 4 2 3
1 4 5 3 2

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 10

S =
2
Y =
3
1
0

X =

1 2 3 4

Bài 3.2 Thӵc hiӋn phép tính:
I=A+D

I =

4 4 7 5 10
4 7 4 6 9
2 7 5 7 9
2 4 8 7 9

J=A*D'

J =

48 44 49 52
48 44 49 52
48 44 49 52
48 44 49 52

%D'la ma tran chuyen vi vi de nhan 2 ma tran thi so hang ma tran A = so cot ma tran D
K=A\D
Warning: Rank deficient, rank = 1, tol = 1.1102e-014.

K =

0 0 0 0 0

0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0.4000 0.7000 0.6000 0.4500 0.8500




L=A/D

L =

-0.0000 0.0000 0.3333 0.6667
-0.0000 0.0000 0.3333 0.6667
-0.0000 0.0000 0.3333 0.6667
-0.0000 0.0000 0.3333 0.6667
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 11

BAI THӴC HÀNH SӔ 4.
XӰ LÝ SӔ LIӊU


Bài 4.1 Ma trұn A như bài 3.1.Xác đӏnh:
-Giá trӏ cӵc đҥi cӫa mӛi cӝt(C1);
-Cӵc tiӇu cӫa mӛi cӝt(C2);
-Cӵc đҥi cӫa mӛi dòng(C3);
-Cӵc tiӇu cӫa mӛi dòng(C4);
-Cӵc đҥi cӫa mӛi cӝt và chӍ sӕ cӫa chúng(C1,I)

A
C1=max(A)
C2=min(A)
C3=max(A,[],2)
C4=min(A,[],2)
[C1,I]=max(A)
KӃt quҧ:
A =

1 2 3 4
2 7 1 0
0 3 0 9


C1 =

2 7 3 9


C2 =

0 2 0 0


C3 =

4
7
9

C4 =

1
0
0



C1 =
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 12


2 7 3 9


I =

2 2 1 3
Bài 4.2 Cho bang sӕ liӋu:

1
2
3
4
5
6
7

8
9
10
Y
74,3
99,05
136,2
191,9
235,2
278,6
315,7
341,8
421,8
478,6
X
1
3,13
7,84
13,06
15,78
18,16
23,34
25,92
27,63
33,53
40,21
X
2
8,33
10,75

12,86
23,85
28,34
33,23
35,67
39,44
44,34
49,45
1.1 Tính giá trӏ trung bình cӫa các đҥi lưӧng x
1
,x
2
,y.
1.2 Phương sai cӫa các đҥi lưӧng Cx
1
,Cx
2
,Cy.
1.3 Đӝ lӋch chuҭn cӫa các đҥi lưӧng.
1.4 Đӝ lӋch trung bình bình phương cӫa các đҥi lưӧng.
1.5 HӋ sӕ tương quan cӫa các đҥi lưӧng quan sát giӳa các đҥi lưӧng Rx1_y và
Rx2_y.
2.Xây dӵng hàm hӗi quy thӵc nghiӋm và đánh giá sai sӕ vӟi các hàm:
2.1 Bұc hai y=f(x1);y=f(x2);
2.2 TuyӃn tính: y=f(x1);y=f(x2);
2.3 Hàm hӗi quy bӝi y=f(x1,x2);
1.
y=[74.3 99.05 136.2 191.9 235.2 278.6 315.7 341.8 421.8 478.6];
x1=[3.13 7.84 13.06 15.78 18.16 23.34 25.92 27.63 33.53 40.21];
x2=[8.33 10.75 12.86 23.85 28.34 33.23 35.67 39.44 44.34 49.45];

%a gia tri trung binh
Xtb1=mean(x1)
Xtb2=mean(x2)
Y_tb=mean(y)
%b phuong sai
Cx1=cov(x1)
Cx2=cov(x2)
Cy=cov(y)
%c Do lech chuan
std(x1,0)
std(x2,0)
std(y,0)
%d Do lech trung binh binh phuong
std(x1,1)
std(x2,1)
std(y,1)
% He so tuong quan
Rx1_y=corrcoef(x1,y)
Rx2_y=corrcoef(x2,y)



TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 13

KӃt quҧ
Xtb1 =

20.8600

Xtb2 =

28.6260
Y_tb =

257.3150
Cx1 =

132.2609
Cx2 =

207.9182

Cy =

1.8333e+004
ans =

11.5005
ans =

14.4194
ans =

135.4012
ans =

10.9103
ans =

13.6794
ans =

128.4528
Rx1_y =
1.0000 0.9938
0.9938 1.0000
Rx2_y =

1.0000 0.9895
0.9895 1.0000
2.
y=[74.3 99.05 136.2 191.9 235.2 278.6 315.7 341.8 421.8 478.6];
x1=[3.13 7.84 13.06 15.78 18.16 23.34 25.92 27.63 33.53 40.21];
x2=[8.33 10.75 12.86 23.85 28.34 33.23 35.67 39.44 44.34 49.45];
%a xay dung ham hoi quy bac 2 va danh gia sai so
c1=polyfit(x1,y,2)
r1=corrcoef(x1,y)
c2=polyfit(x2,y,2)
r2=corrcoef(x2,y)
%b tuyen tinh
c1=polyfit(x1,y,1)
r1=corrcoef(x1,y)
c2=polyfit(x2,y,1)
r2=corrcoef(x2,y)
%c Hoi quy boi
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 14

y=[74.3 99.05 136.2 191.9 235.2 278.6 315.7 341.8 421.8 478.6]';
x1=[3.13 7.84 13.06 15.78 18.16 23.34 25.92 27.63 33.53 40.21]';
x2=[8.33 10.75 12.86 23.85 28.34 33.23 35.67 39.44 44.34 49.45]';
XX=[ones(size(x1)) x1 x2];
a=XX\y
Y=XX*a;
ssb=(Y-y).^2;
er=sqrt(sum(ssb)./size(x1))*100/mean(Y);
ert=er(1)
KӃt quҧ.
c1 =

0.0358 10.1681 25.3768


r1 =

1.0000 0.9938
0.9938 1.0000


c2 =

0.1010 3.6945 49.8446


r2 =

1.0000 0.9895
0.9895 1.0000

c1 =

11.7010 13.2323


r1 =

1.0000 0.9938
0.9938 1.0000


c2 =

9.2918 -8.6735


r2 =

1.0000 0.9895
0.9895 1.0000


a =

2.2332
7.1792
3.6793

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 15

ert =

3.9513
Bài 4.3 Hãy xác đinh các giá trӏ nӝi suy cӫa ma trұn y ӭng vӟi các giá tri cӫa x
Y=[x

x
1,6
x
2,8
];x=3;4,5 và 8,5;
x=[3:9]';
y=[x.^1,x.^1.6,x.^2.8];
xi=[3,4.5,8.5];
yi=interp1(x,y,xi)
KӃt quҧ.
yi =

3.0000 5.7995 21.6740
4.5000 11.1611 69.5502
8.5000 30.7462 403.7786
Bài 4.4 Hãy xây dӵng hàm y trong khoҧn biӃn thiên x
Y=2,1 sinx;x=-5÷5 và các hàm nӝi suy:bұc thang,tuyӃn tính hóa,lұp phương
mô phӓng và hàm lұp phương và đánh giá sai sӕ so vӟi hàm gӕc.

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 16

Chương 5
ĐӖ HӐA


Bai 5.1 VӁ đӗ thӏ hàm sӕ y trong khoҧng biӇn thiên x=-4÷4,5.Khoҧng biӃn
thiên h=0,2.Đӭt/đen X;
࢟ൌࢋ
ି࢞
૛
ሺ૚൅૜࢞െ࢞
૛
ሻ
a)VӁ bҵng lӋnh plot.
x=-4:0.2:4.5;
y=exp(-x.^2).*(1+x*3-x.^2);
plot(x,y,'k ',x,y,'x')

b)VӁ bҵng fplot
fplot('exp(-x.^2).*(1+x*3-x.^2)',[-4,4.5],'k ','x'),grid
c)VӁ bҵng lӋnh ezplot
ezplot(
'exp(-x.^2).*(1+x*3-x.^2)',-4,4.5),grid
d)VӁ đӗ thӏ trên vӟi đưӡng liӅn nét đӓ,điӇm dҩu hoa thӏ màu đen.
x=-4:0.2:4.5;
y=exp(-x.^2).*(1+x*3-x.^2);
plot(x,y,'r-',x,y,'*k')
Bài 5.2 VӁ đӗ thӏ cӫa hàm y=5.e
x
Trong khoҧng 0ൊ ૜૙.

a)Các trөc toҥ đӝ chia theo tӍ lӋ log cơ sӕ 10.
-4 -3 -2 -1 0 1 2 3 4 5
-1.5
-1
-0.5
0
0.5
1
1.5
2
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 17

x=0:30;
y=5*exp(x);
loglog(x,y),gridss
xlabel(
'x,log10');ylabel('Y,log10')
title(
'Do thi y=f(x) lenh loglog(x,y)')
b) Vӟi trөc y chia theo tӍ lӋ log cơ sӕ 10.
x=0:30;
y=5*exp(x);
semilogx(x,y),grid
xlabel(
'x,log10');ylabel('Y')
title(
'Do thi y=f(x) lenh semilogx')
c)Vӟi trөc x chia theo tӍ lӋ log cơ sӕ 10.

x=0:30;
y=5*exp(x);
semilogy(x,y),grid
xlabel(
'x');ylabel('Y,log10')
title(
'Do thi y=f(x) lenh semilogy')
Bài 5.3 Xây dӵng đӗ thӏ hàm d trong khoҧng biӃn thiên cӫa ij trong hӋ tӑa đӝ
cӵc.D=e
cosij
+2cos4ij Vӟi ij=0÷6ʌ.
phi=0:6*pi;
d=exp(cos(phi))+sin(phi/6);
polar(phi,d)
title(
'Do thi ham d=cos(5*phi) trong toa do cuc')




Chương 6
GIҦI PHƯƠNG TRÌNH VÀ TÍNH TOÁN TӔI ƯU


Bài 6.1 Tìm nghiӋm cӫa phương trình y=f(x) vӟi :
y=x
5
+0,7x
4
-1,5x

2
+0,4x-5=0

p=[1 0.7 0 -1.5 0.4 -5];
x=roots(p)
KӃt quҧ:
x =

-1.2652 + 0.9939i
-1.2652 - 0.9939i
1.3688
0.2307 + 1.1653i
0.2307 - 1.1653i
Bài 6.3 Tìm nghiӋm cӫa hӋ phương trình vӟi:
a=[2 6 10 3;11 0.7 6.6 0.8;5 -3.2 5 4;0 4.5 6 -4.2];
b=[6 -5 6 5];
x=a\b'
KӃt quҧ.
x =
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 18


-1.4457
-1.4838
1.8294
-0.1668
Bài 6.4 Giҧi phương trình phi tuyӃn sau:F(x)=x
5

+x
3
+2e
3x
+2;
Tҥo ham ftu2.m.
function F=ftu2(x);
F=x^5+x^3+2*exp(3*x)+2
Thӵc hiӋn lӋnh.
options=optimset(
'fzero');
x0=-1;
x=fzero(
'ftu2',x0,options)
KӃt quҧ
x =

-1.0118
Bài 6.5 Giҧi hӋ phương trình phi tuyӃn sau:
2x
1
+sinx
2
=e
-x1

2x
2
=e
-x2

-cosx1
Tҥo hàm funsolve2 lưu dưӟi dҥng m.file:
function F=funsolve2(x)
F=[2*x(1)+sin(x(2))-exp(-x(1));2*x(2)-exp(-x(2))-cos(x(1))]
Thӵc hiӋn lӋnh:
x0=[1;1];
options=optimset(
'Display','off');
x=fsolve(
'funsolve2',x0,options)
KӃt quҧ.
x =

0.1115
0.73

Bài 6.6 Tìm nghiӋm cӫa phương trình phi tuyӃn.
F(x)=x
6
+x
4
-x
2
+1,5e
2x
-5,4
Thӵc hiӋn câu lӋnh:
x=fzero(
'x^6+x^4-x^2+1.5*exp(2*x)-5.4',1.0)

KӃt quҧ.
x =

0.6556
Bài 6.7 Tìm cӵc trӏ cӫa hàm 2 biӃn:F(x)=e
x1
(4x
1
2
+2x
2
2
+4x
1
x
2
+2x
2
+1);
Vӟi các rҵng buӝc:
x
1
x
2
-x
1
-x
2
-1,5;
x

1
x
2
10;
Tҥo hàm confun2b
function [c,ceq]=confun2b(x)
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 19

c=[x(1)*x(2)-x(1)-x(2)+1.5;-10-x(1)*x(2)];
ceq=[];
Tҥo hàm ftoiu2b
function f=ftoiu2b(x)
f=exp(x(1))*(4*x(1).^2+2*x(2)+4*x(1)*x(2)+2*x(1)*x(2)+1);
Thӵc hiӋn câu lӋnh:
x0=[-0;2];
options=optimset(
'LargeScale','off');
[x,fval]=fmincon(
'ftoiu2b',x0,[],[],[],[],[],[],'confun2b',options)
KӃt quҧ.
x =
-0.7460
13.4041
fval =
-14.2109
Bài 6.8 Giҧi hӋ phương trình phi tuyӃn sau:
ࡲ
ሺ࢞ሻ

ൌቈ
࢞
૚
૝
൅૛ǡ૝࢞
૛
૛
െࢉ࢕࢙
ሺ
࢞
૛
ሻ
ࢋ
ି૚ǡ૞࢞
૚
൅૛ǡૠ࢞
૛
െ૛ǡ૛࢙࢏࢔ሺ࢞
૚
ሻ
቉ൌ૙
Tҥo hàm f2(x).
function F=f2(x)
F=[x(1)^4+2.4x(2)^2-cos(x(2));exp(-1.5*x(1))+2.7*x(2)-2.2*sin(x(1))];
Thӵc hiӋn câu lӋnh.
x0=[1 1];
options=optimset(
'fsolve');
options=optimset(options,
'Display','iter');

[x,y,h]=fsolve(
'f2',x0,options)
KӃt quҧ.
x =
0.7859 0.4625
y =
1.0e-012 *

0.1865
0.0260
h =
1


TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 20

BAI THӴC HÀNH SӔ 9.
GIҦI MҤCH ĐIӊN


Bài 9.1Cho mҥch điӋn 1 chiӅu vӟi:

E
1

E
2

E
3

E
4

R
1

R
2

R
3

R
4

R
5
120
170
0
0
12,2
15,6
9,4
13,2
2,8
Thӵc hiӋn lӋnh.

clear
E=[120 117 0 0 0];
R=[12.2 15.6 9.4 13.2 2,8];
disp(
'a)phuong phap dong dien nhanh');
Z=[1 1 1 1 1;R(1) -R(2) 0 0 0;0 R(2) -R(3) 0 0;0 0 R(3) -R(4) 0;0 0 0 R(4) -
R(5)];
U=[0;E(1)-E(2);E(2)-E(3);E(3)-E(4);E(4)-E(5)];
I=Z\U;
disp(
'I1 I2 I3 I4 I5');
disp([I(1) I(2) I(3) I(4) I(5)]);
KӃt quҧ.
a)phuong phap dong dien nhanh
I1 I2 I3 I4 I5
8.1203 6.1582 -2.2268 -1.5858 -10.4660
Bài 9.2 Giҧi bài toán 9.1 vӟi các than sô phӭc.
E
1

E
2

E
3

E
4

R

1

R
2

R
3

R
4

R
5
145+j88
170+j125
0
150+j110
9+j6,5
5,2+j7
6,2+j5,4
5,4+j7
5,3+j7

Thӵc hiӋn lӋnh.
clear
E=[145+88*i 170+125*i 0 150+110*i 0];
R=[9+6.5*i 5.2+7.0*i 6.2+5.4*i 5.4+7.0*i 5.3+7.0*i];
disp(
'a)phuong phap dong dien nhanh');
Z=[1 1 1 1 1;R(1) -R(2) 0 0 0;0 R(2) -R(3) 0 0;0 0 R(3) -R(4) 0;0 0 0 R(4) -

R(5)];
U=[0;E(1)-E(2);E(2)-E(3);E(3)-E(4);E(4)-E(5)];
I=Z\U;
disp(
'I1 I2 I3 I4 I5');
disp([I(1) I(2) I(3) I(4) I(5)]);
KӃt quҧ.
a)phuong phap dong dien nhanh
I1 I2 I3 I4 I5
Columns 1 through 4

5.3516 - 0.8917i 11.2693 - 2.9086i -13.2417 + 1.6557i 8.4405 - 1.9117i

Column 5

-11.8198 + 4.0563i


TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 21

Bài 9.3 Hãy xác đӏnh dòng điӋn và công suҩt cӫa mҥch điӋn.BiӃt:

U
m
,V

Z
m

127
30
3,49
88
Thӵc hiӋn lӋnh.
Um=127;tetau=30;
Z=3.49;gam=88;
tetai=tetau-gam;
teta=(tetau-tetai)*pi/180;
Im=Um/Z;
wt=0:.05:2*pi;
u=Um*cos(wt+tetau);
i=Im*cos(wt+tetai*pi/180);
p=u.*i;
U=Um/sqrt(2)
I=Im/sqrt(2)
P=U*I*cos(teta)
Q=U*I*sin(teta)
S=P+Q*i
pr=P*(1+cos(2*(wt+tetau)));
px=Q*sin(2*(wt+tetau));
pp=P*ones(1,length(wt));
xline=zeros(1,length(wt));
wt=180/pi*wt;
subplot(2,2,1),plot(wt,u,wt,i,wt,xline),grid
title(
'Dien ap va dong dien')
xlabel(
'wt,do'),ylabel('U,I')

subplot(2,2,2),plot(wt,pr),grid
title('Cong suat toan phan')
xlabel(
'wt,do'),ylabel('S,VA')
subplot(2,2,3),plot(wt,pr),grid
title(
'Cong suat ta dung')
xlabel(
'wt,do'),ylablel('p,W'),grid
subplot(
'Cong suat bieu kien')
xlabel(
'wt,do'),ylabel('q,VAr'),grid

KӃt quҧ.
U =

89.8026
I =

25.7314
P =

80.6438
Q =

2.3093e+003
S =

1.0e+004 *

TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 22



Bài 9.4 Nguӗn điӋn xoay chiӅu có điӋn áp u(t)=U
m
sin(wt) tҫn sӕ 50Hz.Đưӧc
nӕi vӟi điӋn trӣ R và cuӝn dây L,điӋn trӣ shun R
sh
.Xác đӏnh các giá trӏ hiӋu
dөng cӫa điӋn áp U
RL
,các thành phҫn U
R
và U
L
.VӁ đӗ thӏ cӫa các đҥi lưӧng
này.
R,
L, 10
-3
H
R
sh
,
I
sh

,A
4,55
30
0,25
1,5

Thӵc hiӋn lӋnh
disp(
'xac dinh gia tri hieu dung cua dien ap URL. Cac thanh phan Ur va UL');
R=4.55;
L=30*10^-3;
Rsh=0.25;
Ish=1.5;
f=50;
t=0:0.001:0.04;
omega=2*pi*f;
X=omega*L;
Ushr=Ish*R;
Ushx=Ish*X;
UR=Ushr/sqrt(2);
UL=Ushx/sqrt(2);
0 100 200 300 400
-200
-100
0
100
200
Dien ap va dong dien
wt,do
U,I

0 100 200 300 400
0
50
100
150
200
Cong s uat toan phan
wt,do
S ,V A
0 100 200 300 400
0
50
100
150
200
Cong s uat ta dung
wt,do
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 23

URL=UR*sqrt(2);
UshRL=Ushr*sqrt(2);
uR=Ushr*sin(omega*t);
uL=Ushx*sin(omega*t+pi/2);
phi=atan(Ushx/Ushr);
uRL=UshRL*sin(omega*t+phi);
disp(
' UR,V UL,V URL,V')
fprintf(

'%g'),disp([UR UL URL])
plot(t,uR,
'r-',t,uL,'g ',t,uRL,'kp-'),grid;
xlabel(
't');ylabel('U,V');
title(
'Do thi dien ap uR, uL, uRL: Sinh vien ± Bui Bao Linh± Lop D2 ± H3 ')

KӃt quҧ.
xac dinh gia tri hieu dung cua dien ap URL. Cac thanh phan Ur va UL
UR,V UL,V URL,V
4.8260 9.9965 6.8250

Bài 9.5 Cho mҥch điӋn hình sin tҫn sӕ 50Hz:u=125sin(100ʌt).Hãy xác đӏnh các
giá trӏ cӫa điӋn áp và dòng điӋn chҥy trong mҥch,biӃt U
m
;
U
m
=125 V;R=33;C=220 µF

0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
-15
-10
-5
0
5
10
15
U

Do h d en ap uR uL uRL: nh v en ± Bu B ao L nh± Lop D2 ± H3
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 24


Thӵc hiӋn câu lӋnh.

disp(
'Gia tri hieu dung dien ap va dong dien chay trong mach');
Um=125;
R=33;
C=220*10^-6;
f=50;
t=0:0.001:0.04;
omega=2*pi*f;
X=1/(omega*C);
Z=sqrt(R^2+X^2);
Im=Um/Z;
UmR=Im*R;
UmC=Im*X;
phi=atan(UmC/UmR);
i=Im*sin(omega*t+phi);
ur=UmR*sin(omega*t+phi);
uc=UmC*sin(omega*t+phi-pi/2);
I=Im/sqrt(2);
Ur=UmR/sqrt(2);
Uc=UmC/sqrt(2);
disp(
' I,A Ur,V Uc,V')

fprintf(
'%g'),disp([I Ur Uc])
plot(t,ur,
'r-',t,uc,'gp', t,uc,'g-'),grid;
xlabel(
't');ylabel('ur,uc,V');
title(
'Do thi dien ap: Sinh vien ± Bui Bao Linh ± Lop D2 ± H3 ')

KӃt quҧ.
Gia tri hieu dung dien ap va dong dien chay trong mach
I,A Ur,V Uc,V
2.4530 80.9496 35.4918

R
u
u
C
u
R
i
C
TRƯӠNG ĐҤI HӐC ĐIӊN LӴC

SV:Bùi Bҧo Linh 25



Bài tұp 9.6 Hãy xác đӏnh giá trӏ dòng điӋn và công suҩt tác dөng

Chương 10 Tính toán phө tҧi điӋn
Bài tұp 10.1

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à:
Pdl=[25,16,10,7.5,5.6,3];
cosphi=[0.74,0.78,0.72,0.67,0.8,0.78];
Hãy xác đӏnh phө tҧi tính toán theo hai phương pháp:
a) Phương pháp sӕ gia;
b) Phương pháp hӋ sӕ nhu cҫu;
* Phương pháp s͙ gia:
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
-150
-100
-50
0
50
100
150
ur uc
Do h d en ap: nh v en ± B u B ao L nh ± Lop D2 ± H3