HỌC VIỆN CÔNG NGHỆ BƯU CHÍNH VIỄN THÔNG CƠ SỞ TẠI THÀNH PHỐ HỒ CHÍ MINH KHOA KỸ THUẬT ĐIỆN TỬ II BÁO CÁO MÔN HỌC XỬ LÝ TIẾNG NÓI NIÊN KHÓA 2017 – 2022 Nhóm thực hiện Nhóm 8 Giáo viên hướng dẫn ThS H.
Trang 1HỌC VIỆN CÔNG NGHỆ BƯU CHÍNH VIỄN THÔNG
CƠ SỞ TẠI THÀNH PHỐ HỒ CHÍ MINH
KHOA KỸ THUẬT ĐIỆN TỬ II
BÁO CÁO MÔN HỌC
Trang 3PHẦN I: THỰC HIỆN MÔ PHỎNG
PROJECT BẰNG PHẦN MỀM MATLAB
1 Đọc ra một file tiếng nói và bộ lọc tín hiệu tiếng nói với băng tần 5.5kHz, 4 kHz và 3.2 kHz (Tự thiết kế các bộ lọc số với tốc độ lấy mẫu chính xác để thực hiện việc này) Lắng nghe mỗi file tiếng nói sau khi được lọc và mô tả ảnh hưởng của việc lọc thông thấp đối với chất lượng và tính dễ hiểu của tiếng nói
[x,fs]=audioread('demso.wav'); %doc dinh dang file wav
X=fft(x); %bien doi Fourier khong nam o trung tam
X = X(:,2);
N=length(x); %Xac dinh kich thuoc cua x
a = 5500; %tan so cat
w1 = (-N/2+1:(N/2)); % Vector tan so trung tam
w = w1.*fs/N; % LAY MOT TAN SO MAU
H = a./(a + 1i*w); %H nam o trung tam
Hshift = fftshift(H); %H khong nam o trung tam
Y = X *Hshift' ; % loc tin hieu
y = real(ifft(Y));
sound(x,fs);% am thanh goc
sound(y,fs); % am thanh sau khi qua bo loc thong thap
subplot(2,1,1);%vi tri ve
plot(w,abs(fftshift(X))) % abs lay bien do cua so phuc, dich tan so ve trung tam
Trang 4title('Tin hieu ban dau'); %ten tieu de do thi
subplot(2,1,2);%vi tri ve
plot(w,abs(fftshift(Y)))
title('Tin hieu loc') %ten tieu de do thi
2 Với cùng một nội dung thông tin, các thành viên trong nhóm lần lượt phát âm một bản tin ngắn và ghi âm Lưu kết quả ở dạng
*.wav Xác định tần số cơ bản, tần số Formant đầu tiên của mỗi thành viên Lập bản đồ phân bố của các nguyên âm trong tiếng Việt.
% get a section of vowel
% su dung bo loc du doan tuyen tinh
ncoeff=2+fs/1000; % Cong thuc tinh formant
a=lpc(x,ncoeff); % he so bo loc tuyen tinh
% ve tan so dap ung
Trang 5[h,f]=freqz(1,a,512,fs); % dap ung tan so cua bo loc
% Tim tan so bang roots
r=roots(a); % Tim nghiem nguyen
r=r(imag(r)>0.01); % chi tim >0Hz up to fs/2
Trang 6% tinh toan cepstrum phuc và thuc cua các tin hieu
xhat1=cceps(x1,NFFT); %cepstrum phuc cua tin hieu 1
rhat1=rceps(x1); %cepstrum thuc cua tin hieu 1
xhat2=cceps(x2,NFFT); %cepstrum phuc cua tin hieu 2
rhat2=rceps(x2); %cepstrum thuc cua tin hieu 2
xhat3=cceps(x3,NFFT); %cepstrum phuc cua tin hieu 3
rhat3=rceps(x3); %cepstrum thuc cua tin hieu 3
% ve dang song, cepstrum phuc, cepstrum thuc cua tin hieu 1subplot(3,1,1);
plot(n1,x1);
title('dang song')
subplot(3,1,2);
plot(n1,xhat1);
Trang 85 Sử dụng phương pháp tổng hợp trực tiếp đơn giản để xây dựng một hệ thống thông báo điểm đỗ xe bus công cộng
Trang 9clc
stt = input ('Nhap stt: '); play(tb);
Trang 106 Sử dụng phương pháp tổng hợp trực tiếp đơn giản xây dựng một
hệ thống thông báo số thứ tự khách hàng đến lượt được phục vụ tại một điểm giao dịch ngân hàng
Trang 12pause(3); play(y1); case 2
pause(3); play(y2); case 3
pause(3); play(y3); case 4
pause(3); play(y4); case 5
pause(3); play(y5); case 6
pause(3); play(y6); case 7
pause(3); play(y7); case 8
pause(3); play(y8); case 9
pause(3); play(y9);
Trang 13pause(3); play(y1); case 2
pause(3); play(y2); case 3
pause(3); play(y3); case 4
pause(3); play(y4); case 5
pause(3); play(y5); case 6
pause(3); play(y6); case 7
pause(3);
Trang 14play(y7); case 8
pause(3); play(y8); case 9
pause(3); play(y9); end
switch c
case 0
pause(3); play(y0); case 1
pause(3); play(y1); case 2
pause(3); play(y2); case 3
pause(3); play(y3); case 4
pause(3); play(y4); case 5
Trang 15pause(3); play(y5); case 6
pause(3); play(y6); case 7
pause(3); play(y7); case 8
pause(3); play(y8); case 9
pause(3); play(y9); end
switch d
case 0
pause(3); play(y0); case 1
pause(3); play(y1); case 2
pause(3); play(y2); case 3
Trang 16pause(3); play(y3); case 4
pause(3); play(y4); case 5
pause(3); play(y5); case 6
pause(3); play(y6); case 7
pause(3); play(y7); case 8
pause(3); play(y8); case 9
pause(3); play(y9); end
end
Trang 17
khối phi tuyến tạo ra tín hiệu đầu ra y(n) =
a Nếu tín hiệu đầu vào có dạng x1[n] = cos(ω0n), -∞< n <∞ xác định phổ
Trang 18đầu ra Y1(ejω), và vẽ đồ thị biên độ của Y1(ejω) (giả sử ω0 << ωs, tần
b Nếu tín hiệu đầu vào có dạng x2 (n) = cos(n)u(n), | r |<1, vẽ biên độ
của tín hiệu đầu ra Y2() (giả sử << ) và giả sử = 2π.500, r = 0.9 và Fs = 10000Hz
Trang 204 Xem xét một hệ thống bậc nhất: y(n) = αy(n-1) + x(n).
a Tìm hàm H(z) của hệ thống
Biến đổi z 2 vế cho y(n)
⇨ Y(z)(1- α) = X(z)
H(z)=Y(z)/X(z) =1/1-�
b Tìm hàm đáp ứng xung ngẫu nhiên của hệ thống
Đáp ứng xung ngẫu nhiên của hệ thống là :
Trang 216 Việc tính toán năng lượng ngắn hạn thông qua biểu thức sau
sẽ có giá trị hữu hạn vậy sẽ tồn tại năng lượng En
7 Hàm nào dưới đây không phải là hàm tự tương quan? Tại sao?
a R(τ) = 2, -∞ < � < ∞
Có vì R(0) có giá trị lớn nhất
b.R(τ) = |�|, -∞ < � < ∞
Trang 22Không vì R(0) không phải là giá trị lớn nhất
c R(τ) = (sin�� (��))2, -∞ < � < ∞
Có vì R(0) có giá trị lớn nhất
d R(τ) = 2(+4) / +6), -∞ < � < ∞
Trang 23Không vì R(0) không phải là giá trị lớn nhất
Trang 24Vậy ta sẽ không vẽ được phổ chính xác của tín hiệu.
c Tìm và biểu diễn hàm tự tương quan của tín hiệu:
Trang 28b Giả sử thời tiết ngày đầu tiên là nắng, tìm xác suất để thời tiết cho 7 ngày
tiếp theo là “nắng-nắng-mưa-mưa-nắng-mây-nắng
trạng thái đó trong vòng d ngày
Gọi Pi là xác xuất giữ 1 trạng thái thời tiết trong d ngày pn xác xuất củathời tiết đó được lặp lại vào ngày hôm sau ta có :
Và xác suất chuyển bằng nhau bằng 1/3
(Giả sử xác suất trạng thái đầu tiên là 1/3 )
a Quan sát chuỗi sự kiện O=HHHHTHTTT Chuỗi sự kiện trên tương đương
với chuỗi trạng thái nào nhất? Xác suất của chuỗi sự kiện quan sát và
Trang 29trạng thái này là bao nhiêu?
Xác suất của chuỗi sự kiện trên là 1/ = 1/512
Xác suất của chuỗi sự kiện trên gần với chuổi trạng thái
Xác suất của chuỗi sự kiện trên là 1/ = 1/512
Xác suất của chuỗi sự kiện trên gần với chuỗi trạng thái
S2S2S2S2S3S2S3S3S3
Trạng thái = 0.759 *0.1*0.1*0.1*0.45*0.45*0.45*0.1*0.1 = 6.84x1