多速率处理研究性专题.docxVIP

《数字信号处理》课程研究性学习报告 ——多速率处理研究性专题  PAGE15 / NUMPAGES16 例 1 抽样频率转换 试将一段抽样频率为44.1kHz的音乐信号转换为抽样频率为22.05kHz的音乐信号。 解： Matlab仿真 信号x H(z)=1时输出信号y H(z)为Ws=π/2 低通滤波器时输出信号y 直接将抽样频率设置为22.05kHz Matlab程序 1、 [x,fs,Nbits]=wavread(Love.wav); sound(x,fs); plot(x); title(语音信号x的时域波形); ylabel(x); 2、 M=2; [x,fs,Nbits]=wavread(Love.wav); y=x(1:M:end); sound(y,fs); plot(y); title(两倍抽取后语音信号y的时域波形); wavwrite(y,fs/M,Love1_1.wav) ylabel(y); 3、 M=2; [x,fs,Nbits]=wavread(Love.wav); y=x(1:M:end); wavwrite(2*y,fs/M,Love1_2.wav); k=linspace(-fs,fs,length(y)); Y=abs(fft(y)); plot(k,fftshift(Y)); xlabel(k); ylabel(Y); 4、 [x,fs,Nbits]=wavread(Love.wav); wavwrite(x,fs/2,Love1_3.wav); [y,fs,Nbits]=wavread(Love1_3.wav); sound(y,fs); plot(y); title(两倍抽取后语音信号的时域波形y); 例2 抽样频率转换 试将一段抽样频率为44.1kHz的音乐信号转换为抽样频率为11.025kHz的音乐信号。 解： Matlab仿真 信号x H(z)=1时输出信号y H(z)为Ws=π/4 低通滤波器时输出信号y Matlab程序 1、 [x,fs,Nbits]=wavread(Love.wav); sound(x,fs); plot(x); title(语音信号x的时域波形); ylabel(x); 2、 M=4; [x,fs,Nbits]=wavread(Love.wav); y=x(1:M:end); sound(y,fs); plot(y); title(四倍抽取后语音信号y的时域波形); wavwrite(y,fs/M,Love2_1.wav) ylabel(y); 3、 M=4; [x,fs,Nbits]=wavread(Love.wav); y=x(1:M:end); wavwrite(4*y,fs/M,Love2_2.wav); k=linspace(-fs,fs,length(y)); Y=abs(fft(y)); plot(k,fftshift(Y)); xlabel(k); ylabel(Y); 例3 抽样频率转换 试将一段抽样频率为44.1kHz的音乐信号转换为抽样频率为8kHz的音乐信号。 解： 方案1（6倍抽取） Matlab仿真 H(z)=1时输出信号y H(z)为Ws=π/6 ??通滤波器时输出信号y Matlab程序 1、 M=6; [x,fs,Nbits]=wavread(Love.wav); y=x(1:M:end); sound(y,fs); plot(y); title(六倍抽取后语音信号y的时域波形); wavwrite(y,fs/M,Love3_1_1.wav) ylabel(y); 2、 [x,fs,Nbits]=wavread(Love.wav); k=linspace(-fs,fs,length(x)); X=abs(fft(x)); plot(k,fftshift(X)); xlabel(k); ylabel(X); figure M=6; [x,fs,Nbits]=wavread(Love.wav); y=x(1:M:end); wavwrite(6*y,fs/M,Love3_1_2.wav); k=linspace(-fs,fs,length(y)); Y=abs(fft(y)); plot(k,fftshift(Y)); xlabel(k); ylabel(Y); 方案2（5倍抽取） Matlab仿真 H(z)=1时输出信号y H(z)为Ws=π/5 低通滤波器时输出信号