数字信号处理中语音信号的抽取，混叠和数字限带滤波MATLAB

clear;clc; D=30; %% 原声频域时域图像 [y,fs]=read('乐声.mp3'); figure(1); subplot(2,1,1); pinpu(y,fs);title('原音频谱'); subplot(2,1,2);shiyu(y,fs);title('原音时域波形'); %% resample频域时域图像 [y1,fs]=resampleDb('乐声.mp3',D); figure(2); subplot(2,1,1);pinpu(y1,fs);title('resample 抽样的频谱'); subplot(2,1,2);shiyu(y1,fs);title('resample 抽样的时域波形'); %% 直接抽取频域时域图像 [y2,fs]=directDb('乐声.mp3',D); figure(3); subplot(2,1,1);pinpu(y2,fs);title('直接抽样的频谱'); subplot(2,1,2);shiyu(y2,fs);title('直接抽样的时域波形'); %% 卷积滤波后抽样的频域时域图像 tic; [y3,fs]=lvbodirectD('乐声.mp3',D); figure(4); subplot(2,1,1);pinpu(y3,fs);title('卷积滤波后抽样的频谱'); subplot(2,1,2);shiyu(y3,fs);title('卷积滤波后抽样的时域波形'); t1=toc; fprintf('卷积滤波后抽样时间%d\n ',t1); %% 多相滤波后抽样的频域时域图像 tic; [y4,fs]=polyphase_filterD('乐声.mp3',D); figure(5); subplot(2,1,1);pinpu(y4,fs);title('多相滤波抽样的频谱'); subplot(2,1,2);shiyu(y4,fs);title('多相滤波抽样的时域波形'); t2=toc; fprintf('多相滤波后抽样时间%d\n',t2); %% 插值时域频域图像 % [y8,fs]=polyphase_filterI1('乐声.mp3',3); % figure(6); % subplot(2,1,1);pinpu(y8,fs);title('多相插值的频谱'); % subplot(2,1,2);shiyu(y8,fs);title('多相插值时域波形'); %% 多相 分数倍抽样的频域时域图像 tic; [y5,fs] = polyphase_filterID1('乐声.mp3',5/3); figure(7); subplot(2,1,1);pinpu(y5,fs);title('多相分数倍抽样的频谱'); subplot(2,1,2);shiyu(y5,fs);title('多相分数倍抽样的时域波形'); t3=toc; fprintf('分数倍抽样时间%d\n',t3); %% resample 分数倍 抽样的频域时域图像 [y6,fs]=resampleIDb('乐声.mp3',5/3); figure(8); subplot(2,1,1);pinpu(y6,fs);title('resample分数倍的频谱'); subplot(2,1,2);shiyu(y6,fs);title('resample分数倍时域波形'); %% resample 插值 时域频域图像 [y7,fs]=resampleIb('乐声.mp3',3); figure(9); subplot(2,1,1);pinpu(y7,fs);title('resample插值的频谱'); subplot(2,1,2);shiyu(y7,fs);title('resample插值时域波形'); %% 多相 插值 时域频域图像 [y8,fs]=polyphase_filterI1('乐声.mp3',3); figure(10); subplot(2,1,1);pinpu(y8,fs);title('多相插值的频谱'); subplot(2,1,2);shiyu(y8,fs);title('多相插值时域波形'); %% 计算信噪比 s1=power1(y1,0); s2=power1(y2,y1);SNR2=10*log10(s1/s2); s3=power1(y3,y1);SNR3=10*log10(s1/s3); s4=power1(y4,y1);SNR4=10*log10(s1/s4); s5=power1(y6,0); s6=power1(y5,y6);SNR5=10*log10(s5/s6); s7=power1(y7,0); s8=power1(y8,y7);SNR6=10*log10(s7/s8); d1=SNR3-SNR2; d2=abs(SNR3-SNR4); %% 显示信噪比 fprintf('直接抽取信噪比%d\n',SNR2); fprintf('直接滤波后抽取信噪比%d\n',SNR3); fprintf('多相滤波后抽取信噪比%d\n',SNR4); fprintf('卷积滤波与直接抽取信噪比之差%d\n',d1); fprintf('卷积滤波与多相滤波信噪比之差%d\n',d2); fprintf('分数倍的信噪比%d\n',SNR5); fprintf('插值的信噪比%d\n',SNR6); %% 计算乘法次数 N=600; D=30; c1=length(y)*N; c2=length(y)*N/D; bi=c1/c2