Lab 6 xu li tin hieu trong truyen thong

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. Thiết kế bộ lọc thông thấp Butterworth theo các tham số sau: wp=0.2pi, Rp=1dB ws=0.3pi, As=15dB % Butterworth Lowpass Filter Design: % % Digital Filter Specifications: >>wp = 0.2pi; %digital Passband freq in Hz >>ws = 0.3pi; %digital Stopband freq in Hz >>Rp = 1; %Passband ripple in dB >>As = 15; %Stopband attenuation in dB >>ep = sqrt(10(Rp10)1); >>Ripple = sqrt(1(1+epep)); % Passband Ripple >>Attn = 1(10(As20)); % Stopband Attenuation

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);