数字信号处理上机实验 作业结果与说明 实验三、四、五

上机频谱分析过程及结果图 上机实验三：IIR 低通数字滤波器的设计

姓名：赵晓磊 学号：赵晓磊 班级：02311301 科目：数字信号处理B

一、实验目的

1、熟悉冲激响应不变法、双线性变换法设计IIR 数字滤波器的方法。

2、观察对实际正弦组合信号的滤波作用。

二、实验内容及要求

1、分别编制采用冲激响应不变法、双线性变换法设计巴特沃思、切贝雪夫I 型，切贝雪夫II 型低通IIR 数字滤波器的程序。 要求的指标如下：

通带内幅度特性在低于πω3.0=的频率衰减在1dB 内，阻带在πω6.0=到π之间的频率上衰减至少为20dB 。抽样频率为2KHz ，求出滤波器的单位取样响应，幅频和相频响应，绘出它们的图，并比较滤波性能。

（1）巴特沃斯，双线性变换法

Ideal And Designed Lowpass Filter Magnitude Response

frequency in Hz

|H [e x p (j w )]

|

frequency in pi units

|H [e

x p (j w )]|

Designed Lowpass Filter Phase Response in radians frequency in pi units

a r g (H [e x p (j w )]

（2）巴特沃斯，冲激响应不变法

（3）切贝雪夫I 型，双线性变换法

（4）切贝雪夫Ⅱ型，双线性变换法

综合以上实验结果，可以看出，使用不同的模拟滤波器数字化方法时，滤波器的性能可能产生如下差异：

使用冲击响应不变法时，使得数字滤波器的冲激响应完全模仿模拟滤波器的冲激响应，也就是时域逼急良好，而且模拟频率和数字频率之间呈线性关系；但频率响应有混叠效应。

frequency in Hz

|H [e x p (j w )]

|

Designed Lowpass Filter Magnitude Response in dB

frequency in pi units

|H [e x p (j w )]

|

frequency in pi units

a r g (H [e x p (j w )]

Ideal And Designed Lowpass Filter Magnitude Response

frequency in Hz

|H [e x p (j w )]

|

frequency in pi units

|H [e x

p (j w )]|

frequency in pi units

a r g (H [e x p (j w )

]

Ideal And Designed Lowpass Filter Magnitude Response

frequency in Hz

|H [e x p (j w )]

|

frequency in pi units

|H [e

x p (j w )]|

Designed Lowpass Filter Phase Response in radians

frequency in pi units

a r g (H [e x p (j w )]

使用双线性变换法时，克服了多值映射的关系，避免了频率响应的混叠现象；在零频率附近，频率关系接近于线性关系，高频处有较大的非线性失真。

综合以上实验结果，通过下图来说明使用不同的模拟滤波器逼近方法设计滤波器时的性能差异：

巴特沃斯滤波器通频带内的幅频响应曲线最大限度平坦，没有起伏，而在阻频带内则逐渐下降为零。无论在通带内还是阻带内都是频率的单调函数，在通带边缘满足指标时，通带内会有富裕量，也就是会超过指标的要求。

切比雪夫滤波器则将指标精度均匀分布在通带或阻带内，在通带或阻带上频率响应幅度等波纹波动。对于切比雪夫Ⅰ型，在通带中等波纹，在阻带中单调。对于切比雪夫Ⅱ型，在通带内是单调的，在阻带内是等波纹的。

当滤波器具有相同的阶数时：

巴特沃思滤波器通带最平坦，阻带下降慢；切比雪夫滤波器通带等纹波，阻带下降较快。

2、通过上述滤波器对以下正弦组合信号进行滤波，绘出原始信号和滤波后的信号曲线，并比较滤波效果。

X(t)=sin(2*pi*100*t)+sin(2*pi*800*t) （1）巴特沃斯，双线性变换法

（2）巴特沃斯，冲激响应不变法

（3）切贝雪夫I 型，双线性变换法

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s

x Filtered signal:y

time in s

y

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s x

Filtered signal:y

time in s

y

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s x

Filtered signal:y

time in s

y

（4）切贝雪夫Ⅱ型，双线性变换法

巴特沃思滤波器能够较好滤除噪声，原信号几乎没有失真，使用双线性变换法和冲击响应不变法数字化的效果差别不大。

切贝雪夫I 型原信号能够较好滤除噪声，原信号幅度有微小衰减，波形几乎没有失真。 切贝雪夫Ⅱ型相比于前两种方法滤波效果较差，掺杂了一定幅度的高频噪声信号，使信号的波形发生轻微失真。

三、问题与作业：

1、截止频率πω3.0=和πω6.0=对应的模拟频率为多少赫兹？

2、改变输入正弦信号的频率，比较滤波效果。

以巴特沃斯，双线性变换法为例

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s x

Filtered signal:y

time in s

y

100hz 800hz

100hz 600hz

100hz 400hz 100hz 200hz

随着高频混杂信号频率逐渐降低，滤波效果越来越差。

因为对于低通滤波器而言，随着信号频率的逐渐降低，该高频信号对应的幅频响应值越来越大，直至几乎不能滤波，与原信号混叠，使信号波形产生严重失真。

3、程序中，分别使用不同算法得到的数字滤波器系数，进行滤波比较。

可以直接使用[bb1,ab1]=butter(N1,wc1)求解巴特沃思滤波器系数，也可以根据巴特沃思滤波器的设计思想来编制模拟逼近再数字化的方式来求解系数，两种算法得到的滤波器的滤波效果没有明显差别。

4、比较采用冲激响应不变法、双线性变换法设计的巴特沃思、切贝雪夫I 型，切贝雪夫II

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s

x

Filtered signal:y

time in s

y

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s

x

Filtered signal:y

time in s

y

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s x

time in s

y

Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t)

time in s x

time in s

y

型低通IIR 数字滤波器的系统函数的区别。 系统函数阶数不同。

m0301_Butterworth_bilinear N = 4 m0302_Butterworth_impinvar N = 6 m0303_Chebyshev1_bilinear N = 4 m0304_Chebyshev2_bilinear N = 4

5、对输入信号叠加不同幅度的噪声，比较滤波效果。

X1 + 1 * x2

X1 + 10 * x2

time in s

x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

time in s

x

Filtered signal:y

time in s

y

X1 + 100 * x2 X1 + 1000 * x2

随着叠加噪声幅度的逐渐增大，滤波效果越来越差，直至几乎不能完成滤波功能。

原因为信噪比很小时，叠加的噪声幅度远大于滤波器阻带衰减参数20dB，于是非理想的滤波器对于阻带内噪声的衰减不足以消减噪声幅度，就不能完成预期的滤波功能。

上机实验四：FIR 低通数字滤波器的设计

一、实验目的

1、熟悉窗函数法、频率取样法设计IIR 数字滤波器的方法。

2、观察对实际正弦组合信号的滤波作用。

二、实验内容及要求

1、分别编制采用矩形窗、三角窗、汉宁窗、海明窗、布拉克曼窗、凯泽窗设计线性相位低通FIR 数字滤波器的程序。要求的低通FIR 数字滤波器指标如下：

通带截止频率

πω2.0=p ，阻带起始频πω4.0=s ，阻带最小衰减 As=50dB ，求

出滤波器的单位取样响应，幅频和相频响应，绘出它们的图形，并比较各种窗函数和窗口长度对滤波性能的影响。

n h d (n

)

Window

n

w (n )

n

h (n )

n h d (n

)

Window

n

w (n )

n

h (n )

n h d (n

)

Window

n

w (n )

n

h (n )

n h d (n

)

Window

n

w (n )

n

h (n )

n h d (n

)

n

w (n )

n

h (n )

n h d (n )

n

w (n )

n

h (n )

Actual Magnitude Response in dB

frequency in pi units

|H [e x p (j w )] d B |

Actual Magnitude Response

w/rad

|H [e x p (j w )]

|

w/rad

a r g

(H [e x p (j w )]

Actual Magnitude Response in dB

frequency in pi units

|H [e x p (j w )] d B |

Actual Magnitude Response

w/rad

|H [e x p (j w )]

|

w/rad

a r g (H [

e x p (j w )]

Actual Magnitude Response in dB

frequency in pi units |H [e x p (j w )] d B |

Actual Magnitude Response

w/rad

|H [e x p (j w )]

|

w/rad

a r

g (H [e x p (j w )]

Actual Magnitude Response in dB

frequency in pi units |H [e x p (j w )] d B |

Actual Magnitude Response

w/rad

|H [e x p (j w )]|

w/rad

a r g (H [e x p (j w )]

分析比较matlab 运行结果，可知各种窗函数和窗口长度对滤波性能的影响如下： 改变窗口截取长度N ，将会改变窗谱的主瓣宽度、改变w 坐标的比例以及改变以及旁边幅度绝对值的大小。当N 增加时，会减小过渡带的宽度，改善衰减特性；但同时由于吉布斯效应，起伏震荡变密。

改变窗口截取长度N 不会减少带内波动或加大阻带衰减。为了加大阻带衰减，改善滤波器的特性，只能改善窗函数的形状。但主瓣宽度要窄和旁瓣要尽可能小这两点不能兼得，一般总是通过牺牲主瓣宽度来换取对旁瓣的抑制。对于矩形窗、汉宁窗、海明窗。布莱克曼窗，旁瓣的衰减逐步增加，主瓣宽度响应加宽。

用矩形窗设计的滤波器过渡带最窄，但阻带最小衰减也最小，仅-21dB ，布莱克曼窗设计的阻带最小衰减最大，达-74dB ，但过渡带最宽，约为矩形窗的三倍。

矩形窗的优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频干扰和泄漏，甚至出现负谱现象。

三角窗与矩形窗比较，主瓣宽约等于矩形窗的两倍，但旁瓣小，而且无负旁瓣。 汉宁窗可以看作是3个矩形时间窗的频谱之和，它可以使用旁瓣互相抵消，消去高频干扰和漏能。主瓣加宽并降低，旁瓣则显著减小，旁瓣衰减速度也较快。

海明窗是改进的升余弦窗，海明窗与汉宁窗都是余弦窗，只是加权系数不同。海明窗加权的系数能使旁瓣达到更小。

frequency in pi units |H [e x p (j w )] d B

|

w/rad

|H [e x p (j w )]

|

Actual Phase Response in radians w/rad

a r g (H [e x p (j w )

]

frequency in pi units |H [e x p (j w )] d B

|

w/rad

|H [e x p (j w )]

|

Actual Phase Response in radians w/rad

a r g (H [e x p (j w )]

布莱克曼窗主瓣宽，旁瓣比较低，但等效噪声带宽比汉宁窗要大一点，波动却小一点。频率识别精度最低，但幅值识别精度最高，有更好的选择性。

凯泽窗是一种可调窗,可以通过窗函数的形状来控制频谱波纹或阻带衰减指标。定义了一组可调的由零阶贝塞尔Bessel 函数构成的窗函数，通过调整参数β可以在主瓣宽度和旁瓣衰减之间自由选择它们的比重。对于某一长度的Kaiser 窗，给定β，则旁瓣高度也就固定了。

2、编制采用频率取样法设计一个20阶的线性相位低通FIR 数字滤波器的程序。

要求通带截止频率

πω2.0=p ，阻带起始频率πω4.0=s 。求出滤波器的单位取样

响应，幅频和相频响应，绘出它们的图形，并比较阶数对滤波性能的影响。

过渡带抽样点增加 (a)一点：T=0.5； (b)二点：T=0.66,T=0.33；比较对滤波性能的影响。

% wp = 0.3pi,ws = 0.6pi, M =20

Hrs = [1,1,1,1,zeros(1,13),1,1,1]

Hrs=[ones(1,kwp),zeros(1,M-2*kwp+1),ones(1,kwp-1)]

frequency in pi units H r (k )

n

h (n )

frequency in pi units

H r (w )

frequency in pi units

D

e c i b e l s

w/rad

|H [e x p (j w

)]|

w/rad

a r g (H

[e x p (j w )]

time in s x

time in s

y

frequency in pi units H r (k )

n

h (n )

frequency in pi units

H r (w )

frequency in pi units

D

e c i b e l s

w/rad

|H [e x p (j

w )]|

w/rad

a r g (H [

e x p (j w )]

time in s x

time in s

y

% wp = 0.3pi,ws = 0.6pi,T=0.5,N=20

Hrs = [1,1,1,1,0.5,zeros(1,11),0.5,1,1,1]

Hrs=[ones(1,kwp),zeros(1,M-2*kwp+1),ones(1,kwp-1)]

% wp = 0.3pi,ws = 0.6pi,T=0.66,T=0.33,N=20

Hrs = [1,1,1,1,0.66,0.33,zeros(1,9),0.33,0.66,1,1,1]

Hrs=[ones(1,kwp),zeros(1,M-2*kwp+1),ones(1,kwp-1)]

增加了过渡带的抽样点之后，增加了过渡带，减小了频带边缘的突变，也就减小了起伏振荡，但优点是增大了阻带的最小衰减。

增加滤波器的阶数，取样点数N 增加，对给定的幅频曲线频率抽样更密集，拟合效果更好。

frequency in pi units H r (k )

n

h (n )

frequency in pi units

H r (w )

frequency in pi units

D

e c i b e l s

w/rad

|H [e x p (j w

)]|

w/rad

a r g (H [

e x p (j w )]

time in s x

time in s

y

frequency in pi units H r (k )

n

h (n )

frequency in pi units

H r (w )

frequency in pi units

D

e c i b e l s

w/rad

|H [e x p (j w )

]|

w/rad

a r g (H [

e x p (j w )]

time in s x

time in s

y

frequency in pi units H r (k )

n

h (n )

frequency in pi units

H r (w )

frequency in pi units

D e c i b e l s

w/rad

|H [e x p (j w

)]|

w/rad

a r g (H [

e x p (j w )]

time in s x

time in s

y

frequency in pi units H r (k )

n

h (n )

frequency in pi units

H r (w )

frequency in pi units

D e c i b e l s

w/rad

|H [e x p (j w

)]|

w/rad

a r g (H [

e x p (j w )]

time in s x

time in s

y

3、通过上述滤波器对以下正弦组合信号进行滤波，抽样频率为2KHz ，绘出原始信号和滤波后的信号曲线，并比较滤波效果。 x(t)=sin(2*pi*100*t)+sin(2*pi*800*t)

过渡带上增加抽样点数，表现在滤波器对正弦信号组的滤波效果上差别不大，增加了抽样点数后的滤波器对高频噪声的衰减抑制特性更好。

time in s x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

time in s x

Filtered signal:y

time in s

y

上机实验五：等波纹数字滤波器的设计

一、实验目的

熟悉等波纹IIR 、 FIR 数字滤波器的设计方法。

二、实验内容及要求

分别编制采用remez 交换算法设计等波纹FIR 数字滤波器和采用椭圆低通模拟滤波器设计等波纹IIR 的程序。

要求的低通FIR 数字滤波器指标如下：

通带截止频率kHz f c 15=，允许的幅度失真为02.0%,21=±δ，阻带边界频率为

kHz f s 20=，最小衰减为 R S =40dB ，抽样频率为 80KHz ，求出滤波器的单位取样响应，

幅频和相频响应，绘出它们的图形。

2

4

6810

12

Actual Impulse Response

n

h (n )0

50

100

150

200250300350400

450

-80-60-40-200Actual Magnitude Response in dB

frequency in pi units

|H [e x p (j w )]|

50

100

150200250

300

350

-505Actual Phase Response in radians frequency in pi units

a r g (H [e x p (j w )]

2、编制采用椭圆低通模拟滤波器设计低通IIR 数字滤波器的程序。 指标同上，比较滤波性能。

采用remez 交换算法设计等波纹FIR 数字滤波器，得到严格的线性相位。

采用椭圆低通模拟滤波器设计等波纹IIR 相位非线性，过渡带宽度比FIR 略宽，可用比FIR 较少的阶数来满足指标要求。

5

10

15

20253035

40

Actual Impulse Response

n

h (n )

0.10.2

0.30.40.50.60.70.80.91

Actual Magnitude Response in dB

frequency in pi units

|H [e x p (j w )]|

0.1

0.2

0.3

0.40.50.60.7

0.8

0.9

1

-505Actual Phase Response in radians

frequency in pi units

a r g (H [e x p (j w )]

数字信号处理实验报告

实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形

数字信号处理实验一

实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。

（2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结

数字信号处理实验报告

数字信号处理作业提交日期：2016年7月15日

实验一 维纳滤波器的设计 第一部分 设计一维纳滤波器。 (1)产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪(信噪比分别为20,10,6dB dB dB )，得到观测数据123(),(),()x n x n x n 。 (2)估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 1 实验原理 滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 维纳（Wiener ）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。 设一线性系统的单位样本响应为()h n ，当输入以随机信号()x n ，且 ()() () x n s n v n =+，其中()s n 表示原始信号，即期望信号。()v n 表示噪声，则输出()y n 为()=()()m y n h m x n m -∑，我们希望信号()x n 经过线性系统()h n 后得到的()y n 尽可能接近 于()s n ，因此称()y n 为估计值，用?()s n 表示。 则维纳滤波器的输入-输出关系可用下面表示。 设误差信号为()e n ,则?()()()e n s n s n =-,显然)(n e 可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即 它的平方的统计期望最小：222?[|()|][|()()|][|()()|]E e n E s n s n E s n y n =-=-=min 。而要使均方误差最小，则需要满足2[|()|]j E e n h ?=0. 进一步导出维纳-霍夫方程为：()()()()*(),0,1,2...xs xx xx i R m h i R m i R m h m m =-==∑ 写成矩阵形式为：xs xx R R h =，可知：1xs xx h R R -=。表明已知期望信号与观测数据的互相关函数以及观测信号的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的

西电数字信号处理上机实验报告

数字信号处理上机实验报告 14020710021 张吉凯 第一次上机 实验一： 设给定模拟信号()1000t a x t e -=，t 的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的()a x t 进行采样。 ○1()()15000s a f x t x n =以样本秒采样得到。 ()()11j x n X e ω画出及其频谱。 ○2()()11000s a f x t x n =以样本秒采样得到。 ()() 11j x n X e ω画出及其频谱。 比较两种采样率下的信号频谱，并解释。 （1）MATLAB 程序： N=10; Fs=5; T s=1/Fs; n=[-N:T s:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e^(j\omega))');

ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为%fKHZ\n',eband); 运行结果： 等效带宽为12.110000KHZ

数字信号处理实验作业

实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。

数字信号处理实验

实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π=

1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2）

数字信号处理上机实验代码

文件名：tstem.m（实验一、二需要） 程序： f unction tstem(xn,yn) %时域序列绘图函数 %xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称（字符串）n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel('yn'); axis([0,n(end),min(xn),1.2*max(xn)]); 文件名：tplot.m（实验一、四需要） 程序： function tplot(xn,T,yn) %时域序列连续曲线绘图函数 %xn:信号数据序列,yn:绘图信号的纵坐标名称（字符串） %T为采样间隔 n=0;length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]); 文件名：myplot.m（实验一、四需要）

%(1)myplot;计算时域离散系统损耗函数并绘制曲线图。function myplot(B,A) %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000) m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0,1,-80,5]);title('损耗函数曲线'); 文件名：mstem.m（实验一、三需要） 程序： function mstem(Xk) %mstem(Xk)绘制频域采样序列向量Xk的幅频特性图 M=length(Xk); k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率（关于pi归一化值） stem(wk,abs(Xk),'.');box on;%绘制M点DFT的幅频特性图xlabel('w/\pi');ylabel('幅度'); axis([0,2,0,1.2*max(abs(Xk))]); 文件名：mpplot.m（实验一需要）

数字信号处理实验报告一

武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02

一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分

实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列');

实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列');

正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列');

随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列');

数字信号处理实验答案完整版

数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】

实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。 clear all; a=[1 2 3 4]; b=[3 4 5 6]; c=a+b; d=a-b; e=a.*b; f=a./b; g=a.^b; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。

数字信号处理上机实验(第三版)

数字信号处理实验（Matlab） 实验一: 系统响应及系统稳定性 %实验1：系统响应及系统稳定性 close all;clear all %======内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=R8(n) x2n=ones(1,128); %产生信号x2(n)=u(n) hn=impz(B,A,58); %求系统单位脉冲响应h(n) subplot(2,2,1);y='h(n)';tstem(hn,y); %调用函数tstem绘图 title('(a)系统单位脉冲响应h(n)');box on y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) subplot(2,2,2);y='y1(n)';tstem(y1n,y); title('(b)系统对R8(n)的响应y1(n)');box on y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) subplot(2,2,4);y='y2(n)';tstem(y2n,y); title('(c)系统对u(n)的响应y2(n)');box on %===内容2：调用conv函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=R8(n) h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)];

数字信号处理实验报告(实验1_4)

实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B

（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；

（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；

（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。

5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]

（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;

数字信号处理实验三

实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期

x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response

数字信号处理上机实验答案(全)1

第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一 系统响应及系统稳定性。 实验二 时域采样与频域采样。 实验三 用FFT 对信号作频谱分析。 实验四 IIR 数字滤波器设计及软件实现。 实验五 FIR 数字滤波器设计与软件实现 实验六 应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 10.1 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握 求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞→n 时，系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n 的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤 （1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter 函数或conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 （2）给定一个低通滤波器的差分方程为

数字信号处理实验作业

实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：

原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。

数字信号处理实验4

数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线

（2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：

西南交大数字信号处理报告

信息科学与技术学院本科三年级 数字信号处理实验报告 2011 年12 月21日

实验一 序列的傅立叶变换 实验目的 进一步加深理解DFS,DFT 算法的原理;研究补零问题;快速傅立叶变换 （FFT ）的应用。 实验步骤 1. 复习DFS 和DFT 的定义，性质和应用； 2. 熟悉MATLAB 语言的命令窗口、编程窗口和图形窗口的使用；利用提供的 程序例子编写实验用程序；按实验内容上机实验，并进行实验结果分析；写出完整的实验报告，并将程序附在后面。 实验内容 1. 周期方波序列的频谱试画出下面四种情况下的的幅度频谱,并分析补零后，对信号频谱的影响。 实验结果： 60 ,7)4(;60,5)3(; 40,5)2(;20,5)1()] (~[)(~,2,1,01 )1(,01,1)(~=========±±=???-+≤≤+-+≤≤=N L N L N L N L n x DFS k X m N m n L m N L m N n m N n x ) 52.0cos()48.0cos()(n n n x ππ+=

2. 有限长序列x(n)的DFT （1） 取x(n)(n=0:10)时，画出x(n)的频谱X(k) 的幅度； （2） 将（1）中的x(n)以补零的方式，使x(n)加长到(n:0~100)时，画出 x(n)的频谱X(k) 的幅度； （3） 取x(n)(n:0~100)时，画出x(n)的频谱X(k) 的幅度。利用FFT 进行谱分析 已知：模拟信号 以t=0.01n(n=0:N-1)进行采样，求N 点DFT 的幅值谱。 请分别画出N=45; N=50;N=55;N=60时的幅值曲线。 实验结果： ) 8cos(5)4sin(2)(t t t x ππ+=

数字信号处理上机报告-一

数字信号处理上机报告-一

数字信号处理第一次上机实验报告 实验一： 设给定模拟信号()1000t a x t e -=，的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的进行采样。 ○1 。 ○2 。 比较两种采样率下的信号频谱，并解释。 实验一MATLAB 程序： （1） N=10; Fs=5; Ts=1/Fs; n=[-N:Ts:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); t ()a x t ()()15000s a f x t x n =以样本秒采样得到。()()11j x n X e ω画出及其频谱()()11000s a f x t x n =以样本秒采样得到。()() 11j x n X e ω画出及其频谱

axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e ^(j\omega))'); ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为 %fKHZ\n',eband); 运行结果：

数字信号处理第二章上机作业

第二章上机作业 1、ljdt(A,B)函数定义 function ljdt(A,B) p=roots(A); q=roots(B); p=p'; q=q'; x=max(abs([p q 1])); x=x+0.1; y=x; clf hold on axis([-x x -y y]) w=0:pi/300:2*pi; t=exp(i*w); plot(t) axis('square') plot([-x x],[0 0]) plot([0 0],[-y y]) text(0.1,x,'jIm[z]') text(y,1/10,'Re[z]') plot(real(p),imag(p),'x') plot(ral(q),imag(q),'o') title('pole-zero diagram for discrete system') hold off 例2.26 a=[3 -1 0 0 0 1]; b=[1 1]; ljdt(a,b) p=roots(a) q=roots(b) pa=abs(p) 程序运行结果如下： P= 0.7255+0.4633i 0.7255+0.4633i -0.1861+0.7541i -0.1861-0.7541i -0.7455 q=

-1 pa= 0.8608 0.8608 0.7768 0.7768 0.7455 例2.27 b=[0 1 2 1];a=[1 -0.5 -0.005 0.3]; subplot 311 zplane(b,a);xlabel('实部');ylabel('虚部'); num=[0 1 2 1];den=[1 -0.5 -0.005 0.3]; h=impz(num,den); subplot 312

数字信号处理上机实验答案完整版

数字信号处理上机实验 答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】

第十章上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一系统响应及系统稳定性。 实验二时域采样与频域采样。 实验三用FFT对信号作频谱分析。 实验四 IIR数字滤波器设计及软件实现。 实验五 FIR数字滤波器设计与软件实现 实验六应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 functiontstem(xn,yn) %时域序列绘图函数 %xn:信号数据序列，yn:绘图信号的纵坐标名称（字符串） n=0:length(xn)-1; stem(n,xn,'.');boxon xlabel('n');ylabel(yn); axis([0,n(end),min(xn),*max(xn)]) 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可

数字信号处理实验报告

3.(1)用双线性变换法设计一个Chebyshev型高通滤波器程序如下 Rp=1.2;Rs=20;T=0.001;fp=300;fs=200; wp=2*pi*fp*T;ws=2*pi*fs*T; wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); [n,wn]=cheb1ord(wp1,ws1,Rp,Rs,'s'); [b,a]=cheby1(n,Rp,wn,'high','s'); [bz,az]=bilinear(b,a,1/T); [db,mag,pha,grd,w]=freqz_m(bz,az);plot(w/pi,db); axis([0,1,-30,2]); 3.(2) a用双线性变换法设计一个Butterworth型数字低通滤波器程序如下Rp=1;Rs=25;T=0.001;fp=300;fs=200; wp=2*pi*fp*T;ws=2*pi*fs*T; wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); [n,wn]=buttord(wp1,ws1,Rp,Rs,'s'); [b,a]=butter(n,wn,'low','s'); [bz,az]=bilinear(b,a,1/T); [db,mag,pha,grd,w]=freqz_m(bz,az);plot(w/pi,db); axis([0,1,-30,2]); b用脉冲响应不变法设计一个Butterworth数字低通滤波器的程序如下：wp=400*pi;ws=600*pi;Rp=1;Rs=25; [n,wn]=buttord(wp,ws,Rp,Rs,'s') [b,a]=butter(n,wn,'s') [db,mag,pha,w]=freqs_m(b,a,500*2*pi);