Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 12 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
12
Dung lượng
517,5 KB
Nội dung
Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 !"#$ !" #"#"$%&'! !"()*+#,# !")-&./012) #,##/!("3#4+&15!( #,#"4678!(9:;9<=!>?@ r AABACADE:;ℎ=!>?@ D r ACABAAE !% &#$ FG"-"H!(IJK#./012)%"JL!(MN!"4+)%"O!&'! !""-4H"5!(P)0")' ! ≥ =!>Q D =!R>S =!R>?<=!> +> N.I,%3!(<&!(:;I,%3!(M+3!(:T)I)U&V)-!IW&=Q>?=Q>?A V #"" #":;X<=!>? ! &=!>P 1> Y#,#4LIZ"[#")-!"-"H!("\X#,#7]!(#"&^!AP #>F,#I_!"";.M&U!I]#`+"-"H!(Pa"X1)'"-"H!(#bc!I_!""+ V"/!(:;()0)" #":N4+X#bV"/!(P % +#W! !"=!>?<=!>d"=!> 1 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 → =!>?F=!>P=!> +#b1)'!Ic)9#`+<=!>:;"=!>e; F=!>? B C B C Dz z z z z − − − − + + + + =!>? B C D C B z z z z z − − − − + + + + → =!>?Df Qg SCf Qh Sf Q SCf QD SDDf QC SCf QB Sf Q SCf Q SD → =!>?9@=!>E ? @DACAACADDACAACADE r !% =!>Q D =!R>S =!R>?<=!> => 'N.I,%3!(<&!(AI,%3!(M+ (%N.I,%3!(<&!( i)'!Ic)9"+):'+IJK# =f>Q D jf Q =f>S=>kS jf Q =f>Sf Q =Q>S=Q>k?F=f> → =f>PjQ D f Q S f Q k?F=f> → =f>? = > = > Y z X z ? D z z − − − + ? = B>= >z z − − − − ? B z z − − − − − ? B B z z − − − − − → "=!>? BP= > = > P= > = > B n n u n u n − P 2 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 ((%N.I,%3!(M+ i)'!Ic)9"+):'+IJK# =f>Q D jf Q =f>S=>kS jf Q =f>Sf Q =Q>S=Q>k? z z − → =f>PjQ D f Q S f Q k? z z − → =f>? = >= > D zz z z − − − +− ? = >= >= > B z z z − − − − − − ? g D D B z z z − − − + − − − − → =!>? g = > = > = > = > = > D D B n n n u n u n u n+ − P )%+#b=f>j D z z − − − + k?F=f> "+=f>? D = > = > = > X z z Y z z Y z − − + − LIZ#"&^!l#!"J4+& 3 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 LIZ#"&^!l#!"J4+& *%"\X#O&++#b "=!>? BP= > = > P= > = > B n n u n u n − → "=!>?:T).m)!n $"-"H!(Mo!!"O!p&0P 4 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 qrc!I_!""Nsf %A snPtV",#sf %s ?nusf % s?Bn!o!4&M+"- "H!(c!I_!"P % +P (+ (,-./0*1'23-/(42/'-/5+2/67589*./0*1'23-/(42:- ;<*=)>-?@2:-;<-/(,P Q vJ2!("[#")-!1)'!Ic)ww=7[+:;X.+e+1Ir !"X,!>:;1)r& 7)x!%"c!"$!IJK#yww Q qr1)'!Ic)ww+.2Mz!(Iz7;) !")-&";!" :;%"0)#"3+X;! 1z !")-&IW&:;XP Q aX7\ [y, Fs, nbits]=wavread('dsp_test01_01.wav'); L = length(y); % Lấy độ dài chuối tín hiệu T = 1/Fs; t = (0:L-1)*T; NFFT = 2^nextpow2(L); % Hàm nextpow2(L) sẽ trả về số p gần nhất thỏa mãn 2^p >= abs(L) % Với giá trị đó đủ để FFT % Thực hiện FFT với số điểm là NFFT Y = fft(y,NFFT)/L; %%% %%%Biểu diễn phổ 1 phía của tín hiệu đầu vào %%% f = Fs/2*linspace(0,1,NFFT/2+1); %Hàm linspace(a,b,c) tạo ra vecto 1*n gồm các số cách đều nhau từ a -> b %Ví dụ: linspace(0,1,5) = [0 0.25 0.5 0.75 1] %Vẽ phổ 1 phía plot(f,2*abs(Y(1:NFFT/2+1)),'y','linewidth',2); title('Single-Sided Amplitude Spectrum of Input Signal'); xlabel('f (Hz)'); ylabel('|Y(f)|'); set(gca,'Color',[0.5 0.5 0.5]); Q {'p&0 5 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 )% /(A2BA)CDE*F Q X])1zem#ww=N!")x&2";!"%"W!W!4H#+X>P Q "JL!(%",%")'V'[+:;X";.#6+4c MX!(Ib :T)w Q a"|)o&V}"&$=[#"m!S7[+:;X1)r&7)x!%"c2Mo!> o =#"|)o&!;#b"r()*e])";!"%"W!W!4H#L 10!nDDf> o =#"|)o&!;#b"r()*e])";!"%"W!W!4H#L10! nDDf> o "+ 6 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 o "+ o =7*e)-&y.+e+1> T) !"JMo!+4Y7~!(1zem#+ )!(P !"X,!"/!(4H • qzMz!(70)#"&r!)'% • qz7;)#`+1zem# +e•?D=!?€D>P • W!4H#l • F,#I_!"#,#"-4H#`+1zem# $ X"=!>#b !"IH)<3!(!o!+#"| !"#,#(),M_y"=> I'!"=D>P 7 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 +&Ib7_#"#,#(),M_"=!>:y+ !"IJK#4+!(%"0)DIL! :_Ae•IH)<3!(=IrI0.10X1zem##b !"!"O!p&0>P !"X,!7[+:;X.+e+1+#b • aX7\ fc = 0.10204; n = 25; Hd_n = zeros(1,26); W_n = zeros(1,26); Hn = zeros(1,26); N = 51; for i = 0:25 Hd_n(i+1) = 2*fc*sin(2*i*pi*fc)/(2*i*pi*fc); W_n(i+1) = 0.54+0.46*cos(2*pi*i/(N-1)); Hn(i+1) = Hd_n(i+1)*W_n(i+1); end; %Chỉnh lại Hn(1) = 2fc Hn(1) = 2*fc; %Dịch và lấy đối xứng Hn = [zeros(1,25),Hn]; %Lấy đối xứng for i = 0:25 Hn(i+1) = Hn(N-i); end; {'p&0+"&IJK##,#"-4H1zem# Hn = Columns 1 through 9 -0.0003 0.0003 0.0011 0.0016 0.0016 0.0007 -0.0013 -0.0037 -0.0055 Columns 10 through 18 -0.0051 -0.0016 0.0045 0.0111 0.0149 0.0126 0.0028 -0.0128 -0.0286 Columns 19 through 27 -0.0369 -0.0301 -0.0037 0.0409 0.0964 0.1504 0.1897 0.2041 0.1897 Columns 28 through 36 0.1504 0.0964 0.0409 -0.0037 -0.0301 -0.0369 -0.0286 -0.0128 0.0028 Columns 37 through 45 8 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 0.0126 0.0149 0.0111 0.0045 -0.0016 -0.0051 -0.0055 -0.0037 -0.0013 Columns 46 through 51 0.0007 0.0016 0.0016 0.0011 0.0003 -0.0003 *% G-/H(G)CDE* Q "[#")-!e•I,%3!(M+:T)V #"" #":;Xe; !")-&O."+!":;I,% 3!(<&!(#b#,#"-4H!"J%"W!1Pab Yn = conv(Hn,y); Q i)r&7)x!%"c#`+";.M&U!I] Q ("\"6 !")-&IW&M+:;4X4,!":T)IW&:;X %Đầu ra sound(Yn,Fs); %Đầu vào sound(y,Fs); Q "[#")-!ww !")-&"&IJK#Mo!IrN.M+%"cP L = length(Yn); T = 1/Fs; NFFT = 2^nextpow2(L); Yn_fft = fft(Yn,NFFT)/L; 9 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 Q i)r&7)x! Q X4,!":T)%"c !")-&IW&:;X 10 [...]... home exam DSP D10VT6Trịnh Phú Sơn – Mã SV: 1021010291 - Nhận xét: o Các thành phần tần số bậc cao đã được loại bỏ khi đi qua bộ lọc được thiết kế ở trên o Bộ lọc tuy có dải chuyển tiếp khá rộng nhưng vẫn có thể thực hiện lọc khá tốt Tuy nhiên nếu các tần số nhiễu nằm trong dải chuyển tiếp này thì hiệu quả của bộ lọc sẽ giảm đi khá nhiều Code toàn bộ bài 3 – Không có chú giải: %////////////////... plot(f,2*abs(Y(1:NFFT/2+1)),'y','linewidth',2); title('Single-Sided Amplitude Spectrum of Input Signal'); xlabel('f (Hz)'); ylabel('|Y(f)|'); set(gca,'Color',[0.5 0.5 0.5]); 11 Take home exam DSP D10VT6Trịnh Phú Sơn – Mã SV: 1021010291 %//////////////// %////Câu b %//////////////// fc = 0.10204; n = 25; N = 51; Hd_n = zeros(1,26); W_n = zeros(1,26); Hn = zeros(1,26); for i = 0:25 Hd_n(i+1) = 2*fc*sin(2*i*pi*fc)/(2*i*pi*fc); . dài chuối tín hiệu T = 1/Fs; t = (0:L-1)*T; NFFT = 2^nextpow2(L); % Hàm nextpow2(L) sẽ trả về số p gần nhất thỏa mãn 2^p >= abs(L) % Với giá trị đó đủ để FFT % Thực hiện FFT với số điểm là. fft(y,NFFT)/L; %%% %%%Biểu diễn phổ 1 phía của tín hiệu đầu vào %%% f = Fs/2*linspace(0,1,NFFT/2+1); %Hàm linspace(a,b,c) tạo ra vecto 1*n gồm các số cách đều nhau từ a -> b %Ví dụ: linspace(0,1,5). exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 Q i)r&7)x! Q X4,!":T)%"c !")-&IW&:;X 10 Take home exam. DSP D10VT6Trnh Ph Sơn – Mã SV: 1021010291 Q