DSP数字信号处理习题集及matlab编程.doc

P1．已知序列绘制。 for i=1:5 x(i)=0.8^(i-1); end h=[1 1 1 1 1]; y=conv(x,h) m=[0 1 2 3 4 5 6 7 8]; stem(m,y,filled) P2．已知复指数序列绘制0点该序列的实部和虚部。 所以 for n=1:20 Re(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)); end for i=1:20 x(i)=i-1; end subplot(2,1,1) stem(x,Re,filled); lab1=\rightarrowRe(x[n]);text(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); P3．编写长度为5的中值滤波器程序。原始未受干扰的序列为： s[n]=3[n(0.) ]n，加性噪声信号d[n]为随机序列，幅度0.，分别绘制长度为40的受干扰序列，以及中值滤波器的输出。for n=1:40 s(n)=3*(n-1)*0.5^(n-1); end a=(rand(1,40)-0.5)*4/5; for i=1:40 x(i)=s(i)+a(1,i); end for i=1:40 b(i)=i-1; end subplot(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,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; for i=5:40 y(i)=1/5*(x(i)+x(i-1)+x(i-2)+x(i-3)+x(i-4)); end subplot(2,2,4); stem(b,y,r--,filled); subplot(2,2,3); stem(b,x,b,filled); P4. 已知序列xn]={2.2,3,1.5,4.2,1.8}, x2[n]= {0.8,1,1.6,0.8}，x[n]=x1[n]?x2[n] （卷积），分别绘制序列xn] ,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); for i=1:5 m(i)=i-1; end for i=1:4 n(i)=i-1; end for i=1:8 p(i)=i-1; end subplot(3,1,1); stem(m,x1,filled); subplot(3,1,2); stem(n,x2,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的原始未受干扰的序列，噪声序列和受干扰序列，以及滑动平均滤波器的输出。for n=1:40 s(n)=3*(n-1)*0.8^(n-1); end a=(rand(1,40)-0.5)*6/5; for i=1:40 x(i)=s(i)+a(1,i); end for i=1:40 b(i)=i-1; end subplot(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,40); y(1)=x(1)/4;y(2)=(x(1)+x(2))/4;y(3)=(x(1)+x(2)+x(3))/4; for i=4:40 y(i)=1/4*(x(i)+x(i-1)+x(i-2)+x(i-3)); end subplot(2,2,4); stem(b,y,r--,filled); subplot(2,2,3); stem(b,x,b,filled); P6． clear for i=1:8 x(i)=0.8^(i-1)*cos(0.75*pi*(i-1)); end for j=1:8 y(j)=x(-j+9); en