語音信號處理與濾波課程設(shè)計-畢業(yè)設(shè)計-資料下載頁
2025-01-19 02:12本頁面
【導(dǎo)讀】所學(xué)知識的理解。熟悉離散信號和系統(tǒng)的時域特性。掌握序列快速傅里葉變換FFT方法。學(xué)會MATLAB的使用，掌握MATLAB的程序設(shè)計方法。利用MATLAB對語音信號進(jìn)行頻譜分析。掌握MATLAB設(shè)計FIR和IIR數(shù)字濾波器的方法。y=wavread，讀取file所規(guī)定的wav文件，返回采樣值放在向量y中。[y,fs,nbits]=wavread，采樣值放在向量y中，fs表示采樣頻率，nbits表示采樣位數(shù)。②高通濾波器性能指標(biāo)：fc=4800Hz，fp=5000Hz，As=100dB，Ap=1dB；③帶通濾波器性能指標(biāo)：fp1=1200Hz，fp2=3000Hz，fc1=1000Hz，首先用窗函數(shù)法（矩形窗(Rectangularwindow)、三角窗(Triangularwindow)、漢。寧窗、海明窗、布拉克曼窗(Blackman. window)、切比雪夫窗、巴特里特窗及凱塞。）設(shè)計上面要求的三種濾波器。
【正文】 lot(2,1,1)。plot(f,abs(y2(1。512)))。title(39。IIR 帶通濾波器濾波前的頻譜 39。)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 subplot(2,1,2)。plot(f,abs(y1(1。512)))。title(39。IIR 帶通濾波器濾波后的頻譜 39。)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 fs=22050。 x1=wavread(39。D:\39。)。 Ts=1/fs。R1=50。 wp=2*pi*5000/fs。 ws=2*pi*4800/fs。 Rp=1。 Rl=100。 wp1=2/Ts*tan(wp/2)。ws1=2/Ts*tan(ws/2)。 [N,wn]=chebyord(wp1,ws1,Rp,Rl,39。s39。)。 [Z,P,K]=chebyap(N,R)。 [Bap,Aap]=zpwtf(Z,P,K)。 [b,a]=Ip2bp(Bap,Aap,2100*2*pi,1800*2*pi)。 [bz,az]=bilinear(b,a,fs)。 [H,W]=freqz(bz,az)。 figure(1) plot(w*fs/(2*pi),abs(H)) grid xlabel(39。頻率 /Hz39。)。 ylabel(39。頻率響應(yīng)幅度 39。)。 title(39。IIR高通濾波器 39。) fl=filter(bz,az,x1)。 figure(2) y1=fft(f1,1024)。 y2=fft(x1,1024)。 subplot(2,1,1)。plot(f,abs(y2(1。512)))。title(39。IIR 高通濾波器濾波前的頻譜 39。)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 subplot(2,1,2)。plot(f,abs(y1(1。512)))。title(39。IIR 高通濾波器濾波后的頻譜 39。)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 fs=22050。 x1=wavread(39。D:\39。)。 Ts=1/fs。R1=10。 wp=2*pi*1000/fs。 ws=2*pi*1200/fs。 Rp=1。 Rl=100。 wp1=2/Ts*tan(wp/2)。ws1=2/Ts*tan(ws/2)。 [N,wn]=buttord(wp1,ws1,Rp,Rl,39。s39。)。 [Z,P,K]=buttap(N,R)。 [Bap,Aap]=zpwtf(Z,P,K)。 [b,a]=Ip2bp(Bap,Aap,2100*2*pi,1800*2*pi)。 [bz,az]=bilinear(b,a,fs)。 [H,W]=freqz(bz,az)。 figure(1) plot(w*fs/(2*pi),abs(H)) grid xlabel(39。頻率 /Hz39。)。 ylabel(39。頻率響應(yīng)幅度 39。)。 title(39。IIR低通濾波器 39。) fl=filter(bz,az,x1)。 figure(2) y1=fft(f1,1024)。 y2=fft(x1,1024)。 subplot(2,1,1)。plot(f,abs(y2(1。512)))。title(39。IIR 低通濾波器濾波前的頻譜 39。)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 subplot(2,1,2)。plot(f,abs(y1(1。512)))。title(39。IIR 低通濾波器濾波后的頻譜 39。)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 [x1,fs,bits]=wavread(‘D:\’)。 y1=fft(x1,1024)。 f=fs*(0:511)/1024。 figure(1) subplot(2,1,1)。 plot(f,abs(1:512)))。 title(’原始語音信號頻譜 ’)。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 subplot(2,1,2)。 plot(abs(y1(1:1024)))。 title(‘原始語音信號 FFT頻譜 ’)。 xlabel(39。點數(shù) N39。)。 ylabel(39。幅值 39。)。 [x1,fs,bits]=wavread(‘D:\’)。 sound(x1,fs)。 plot(x1)。title(39。原始語音采樣后時域信號 39。)。 xlabel(39。時間軸 n39。)。 ylabel(39。幅值 A39。)。
試題試卷相關(guān)推薦
鄂ICP備17016276號-1