数字音乐采样报告.docxVIP

音乐采样实验 姓名： 罗金桃学号：专业：集成电路设计与集成系统学院： 通信工程学院 实验背景： 声音信号的采集与分析处理在工程应用中是经常需要解决的题,如何实时采集声音信号并对其分析处理,找出声音信号的特征在科学研究中是一项非常有意义的工作。采样定理是信号处理中最重要的定理之一，通过采样定理，可以确定对原始信号的采样频率及采样周期，确保在对已采样信号进行恢复时不失真,又不浪费频带。 实验目的： 本实验旨在对采样定理进行初步验证，体会频谱混叠现象，并大致确定音频信号的最低采样频率。 实验内容： 本实验通过 MATLAB 软件，完成以下三项任务 1、用 fs=44100HZ 采集一段音乐 2、改变采样频率，用fs=5512HZ 采集一段音乐，体会混叠现象 3、录制一段自己的声音，试验当 fs=？时，发生混叠实验内容如下： 一、 用 44100HZ 的频率对音乐进行采样，用如下程序对该音乐进行 频率为 44100HZ 的采样： clc; [y,fs,nbits]=wavread(ww.wav); %采样频率 44100HZ NFFT=2^nextpow2(length(y)); fy=fft(y,NFFT)/length(y); %fft 变换， p=2*abs(fy(1:NFFT)); %abs 取幅度b=fs/2*linspace(0,1,NFFT); %b 定义一个数组，作为图像的横坐标plot(b,p); sound(y,44100); 图 1,采样前 二、用 5512HZ 的频率对音乐进行采样 用实验一的程序对同一段音乐进行采样，将采样频率改为 5512HZ ， 采样程序如下： clc; [y,fs,nbits]=wavread(ww.wav); y=y(1:1000000); y1=y(1:8:end); figure; NFFT=2^nextpow2(length(y1)); fy=fft(y1,NFFT)/length(y1); p=2*abs(fy(1:NFFT)); b=5512/2*linspace(0,1,NFFT); plot(b,p); wavwrite(y1,5512,ww_0.wav); %fft 变换， %abs 取幅度 %b 定义一个数组，作为图像的横坐标 sound(y1,5512); 图 2，5512HZ 采样 播放音乐，与原始音乐相比，用 5512HZ 频率采样后的信号高频处有较多的杂音。 人耳能听到的音频信号的频率介于 20——20000HZ 之间，根据采样定理，原始信号频率 fm 与采样频率 fs 之间如满足 fs=2fm,则对已采样信号进行恢复时不会产生失真，实验一中的采样频率44100HZ2*20000HZ，因此采样后恢复的信号不失真，几乎和原始信号一致；B 实验中采样频率 5512HZ 不满足与原始信号频率的 2 倍关系，因此采样后恢复的信号在高频处有失真。 三、录制一段自己的声音并对其进行采样，测出大致的采样最低频率 clc; [y,fs,nbits]=wavread(recording.wav); sound(y,fs); %以 44100HZ 采样播放 figure; NFFT=2^nextpow2(length(y)); fy=fft(y,NFFT)/length(y); p=2*abs(fy(1:NFFT)); b=fs/2*linspace(0,1,NFFT); plot(b,p); figure; y1=y(1:5:end); NFFT=2^nextpow2(length(y1)); fy=fft(y1,NFFT)/length(y1); p=2*abs(fy(1:NFFT)); b=8820/2*linspace(0,1,NFFT); plot(b,p); figure; y2=y(1:10:end); NFFT=2^nextpow2(length(y2)); fy=fft(y2,NFFT)/length(y2); %fft 变换， %abs 取幅度 %b 定义一个数组，作为图像的横坐标 %fft 变换， %abs 取幅度 %fft 变换， p=2*abs(fy(1:NFFT)); %abs 取幅度 b=4410/2*linspace(0,1,NFFT); plot(b,p); y3=y1(1,1:end); y4=y2(1,1:end); sound(y3,8820); %以 8820HZ 采样播放 sound(y4,4410); %以 4410HZ 采样播放 wavwrite(y1,8820,recording0.wav); %以 8820 采样存储wavwrite(y2,4410,recording1.wav); %以 4410 采