数字信号处理实验（西电高西全教材）.pdfVIP

百学须先立志。——朱熹

《数字信号处理》第三版，高西全、丁玉美，

实验程序清单：

实验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）figure

2、常用序列的产生及其频谱分析

（1）单位脉冲序列及其频谱

n=0:50;%定义序列的长度是50

x=[1,zeros(1,50)];%注意：MATLAB中数组下标从1开始

%x(1)=1;

closeall;

subplot(3,1,1);stem(x);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:50

x=sign(sign(10-n)+1);

closeall;

subplot(3,1,1);stem(x);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(„矩形序列的相位谱‟)

(3)特定冲击串及其频谱

%h(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)

2

n=1:50;%定义序列的长度是50

x=zeros(1,50);%注意：MATLAB中数组下标从1开始

x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;

closeall;

subplot(3,1,1);stem(x);title(„特定冲击串‟);

k=-25:25;

0

百学须先立志。——朱熹

X=x*(exp(-j*pi/12.5)).^(n‟*k);

magX=abs(X);%绘制x(n)的幅度谱

subplot(3,1,2);stem(magX);title(„特定冲击串的幅度谱‟);