Question

matlab对语音信号采样频率怎么改变

Answer 1

（1）做波形显示以及fft变换，程序如下:

[y,fs]=wavread('E:\MATLAB6p5\work\3.wav');%读出信号，采样率。
y=y(:,1);%取单声道。
sigLength=length(y);
Y = fft(y,sigLength);
Pyy = Y.* conj(Y) / sigLength;
halflength=floor(sigLength/2);
f=Fs*(0:halflength)/sigLength;
figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)'); %画频域波形

t=(0:sigLength-1)/Fs;
figure;plot(t,y);xlabel('Time(s)'); %画时域波形

(2)关于滤波
声音频率主要集中在0~1KHZ，我想虑掉500hz以下的频率，因此采用一个高通滤波器
这里我使用了一个10阶butterworth高通滤波器，边带是500hz，但是这不能直接用，因为声音文件的采样率是44k，500hz相对于44k来说太小了。所以我得先把我的声音欠采样，然后再滤波。程序如下:
[k,Fs,bits]=wavread('E:\MATLAB6p5\work\3.wav');
k=k(:,1);
y_temp=k(1:90000);

dfactor=3;
y=decimate(y_temp,dfactor);
[b,a] = butter(10,500/(Fs/(dfactor*2)),'high');
y=filter(b,a,y);
y=interp(y,dfactor);

sigLength=length(y);
Y = fft(y,sigLength);
Pyy = Y.* conj(Y) / sigLength;
halflength=floor(sigLength/2);
f=Fs*(0:halflength)/sigLength;
figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)'); %滤波后的频域波形

t=(0:sigLength-1)/Fs;
figure;plot(t,y,t,y_temp);xlabel('Time(s)'); %滤波后的时域波形

wavwrite(y,Fs,'3_'); %保存处理后的声音文件，文件名为”3_”.