资源描述：
《音乐信号的处理 matlab.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、1、音乐信号的音谱和频谱观察[y,fs,bit]=wavread('2');sound(y,fs);plot(y);y1=y(:,1);M=512;fh1=fft(y1,M);w=2/M*[0:M-1];figure;plot(w,abs(fh1));2、音乐信号的抽取（减抽样）[y,fs,bit]=wavread('2');sound(y,fs);figure;subplot(1,2,1);plot(y);title('原音乐信号波形');y1=y(:,1);M=512;fh1=fft(y1,M);w=2/M*[0:M-1];subplot(1,2,2);p
2、lot(w,fh1);title('原音乐信号频谱');gridon;j=0;d=10;%减抽样间隔取3和10代表混叠与非混叠y3=y1(1:100000);fori=1:d:length(y3);j=j+1;yd(j)=y1(i);endsound(yd,fs/d);figure;subplot(1,2,1);plot(yd);title('减抽样时域');fh2=fft(yd);l1=length(fh2);w2=2/l1*(0:l1-1);subplot(1,2,2);plot(w2,abs(fh2));title('减抽样频域');3、音乐信号的AM调
3、制和巴特沃斯IIR滤波器滤波[w,fs,bit]=wavread('2');w1=w(:,1);l=length(w1);ww=2/l*(0:l-1);n=0:length(w1)-1;t1=w1.*cos(0.6*pi*n)';T1=fft(t1);plot(t1);title('AM低/高频调制后波形');figure;plot(ww,abs(T1));title('AM低/高频调制后频谱');sound(t1,fs);t2=t1.*cos(0.6*pi*n)';T2=fft(t2);figure;plot(t2);title('AM低/高频解调后波形')
4、;figure;plot(ww,abs(T2));title('AM低/高频解调后频谱');sound(t2,fs);[N,wc]=buttord(0.2,0.3,2,15);[b1,a]=butter(N,wc);[h,wi]=freqz(b1,a);figure;plot(wi/pi,abs(h));title('iir滤波器频谱');%用iir对信号滤波y=filter(b1,a,t2);Y=fft(y);figure;plot(y);title('低/高频解调滤波后时域');figure;plot(ww,abs(Y));title('低/高频解调滤波后
5、的频域');sound(y,fs);用窗函数法设计FIR滤波器（矩形窗和布莱克曼窗）[w,fs,bit]=wavread('2');w1=w(:,1);n=0:length(w1)-1;t1=w1.*cos(0.6*pi*n)';j1=t1.*cos(0.6*pi*n)';%用FIR对信号同步解调N=55;wc=pi/4;%boxcarN取9hd=ideal(N,wc);w=blackman(N);%blackman和boxcarh=hd.*w';M=512;fh=fft(h,M);ww=2/M*(0:M-1);plot(ww,abs(fh));title('
6、矩形窗fir滤波器频率响应');y=conv(j1,h);Y=fft(y,M);figureplot(y);title('矩形窗fir滤波后时域');figureplot(ww,abs(Y));title('矩形窗fir滤波后频域');sound(y,fs);三余弦和白噪声的实验[y,fs,bit]=wavread('2');y1=y(:,1);f=fft(y1);l=length(f);ww=2/l*(0:l-1);n=0:length(y1)-1;x=0.05*(cos(2*pi*3000*n/fs)+cos(2*pi*5000*n/fs)+cos(2*p
7、i*8000*n/fs));fx1=fft(x);plot(ww,abs(fx1));title('三余弦噪声频谱');r=rand(size(y1),1)-0.5;fx2=fft(r);figureplot(ww,abs(fx2));title('随机白噪声频谱');x1=y1'+x;figureplot(x1);title('加三余弦噪声的信号音谱');fx3=fft(x1);figureplot(ww,abs(fx3));title('加三余弦噪声的信号频谱');x2=y1+r;figureplot(x2);title('加白噪声的信号音谱');fx4=
8、fft(x2);figureplot(
