Đang tải... (xem toàn văn)
Bài 3.3 Cho phương trình sai phân: y(n)0.6(n1)+0.9y(n2)=x(n)+0.8x(n1)FIR 1: cho đáp ứng xung bộ lọc fir như sau h1=1 2 1.3 2.2 0.6 3 0.6 2.2 1.3 2 1;Xác định đáp ứng tần số và phân bố cực không của bộ lọcHàmcho đáp ứng xung bộ lọc fir 4cho đáp ứng xung bộ lọc fir 2cho đáp ứng xung bộ lọc fir 3Thiết kế bộ lọc thông thấp dùng cửa sổ blackmanThiết kế bộ lọc thông cao dùng cửa sổ blackmanThiết kế bộ lọc thông thấp dùng cửa sổ kasiser
Bài 3.3 Cho phương trình sai phân: y(n)-0.6(n-1)+0.9y(n-2)=x(n)+0.8x(n-1) Y ( Z ) − 0.6Z −1Y ( Z ) + 0.9Z −2Y ( Z ) = X ( Z ) + 0.8Z −1 X ( Z ) Y (Z ) + 0.8Z −1 H (Z ) = = X ( Z ) − 0.6 Z −1 + 0.9 Z −2 H(e => jω ) Y (e jω ) + 0.8e − jω = H (Z) |z =e jω = = − jω −2 jω jω X ( e ) − 0.6e + 0.9e Chương trình mô sau: num=[1 0.8]; den=[1 -0.6 0.9]; w=-4*pi:8*pi/511:4*pi h=freqz(num,den,w); subplot(231) plot(w/pi,real(h));grid on;title('Phan thuc cua X(e^{j\omega})'); xlabel('\omega/\pi'); ylabel('bien do'); subplot(232) plot(w/pi,imag(h));grid on; title('Phan ao cua X(e^{j\omega})'); xlabel('\omega/\pi'); ylabel('bien do'); subplot(233); plot(w/pi,abs(h));grid on; title('Pho bien X(e^{j\omega})'); xlabel('\omega/\pi'); ylabel('bien do'); subplot(234); plot(w/pi,angle(h));grid on; title('pho phase X(e^{j\omega})'); xlabel('\omega/\pi'); ylabel('pha bang radian'); subplot(235) h=impz(num,den); stem(h); grid on; title ('dap ung xung don vi'); xlabel('chi so thoi gian'); ylabel('bien do'); subplot(236) zplane(num,den);grid on; title('Bieu diem cuc- diem khong'); FIR 1: cho đáp ứng xung lọc fir sau h1=[-1 1.3 -2.2 0.6 0.6 -2.2 1.3 -1]; Xác định đáp ứng tần số phân bố cực không lọc *Hàm function [a,w,L,Hr]=loc1(h) M=length(h); L=(M-1)/2; a=[h(L+1) 2*h(L:-1:1)] n=[0:1:L]; w=linspace(0,2*pi,100)'; Hr=cos(w*n)*a'; stem(Hr); end *tạo file gõ h1=[-1 1.3 -2.2 0.6 0.6 -2.2 1.3 -1]; M=length(h1); n=0:M-1; subplot(221) stem(n,h1); [a,w,L,Hr]=loc1(h1); title(' dap ung xung '); xlabel('n'); ylabel('h(n)'); subplot(222) stem(0:L,a); title('He so a(n)'); xlabel('n'); ylabel('a(n)'); subplot(223) w=linspace(0,2*pi,100); plot(w,Hr); title('dap ung tan so'); xlabel('\omega'); ylabel('H(\omega'); subplot(224) zplane(h1,1); title('Bieu diem cuc _ diem khong'); xlabel('thuc’); ylabel('ao'); FIR 2: cho đáp ứng xung lọc FIR sau h2=[1.5 -2.5 3 -2.5 1.5]; Xác định đáp ứng tần số phân bố cực không lọc *Hàm function [b,w,L,Hr]=loc2(h) M=length(h); L=M/2; b=[2*h(L:-1:1)]; n=[1:1:L]; w=linspace(0,2*pi,100)'; Hr=cos(w*(n-0.5))*b'; Stem(Hr); end *tạo file gõ: h2=[1.5 -2.5 3 -2.5 1.5]; M=length(h2); n=0:M-1; subplot(221) stem(n,h2); [b,w,L,Hr]=loc2(h2); title('Dap ung xung'); xlabel('n'); ylabel('h(n)'); subplot(222) stem(1:L,b); title('Cac he so b(n)'); xlabel('n'); ylabel('b(n)'); subplot(223) w=linspace(0,2*pi,100); plot(w,Hr); title('Dap ung tan so'); xlabel('\omega'); ylabel('H(\omega)'); subplot(224) zplane(h2,1); title('Bieu diem cuc diem 0'); xlabel('Thuc'); ylabel('ao'); FIR 3: cho đáp ứng xung lọc FIR sau h3=[1.5 -2.5 2.5 -1.5]; Xác định đáp ứng tần số phân bố cực không lọc *hàm function [c,w,L,Hr]=loc3(h) M=length(h); L=(M-1)/2; c=[2*h(L:-1:1)]; n=[1:1:L]; w=linspace(0,2*pi,100)'; Hr=sin(w*n)*c'; stem(Hr); end *tạo file mới: h3=[1.5 -2.5 2.5 -1.5] M=length(h3); n=0:M-1; subplot(221) stem(n,h3); [c,w,L,Hr]=loc3(h3) title(' Dap ung xung'); xlabel('n'); ylabel('h(n)'); subplot(222) stem(1:L,c); title('Cac he so c(n)'); xlabel('n'); ylabel('c(n)'); subplot(223) w=linspace(0,2*pi,100); plot(w,Hr); title('dap úng tan so'); xlabel('\omega'); ylabel('H(\omega)'); subplot(224) zplane(h3,1); title('Bieu diem cuc - diem khong'); xlabel('Thuc'); ylabel('Ao'); FIR 4: cho đáp ứng xung lọc FIR sau h4=[1.5 -2.5 -3 2.5 -1.5]; Xác định đáp ứng tần số phân bố cực không lọc *Hàm: function [d,w,L,Hr]=loc4(h) M=length(h); L=M/2 d=[2*h(L:-1:1)]; n=[1:1:L]; w=linspace(0,2*pi,100)'; Hr=sin(w*(n-0.5))*d'; stem(Hr); end *tạo file mới: h4=[1.5 -2.5 -3 2.5 -1.5] M=length(h4); n=0:M-1; subplot(221) stem(n,h4); title('dap ung xung'); xlabel('n'); ylabel('h(n)'); [d,w,L,Hr]=loc4(h4); subplot(222) stem(1:L,d); title('cac he so d(n)'); xlabel('n'); ylabel('d(n)'); subplot(223) w=linspace(0,2*pi,100); plot(w,Hr); title('dap ung tan so'); xlabel('\omega'); ylabel('H(\omega)'); subplot(224) zplane(h4,1); title('Bieu diem cuc-diem khong'); xlabel('thuc'); ylabel('ao'); Chương 2: BÀI 1:Viết chương trình vẽ tín hiệu dạng hàm u(n-3),u(n+2) *Hàm stepseq function [ x,n ] = stepseq( n0,n1,n2 ) if((n0>n2)|(n0n2)|(n0n2)) error('arguments must satisfy n1[...]... X[k]'); xlabel('chi so tan so k'); ylabel('bien do'); subplot(2,1,2); plot(k,angle(X)); grid on; title('pha cua DFT X[k]'); xlabel('chi so tan so k'); ylabel('pha cua X(k)'); Bài 5c: Viết chương trình tính và thể hiện trên đồ thị fourier rời rạc của các dãy sau: N=60 và 1 0 ≤ n ... title('Bieu diem cuc-diem khong'); xlabel('thuc'); ylabel('ao'); Chương 2: BÀI 1:Viết chương trình vẽ tín hiệu dạng hàm u(n-3),u(n+2) *Hàm stepseq function [ x,n ] = stepseq( n0,n1,n2 ) if((n0>n2)|(n0