数字信号处理课程设计（对音乐信号的各种处理)

《数字信号处理课程设计（对音乐信号的各种处理)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、WORD格式整理实验11、音乐信号的音谱和频谱观察使用wavread语句读取音乐信号，获取抽样率；输出音乐信号的波形和频谱，观察现象；使用sound语句播放音乐信号，注意不同抽样率下的音调变化，解释现象。clearall;closeall;clc[a,fs,bit]=wavread('c:MATLAB6p5work陪你一起看草原.wav');size(a);y1=a(:,1);a1=y1(10000:60000)figure;subplot(2,1,1),plot(a);subplot(2,1,2

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('时域波形');%

3、时域波形;学习参考资料分享WORD格式整理学习参考资料分享WORD格式整理1，以二倍的抽样率听声音信号时，音乐播放的特别快，像被压缩了，播放的时间比原信号短。2，以二分之一的抽样率听声音信号时，音乐播放的特别慢，像被拉长了，播放的时间比原信号长。3，原信号频谱截止频率为0.5*pi实验21、音乐信号的抽取（减抽样）观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔，代表混叠和非混叠）；输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释；播放减抽样音乐信号，注意抽样率的改变，比

4、较不同抽取间隔下的声音，解释现象。clearall;closeall;clc[a,fs,bit]=wavread('c:MATLAB6p5work陪你一起看草原.wav');size(a);y1=a(:,1);a1=y1(10000:60000)D=2;%减抽样；l=length(a1);学习参考资料分享WORD格式整理yd=a1(1:D:l);sound(yd,fs/D);N3=length(yd);t=0:1/N3:1/N3*(N3-1);%横坐标设置figure;plot(yd);title

5、('减抽样时域波形');%时域波形;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,减抽样学习参考资料分享WORD格式整理D=4,减抽样学习参考资料分享WORD格式整理1，原信号频谱截止频率为0.5*pi，当D=2时，频谱刚好不混叠，当D>2时，频谱就会混叠。2，减抽样后

6、的音乐信号听起来变得尖锐，有失真。3，抽样率随着抽样间隔的增大而逐渐变小，声音越来越失真，音调变得急促，而尖锐，信号产生混叠实验31、音乐信号的AM调制观察音乐信号的频率上限，选择适当调制频率对信号进行调制（给出高、低两种调制频率）；输出调制信号的波形和频谱，观察现象，给出理论解释；播放调制音乐信号，注意不同调制频率下的声音，解释现象。clearall;closeall;clc[a,fs,bit]=wavread('c:MATLAB6p5work陪你一起看草原.wav');size(a);y1=a

7、(:,1);a1=y1(10000:100000);%sound(a1,fs);N=length(a1);n1=0:[N-1];学习参考资料分享WORD格式整理y=cos(0.5*pi*n1);%产生余弦信号N1=length(a1);F1=fft(y,N1);w=2/N1*[0:N1-1];figure;plot(w,abs(F1));title('cos(wt)频谱');%余弦频谱xlabel('f');ylabel('幅度');N2=length(a1);F2=fft(a1,N2);w=2/N2*

8、[0:N2-1];figure;plot(w,abs(F2));title('yinyue频谱');%原信号频谱xlabel('f');ylabel('幅度');x1=a1.*y';%音乐信号与余弦信号点乘进行调制figureplot(x1);N3=length(a1);F3=fft(x1,N3);%调制信号傅立叶变换w=2/N3*[0:N3-1];figure;plot(w,abs(F3));title('调制频谱');%调制信号频谱xl