数字信号处理附录pdf.docVIP

实验一 MATLAB 数学软件的使用 1、绘制单位脉冲序列 clf; n=-10:20; u=[zeros(1,10) 1 zeros(1,20)]; stem(n,u); xlabel(时间序号 n); ylabel(振幅); title(单位脉冲序列); axis([-10 20 0 1.2]); 2、绘制长度为 20 的单位阶跃序列 clf; n=1:20; u=[ones(1,20)]; stem(n,u); xlabel(时间序号 n); ylabel(振幅); title(单位阶跃序列); axis([-10 20 0 1.2]); 3、产生一个正弦序列 clf; n=0:40; f=0.1; phase=0; A=1.5; arg=2*pi*f*n-phase; x=A*cos(arg); stem(n,x); axis([0 40 -2 2]); grid; title(正弦序列); xlabel(时间序号 n);ylabel(振幅); 4、绘制复指数序列 clf; c=-(1/12)+(pi/6)*i; k=2; n=0:40; x=k*exp(c*n); subplot(2,1,1); stem(n,real(x)); xlabel(时间序号 n);ylabel(振幅);  title(实部); subplot(2,1,2); stem(n,imag(x)); xlabel(时间序号 n);ylabel(振幅); title(虚部); 实验二 数字信号处理中常见算法的 MATLAB 实现 1、序列移位与周期延拓运算 N=24;M=8;m=3; n=0:N-1; x1=(0.8).^n; %生成指数序列 x2=[(n=0)(nM)]; %形成指数序列 x=x1.*x2; %截取操作形成新序列 x(n) xm=zeros(1,N); for k=m+1:m+M xm(k)=x(k-m); %产生序列移位 x(n-3) end; xc=x(mod(n,M)+1); %产生 x(n)的周期延拓 x((n))8 xcm=x(mod(n-m,M)+1); %产生移位序列 x(n-3)的周期延拓 x((n-3))8 subplot(4,1,1),stem(n,x,.);ylabel(x(n)); subplot(4,1,2),stem(n,xm,.);ylabel(x(n-3)); subplot(4,1,3),stem(n,xc,.);ylabel(x(n)_8); subplot(4,1,4),stem(n,xcm,.);ylabel(x(n-3)_8); 2、序列翻褶与叠加 n=0:N-1; %x(n)序列的时间序列 x=3*exp(-0.2*n); %x(n)序列大小 y=fliplr(x); %x(n)序列翻褶 n1=-fliplr(n); %时间序列的翻褶（翻褶点为原点） n2=fliplr(-(n-3)); %在指定位置 m=3 处的时间序列的翻褶 subplot(2,1,1);stem(n,x); %绘制 x(n)序列 xlabel(n),ylabel(x(n)); subplot(2,1,2);stem(n1,y); %绘制 y(n)=x(-n)序列 xlabel(n),ylabel(y(n)=x(-n)); s=cumsum(x); %求累加序列 figure(2); subplot(2,1,1);stem(n2,y); %绘制点为 3 处的翻褶序列 y(n)=x(-n+3) xlabel(n),ylabel(y(n)=x(-n+3)); subplot(2,1,2);stem(n,s); %绘制 s(n)序列 xlabel(n),ylabel(s(n));  2、 序列的卷积与相关运算 （1） 卷积程序： Nx=20;Nh=10;m=5; %设定 Nx,Nh 和位移值 m n=0:Nx-1; x1=(0.9).^n; %产生 x1(n) x2=zeros(1,Nx+m); for k=m+1:m+Nx %产生 x2(n)=x1(n-m) x2(k)=x1(k-m); %完成位移 end nh=0:Nh-1;h1=ones(1,Nh); %产生 h1(n) h2=h1; y1=conv(x1,h1); %计算 y1(n)=x1(n)*h1(n) subplot(3,2,1);stem(n,x1); %绘制 x1(n) xlabel(n),ylabel(x1(n)); subplot(3,2,3);stem(nh,h1); %绘制 h1(n) xlabel(n),ylabel(h1(n)); subplot(3,2,5);stem(y1);