NHAN XET CUA GIAO VIEN HUONG DAN Diém lab report:.
Trang 1VỆ INDUSTRIAL
TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP TP HÒ CHÍ MINH
Khoa: Công Nghệ Thông Tin
TY
ol OF 46
VỆ ©
Gy
LAB REPORT
Student’s ID
Student’s name
Subject
Instructor
Faculty
Completed Date
: 22001985 : Vt Phi Long : PTHTDPT
: Nguyễn Thành Thái
: Công Nghệ Thông Tin
28 /10 /2024
Trang 2
NHAN XET CUA GIAO VIEN HUONG DAN
Diém lab report:
Trang 3
+ Sine Wave (tạo, điêù chỉnh biên độ, tân sôƒ
% tạo sóng sin đơn giản
i = 0:0.2:2*pi; % Tao mot mang i từ 6 đêñ 2m với khoảng cách 9.2
y = sin(i); % Tính giá trị sine cho môi phâh tư trong mảng i
figure(1);
ptot(y);
title('Simple 1 Period Sine Wave');
% thay đôi biên độ sóng sin
amp = 2.6; % Đặt biên độ mới Là 2
y = amp * sin(1); % Tính sóng sine với biên độ 2
figure(2);
ptot(y);
title('Simple 1 Period Sine Wave Modified Amplitude');
% điêù chỉnh tâà sô “sóng sin
Fs = 11925; % Tâà sô Lâý mâũ (Hz)
nsec = 2; % Thời tượng sóng (s)
dur = nsec * Fs; % Sô mâũ câà thiêft trong thời gian nsec
n = O:dur; % Tao mang mau
y = amp * sin(2 * pi * n * Fw / Fs); % Tính sóng sine voi tan sô“Fw và biên độ amp figure(3);
plot(y(1:500)); % Hién thi 500 mai dau tién
title('N Second Duration Sine Wave');
% tạo N chu kì sóng sin
ncyc = 2; % Sô chu kỳ mong muôá
n = 6:fLoor(ncyc * Fs / Fw); % Sô mâũ đê tạo `n` chu kỳ
y = amp * sin(2 * pi * n * Fw / Fs);
figure( 7);
ptot(y);
title('N Cycle Duration Sine Wave');
Trang 4
New Open Seve fay print ~ GoTo Afind ~ Refactor 2 281 ‘Analyze Pun Fah unand Advance pun Step Stop
~~ oe x + — Bacckmark © ~ te Section Qj Run to End >
RE NenGaTE ccoE ANALYTE SECTION RUN
weer
Current Folder itor - F\2024- er thdeg da prxsong tem THUC HANH\matlab\untitled.m
Name ^ untl0edm ‘untitied2 * unDed3* X| uslled4* %] wnUWed5* | untiveat2 * Ì Name^+ Value
© & mavab 1 % tạo sóng sin đơn giải + ®1; (1-2891929484 nhe hạn toàn thêng ti 2 i = 0:0.2:2"pi; Tạo nit mang i từ Ø đến 2m với khoảng cách 6.2 >
3 y = sin(i); % Tinh giá trị sine cho mỗi phần tử trong lampimơe — +1
& Oe THUC HANH 7 to s bie in sn Wave"); Bile Edit View Insert Tools Desktop Window Hep
# %® AUDIO 8 ” ¡ biên độ aM sin
16 ye co * T sinc; xì Lên tổng sine với biên để 561256 double
ă 1 8imple 1 Period Sine Wa.£, 3 =) QQ) 100
14 title('Sixple 1 Peried Sine Meve Hedified AnpÏ ị \ oy fhx32 doutie
16 % điều chính tần số sóng s ị \ 1500
1 Fs = 11025; % Tân số 1 mẫu (Hz) 04
18 ‘Fu = 440; it số sóng (Hz) ˆ 5! double
19 nsec = 2; % Thoi_lugne sone (s) 02 256
20 aur nsec * Fs; % so nău cần thiết trong thời 2
23 y= anp * sỉn(2 * pÌ *n * Fw / Fs); % Tinh sé
3 Thống bê máy tính và ứng đụng ‘Command Window | / Jhx7500 double
Tương tác người máy
sinwaw (Têp WAV) i
itude:
Phase angle (radians): 0.7954
No details available Resultant Amplitude: 2.4183
Resultant Phase (radians): 0.63094
Resultant Phase (degrees): 46.0675
»
= 2 @& %'@#2 E Section Break
qP (HE öemev SỈ € g _ *?? na„a„ TE] a b ®
New Open Save ja print ~ GoTo Afind ~ Refactor Od 2H ie [o> Analyze | Fun Runand AđGC pun sụp Stop
‘Current Folder (Bi Eottor - F\2024- Do nh han tiện THỰC HÀN##/natlaÐ\ưrttled.en @ * Workspace 2 Name ~ u ‘untitled * ntiveds* X | untled4° €3] ưn0eds* X] unHedl2^ X| + | Name Value
I8 Ê9 matab 1 "x tạo sóng sẵn ôm giản {2 Figure 2 ¬ 5B Nhập môn an toàn thông từ 2 9:0.2:2*pi; % Tạo một nằng ¡ từ Ø đến 2m với khoả
8 3 y = sin(i); Tinh gid tri sine cho mdi phan ti {Ble Ed Wew Insert Tools Desktop Window Help »
5 Figure(1);
6 plot(y); Simple 1 Period Sine Wave Modified Amplitude ie
7 title('Sisple 1 Peried Sine Neve"); 2 at Inpiex
8 * they ot biến độ sine, se / \
9 anp % Đặt biên độ nới a
10 y= co * T sinQ); * Tính sống cine vũ “hiên độ 2 / \ buble
1
14 title("Simple 1 Period Sine Weve Modified Amplitude’)} 05 / 2
15 / ubie
16 % di8u chinh tan e Sóng s
17 Fs = 12025; số ty miu (Hz) 9
ine (s) thiết trong thời gian nệ ao f
X Tạo mảng „ `" mẫu 4 \ / amp * sin(2 * pi * n * Fw / Fs); % Tinh séne sine)
De Tượng tác người máy
- 0.7071 + 0.70724 | 5 10 15 20 2 a 35 50001
Anplitude: 1
Phase angle (degrees): 45 rotated phasor O7071 + O70Tti
No detaiis available Resultant Amplitude: 2.4183 is 11500 double
Resultant Phase (radians): 0.63094 signal 88200! double Resultant Phase (degrees): 46.0675 1 it 1x2007 double
>> T
Trang 5
©
oF Eb come + i ®) Some El a $
New Open Save (i print ~ GoTo AFind © Refactor 2) 2 lit Ly Analyze Pun HE Ranand Advance Rụn Step Stop
~~ oe x + — Bacckmark © ~ We Section Qi} Run to End x
RE NenGaTE cơø snare SECTION RUN
6» 8 aIZl = › rỉ
|Current Foider 1d“
Current Folder @ (Bh Eortor - Fx2024-2025\Pnat iến hệ thống đa phương tiệnTHỰC HÀNt#ynarla/rrttled.m _ Saue3 = n
Name ~ untiveam tưfftled2 + untitled: untitedt2 Bile Edit Vi “ wi Hel
or — , — Bile Edit View Insert Tools Desktop Window Help
® % matiat 7 title('Simle 1 Period Sine Wave"); = =
tien an tồn thơng tự 8 % thay đối biến độ sĩng sin 36 3.đE1 hfã
a " 9 anp = 2.0; Đặt biên độ nởi là 2
1Ø y = amp * sin(i); % Tinh sĩng sine với biên độ 2 2 N Second “1= Sine Wave
14 titleC “Sipe 1 Period Sine Wave Nodified Amplitude’); | | | | | | |
5 | | | HHI
7 Fs = 11025; Ä Tần số lấy mẫu (
18 M = 440; Tan 38 s6ny is) 05 Soman 20 dur = nsec * Fs; ặ _ oe <n n thiết trong thời gian nsec 0
ena so | | | Ì | | | | |
THIGK 2 23 anp * sin(2 * pi * n * Fw / Fs); % Tinh sĩng sine voi tan sf 05 \| | | | | | | |
24
26 pletG( 500); * Hiến thị 500 mẫu đầu tiên z7 itle('N Second Duration Sine Wave’); | Ì | |
28 x tye chu As ning si a5 | || | | | |
Comenand Window 50 100 we 200 250 300 xã 400 450 600
os 0.7071 + 0.7014
sinwaw (Tép WAV) Amplitude: 1 0.7500 + 1.29901
Resultant Amplitude: 2.4163 LH rotsted_phasor oven + 07071¡
No detaits available Resultant Phase (radians): 0.83894 Lis
Resultant Phase (degrees): 48.0675 Hd signal
©
qP (HE ưemev SỈ € BB B® nay | CS PGS
New Open Save ay print ~ GoTo Afind = Refactor [J 2 lit Le Analyze Fun PB Runand Adore — pun Step Stop
~ oe x ~ Bi acckmark © ~ ae Section Qj Run to End x
ne hasncare còt ANaUzE section RUN
6> mi ểl m › rỉ
Current Folder © [r6 -rv2024-2025 mát iếnhệ Đếng dạghương tem THỰC HÀygsynetstoroueoet Sreue + - o x
7 tat Simple 1 Period Sine Wave"); ¥ ï
9 amp = 3; 6; X Đặt biên độ mới là 2 1Ø y = emp * sin(i); % Tinh sĩng sine với biên độ 2 2 N Cycle Duration Sine Wave 2
11 \ \
14 title "Siepie 1 Period Sine Weve Modified Amplitude’); i
GK 16 % didu chinh tan số sĩng sỉin /
matiad 17 Fs = 11025; X Tan 36 léy miu an / \
mm "¬ 19 25 % Thai luong séng (s) } \ \
#) umic 20 nsec * Fs; X Số nầu cần thiết trong thời gian nsec 0 \ \
_ men 3 n = 9:dur, ®% Tạo mảng mẫu \ \ Ị
+ K 2 \ \ Ị
TS 23 anp * sin(2 * pi * n * Fw / Fs); % Tính sĩng sine với 05 \ \ }
' 26 plot(y(1 500)); X Hiến thị S00 mẫu đầu tiên
: 27 title("N Second Duration Sine Wave’); |
29 % tạo mi chu ki séng sin \ \
sen tem & h \
sinwav (Tep WAY) Phase angle (radians): 0.7854 He
Phase angle (degrees): 45
Resultant Amplitude: 2.4162
Resultant Phase (radians): 0.63094 br
No detaits available Resultant frase (degrees): 48.0675 Hs
+ Waveform (Square/Sawtooth)
Square Wave (Sĩng Vuơng)
Fs = 1000; % Tâà sơ Lâý mati (Hz)
Fw = 5; % Tâh sơ sĩng (Hz)
n = 0:1/Fs:1; % Thời gian từ 60 đêđ 1 giây với bước nhảy 1/Fs
% Biện độ sĩng
Trang 6
y_square = amp * square(2 * pi * Fw * n);
figure;
plot(n, y_square);
title('Square Wave');
xLabeL('Time (s)');
ylabel('Amplitude');
% Sawtooth Wave (Sóng Răng Cưa)
y_sawtooth = amp * sawtooth(2 * pi * Fw * n);
figure;
title('Sawtooth Wave');
xlabel('Time (s)');
ylabel('Amplitude');
% Hiêñ Thị Sóng Vuông và Sóng Răng Cura Cùng Lúc
figure;
hold on;
title('Square (Blue) and Sawtooth (Red) Waveforms');
xlabel('Time (s)');
ylabel('Amplitude');
hold off;
a K
New Oph Sore gee pate Bind — peace ab PO Rab mmand Ate mm sp sọ
+ xi Nhedmav v Fr Section Rj Run to End x
Current Falde @ F\2024-2025\Phát triển hệ th ng eN\THUC HANH\matlab\unttied2 n Bile Edit View Insert Tools Desktop Window Help » Nome + Untitiedm XÍ ưwHledzm untitleds * WuWeds" X| WMHEr) c3 gị là @ | [E3 bíE
2 % Square Wave (Sỏng Vì
Fs = 1000; % Ta miu (Hi Square Wave
® " 1; X Th ừ 6 đến 1 giây với bước nhảy 1/
* unitledrm th Wave (Séng Ran « -02+
# untiuled2m w Sống Răng Cus)
¥ untitled th = emp * sawtooth(2 * pi * Fw * np mat
a ú ("Sawtooth Mave")
B va Cùne Lúc 0 01 02 03 04 OS 06 07 08 09 1
= 078%
phasor 1
LE resstant_amp 24183
bi reostant_phase 0.8289
Hi resutant_ pha 1.6160 + 1.7990
Ei rotated phasor 0.71
Zoom: 100% | UTF-B CƯ [script na Col 24
Trang 7
New Đen Sen i Print GoTo Afind ~ Refactor 2 281 Analyze fun HÀRaEMAđAPE pm sụp Stop
Ñbcemav v ry Section Qj Run to End
'* 8 aiðl = › r>
Current Folder đeo F2024-2025\Phát triển hệ thống đa phương tiệt \THỰC a led2.m # Figure 2 = D x Name ^ uNWf@đm Xj ưwiledzm uHÚ@đ3* X | untitied4 * Untiteds =H] untitledt2 File Est View Insert Tools Desktop Window Help »
© & mavab % Saượcc Mave (Sông Vuông)
Fu = 5; số sóng (Hz)
n~6:1/Fs:i; X Thời gian | từ 6 đến 1 giây với bước nhảy 1/Fs 1 Sawtooth Wave amp = 1.0; % Biên độ sống
De THUC HANH Lý tuy yusquare = amp * square(2 * pi * Fw * n);
BS AUDIO
® % Kick ¬ title(’square Wave"); Z
thew sinwav ylabel("Arplitude");
5
Sioned % Sawtooth Wave (Séng Răng Cưa)
umnledem y_sawtooth = amp " sautooth(2 * pí * Fw * n);
@ Se THIGK figu
yo qua plot(n, y_savtooth
‘ ¬- title: Seatooth lave" ;
Teens z1 xÌabel("Tine (s)*);
8 35 TUAN4 2 yÌabel( "Asplitude *);
D Oh TUANS bs
8Ù 33 TUANG 24 # Hiến Thì Sóng Vuông và Són” Răng ưa Cùnø Lúc i
* ee Command Window
+ 19 he may tinh v8 ứng đựng apeaeiee x
Ấ ®S Thư mục més Phase angle (radians): 0.7954
nu Phase angle (degrees): 45 jphasoet “08660 + 050001
#nwar (êp WAV) Resultant Amplitude: 2.4163 Resultant Phase (radians): 0.63094 ae —
Resultant Phase (degrees): 4€.0675 itent_ pha 1.6160 — + 1.7990
>> untitled H eet | phazoc san +07071i
= 2 ( % @N Section Break
ä (Ữ@mpaev WI © g bì ` ©) profter & iS b ®
Sove fab Print ~ «GoTo AFind © Refactor DJ 9Ï lie [Aye un Et Runand Adare — un sụp Stop
v x Ñ8osemank v ~ ae Section ’
ne NawiGare ccoE anauvte SECTION RUN
jiếi =» +
Current Folder 2 E@tor - FA2024-2025\Phát triển hệ thống đa phương tiệt) THUC noe led2m J4 ngưe 3 - D x Name ~ J untitecin Šj ưnmedzm >[ HNUf@đ3' Ý] HWNed4* Se] Unttiess* SY Untitled Fie cor view Insert Too's Desktop window Help + f8 3$ matlab ĩ ¬ Wave (Song Vuong)
® Nhập món an toàn thông tự 1000; % Tan sd ly mẫu (Hz) 1ø dè 3 1681| h5
i n= O:1/Fs:1; mn ti 0 đến 2 giây với bước nhảy 1/E amp = 1.0; % Biên độ sống
y_zquare = amp * square(2 * pi * Fw * n);
Figure;
plot(n, y_square);
title(‘ Square Wave");
xlabel ("Tine (s)");
ylabel( "Amplitude" );
% Sawtooth Mave (Sing Răng Cưa) y_sawtooth = emp * sawtooth(2 * pi * Fw * n);
đang;
plot
tu sawtooth | me 3;
xabel(*Tine (s)*);
ylabel(‘Arp]itude");
2 Tal TUANG 24 #% Hiến Thì Sóng Vuông và Sónr Răng Cưa Cùnz Lúc Ẫ
5 vom 9 kể máy tính vô ứng đựng ‘Command Window U22 GUJAE tiauipiia,- 021024 | B5 Thư mục mã Phase angle (degrees): 45
Resultant Amplitude: 2.4163 Resultant Phase (radians): 0.83094
>> untitled
> untitled
> untitled?
No details available
vyy vy untitied2
fe >
Trang 8
2 B Section Break €
alk Ga Runond Advance op, Step Stop
lyze _ Pun Section P} Run to End
ON
© |B toto - Frzoes Sp file Est View Inert Tools Desktop Window Help
OUces BoB ra Name
9.29404, Sawtooth Wave
10.1582
“0 Of 02 03 04 05 08 07 08 09 1
No details available
fat 1.0060e-03 Zoom: 100% — 'UTF-B CRƯ script ln 20 G24
+ Filtering (các bộ Lọc IIR/FIR)
% Thiêf Lập Bộ Lọc Cơ Bản với filter()
Fs = 1000; % Tan sd’lay maui (Hz)
Fe = 100; % Tan sô căt của bộ Lọc (Hz)
Tạo tín hiệu đâù vào - một sóng sine + nhiêũ
0:1/FS:1-1/FS; % Thời gian từ @ dén 1 giây
sin(2 * pi * 50 * t) + 0.5 * randn(size(t)); % Séng sine 50 Hz với nhiêũ
oe Áp dụng bộ Lọc
= filter(b, a, x);
figure;
subplot(2,1,1);
plot(t, x);
title('Tin hiéu géé vai nhiéii');
xLabeL('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
plot(t, y);
title('Tin hiéu sau khi loc théng thap');
xlabel('Time (s)');
ylabel('Amplitude');
% Bộ Lọc Thông Cao (High-Pass Filter)
y_high = filter(b, a, x); % Ap dung bộ loc
figure;
Trang 9
subplot(2,1,1);
plot(t, x);
titLe('Tín hiệu gôé với nhiêñ');
xLabeL('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
plot(t, y_high);
title('Tin hiệu sau khi Lọc thông cao');
xlabel('Time (s)');
ylabel('Amplitude');
% Loc Bang Tan (Bandpass Filter)
f2 = 159; % Tâh sô cắt trên
y_band = filter(b, a, x);
figure;
subplot(2,1,1);
plot(t, x);
title('Tin hiệu gôé với nhiêñ");
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
plot(t, y_band);
title('Tin hiéu sau khi loc bang tan');
xlabel('Time (s)');
ylabel('Amplitude');
‘bandpass'); % Thiét ké’bé loc bang tan bac 2
New Open Save iy p GoTo find © ato SP analyze Pu
v.v v h§esmav v
untiveds *
Tin hiệu gẻc vời nhiễu
+ ba i Hi ah
S4 matlao ;
#3 umid vn hiệu Thi se=
®\ uniul g
# Se TUANZ —- Tme 6)
31 figure;
2 subplot(2,1,1 sinway (Tép WAV)
No details available
Eile Edt View Inset Joo’ Destop Window Help Oé6a8 3/08/ka
(1-26919,29484,
1 [00791,0,-0.1582
2050
file Edit View keed Tools Desktop Window Help
¬ø dã n8»
bai th lì i iN Ai i HA iM I
Tín hiệu sen tạo tông ân
Ìj x Thiết kế bộ lọc thôn
X Áp dụng bộ lọc
an Il
01 02 043 04 "6 07 08 09
Trang 10
% Biêñ Đôi Fourier Một Chiéu (1D Fourier Transform)
Fs = 1000; % Tâà sô Lâý maui (Hz)
T = 1/Fs; % Chu ky mat (s)
L = 1599; % Sô mâũ
t = (@:L-1) * T; % Vector thời gian
% Tạo một tín hiệu gôh hai tan sd! 50 Hz va 120 Hz
S = 0.7 * sin(2 * pi * 50 * t) + sin(2 * pi * 120 * t);
Thực hiện phép biêñ đôi Fourier
= fft(S);
P2 = abs(Y / L); % Phô hai phía
Pl = P2(1:L/2+1); % Chi Lay phô một phía
P1(2:end-1) = 2 * PI(2:end-1); % Nhân đôi các giá trị giữa
% Tạo trục tan sé’
f = Fs * (@:(L/2)) / L;
% Hiêñ thị phô biên độ
figure;
pLot(f, P1);
title('Single-Sided Amplitude Spectrum of S(t)');
xLabet('f (Hz)');
ylabel('|P1(f)|');
% Biêñ Đôi Fourier Hai Chiéu (2D Fourier Transform)
% Tao anh mâũ đơn giản (256x256)
M = 256;
N = 256;
image = zeros(M, N);
image(97:160, 97:160) = 1; % Hinh vuéng trang trén nén den
% Thực hiện biêñ đôi Fourier 2D
F = fft2(image);
% Chuyên đôi thành phô biên độ
F_magnitude = abs(fftshift(F)); % fftshift đê đặt tân sô “thâ6 ở trung tam
% Hiêñ thị phô biên độ
figure;
1magesc (Log(1 + F_magnitude)); % Dùng log dé’hiénh thi ré hon
colormap(jet);
colorbar;
title('Magnitude Spectrum of 2D Fourier Transform');
Bién Đôi Fourier Nghich (Inverse Fourier Transform)
Thực hiện phép biêñ đôi Fourier ngược
S_reconstructed = ifft(Y);
% Hiêñ thị tín hiệu ban đâù và tín hiệu khôi phục
figure;
subplot(2, 1, 1);
ptot(t, S);
titLe('Tín hiệu ban đâù');
xlabel('Time (s)');