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 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ộ