BÀI TẬP THỰC HÀNH
XỬ LÝ TÍN HIỆU TRONG TRUYỀN THÔNG
BỘ LỌC SỐ IIR
Bài 6.1:
Kết quả sau khi chạy:
Bài 6.2
a. Bộ lọc lowpass
%Doc file am thanh file_goc.wav (duoc thu am o che do mono, tan so lay mau 44100 Hz)
>> x = audioread('file_goc.wav');
%Ve am thanh tren do thi
>> plot(x);
%Phat am thanh
>> sound(x);
%Tao cac diem tan so
>> f = [0 0.6 0.6 1];
%Dap ung bien do tuong ung
>> m = [1 1 0 0];
%Tao bo loc lowpass có bac la 40, cutoff frequency là 0.6
>> [b,a] = yulewalk(40,f,m);
%Xem dap ung tan so cua bo loc
>> freqz(b,a);
%Ap dung bo loc doi voi tin hieu x
>> y = filter(b,a,x);
%Phat am thanh sau khi loc
>> sound(y);
%Vẽ đồ thị âm thanh sau khi loc
>> plot(y);
b.Bộ lọc highpass
%Doc file am thanh file_goc.wav (duoc thu am o che do mono, tan so lay mau 44100 Hz)
>> x = audioread('file_goc.wav');
%Ve am thanh tren do thi
>> plot(x);
%Phat am thanh
>> sound(x);
%Tao cac diem tan so
>> f = [0 0.9 0.9 1];
%Dap ung bien do tuong ung
>> m = [0 0 1 1];
%Tao bo loc hightpass có bac la 40, cutoff frequency là 0.9
>> [b,a] = yulewalk(40,f,m);
%Xem dap ung tan so cua bo loc
>> freqz(b,a);
%Ap dung bo loc doi voi tin hieu x
>> y = filter(b,a,x);
%Phat am thanh sau khi loc
>> sound(y);
%Vẽ đồ thị âm thanh sau khi loc
>> plot(y);
c.Bộ lọc bandpass
%Doc file am thanh file_goc.wav (duoc thu am o che do mono, tan so lay mau 44100 Hz)
>> x = audioread('file_goc.wav');
%Ve am thanh tren do thi
>> plot(x);
%Phat am thanh
>> sound(x);
%Tao cac diem tan so
>> f = [0 0.6 0.9 1];
%Dap ung bien do tuong ung
>> m = [0 1 1 0];
%Tao bo loc lowpass có bac la 40, cutoff frequency là 0.6-0.9
>> [b,a] = yulewalk(40,f,m);
%Xem dap ung tan so cua bo loc
>> freqz(b,a);
%Ap dung bo loc doi voi tin hieu x
>> y = filter(b,a,x);
%Phat am thanh sau khi loc
>> sound(y);
%Vẽ đồ thị âm thanh sau khi loc
>> plot(y);