资源描述:
《dsp数字信号处理习题集及matlab编程》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、DSP数字信号处理习题集及matlab编程P1.已知两序列x[n]?0.8{u[n]?u[n?5]},h[n]?{1,1,1,1,1}计算两序列的卷积并绘制其波形。解:fori=1:5x(i)=0.8^(i-1);endh=[11111];y=conv(x,h)m=[012345678];stem(m,y,'filled')P2.已知复指数序列x[n]?1.2e(1.5?j2)n,绘制20点该序列的实部和虚部。解:x[n]?1.2e(1.5?j2)n?1.2e1.5nej2n?1.2e1.5n(cos2n?jsin
2、2n)?1.2e1.5ncos2n?j1.2e1.5nsin2n所以Re(x[n])?1.2e1.5ncos2n,Im(x[n])?1.2e1.5nsin2nforn=1:20Re(n)=1.2*exp(1.5*(n-1))*cos(2*(n-1));Im(n)=1.2*exp(1.5*(n-1))*sin(2*(n-1));endfori=1:20x(i)=i-1;endsubplot(2,1,1)stem(x,Re,'filled');lab1='rightarrowRe(x[n])';tex
3、t(14,1.2e+12,lab1,'Fontsize',18);subplot(2,1,2)stem(x,Im,'filled');lab2='rightarrowIm(x[n])';text(14,4e+11,lab2,'Fontsize',18);nP3.编写长度为5的中值滤波器程序。原始未受干扰的序列为:s[n]=3[n(0.5)]n,加性噪声信号d[n]为随机序列,幅度0.4,分别绘制长度为40的受干扰序列,以及中值滤波器的输出。解:forn=1:40s(
4、n)=3*(n-1)*0.5^(n-1);enda=(rand(1,40)-0.5)*4/5;fori=1:40x(i)=s(i)+a(1,i);endfori=1:40b(i)=i-1;endsubplot(2,2,1);stem(b,s,'r','filled');subplot(2,2,2);stem(b,a,'g','filled');subplot(2,2,3);stem(b,x,'b','filled');y=zeros(1
5、,40);y(1)=x(1)/5;y(2)=(x(1)+x(2))/5;y(3)=(x(1)+x(2)+x(3))/5;y(4)=(x(1)+x(2)+x(3)+x(4))/5;fori=5:40y(i)=1/5*(x(i)+x(i-1)+x(i-2)+x(i-3)+x(i-4));endsubplot(2,2,4);stem(b,y,'r--','filled');subplot(2,2,3);stem(b,x,'b','filled');P4.已知序列x1[n]=
6、{2.2,3,1.5,4.2,1.8},x2[n]={0.8,1,1.6,0.8},x[n]=x1[n]?x2[n]分别绘制序列x1[n],x2[n]和x[n]的波形。解:x1=[2.2,3,1.5,4.2,1.8];x2=[0.8,1,1.6,0.8];y=conv(x1,x2);fori=1:5m(i)=i-1;endfori=1:4n(i)=i-1;endfori=1:8p(i)=i-1;endsubplot(3,1,1);stem(m,x1,'filled');subplot(3,1,2);stem(n,x
7、2,'filled');subplot(3,1,3);(卷积),stem(p,y,'filled');P5.编写4点滑动平均滤波器程序。原始未受干扰的序列为:s[n]=3[n(0.8)]n,加性噪声信号d[n]为随机序列,幅度0.6,受干扰的序列为:x[n]=s[n]+d[n],分别绘制长度为40的原始未受干扰的序列,噪声序列和受干扰序列,以及滑动平均滤波器的输出。解:forn=1:40s(n)=3*(n-1)*0.8^(n-1);enda=(rand(1,40)-0.5)*6/5;fori=1:4
8、0x(i)=s(i)+a(1,i);endfori=1:40b(i)=i-1;endsubplot(2,2,1);stem(b,s,'r','filled');subplot(2,2,2);stem(b,a,'