读源文件.m

输入自己事先录制好的语音信号

[y,fs,bits]=wavread('d:\test.wav');  %将文件读进来
subplot(3,1,1);
plot(y,'LineWidth',2);
xlabel('时间/t')  ; ylabel('幅度/A');
title('原信号波形');   
spectrum=fft(y);                %对采样点进行fft变换
n2=length(spectrum);             %计算fft变换后向量的长度
t2=0:n2-1;   
subplot(3,1,2);
spreal=abs(spectrum)   ;        %计算fft变换后序列的实部
spimag=angle(spectrum)  ;         %计算fft变换后序列的虚部
plot(t2,spreal) ;
xlabel('k')                     %把x轴定义为w
title('原信号幅度')                   %标明此图形画的是幅度
subplot(3,1,3)
plot(t2,spimag)  ;  
xlabel('k')                     %把x轴定义为虚部
title('原信号相位')                   %把x轴定义为虚部