实验七多采样率数字信号处理

《实验七多采样率数字信号处理》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、⑻⑹I-3000200010000原始信号频域阁0200040006000⑼抽取后的信号频域图实验七多样采样率数字信号处理7.1实验目的1.掌握信号抽取和插值的基本原理和实现：2.掌握信号的有理数倍率转换。7.2实验原理例7.1对信号进行抽取，使采样频率为原来的1/4倍解：t=0:.00025:1;x二sin(2*pi*30*t)+sin(2*pi*60*t);y=decimate(x,4);figure;subplot(2,2,1),stem(x(l:120));titleC原始信号时域图’)；xlabelC(a)’)；subplot(2,2,2);plot(abs

2、(fft(x)));titleC原始信号频域图’)；xlabel(’(b)’);subplot(2,2,3);stem(y(1:30));titleC抽取后的信号时域图’);xlabelC(c)’)；subplot(2,2,4);plot(abs(fft(y)));titleC抽取后的信号频域图’)；xlabelC(d)’)；6004002000050010001500(d)说明：decimate用于信号整数倍抽取，其调用形式为pdecimate(x,M)，其是将x的采样率降低为原來的1/M倍，而抽取前默认采川8阶ChebyshenT型低通滤波器压缩频带：plot(a

3、bs(fft(x))):指画出x的幅度谱；figure：指产生一个空白的图形窗口。例7.2信号，采样频率f/fs=l/16,现将采样频率提高为原来的四倍。解：n二0:30;x=cos(3*pi*n/16);y=interp(x,4);figure;subplot(2,2,1),stem(x);titleC原始信号吋域图’)；xlabel(’(a)’)；subplot(2,2,2);plot(abs(fft(x)));title(’原始信号频域图’);xlabel(’(b)’);subplot(2,2,3);stem(y);title(’插值后的信号时域图’);xlab

4、elC(c)>);subplot(2,2,4);plot(abs(fft(y)));title(’插值后的信号频域图’)；xlabelC(d)’);原始信号时域图(b)插值后的信号频域图10-1-2⑻插值后的信号时域图⑹80050100150604020ok」050100150(d)说明：interp用于信号整数倍插值，其调用格式为y=interp(x,L)，是将信号采样率提高L倍。figure：指产生一个空白的图形窗口。例7.3序列x(n)=sin(0.5pin)+cos(l.2pin)+l.2,调用resample函数对x(n)按因子3/7进行采样率转换，并绘出图

5、形。解：n=0:34;x=sin(0.5*pi*n)+cos(1.2*pi*n)+l.2;[y,b]二resample(x,3,7);figure;subplot(2,2,1),stem(x);titleC原始信号吋域图’)；xlabel(’U)’)；subplot(2,2,2);plot(b,’•’)；title(’滤波器的单位脉冲响应’)；xlabel(’(b)’)；subplot(2,2,3);stem(y);title('变换后的信号时域图’)；xlabelC(c)’)；w=(0:1023)*2/1024;subplot(2,2,4);plot(w,20*lo

6、gl0(abs(fft(b,1024))));axis([0,1/4,-100,20]);gridontitleC滤波器的频率响应’)；xlabel(’(d)’)；4原始信号吋域图010203040滤波器的单位脉冲响应⑼⑻变换后的信号时域图⑹-5010°00.050.10.150.2说明：resample用于有理倍数率转换，其调用格式为y=resample(x,L,M)，是将信号x的采样率转换为原来的L/M倍。lot(w，20*logl0(abs(fft(b,1024))))是求滤波器的单位脉冲响应函数7.3实验内容1.令x(n)=cos(2nf/fs),f/fs=l

7、/12,实现以下采样率的转换:(1)作＞2倍的插值；(2)作M=3倍的抽取；(3)作L/M-2/3倍的采样率转换。解(1)n=0:30;x=cos(2*pi*n/12);y=interp(x,3);figure;subplot(2,2,1),stem(x);title('原始信号时域图’);xlabelC(a)’)；subplot(2,2,2):plot(abs(fft(x)));title(’原始信号频域图’);xlabel(’(b)’);subplot(2,2,3);stem(y);title(’插值后的信号时域图’);xlabelC(c)>);