语音怎么进入滤波器matlab,基于matlab的语音信号滤波器的设计与实现 -尊龙游戏旗舰厅官网

(4)双线性变换法设计butterworth滤波器

程序如下：

fs=22050;

x1=wavread('wo.wav');

t=0:1/22050:(size(x1)-1)/22050;

au=0.03;

d=[au*cos(2*pi*5000*t)]';

x2=x1 d;

wp=0.25*pi;

ws=0.3*pi;

rp=1;

rs=15;

fs=22050;

ts=1/fs;

wp1=2/ts*tan(wp/2);             %将模拟指标转换成数字指标

ws1=2/ts*tan(ws/2);

[n,wn]=buttord(wp1,ws1,rp,rs,'s');%选择滤波器的最小阶数

[z,p,k]=buttap(n);                 %创建butterworth模拟滤波器

[bap,aap]=zp2tf(z,p,k);

[b,a]=lp2lp(bap,aap,wn);

[bz,az]=bilinear(b,a,fs);          %用双线性变换法实现模拟滤波器到数字滤波器的转换

[h,w]=freqz(bz,az);                %绘制频率响应曲线

figure(1)

plot(w*fs/(2*pi),abs(h))

grid

xlabel('频率／hz')

ylabel('频率响应幅度')

title('butterworth')

f1=filter(bz,az,x2);

figure(2)

subplot(2,1,1)

plot(t,x2)                %画出滤波前的时域图

title('滤波前的时域波形');

subplot(2,1,2)

plot(t,f1);               %画出滤波后的时域图

title('滤波后的时域波形');

sound(f1,22050);        %播放滤波后的信号

f0=fft(f1,1024);

f=fs*(0:511)/1024;

figure(3)

y2=fft(x2,1024);

subplot(2,1,1);

plot(f,abs(y2(1:512)));          %画出滤波前的频谱图

title('滤波前的频谱')

xlabel('hz');

ylabel('fuzhi');

subplot(2,1,2)

f1=plot(f,abs(f0(1:512)));     %画出滤波后的频谱图

title('滤波后的频谱')

xlabel('hz');

ylabel('fuzhi');

波形如下：

(5)窗函数法设计滤波器：

程序如下:

fs=22050;

x1=wavread('wo.wav');

t=0:1/22050:(size(x1)-1)/22050;

au=0.03;

d=[au*cos(2*pi*5000*t)]';

x2=x1 d;

wp=0.25*pi;

ws=0.3*pi;

wdelta=ws-wp;

n=ceil(6.6*pi/wdelta);                  %取整

wn=(0.2 0.3)*pi/2;

b=fir1(n,wn/pi,hamming(n 1));           %选择窗函数，并归一化截止频率

figure(1)

freqz(b,1,512)

f2=filter(bz,az,x2)

figure(2)

subplot(2,1,1)

plot(t,x2)

title('滤波前的时域波形');

subplot(2,1,2)

plot(t,f2);

title('滤波后的时域波形');

sound(f2,22050);               %播放滤波后的语音信号

f0=fft(f2,1024);

f=fs*(0:511)/1024;

figure(3)

y2=fft(x2,1024);

subplot(2,1,1);

plot(f,abs(y2(1:512)));

title('滤波前的频谱')

xlabel('hz');

ylabel('fuzhi');

subplot(2,1,2)

f2=plot(f,abs(f0(1:512)));

title('滤波后的频谱')

xlabel('hz');

ylabel('fuzhi');

波形如下：

(6)回放语音信号

在matlab中，函数sound可以对声音进行回放。其调用格式：sound(x，fs，bits)，x为要播放的音频信号，fs采样频率，bits采样位。在运行matlab程序时，可以听到回放的声音，经过比较，明显感觉滤波前后的声音有变化。

3.结语

遇到的问题：出现错误提示：

??? error using ==> wavread

error using ==> wavread

data compression format (ima adpcm) is not supported.

解决办法：wave格式的音频分为pcm和ima adpcm两种格式，matlab中用waveread函数做音频处理时，只能对pcm格式的.wav音频进行处理，因此需要将要处理的格式事先转换成pcm格式的.wav音频。可采用adensoft audio mp3 converter音频转换工具进行转换。

11

评分

威望 2

5

积分 2

收起

理由

2

5

2

您的付出是论坛的东西,感谢您一直支持!.

总结

以上是尊龙游戏旗舰厅官网为你收集整理的语音怎么进入滤波器matlab,基于matlab的语音信号滤波器的设计与实现的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得尊龙游戏旗舰厅官网网站内容还不错，欢迎将尊龙游戏旗舰厅官网推荐给好友。