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 sau chạy: Bài 6.2 a Bộ lọc lowpass %Doc file am file_goc.wav (duoc thu am o che mono, tan so lay mau 44100 Hz) >> x = audioread('file_goc.wav'); %Ve am tren thi >> plot(x); %Phat am >> sound(x); %Tao cac diem tan so >> f = [0 0.6 0.6 1]; %Dap ung bien tuong ung >> m = [1 0]; %Tao bo loc lowpass có bac la 40, cutoff frequency 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 sau loc >> sound(y); %Vẽ đồ thị âm sau loc >> plot(y); b.Bộ lọc highpass %Doc file am file_goc.wav (duoc thu am o che mono, tan so lay mau 44100 Hz) >> x = audioread('file_goc.wav'); %Ve am tren thi >> plot(x); %Phat am >> sound(x); %Tao cac diem tan so >> f = [0 0.9 0.9 1]; %Dap ung bien tuong ung >> m = [0 1]; %Tao bo loc hightpass có bac la 40, cutoff frequency 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 sau loc >> sound(y); %Vẽ đồ thị âm sau loc >> plot(y); c.Bộ lọc bandpass %Doc file am file_goc.wav (duoc thu am o che mono, tan so lay mau 44100 Hz) >> x = audioread('file_goc.wav'); %Ve am tren thi >> plot(x); %Phat am >> sound(x); %Tao cac diem tan so >> f = [0 0.6 0.9 1]; %Dap ung bien tuong ung >> m = [0 1 0]; %Tao bo loc lowpass có bac la 40, cutoff frequency 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 sau loc >> sound(y); %Vẽ đồ thị âm sau loc >> plot(y); ... >> sound(x); %Tao cac diem tan so >> f = [0 0 .6 0 .6 1]; %Dap ung bien tuong ung >> m = [1 0]; %Tao bo loc lowpass có bac la 40, cutoff frequency 0 .6 >> [b,a] = yulewalk(40,f,m); %Xem dap ung tan... [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 sau loc >> sound(y); %Vẽ đồ thị âm sau loc >> plot(y); b.Bộ... [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 sau loc >> sound(y); %Vẽ đồ thị âm sau loc >> plot(y); c.Bộ