实验1
1、音乐信号的音谱和频谱观察
○
1使用wavread语句读取音乐信号，获取抽样率； ○
2输出音乐信号的波形和频谱，观察现象； ○
3使用sound语句播放音乐信号，注意不同抽样率下的音调变化，解释现象。 clear all;close all;clc
[a,fs,bit]=wavread('c:\\MATLAB6p5\\work\\陪你一起看草原.wav'); size(a); y1=a(:,1);
a1=y1(10000:60000) figure;
subplot(2,1,1),plot(a); subplot(2,1,2),plot(a1);
x1=resample(a1,2,1); %y=resample(x,p,q)返回量的长度是向量x的p/q倍 sound(x1,fs); %sound(a,fs); N1=length(a1); F1=fft(a1,N1);
w=2/N1*[0:N1-1]; %频谱图横坐标设置 figure;
plot(w,abs(F1)); N2=length(a1);
t=0:1/N2:1/N2*(N2-1);
title('傅利叶变换'); %傅利叶变换; figure; plot(a1);
title('时域波形'); %时域波形;
1
2
1， 以二倍的抽样率听声音信号时，音乐播放的特别快，像被压缩了，播放的时间比原信号短。
2， 以二分之一的抽样率听声音信号时，音乐播放的特别慢，像被拉长了，播放的时间比原信号长。
3， 原信号频谱截止频率为0.5*pi
2、音乐信号的抽取（减抽样）
1观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔，○
代表混叠和非混叠）；
2输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释； ○
3播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音，解释现象。 ○
clear all;close all;clc
[a,fs,bit]=wavread('c:\\MATLAB6p5\\work\\陪你一起看草原.wav'); size(a); y1=a(:,1);
a1=y1(10000:60000)
3
实验2
D=2; %减抽样； l=length(a1); yd=a1(1:D:l); sound(yd,fs/D);
N3=length(yd);
t=0:1/N3:1/N3*(N3-1); %横坐标设置 figure; plot(yd);
title('减抽样时域波形'); %时域波形; xlabel('t'); ylabel('幅度'); N4=length(yd); F2=fft(yd,N4); w=2/N4*[0:N4-1]; figure;
plot(w,abs(F2));
title('减抽样频谱'); %减抽样频谱 xlabel('f'); ylabel('幅度');
D=2,减抽样
4
D=4,减抽样
5