基于MATLAB的音乐信号频谱分析和滤波.docxVIP

基于MATLAB的音乐信号频谱分析和滤波摘要：从20世纪60年代以来，随着信息学科和计算机学科的高速发展而发展起来的数字信号处理，在各个领域显示出重要的作用。本课程设计通过对音乐信号的采样、抽取、调制、解调等多种处理过程的理论分析和MATLAB实现，进一步强化对于《数字信号处理》课程的认识与理解。关键字：数字信号处理 MATLAB 音乐信号Music signal processing and analysis based on MATLABJi QiupingCollege of Physics and Electronic Engineering InformationWenzhou universityAbstract:Since the 1960s, with the rapid development of information science and computer science ,digital signal is developing quickly and showing an important role in various fields.This course is designed for music signal by sampling, theoretical analysis extraction, modulation, demodulation and other processing and MATLAB to achieve,for further strengthen the course Digital Signal Processing knowledge and understanding.Keywords:digital signal processing MATLAB music signal一、引言1.1 课程设计的目的理论与实际相结合，进一步扩展《数字信号处理》课程中所学的基本理论知识和分析方法，增强软件编程能力，熟悉 matlab 的操作界面，提高解决实际问题的能力。1.2 课程设计的基本要求学会 matlab的使用，掌握基本的编程语句掌握在windows下音乐信号采集的方法掌握数字信号处理的基本概念，基本理论和方法掌握使用matlab 处理数字信号、进行频谱分析、设计数字滤波器的编程方法二、课程设计内容2.1 音乐信号的音谱和频谱观察 截取一段音乐信号（要求：时间不超过5s，文件格式为wav文件），在matlab中用y=wavread(‘YY.wav’)命令读入原始音乐文件，fs为采样频率。【matlab程序如下：】close all;clear all;clc;[w,fs]=audioread(D:\FFOutput\音乐.wav);%读取一段音乐信号w1=w(:,2);%抽样率为fssubplot(2,2,1);%分块图plot(w1);title(音乐信号时域图);grid on;%sound(w1,fs);显示网格f1=fft(w1);%输入向量函数I1=length(f1);%向量长度；ww1=2/I1*(0:I1-1);%抽样率为fs/2 sound(w1,fs/2);subplot(2,2,3);%抽样率为2fs sound(w1,2*fs);plot(ww1,abs(f1));title(音乐信号频谱图);grid on;【程序运行结果如下图】分析：通过观察音乐信号的波形和频谱可知，所选取的音乐信号集中在0~0.2pi之间.2音乐信号的抽取（减抽样）观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音 理论基础：时域抽样定理：一个频谱受限的信号，如果频谱只占据-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;%间隔为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);I2=length(f2);ww2=2/I2*(0:I2-1);subplot(2,2,4);p