吹管音樂濾波去噪基于漢寧創(chuàng)的fir濾波器-資料下載頁
2025-06-27 23:09本頁面
【正文】 。title(39。以db為單位的幅度特性39。)。xlabel(39。w/pi39。)。ylabel(39。db39。)。 grid on。subplot(2,2,2)。plot(w/pi,mag)。axis([0,])。title(39。以線性為單位的幅度特性39。)。xlabel(39。w/pi39。)。ylabel(39。mag39。)。grid on。subplot(2,2,3)。plot(w,pha)。title(39。濾波器相位響應圖39。)。xlabel(39。w/pi39。)。ylabel(39。相位（pha）39。)。axis([0,3,4,4])。grid on。subplot(2,2,4)。plot(h_bs)。axis([0,800,1])。title(39。濾波器脈沖響應圖39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。grid on。 y_fil=filter(h_bs,1,y)。 % 用設計好的濾波器對y進行濾波Y_fil=fft(y_fil)。Y_fil=Y_fil(1:N/2)。 % 計算頻譜取前一半figure(3)subplot(3,2,1)。plot(t,x)。 axis([0,2,])。title(39。原始音樂信號時間x39。)。xlabel(39。時間（t）39。)。ylabel(39。幅度39。)。grid on。subplot(3,2,2)。plot(f,X)。axis([0,10000,0,1000])。title(39。原始音樂信號幅度譜X39。)。xlabel(39。頻率（f）39。)。ylabel(39。幅度39。)。grid on。subplot(3,2,3)。plot(t,y)。axis([0,2,])。title(39。加干擾音樂信號時間x139。)。xlabel(39。時間（t）39。)。ylabel(39。幅度39。)。grid on。subplot(3,2,4)。plot(f,Y)。axis([0,10000,0,1000])。title(39。加干擾音樂信號幅度譜X139。)。xlabel(39。頻率（f）39。)。ylabel(39。幅度39。)。grid on。subplot(3,2,5)。plot(t,y_fil)。axis([0,2,])。title(39。濾波后音樂信號時間y39。)。xlabel(39。時間（t）39。)。ylabel(39。幅度39。)。grid on。subplot(3,2,6)。plot(f,Y_fil)。axis([0,10000,0,1000])。title(39。濾波后音樂信號幅度譜Y39。)。xlabel(39。頻率（f）39。)。ylabel(39。幅度39。)。grid on。sound (y_fil,fs,bits)。附錄2：function hd = ideal_lp(wc,M)。% 理想低通濾波器計算% % [hd] = ideal_lp(wc,M)% hd = 0 to M1之間的理想脈沖響應% wc = 截止頻率(弧度) % M = 理想濾波器的長度alpha = (M1)/2。n = [0:1:(M1)]。m = n alpha + eps。hd = sin(wc*m) ./ (pi*m)。附錄3：function [db,mag,pha,grd,w] = freqz_m(b,a)。% freqz 子程序的改進版本% % [db,mag,pha,grd,w] = freqz_m(b,a)。% db = [0 到pi弧度]區(qū)間內的相對振幅(db)% mag = [0 到pi弧度]區(qū)間內的絕對振幅% pha = [0 到pi弧度]區(qū)間內的相位響應% grd = [0 到pi弧度]區(qū)間內的群遲延% w = [0 到pi弧度]區(qū)間內的501個頻率樣本向量% b = Ha(z)的分子多項式系數(shù)(對FIR b=h)% a = Ha(z)的分母多項式系數(shù)(對 FIR: a=[1])% [H,w] = freqz(b,a,1000,39。whole39。)。H = (H(1:1:501))39。 w = (w(1:1:501))39。mag = abs(H)。db = 20*log10((mag+eps)/max(mag))。pha = angle(H)。% pha = unwrap(angle(H))。grd = grpdelay(b,a,w)。% grd = diff(pha)。% grd = [grd(1) grd]。% grd = [0 grd(1:1:500)。 grd。 grd(2:1:501) 0]。% grd = median(grd)*500/pi。
環(huán)評公示相關推薦
鄂ICP備17016276號-1