《数字信号处理》课程设计报告--基于MATLAB的音乐信号处理和分析.doc

一、课程设计的目的 2 二、课程设计基本要求 2 三、课程设计内容 2 1、音乐信号的音谱和频谱观察 2 2、音乐信号的抽取（减抽样） 4 3、音乐信号的AM调制 5 4、AM调制音乐信号的同步解调 7 5、音乐信号的滤波去噪 14 四、问题讨论 19 1、IIR数字滤波器和FIR数字滤波器的比较 19 2、音乐信号的音调与信号的什么特征有关？ 20 3、音乐信号的音色与信号的什么特征有关？ 20 4、两种不同音色的音乐信号叠加混叠后，为何人耳还可以分辨？ 20 5、音乐信号的幅度与相位特征对信号有哪些影响？ 21 五、心得体会 21 六、参考文献 21 《数字信号处理》课程设计 设计题目：基于 MATLAB 的信号 学会 MATLAB 的使用， 掌握在 Windows 环境下信号采集的方法。 掌握数字信号处理的基本概念、基本理论和基本方法。 掌握 MATLAB 设计 FIR 和 R 数字滤波器的方法。 掌握MATLAB处理数字信号、进行频谱分析、设计数字滤波器的编程方法。 Wavread格式说明： [w,fs,b]=wavread(‘语音信号’)，采样值放在向量w中，fs表示采样频率（hz），b表示采样位数。 【Matlab程序如下】： close all; clear all; clc; [w,fs,b]=wavread(C:\Users\lenovo\Desktop\matlab设计\梦里水乡 [移动影音转换专家].wav); w1=w(:,2); %抽样率为fs subplot(2,2,1); plot(w1);title(音乐信号时域图);grid on; %sound(w1,fs); f1=fft(w1); l1=length(f1); ww1=2/l1*(0:l1-1); subplot(2,2,3); plot(ww1,abs(f1));title(音乐信号频谱图);grid on; %抽样率为fs/2 sound(w1,fs/2); %抽样率为2fs sound(w1,2*fs); 【程序运行结果如下图】： 分析：通过观察音乐信号的波形和频谱可知所选取的音乐信号频谱集中在0~0.5pi之间，同时抽样点数fs=44000。 2、音乐信号的抽取（减抽样） ① 观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔，代表混叠与非混叠）； ② 输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释； ③ 播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音，解释现象。 理论基础： 时域抽样定理：一个频谱受限的信号f(t)，如果频谱只占据-wm~+wm的范围， 则信号f(t)可以用等间隔的抽样值唯一的表示。而抽样间隔必须不大于1/（2*fm）. 频域抽样定理：一个频谱受限的信号f(t)，它集中在-tm~+tm的时间范围内， 若在频域中以不大于1/（2*tm）的频率间隔对f(t)的频谱F(w)进行抽样，则抽样后 频谱F1(w)可以唯一的表示原信号。 【Matlab程序如下】： %减抽样小间隔 j=0; d=2; wz=w1(1:10000); for i=1:d:length(wz) j=j+1; dwav(j)=w1(i); end %sound(dwav,fs/d); subplot(2,2,2); plot(dwav);title(减抽样时域图);grid on; f2=fft(dwav); l2=length(f2); ww2=2/l2*(0:l2-1); subplot(2,2,4); plot(ww2,abs(f2));title(减抽样频谱图);grid on; %减样大间隔 j=0;d=4;wz=w1(1:10000);for i=1:d:length(wz)j=j+1;dwav(j)=w1(i);end %sound(dwav,fs/d);figure;subplot(2,1,1);plot(dwav);f2=fft(dwav);l2=length(f2);ww2=2/l2*(0:l2-1);subplot(2,1,2);plot(ww2,abs(f2)); 【程序运行结果如下图】： 分析：通过观察两种不同抽样间隔（2倍频和四倍频）下的音乐信号可知，当采用较大的抽样间隔对音乐信号进行抽样时，频谱发生了混叠，而采用较小的抽样间隔对音乐信号进行抽样时，频谱并未发生混叠。这是因为，抽样时频谱发生混叠的条件是fs2fh，即抽样频率小于信号频谱的最高频率。当采用较大的抽样间隔时抽样频率时fs2fh，所以发生混叠，而采用较小的抽样间隔时抽样频率时fs=2fh，则不会发生混叠。 当我们播放不同抽样间隔下的音乐信号时，会发现大抽样的音乐信