信号的波形分析与合成

武汉大学教学实验报告电子信息学院通信工程专业 2015年10月7日

实验名称 Wav信号的波形分析与合成指导教师

姓名年级学号成绩

3.时变音的短时频谱分析

代码如下：

load chirp

len=length(y);

fori=1:floor(len/1024)

seg=y((i1)*1024+1:i*1024);

figure('numbertitle','off','name',['frame',num2str(i)]); plot(linspace(-pi,pi,1024),fftshift(abs(fft(seg))));

end

实验结果如下：

4.分别录制男生和女生发元音“a”的声音，通过对音频文件的频谱分析对比两者的差异，并进行合理的解释。

⑴.通过以下程序录制男生女生声音

%录制声音并播放

title('FFT的幅频特性')

subplot(2,1,2)

plot(linspace(-pi,pi,2048),fftshift(abs(YR2048))) %FFT的幅频特性FFTSHIFT

title('FFT后幅频特性的fftshift')

%fs=1024

YR1024=fft(yr,1024);

figure('numbertitle','off','name','1024点FFT');

subplot(2,1,1)

plot(linspace(-pi,pi,1024),abs(YR1024)) %FFT的幅频特性

title('FFT的幅频特性')

subplot(2,1,2)

plot(linspace(-pi,pi,1024),fftshift(abs(YR1024))) %FFT的幅频特性的FFTSHIFT

title('FFT后幅频特性的fftshift')

图形如下：

①．1024点FFT

②．2048点FFT

%反变换，结果应与原始信号相同

yr1024=100*real(ifft(YR1024));

figure('numbertitle','off','name',' 1024点fft的反变换');

图形如下：

plot(linspace(-pi,pi,1024),abs(YR1024)) %FFT的幅频特性

title('FFT的幅频特性')

subplot(2,1,2)

plot(linspace(-pi,pi,1024),fftshift(abs(YR1024))) %FFT的幅频特性的FFTSHIFT

title('FFT后幅频特性的fftshift')

图形如下：

①．1024点FFT

②.2048点FFT

%反变换，结果应与原始信号相同

yr1024=100*real(ifft(YR1024));

figure('numbertitle','off','name',' 1024点fft的反变换'); subplot(2,1,1);

plot(yr1024)

title(' 1024点fft的反变换')

subplot(2,1,2);

plot(yr);

title('原信号');%重构寻找幅度最大的正弦分量

[maxpeak,peaki]=max(abs(YR1024(1:512)));

MAXSIN=zeros(1,1024);

MAXSIN(peaki)=maxpeak;

MAXSIN(1026-peaki)=maxpeak;

maxsin=100*ifft(MAXSIN);

figure('numbertitle','off','name',’用最大正弦分量重构信号’); subplot(2,1,1)

plot(maxsin);

title('重构的信号');

subplot(2,1,2);

plot(yr1024);

title('原始信号');%试听

wavplay(yr1024,fs);

wavplay(maxsin,fs);

图形如下：

5. 实验结论

⑴.女声比男声听上去高，因为女性的声带比较细，女声比男声的基频高，基频

决定了你听到的声音音调的高低，男声基频在50HZ到180HZ之间，女声基频在160HZ到380HZ之间。