数字信号处理
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字信号处理实验报告（1）
班级：通信09-1
姓名：陈阳
学号：3号
实验7 z 变换及其应用
一．实验目的
（1） 加深对离散系统变换域分析——z 变换的理解。
（2） 掌握进行z 变换和z 反变换的基本方法，了解部分分式法在z 反变换中的应用。 （3） 掌握使用MATLAB 语言进行z 变换和z 反变换的常用子函数。
二．实验内容
(4)用部分分式法求解下列系统函数的z 反变换，写出x(n)的表示式，并用图形与impz 求得的结果相比较，取前10
个点作图。
b=[10,20,0,0];a=[1,8,19,12]; [r p c]=residuez(b,a) N=10;n=0:N-1;
x=r(1)*p(1).^n+r(2)*p(2).^n+r(3)*p(3).^n; subplot(1,2,1),stem(n,x);
title('用部分分式法求反变换x(n)'); x2=impz(b,a,N);
subplot(1,2,2);stem(n,x2); title('用impz 求反变换x(n)');
0510-7
-6-5-4-3-2-1012x 10
6用部分分式法求反变换x(n)0510
-7
-6-5-4-3-2-1012x 106
用impz 求反变换x(n)
b=[0,0,5];a=[1,1,-6]; [r p c]=residuez(b,a); N=10;n=0:N-1;
x=r(1)*p(1).^n+r(2)*p(2).^n; subplot(1,2,1),stem(n,x);
title('用部分分式法求反变换x(n)'); x2=impz(b,a,N);
subplot(1,2,2);stem(n,x2); title('用impz 求反变换x(n)');
5
10
-7000
-6000-5000-4000-3000-2000-100001000
20003000用部分分式法求反变换x(n)0
5
10
-7000
-6000-5000-4000
-3000-2000-100001000
20003000用impz 求反变换x(n)
b=[1,0,0,0];a=[1,-0.9,-0.18,0.729]; [r p c]=residuez(b,a); N=10;n=0:N-1;
x=r(1)*p(1).^n+r(2)*p(2).^n; subplot(1,2,1),stem(n,x);
title('用部分分式法求反变换x(n)'); x2=impz(b,a,N);
subplot(1,2,2);stem(n,x2); title('用impz 求反变换x(n)');
0510
-1.5
-1
-0.5
0.5
1
1.5
用部分分式法求反变换x(n)0510
-1.5
-1
-0.5
0.5
1
用impz 求反变换x(n)
实验9 离散系统的零极点分析
一、实验目的
(1)了解离散系统的零极点与系统因果性和稳定性的关系。 (2)观察离散系统零极点对系统冲激响应的影响。
(3)熟悉MATLAB 中进行离散系统零极点分析的常用子函数。
二、实验内容
(3)已知离散时间系统函数分别为
求该系统的零极点及零极点分布图，并判断系统的因果稳定性。
b=[5,10,-15];a=[1,2,-8]; rz=roots(b) rp=roots(a)
subplot(1,2,1),zplane(b,a); title('系统的零极点分布图'); subplot(1,2,2),impz(b,a,20); title('系统的冲激响应')； xlabel('n');ylabel('x(n)');
-4
-202
-6
-4
-2
2
4
6
Real Part
I m a g i n a r y P a r t
系统的零极点分布图05
1015
-3
-2.5
-2
-1.5
-1-0.5
0.5
1x 1011
n (samples)
A m p l i t u d e
Impulse Response
b=[4,-1.6,-1.6,4];a=[1,0.4,0.35,-0.4]; rz=roots(b) rp=roots(a)
subplot(1,2,1),zplane(b,a); title('系统的零极点分布图'); subplot(1,2,2),impz(b,a,20); title('系统的冲激响应')； xlabel('n');ylabel('x(n)');
-1
-0.5
00.51
-2
-1.5-1-0.5
00.51
1.5
2Real Part
I m a g i n a r y P a r t
系统的零极点分布图05
1015
-4
-20
2
4
6
8
n (samples)
A m p l i t u d e
Impulse Response
b=[2,-1,-0.5];a=[1,-1,-0.25,0.25]; rz=roots(b) rp=roots(a)
subplot(1,2,1),zplane(b,a); title('系统的零极点分布图'); subplot(1,2,2),impz(b,a,20); title('系统的冲激响应')； xlabel('n');ylabel('x(n)');