AEA BO CONG THUONG TP.HCM KHOA CÔNG NGHỆ ĐIỆN Wu BAO CAO CUOI KY MON: MATLAB VA UNG DUNG TRONG KY THUAT DIEN GVHD: Nguyễn Ngọc Anh Tp.
Trang 1
AEA
BO CONG THUONG
TP.HCM KHOA CÔNG NGHỆ ĐIỆN
Wu
BAO CAO CUOI KY
MON: MATLAB VA UNG DUNG TRONG
KY THUAT DIEN
GVHD: Nguyễn Ngọc Anh
Tp Hồ Chí Minh , Ngày tháng năm 2023
Trang 2
BÀI TẬP CHƯƠNG 2
Bài tập 2.1: Xác định điện áp, dòng điện, công suất theo giá trị thời
gian:
2
4 if te2
6 v = 2*t-5
10 v = (t82/4)4t-4 4
12 elseif te=14
16 else
18 v = 67-4*t
19 p= vii pe
21 disp([‘Gia tri dong dien: *,numastr(i)]) 20
3 disp(['Gia tri cong suat: ',numstr(p)]) Gia tri dong dien: 4
Bai 2.2 Giai phwong trNnh phi tuyén bang phương pháp lặp Newton — Rapshon:
Trang 3
NNNNNNNNRP
28
29
30
31
32
33
35
36
37
38
39
41
42
43
45
47
clear all;clc;
n=input(‘Nhap gia tri x01 :');
m=input( 'Nhap gia tri x0@2:');
k=input('Nhap gia tri x@3 :');
x=n;
y=m;
z=k;
f=69;
g=69;
h=69;
s111=0;
s112=0;
s113=0;
x1=0;
y1=0;
z1=0;
while f>1e-15
x1=x;
f1=x^3-x-1;
đf1=3*x^2-1;
x=x-f1/df1;
f1=x^3-x-1;
sS111=s111+1;
=abs(x1-x);
fprintf('%3i %12.8f %12.8f %12.8f %20.15f\n"',s111,x,f1,df1,£)|
end
while g>1e-15
y1=y;
f2=(exp(-y))-Y;
df2=(-exp(-y))-1;
y=y-f2/df2;
f2=(exp(-y))-Y;
s112=s112+1;
g=abs(y1-y);
fprintf('43i %X12.8f %12.8f %12.8f %2Ø.15f\n',s112,y,f2,df2,g)
end
while h>1e-15
z1=z;
f3=z+sin(z)-(p1/2);
df3=1+cos(z);
z=z-f3/df3;
f3=z+sin(z)-(p1/2);
s113=s113+1;
h=abs(z1-z);
fprintf('%3i %12.8f %12.8f %12.8f %2Ø.15f\n',s113,z,f3,df3,h) end
Trang 4Nhap gia tri x01 :1.5
Nhap gia tri x02:0
Nhap gia tri x03 :0
Bài tập 2.3 Xác định số tiền điện theo điện năng tiêu thụ:
Bài tập 2.4 Tính giá trị điện trở của dây dẫn theo loại dây dẫn được nhập vào:
Trang 5(\Ð
1-16; Nhap chat lieu:cu
n=input('Nhap chat lieu:','s');
if n=='Cu!
R=rCu*1/F;
again='no';
break
elseif n=='Al’
R=rAl*1/F;
again='no';
break
else
disp('Xin vui long nhap lai chat lieu')
again='yes';
end
Bài tap 2.5 Khao sat tong tré mạch RLC nối tiếp:
clc;clear all; & Command Window
R = 10;
L = 10e-3; Zmin=10
fmin = 20; fx >>
0 = (1/(2*pi*sqrt(L*C)));
for f=1:1000
w = 2*pi*f;
Z=sart(R^2+(w*L-1/(w%*C))^2);
if Z<Zmin
Zmin=Z;
fmin=f;
end
end
Trang 6Bài tập 2.6 Khảo sát đáp ứng mạch RC:
PRPRR 2
14
clc;clear all; @& Command Window
R=1e3;
C=10e-6; _nn
ucb=E;
A=utd/ucb;
while A>0.05
t1=t1+0.001;
A=utd/ucb;
end
uc=E-utd;
disp([‘uc=', num2str(uc)]);
Bai tap 2.7 Đọc giá trị điện trở theo vạch màu:
clc;clear all;
Den =
Do =
Cam =
Vang
Lục
Tim
Xam
Trang
kta =
while kta==
a=input('Nhap gia tri wach mau thu 1:');
\O “`
kta=a==Den | |a==Nau | |a==Do | |a==Cam | |a==Vang | |a==Luc | |a==Lam| | a==Tim| |a==Xam | |a==Trang;
end
Den = 0;
Do = 2;
Vang = 4;
Lam = 6;
Tim = 7;
Xam = 8;
Trang 7
32
33
35
36
37
38
39
41
42
43
45
47
49
59
51
52
53
55
56
57
58
59
61
62
63
65
67
69
79
71
72
73
74
75
76
77
78
while ktb==
b=input('Nhap gia tri vach mau thu 2:');
ktb=b==Den| | b==Nau| | b==Do| | b==Cam| | b==Vang| | b==Luc| |b==Lam] |
end
Den = 0;
Nau = 1;
Do = 2;
Cam = 3;
Vang = 4;
Lục = 5;
Lam = 6;
Tim = 7;
Xam = 8;
Trang = 9;
Hoangkim = -1;
Bac = -2;
ktc=0;
while ktc==0
c=input('Nhap gia tri vach mau thu 3:');
ktc=c==Den | |c==Nau | |c==Do | |c==Cam | |c==Vang | |c==Luc | |c==Lam| |
end
Den = 0;
Nau = 1;
Cam = 0;
Vang = 0;
Trang = 0;
Do = 2;
Luc = 0.5;
Lam = 0.25;
Tim = 0.1;
Xam = 0.05;
Hoangkim = 5;
Bac = 10;
ktd=0;
while ktd==0
d=input('Nhap gia tri vach mau thu 4:');
d==Tim| |d==Xam | |d==Den | |d==Hoangkim | |d==Trang | |d==Bac ;
end
R=(a*10+b)*10*c;
disp(['Gia tri dien tro theo don vi Ohm: ',num2str(R)])
Trang 8
Nhap gia tri vach mau thu 1:Vang
Nhap gia tri vach mau thu 2:Xam
Nhap gia tri vach mau thu 3:Do
Nhap gia tri vach mau thu 4:Hoangkim
Gia tri dien tro theo don vi Ohm:4800
bt2_7
Bài 2.8 Tính đô hụt áp trong mạch (Tính gần đúng):
klear all; clc;
I=1;
L = input('nhap chieu dai day dan L:'); s=input('Nhap tiet dien S:');
if s<10
X=0.15;
else
X=0.1;
end
again='yes';
while again=='yes'
VL = input('Nhap chat lieu:','s');
if VL=='Cu'
R=22.5*L/s;
break
elseif VL=='Al'
R=36*L/s;
again='no';
break
else
disp('Nhap lai chat lieu’)
again='yes';
end
end
disp(' 1 = Chieu sang,den soi dot’)
disp(' 2 = Chieu sang,den huynh quang')
disp(' 3 = Dong co,khi khoi dong")
disp(' 4 = Dong co,khi van hanh binh thuong') load=input('Nhap loai tai:');
Trang 935
36
37
38
41
42
43
45
47
49
50
51
52
53
54
55
56
57
58
61
case 1
pf=1;
case 2
pf=0.85;
case 3
pf=0.35;
case 4
pf=0.8;
end
n=input(‘Nhap so pha:");
if n==
deltaU = 2*1*(R*pf+X*sqrt(1-pf*2));
again2='no';
break
elseif n==3
deltaU = sqrt(3)*1*(R*pf+X*sqrt(1-pf*2));
again2='no';
break
else
disp('Nhap lai so pha‘)
again2='yes';
end
end
<< Command Window
>> bt2 8
nhap
Nhap
Nhap
ua
Nha
Nhap
chieu dai day dan L:1
tiet dien S:€
Dong co,khi wan hanh binh thuong
Se >>
Trang 10CHƯƠNG 4 VẼ ĐỎ THỊ
Bài tập 4.1: Vẽ đường đặc tuyến của Diode:
6 Vt=(k*T)/q
Đồ thị:
Dac tuyen Diode
0.35 T T T T T T T T T
8
c
„Ð
c
S
a
0 1 L L L 1 L L
Dien ap Diode V
Bài tập 4.2: Khảo sát lực tĩnh điện:
Trang 11POUWUANDUAWN
clear all; clc;
k = 9*10^9;
qi=1.6*10*-19;
q2=q1;
r = 1.0*10*-15:1.0*10*-16:1.0*10*-14;
F = k*((q1*q2)./r.*2)
plot(r,F,'*");
legend( 'F') ;
Đồ thị:
250
200
bk
= an oO
= ° oO
50
Dac tuyen luc tinh dien
Bài tập 4.3: Khảo sát đồ thị nhiệt độ trong năm:
Trang 12thang = [1 2 3 4 5 6 7 8 9 10 11
mt = [thang;avg;min;max];
hold on
Đồ thị:
30 - Temperature in City
Months in year
0 1 L
12];
4);
3];
8];
Min T
Bài tap 4.4: Khao sat tNnh bNnh van hành hệ thống điện quốc gia:
Đồ thị:
Trang 13\O
Bài tập 4.5: Tông hợp chuỗi Fourier:
clc; clear all];
T=4;
N=500;
f=1/T;
w=2*pi*f;
F = 1/2;
k=0;
for i =-T/2:T/10000:T/2
k=k+1;
F(k) = F(k) + 2/pi/n*sin(pi*n/2)*cos(n*w*t(k)); end
end
plot(t,F)
Đồ thị trường hợp N=5
0.8
0.6
0.4
0.2
chuoi Fourier
Trang 14
Đồ thị trường hợp N=50
0.6 F
0.4 -
0.8 5
chuoi Fourier