1、1实验 11、音乐信号的音谱和频谱观察使用 wavread 语句读取音乐信号，获取抽样率； 1输出音乐信号的波形和频谱，观察现象； 2使用 sound 语句播放音乐信号，注意不同抽样率下的音调变化，解释现象。 3clear all;close all;clca,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),plot(a1);x1=resample(a1,2,1); %y=resample(x,p
2、,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(时域波形 ); %时域波形;231，以二倍的抽样率听声音信号时，音乐播放的特别快，像被压缩了，播放的时间比原信号短。2，以二分之一的抽样率听声音信号时，音乐播放的特别慢，像被拉长了，播放的时间比原信号
3、长。3，原信号频谱截止频率为 0.5*pi实验 22、音乐信号的抽取（减抽样）观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔， 1代表混叠和非混叠） ；输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释； 2播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音，解释现象。 3clear all;close all;clca,fs,bit=wavread(c:MATLAB6p5work陪你一起看草原.wav);size(a);y1=a(:,1);a1=y1(10000:60000)D=2; %减抽样；l=length(a1);4yd=a1(1:D:l);
4、sound(yd,fs/D);N3=length(yd);t=0:1/N3:1/N3*(N3-1); %横坐标设置figure;plot(yd);title(减抽样时域波形); %时域波形;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,减抽样5D=4,减抽样61，原信号频谱截止频率为 0.5*pi，当 D=2 时，频谱刚好不混叠，当 D2 时，频谱就会混叠。2，减抽样后的音乐信号听起
5、来变得尖锐，有失真。3，抽样率随着抽样间隔的增大而逐渐变小，声音越来越失真，音调变得急促，而尖锐，信号产生混叠实验 33、音乐信号的 AM 调制观察音乐信号的频率上限，选择适当调制频率对信号进行调制（给出高、低两种调制频 1率） ；输出调制信号的波形和频谱，观察现象，给出理论解释； 2播放调制音乐信号，注意不同调制频率下的声音，解释现象。 3clear all;close all;clca,fs,bit=wavread(c:MATLAB6p5work陪你一起看草原.wav);size(a);y1=a(:,1);a1=y1(10000:100000);%sound(a1,fs);N=length
6、(a1);n1=0:N-1;7y=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*0:N2-1;figure;plot(w,abs(F2);title(yinyue 频谱); %原信号频谱xlabel(f);ylabel(幅度);x1=a1.*y; % 音乐信号与余弦信号点乘进行调制figureplot(x1);N3
7、=length(a1);F3=fft(x1,N3); %调制信号傅立叶变换w=2/N3*0:N3-1;figure;plot(w,abs(F3);title(调制频谱 ); %调制信号频谱xlabel(f);ylabel(幅度);sound(x1,fs);w=0.5*pi 时.8910w=0.3*pi 时.11w=0.7*pi 时.121，由原信号频谱知，信号截止频率约为 0.4pi，则产生混叠的阈值为 0.6pi。2，原信号的调制相当于频谱搬移，左移一个右移一个，当调制频率（余弦频率）小于 0.4pi 或大于 0.6pi 时就会产生混叠或丢失一部分信息。3，当余弦点数取得少时，余弦频谱会产生
8、泄漏。4，当调制频率较高时（发生混叠） ，声音响度低，几乎只能听见兹兹的声音，信号几乎完全失真，当调制频率较低时（未发生混叠） ，声音很尖锐，响度较大，稍微能听出一点调子，但也有兹兹的声音。实验 44、AM 调制音乐信号的同步解调设计巴特沃斯 IIR 滤波器完成同步解调；观察滤波器频率响应曲线； 1用窗函数设计 FIR 滤波器完成同步解调，观察滤波器频率响应曲线；（要求：分别使用 2矩形窗和布莱克曼窗，进行比较） ；输出解调音乐信号，比较不同滤波器下的声音，解释现象。 3clear all;close all;clcfunction hd=ideal(N,wc)for n=0:N-1if n=
9、(N-1)/2hd(n+1)=wc/pi;else hd(n+1)=sin(wc*(n-(N-1)/2)/(pi*(n-(N-1)/2);endenda,fs,bit=wavread(c:MATLAB6p5work陪你一起看草原.wav);size(a);y1=a(:,1);a1=y1(10000:100000);%sound(a1,fs); %播放 yuan 的信号N=length(a1);n1=0:N-1;y=cos(0.5*pi*n1);x1=a1.*y;%点乘figureplot(x1);N1=length(a1);F1=fft(x1,N1); %调制信号傅立叶变换w=2/N1*0:N
10、1-1;13figureplot(w,abs(F1);title(调制频谱 ); %调制信号频谱xlabel(f);ylabel(幅度);%sound(x1,fs); %播放调制后的信号x2=x1.*y;figureplot(x2); %解调后信号F2=fft(x2,N1); %解调信号傅立叶变换w=2/N1*0:N1-1;figureplot(w,abs(F2);title(解调频谱 ); %解调频谱xlabel(f);ylabel(幅度);%sound(x2,fs); %播放解调后的信号N,Wc=buttord(0.4,0.5,1,15);B,A=butter(N,Wc);H,W=freq
11、z(B,A);figureplot(W/pi,abs(H);title(数字巴特沃斯滤波器); %数字巴特沃斯滤波器x3=filter(B,A,x2); %滤波后信号w=2/N1*0:N1-1;fx=fft(x3,N1);fa=fft(a1,N1);figuresubplot(2,1,1),plot(w,abs(fa);title(yuan xin hao pin pu);subplot(2,1,2),plot(w,abs(fx);title(数字巴特沃斯滤波器滤波频谱); %滤波后频谱sound(x3,fs); %播放巴特沃斯滤波器滤波后信号 sheng yinN=33;wc=0.4*pi;
12、 hd=ideal(N,wc); w1=boxcar(N); %矩形窗w2=blackman(N); %布莱克曼窗14h1=hd.*w1;h2=hd.*w2;y3=conv(x2,h1); %解调后信号与矩形窗函数卷积y4=conv(x2,h2); %解调后信号与布莱克曼窗函数卷积%sound(y3,fs); %播放矩形窗滤波后信号 sheng yin%sound(y4,fs); %播放布莱克曼窗滤波后信号 sheng yin figure;subplot(2,1,1),plot(y3);title(矩形窗滤波后信号);subplot(2,1,2);plot(y3);title(布莱克曼窗滤波
13、后信号);fh1=fft(h1,N1);db1=-20*log10(abs(fh1(1)./(abs(fh1)+eps); %理想低通滤波器加窗后幅度响应fh2=fft(h2,N1);db2=-20*log10(abs(fh2(1)./(abs(fh2)+eps);w=2/N1*0:N1-1;figure;subplot(3,1,1),stem(h1); %矩形窗函数grid on;title(矩形窗 );xlabel(n);ylabel(h(n);subplot(3,1,2),plot(w,abs(fh1);grid on;title(矩形窗 );xlabel(w);ylabel(H(k);
14、subplot(3,1,3);plot(w,db1); %矩形窗函数分贝图grid on;figure;subplot(3,1,1);stem(h2);grid on;title(布莱克曼窗);xlabel(n);ylabel(h(n);subplot(3,1,2);plot(w,abs(fh2);grid on;title(布莱克曼窗);xlabel(w);ylabel(H(k);subplot(3,1,3);plot(w,db2); %布莱克曼窗窗函数分贝图15grid on;w=2/N1*0:N1-1;Fy1=fft(y3,N1);Fy2=fft(y4,N1);figure;subplo
15、t(2,1,1);plot(w,abs(Fy1);title(矩形窗滤波后频谱);subplot(2,1,2);plot(w,abs(Fy2);title(布莱克曼窗滤波后频谱);16171819201，解调后信号频谱在高频和低频处均有一部分，且成对称分布，需要滤掉高频才可大致还原原信号。2，原信号的截止频率为 0.4pi，使用数字巴特沃斯滤波器滤波器滤波参数通带截止频率 0.4pi，阻带开始频率 0.5pi，阻带衰减 15db。滤波效果很好，基本还原了原信号。3，使用窗函数滤波要根据过渡带宽算阶数 N，选截止频率为 0.4pi。4，使用矩形窗滤波，矩形窗过渡带窄，但是阻带有波纹，高频部分有小
16、部分未滤掉。5，使用布莱克曼窗滤波，布莱克曼窗过渡带宽，但是阻带较好。6，使用矩形窗和布莱克曼窗滤波，效果都行，基本都能还原原信号实验 55、音乐信号的滤波去噪给原始音乐信号叠加幅度为 0.05，频率为 3kHz、5kHz、8kHz 的三余弦混合噪声，观察 1噪声频谱以及加噪后音乐信号的音谱和频谱，并播放音乐，感受噪声对音乐信号的影响；给原始音乐信号叠加幅度为 0.5 的随机白噪声（可用 rand 语句产生） ，观察噪声频谱以 2及加噪后音乐信号的音谱和频谱，并播放音乐，感受噪声对音乐信号的影响；根据步骤 、 观察到的频谱，选择合适指标设计滤波器进行滤波去噪，观察去噪后 3 1 2信号音谱和频
17、谱，并播放音乐，解释现象。21clear all;close all;clca,fs,bit=wavread(c:MATLAB6p5work陪你一起看草原.wav);size(a);y1=a(:,1);a1=y1(10000:100000);%sound(a1,fs); %播放 yuan 的信号N1=length(a1);n1=0:N1-1;x1=0.05*cos(2*pi*3000*n1/fs);x2=0.05*cos(2*pi*5000*n1/fs);x3=0.05*cos(2*pi*8000*n1/fs);x4=x1+x2+x3;figure;plot(x4);%叠加余弦信号F1=fft
18、(x4,N1);%叠加余弦信号傅立叶变换w=2/N1*0:N1-1;figureplot(w,abs(F1);title(叠加余弦信号频谱);% 叠加余弦信号频谱xlabel(f);ylabel(幅度);y2=a1+x4;%sound(y2,fs); %播放叠加余弦音乐信号figure;plot(y2); %叠加余弦音乐信号F2=fft(y2,N1); %叠加余弦音乐信号傅立叶变换figureplot(w,abs(F2);title(叠加余弦噪声音乐信号频谱); %叠加余弦音乐信号频谱xlabel(f);ylabel(幅度);N,Wc=buttord(0.06,0.3,1,50);B,A=bu
19、tter(N,Wc);H,W=freqz(B,A);figureplot(W/pi,abs(H);title(数字巴特沃斯滤波器); %数字巴特沃斯滤波器22y3=filter(B,A,y2); %滤波后信号fa=fft(a1,N1);fy=fft(y3,N1);figuresubplot(3,1,1),plot(w,abs(fa);title(yuan xin hao pin pu);subplot(3,1,2),plot(w,abs(F2);title(叠加余弦噪声音乐信号频谱); %叠加音乐信号频谱subplot(3,1,3),plot(w,abs(fy);title(数字巴特沃斯滤波器
20、滤波频谱); %滤波后频谱%sound(y3,fs); %播放滤波后音乐信号r=rand(N1,1)-0.5;yr=a1+r;figure;plot(yr);%sound(yr,fs); %播放叠加随机噪声音乐信号N,Wc=buttord(0.05,0.4,1,50);B,A=butter(N,Wc);H,W=freqz(B,A);figure;plot(W/pi,abs(H);title(数字巴特沃斯滤波器 2); %数字巴特沃斯滤波器yf=filter(B,A,yr); %滤波后信号Fz=fft(r,N1);Fr=fft(yr,N1);Ff=fft(yf,N1);figure;subplo
21、t(2,1,1),plot(w,abs(fa);title(yuan xin hao pin pu);subplot(2,1,2),plot(w,abs(Fz);title(随机噪声频谱);figure;subplot(2,1,1),plot(w,abs(Fr);title(叠加随机噪声音乐信号频谱); %叠加音乐信号频谱subplot(2,1,2),plot(w,abs(Ff);title(数字巴特沃斯滤波器滤波频谱); %滤波后频谱sound(yf,fs); %播放滤波后音乐信号2324252627281，三余弦信号的频谱为不同频率处得三根线，加噪声后的音乐信号频谱是在原信号频谱上加了三条
22、不同频率的线。加噪声后音乐信号能听到原有的音调，但里面有非常大的杂音，兹兹的噪声。2，对加余弦噪声的信号进行滤波，用巴特沃斯滤波器滤波，参数为：通带截止频率 0.06pi，阻带开始频率 0.3pi，阻带衰减 50db，滤波后滤掉了一部分高频信息，三根余弦基本滤除。滤波后信号听起来还稍微有点杂音，有稍微的失真。3，对原信号加随机白噪声，白噪声均匀分布，对其用巴特沃斯滤波器滤波，参数为：通带截止频率 0.05pi，阻带开始频率 0.4pi，阻带衰减 50db，滤波后滤掉了一部分高频信息，阻带衰减快，到 0.4pi 就衰减到 0，但还保留了低频处得噪声，滤波后音乐信号听起来能听出原调，但有失真，且伴
23、随有较大声的随机白噪声。实验 66、音乐信号的幅频滤波及相频分析设计低通滤波器（可自行选取不同的截止频率） ，滤除原始音乐信号的高频信息，观察 1滤波前后的幅度频谱，并比较滤波前后的音乐效果，感受高频信息对音乐信号的影响；设计高通滤波器（可自行选取不同的截止频率） ，滤除原始音乐信号的低频信息，观察 2滤波前后的幅度频谱，并比较滤波前后的音乐效果，感受高频信息对音乐信号的影响；选取两端不同的音乐信号，分别将其幅度谱和相位谱交叉组合构成新的音乐信号，播放 3并比较组合后的音乐与原始音乐，感受相频信息对音乐信号的影响。clear all;close all;clca,fs,bit=wavread(
24、c:MATLAB6p5work陪你一起看草原.wav);size(a);y1=a(:,1);a1=y1(10000:100000);%sound(a1,fs); %播放原的信号 1N1=length(a1);F1=fft(a1,N1); %原信号 1 傅立叶变换w=2/N1*0:N1-1;29N,Wc=buttord(0.1,0.2,1,30);B,A=butter(N,Wc);H,W=freqz(B,A);figureplot(W/pi,abs(H);title(数字巴特沃斯滤波器); %数字巴特沃斯滤波器a2=filter(B,A,a1); %滤波后信号fd=fft(a2,N1);figu
25、re;subplot(211),plot(w,abs(F1);title(原信号 1 频谱); % 原信号 1 频谱xlabel(f);ylabel(幅度);subplot(212),plot(w,abs(fd);title(数字巴特沃斯滤波器滤波频谱); %滤波后频谱%sound(a2,fs); %播放滤波后音乐信号N,Wc=buttord(0.1,0.2,1,30);B,A=butter(N,Wc,high);H,W=freqz(B,A);figureplot(W/pi,abs(H);title(数字巴特沃斯滤波器); %数字高通巴特沃斯滤波器a3=filter(B,A,a1); %滤波后
26、信号fg=fft(a3,N1);figureplot(w,abs(fg);title(数字巴特沃斯滤波器滤波频谱); %滤波后频谱%sound(a3,fs); %播放滤波后音乐信号x,fs,bit=wavread(C:MATLAB6p5work刘欢 - 在路上.wav);size(x);x1=x(:,1);x2=x1(10000:100000);%sound(x2,fs); %播放原的信号 2figuresubplot(211),plot(x2)F2=fft(x2,N1); %原信号 2 傅立叶变换w=2/N1*0:N1-1;30subplot(212),plot(w,abs(F2);titl
27、e(原信号 2 频谱); % 原信号 2 频谱xlabel(f);ylabel(幅度);Fa1=abs(F1);Fa2=abs(F2);Fx1=angle(F1);Fx2=angle(F2);f1=Fa1.*exp(j*Fx2); %相位谱和幅度谱交叉组合f2=Fa2.*exp(j*Fx1);figuresubplot(211),plot(w,abs(f1);title(原信号 1 幅度谱，原信号 2 相位谱);subplot(212),plot(w,abs(f2);title(原信号 2 幅度谱，原信号 1 相位谱);b1=ifft(f1);b2=ifft(f2);%sound(real(b1),fs);%sound(real(b2),fs);
