欢迎来到天天文库
浏览记录
ID:34312843
大小:153.30 KB
页数:4页
时间:2019-03-05
《音乐频谱分析(18103317)》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、%y=wavrecord(110250,11025,'int16');%录音;%wavwrite(y,11025,'e:MatlabSPDMTboxingchulixy.wav');%将录制的音乐存盘;[x,Fs,bit,size]=wavread('e:MatlabSPDMTboxingchulixy.wav');%读取文件;wavplay(10*x,Fs);%滤波前音乐效果演示n=length(x);t=(0:n-1)/Fs;figure(1),plot(t,x);xlabel('时间T(s)');ylabel('幅值');title(
2、'音乐波形信号');gridon;%频谱分析y=fft(x);w=2/n*(0:n-1);figure(2),plot(w,abs(y));title('信号频谱');xlabel('数字角频率W');ylabel('幅度');gridon;%巴特沃斯低通滤波器,滤除噪声干扰llllllllllllllllllllllllllllllllllll[N,Wc]=buttord(0.06,0.10,1,15);%确定低通滤波器的阶数和截止频率[b,a]=butter(N,Wc);%确定低通滤波器分子分母系数[H,W]=freqz(b,a);%低通滤波器波形fi
3、gure(3),plot(W,abs(H));xlabel('W');ylabel('H');title('低通滤波器');m=filter(b,a,x);wavplay(10*m,Fs);%经过滤波器后音乐播放效果演示s=length(m);%滤波后信号长度tm=(0:s-1)/Fs;Wm=2/s*(0:s-1);M=fft(m);figure(4),plot(tm,m);gridon;title('滤波后波形');xlabel('T');ylabel('幅度');figure(5),plot(Wm,abs(m));title('滤波后波形频谱图');x
4、label('w/pi');ylabel('幅度');%滤波后频谱分析k=fft(m);z=length(k);w1=2*(0:z-1)/z;figure(6),plot(w1,abs(k));xlabel('数字角频率');ylabel('幅度');title('低通滤波后的频谱');图形
此文档下载收益归作者所有