数字语音处理B第二次实验教学内容Word文档格式.docx
- 文档编号:13710590
- 上传时间:2022-10-13
- 格式:DOCX
- 页数:16
- 大小:475.24KB
数字语音处理B第二次实验教学内容Word文档格式.docx
《数字语音处理B第二次实验教学内容Word文档格式.docx》由会员分享,可在线阅读,更多相关《数字语音处理B第二次实验教学内容Word文档格式.docx(16页珍藏版)》请在冰豆网上搜索。
));
frameLength=fs*timeperFrame*0.001;
nframeLength=length(winfunc(frameLength));
%取窗长
overlap=frameLength*overlapRate;
inc=frameLength-overlap;
if(nframeLength==1)
len=frameLength;
else
len=nframeLength;
end
frameCount=fix((ndata-overlap)/inc);
%计算帧数
frameSet=zeros(len,frameCount);
startIndexf=(0:
(frameCount-1))*inc;
%每帧在数据data中开始的位置
startIndexs=(1:
len);
startIndexsT=startIndexs(:
);
frameSet(:
)=data(startIndexsT(:
ones(frameCount,1))+startIndexf(ones(1,len),:
w=winfunc(frameLength);
if(nframeLength>
1)
frameSet=frameSet.*w(:
ones(frameCount,1));
timeAxis=(((1:
frameCount)-1)*inc+frameLength/2)/fs;
%求出每帧对应的时间
主函数:
clc
clearall
closeall
%----------------------------------------
loadah;
x1=enframe(data,fs,20,0.5,'
boxcar'
x2=enframe(data,fs,20,0.5,'
hamming'
plot(x1(:
20),'
g'
'
LineWidth'
2);
holdon;
plot(x2(:
b'
legend('
矩形窗信号'
汉明窗信号'
xlabel('
时间(采样点数)'
ylabel('
幅度'
suptitle('
电信1405龙为花0121409341029'
)
图1:
语音ah第20帧时域波形图与加hammming窗10k采样
2.语音单帧频谱图
fs=10000;
u1=x1(:
20);
%取得一帧数据
u2=x2(:
N=1024;
t=(0:
N-1)/fs;
y1=20*log10(fft(u1,N));
y2=20*log10(fft(u2,N));
f=(0:
511)*fs/1024;
figure
(1);
Hz'
ylabel('
db'
plot(f,y1(1:
512),'
r'
1);
%做原始语音信号的FFT频谱
汉明窗的频谱'
figure
(2);
plot(f,y2(1:
矩形窗的频谱'
)
图2:
ah第20帧单帧语音频谱图,窗长20ms,帧移10ms,矩形窗和hamming窗,10k采样
3.语音频域功率谱图
y2=20*log10(fft(u1,N).^2);
1.5);
%做原始语音信号的频谱
%做原始语音信号的功率谱
汉明窗功率谱'
z1=20*log10(fft(u2,N));
z2=20*log10(fft(u2,N).^2);
plot(f,z1(1:
plot(f,z2(1:
矩形窗功率谱'
图3:
Ah第20帧时域波形图,10k采样,帧长20ms,hamming窗
图4:
Ah第20帧时域波形图,10k采样,帧长20ms,矩形窗
4.语谱图
(1)
clc
clear
closeall
[x,fs]=audioread('
E:
\1.wav'
m1=42230;
m2=138000;
x=x(m1:
m2);
win=0.01;
inc=win/2;
winlen=win.*fs;
inclen=inc.*fs;
amp1=enframe(filter([1-0.9375],1,x),hamming(winlen),inclen)'
;
n=fix((length(x)-winlen+inclen)/inclen);
w=winlen/2+1;
n2=1:
w;
freq=(n2-1)*fs/winlen;
Y=fft(amp1);
Yn=Y(n2,:
clf;
frameTime=(1:
n)*inclen;
imagesc(frameTime,freq,20*log10(abs(Yn)+eps));
axisxy;
时间/s)'
频率/Hz'
colormap(jet);
图5:
语音WHUT语谱图,10k采样
(2)
win=0.04;
时间/s'
图6
4.2
load('
should.mat'
amp1=enframe(filter([1-0.9375],1,data),hamming(winlen),inclen)'
n=fix((length(data)-winlen+inclen)/inclen);
图7
amp1=enframe(filter([1-0.9375],1,data),hamming(winlen),inclen)'
imagesc(frameTime,freq,20*log10(abs(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字 语音 处理 第二次 实验教学 内容