数字信号处理实验讲义(1-4).doc

数字信号处理实验 目 录 实验一 常见离散信号产生和实现 3 实验二 信号、系统及系统响应 10 实验三 应用FFT对信号进行频谱分析 23 实验四 离散系统的变换域分析 35 实验五 FIR数字滤波器的设计 41 实验六 IIR数字滤波器的设计 52 实验一 常见离散信号产生和实现 1.单位抽样序列 在MATLAB中可以利用zeros()函数实现。 如果在时间轴上延迟了k个单位，得到即： 2．单位阶越序列 在MATLAB中可以利用ones()函数实现。 3．正弦序列 在MATLAB中 4．复指数序列 在MATLAB中 5．指数序列 在MATLAB中 五、实验内容 1．五种基本函数的生成 程序如下： (1)单位抽样序列 clf; n=0:10; x1=[1 zeros(1,10)];x2=[zeros(1,5) 1 zeros(1,5)]; subplot(1,2,1); stem(n,x1);xlabel (时间序列n);ylabel(振幅);title(单位抽样序列x1); subplot(1,2,2); stem(n,x2); xlabel(时间序列n);ylabel(振幅);title(延时了5的单位抽样序列); （2）单位阶越序列 clf; n=0:10; u=[ones(1,11)]; stem(n,u);xlabel (时间序列n);ylabel(振幅);title(单位阶越序列); 所得的图形如下所示： （3）正弦函数 clf; n=1:30; x=2*sin(pi*n/6+pi/3); stem(n,x); xlabel (时间序列n);ylabel(振幅);title(正弦函数序列x=2*sin(pi*n/6+pi/3)); (4)复指数序列 clf; n=1:30; x=2*exp(j*3*n); stem(n,x); xlabel (时间序列n);ylabel(振幅);title(复指数序列x=2*exp(j*3*n)); 图形如下： （5）指数序列 clf; n=1:30; x=1.2.^n; stem(n,x); xlabel (时间序列n);ylabel(振幅);title(指数序列x=1.2.^n); 2．绘出信号，当、时、、、时的信号实部和虚部图； clf; z1=-1/12+j*pi/6;z2=1/12+j*pi/6;z3=1/12;z4=2+j*pi/6;z5=j*pi/6; n=0:20; x1=exp(z1*n);x2=exp(z2*n); x3=exp(z3*n);x4=exp(z4*n); x5=exp(z5*n); subplot(5,2,1); stem(n,real(x1)); xlabel (时间序列n);ylabel(实部);title(复指数z1=-1/12+j*pi/6时序列实部); subplot(5,2,2); stem(n,imag(x1)); xlabel (时间序列n);ylabel(虚部);title(复指数z1=-1/12+j*pi/6时序列虚部); subplot(5,2,3); stem(n,real(x2)); xlabel (时间序列n);ylabel(实部);title(复指数z2=1/12+j*pi/6时序列实部); subplot(5,2,4); stem(n,imag(x2)); xlabel (时间序列n);ylabel(虚部);title(复指数z2=1/12+j*pi/6时序列虚部); subplot(5,2,5); stem(n,real(x3)); xlabel (时间序列n);ylabel(实部);title(复指数z3=1/12时序列实部); subplot(5,2,6); stem(n,imag(x3)); xlabel (时间序列n);ylabel(虚部);title(复指数z3=1/12时序列虚部); subplot(5,2,7); stem(n,real(x4)); xlabel (时间序列n);ylabel(实部);title(复指数z4=2+j*pi/6时序列实部); subplot(5,2,8); stem(n,imag(x4)); xlabel (时间序列n);ylabel(虚部);title(复指数z4=2+j*pi/6时序列虚部); subplot(5,2,9); stem(n,real(x5)); xlabel (时间序列n);ylabel(实部);title(复指数z5=j*pi/6时序列实部); subplot(5,2,10); stem(n