语音信号处理实验报告集锦
【实验一】
一、实验题目
Short time analysis(短时间分析)
二、实验要求
Write a MA TLAB program to analyze a speech and simultaneously, on a single page, plot the following measurements:
1. the entire speech waveform
2. the short-time energy, En
3. the short-time magnitude, Mn
4. the short-time zero-crossing, Zn
5. the narrowband spectrogram
6. the wideband spectrogram
Use both the speech waveforms in the wznjdx_normal.wav. Choose appropriate window sizes, window shifts, and window for the analysis. Explain your choice of these parameters.
写一个程序来分析语音的同时，在一个单页，情节如下措施：
1。整个语音波形
2。短时间能量
3。短时幅度，锰
4。短时过零，锌
5。窄带语谱图
6。宽带语谱图
在wznjdx_normal.wav使用语音波形。选择合适的窗口大小，窗口移动和窗口进行分析。解释你对这些参数的选择。
三、实验程序
clear
[x,fs]=wavread('wznjdx_normal.wav');
n=length(x);
N=320;
subplot(4,1,1);plot(x);
h=linspace(1,1,N);
En=conv(h,x.*x);
subplot(4,1,2);plot(En);
Mn=conv(h,abs(x));
subplot(4,1,3);plot(Mn);
for i=1:n-1
if x(i)>=0 y(i)=1;
else y(i)=-1;
end
if x(i+1)>=0 y(i+1)=1;
else y(i+1)=-1;
end
w(i)=abs(y(i+1)-y(i));
end
k=1;
j=0;
while (k+N-1)Zm(k)=0;
for i=0:N-1
Zm(k)=Zm(k)+w(k+i);
end
j=j+1;
k=k+N/2;
end
for w=1:j
Q(w)=Zm(160*(w-1)+1)/(2*N);
end
subplot(4,1,4);plot(Q);grid;
figure(2);
subplot(2,1,1);spectrogram(x,h,256,200,0.0424*fs); subplot(2,1,2);spectrogram(x,h,256,200,0.0064*fs);
四、实验结果
语谱图：（Matlab 7.0 用不了spectrogram）
【实验二】
一、实验题目
Homomorphic analysis同态分析
二、实验要求
Write a MATLAB program to compute the real cepstrums of a section of voiced speech and unvoiced speech.
Plot the signal, the log magnitude spectrum, the real cepstrum, and the lowpass liftered log magnitude spectrum.
写一个程序来计算实际倒谱一段浊音和清音。
小区的信号，对数幅度谱，真正的倒谱，和低通liftered对数幅度谱。三、实验程序
nfft=256;
[x,fs] = wavread('wznjdx_normal.wav');
fx=x;
Xvm=log(abs(fft(fx,nfft)));
xhv=real(ifft(Xvm,nfft));
lifter=zeros(1,nfft);
lifter(1:30)=1;
lifter(nfft-28:nfft)=1;
fnlen=0.02*fs; % 20ms
win=hamming(fnlen);%加窗
n=fnlen;%窗宽度赋给循环自变量n noverlap=0.5*fnlen;
while(n<=length(x)-1)
fx=x(n-fnlen+1:n).*win;
n=n+noverlap;
end
xhvp=xhv.*lifter';
figure;
subplot(4,1,1)
plot(lifter);
title('倒谱滤波器');
subplot(4,1,2)
plot(x);
title('语音信号波形');
subplot(4,1,3)
plot(Xvm);
title('Xvm');
subplot(4,1,4)
plot(xhv);
title('xhv');
四、实验结果
