东南大学《数字信号处理》实验报告二

《东南大学《数字信号处理》实验报告二》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验二一，实验目的1、二、实验内容1fs=8000t=-0.02:l/fs:0.02x=sin(50*2*pi*t)+sin(500*2*pi*t)+sin(1000*2*pi*t);plol(l,x);fs1=800t=-0.05:l/fsl:0.05x=sin(50*2*pi*t)+sin(500*2*pi*t)+sin(1000*2*pi*t);plot(t,x);3-1N=64;M=4;n=O:N-l;x=cos(0.125*pi*n);y=decimate(x,M,'fir');n=0:N/4-l;stem(n,y);0.5-0.5-110153-2N=48;L=3;n=0:

2、N-l;x=cos(0.125*pi*n);y=interp(x,L)；n=0:N*L-l;stem(n,y);100150-1.51L0504N=50;M=3;L=5;n=O:N-l;x=sin(0.04*2*pi*n)+sin(0.3*2*pi*n);yl=interp(x,L);n=0:N*L-l;y2=decimate(y1,M);n=O:N*L/M-l;stem(n,y2(1:N*L/M));51)[y,fs]=wavread('F:铃声.wav');sound(y,fs);m=fft(y);time=(l:length(y))/fs;figure;subplot(211)

3、;plot(time,y);titleC原始音频信号时域图J;subplot(212);plot(abs(m));titleC原始咅频信号频域图)；原始音频信号时域图原始音频信号频域團2)[y,fs]=wavread('F:铃声.wav');t=(l:2:length(y))/fs;x=decimate(y,2);sound(y,l1025/2);ml=fft(x);figure;subplot(211);plot(t,x)；titleC＜样率减一倍的时域图);subplot(212);plot(abs(ml));titleC釆样率减一倍的频域图);采样率减一倍的时域图采样率减一倍的

4、频域图2000III!!I1500-・1000--500-・020004000600080001000012000140003)[y,fs]=wavread('F:铃声.wav');t2=(0.25:0.25:length(y))/fs;b=resample(y,4,l);sound(y,44100);m3=fft(b);figure;subplot(211);plot(t2,b);titleC＜样率4倍的时域图?;subplot(212);plot(abs(m3));titleCM样率4倍的频域图J;采样率4倍的时域图150000123456789104x10采样率4倍的频域團10

5、000500004)[y,fs]=wavread('F:铃声.wav');tl=(0.4:0.4:length(y))/fs;c=resample(y,5,2);sound(y,32000);m2=fft(c);figure;subplot(211);plot(tl,c);titleC釆样率5/2倍的时域图匕subplot(212);plot(abs(m2));title。采样率5/2倍的频域图);采样率5/2倍的时域图8000采样率5/2倍的频域團6000400020002345°0674x10解释：信号经过采样率提高转换后，数据量为原来的倍数。在时域上，采样率提高的时域图与原来的

6、音频信号的时域图特性保持一样且排列更加紧密，在频域上，坐标点对应原始音频信号频谱坐标也提高了数倍。采样率提高的频谱应是原始音频信号的压缩。三、思考题