资源描述：
《数字信号处理实验(西电 高西全教材)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、《数字信号处理》第三版，高西全、丁玉美，实验程序清单：实验0实验内容：1、实验用Matlab工具箱函数简介熟悉以下函数的功能、格式（1）abs（2）angle（3）conv（4）filter（5）freqz（6）impz（7）fft（8）ifft（9）plot（10）stem（11）subplot（12）figure2、常用序列的产生及其频谱分析（1）单位脉冲序列及其频谱n=0:50;%定义序列的长度是50x=[1,zeros(1,50)];%注意：MATLAB中数组下标从1开始%x(1)=1;closeall;subplot(3,1,1);stem(x
2、);title('单位冲击信号序列');k=-25:25;X=x*(exp(-j*pi/25)).^(n'*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('单位冲击信号的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title('单位冲击信号的相位谱')（2）矩形序列及其频谱n=1:50x=sign(sign(10-n)+1);closeall;subplot(3,1,1);stem(x);title(‘矩形序列’);
3、k=-25:25;X=x*(exp(-j*pi/25)).^(n’*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘矩形序列的幅度谱’);angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(‘矩形序列的相位谱’)(3)特定冲击串及其频谱%h2(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)n=1:50;%定义序列的长度是50x=zeros(1,50);%注意：MATLAB中数组下标从1开始x(1)=1;
4、x(2)=2.5;x(3)=2.5;x(4)=1;closeall;subplot(3,1,1);stem(x);title(‘特定冲击串’);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘特定冲击串的幅度谱’);angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(‘特定冲击串的相位谱’)（4）采样信号序列及其频谱%x(n)=Ae-anTsin(w0n
5、T)u(nT)（0<=n<=50）;%A=444.128;a=50*sqrt(2.0)*pi;fs=1kHz(T=0.001);w0=50*sqrt(2.0)*pin=0:50;%定义序列的长度是50A=444.128;%设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001;%采样率，fs=1kHzw0=50*sqrt(2.0)*pi;x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”closeall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x
6、(n)的图形title(‘理想采样信号序列’);k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(‘理想采样信号序列的相位谱’)实验一 系统响应及系统稳定性%====内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====cl
7、oseall;clearallA=[1,-0.9];B=[0.05,0.05];%系统差分方程系数向量B和Ax1n=[11111111zeros(1,50)];%产生信号x1(n)=R8(n)x2n=ones(1,128);%产生信号x2(n)=u(n)hn=impz(B,A,58);%求系统单位脉冲响应h(n)subplot(3,1,1);stem(hn);title('(a)系统单位脉冲响应h(n)');y1n=filter(B,A,x1n);%求系统对x1(n)的响应y1(n)subplot(3,1,2);stem(y1n);title('(b)系
8、统对R8(n)的响应y1(n)');y2n=filter(B,A,x2n);%求
