数字信号处理实验四
更新时间:2024-02-27 16:18:01 阅读量: 综合文库 文档下载
实 验 报 告
实验名称_____离散系统分析______
课程名称____数字信号处理________
院 系 部:电气与电子工程 专业班级:信息1002 学生姓名:王萌 学 号: 11012000219 同 组 人: 实验台号: 指导教师:范杰清 成 绩: 实验日期:
华北电力大学
实验四 离散系统分析
一、实验目的
深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义,熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 二、 实验原理
MATLAB提供了许多可用于分析线性时不变连续系统的函数,主要包含有系统函数、系统时域响应、系统频域响应等分析函数。
1. 离散系统的时域响应
在调用MATLAB函数时,需要利用描述该离散系统的系数函数。对差分方程进行Z变换即可得系统函数:
Y(z)b0?b1z?1???bM?1z?(M?1)?bMz?Mb(z)?? H(z)??1?(N?1)?NX(z)1?a1z???aN?1z?aNza(z)在MATLAB中可使用向量a和向量b分别保存分母多项式和分子多项式的系数:
a ?[1,a1,?,aN?1,aN]
b?[b0,b1,?,bM?1,bM]这些系数均从z0按z的降幂排列。
2.离散系统的系统函数零极点分析
离散LTI系统的系统函数H(z)可以表示为零极点形式:
Y(z)b0?b1z?1???bM?1z?(M?1)?bMz?M(z?z1)(z?z2)...(z?zM)H(z )???kX(z)1?a1z?1???aN?1z?(N?1)?aNz?N(z?p1)(z?p2)...(z?pN)使用MATLAB提供的roots函数计算离散系统的零极点;
使用zplane函数绘制离散系统的零极点分布图。 注意:在利用这些函数时,要求H(z)的分子多项式和分母多项式的系数的个数相等,若不等则需要补零。
3.离散系统的频率响应
当离散因果LTI系统的系统函数H(z)的极点全部位于z平面单位圆内时,系统的频率响应可由H(z)求出,即
H(ej?)?H(z)z?ej??H(ej?)ej?(?) [H, w]=freqz(b, a, n):
计算系统的n点频率响应H,w为频率点向量。 H=freqz(b, a, w) :计算系统在指定频率点向量w上的频响;
freqz(b,a): 绘制频率响应曲线。
其中:b和a分别为系统函数H(z)的分子分母系数矩阵; 4.利用DTFT和DFT确定离散系统的特性
在很多情况下,需要根据LTI系统的输入和输出对系统进行辨识,即通过测量系统在已知输入x[k]激励下的响应y[k]来确定系统的特性。
若系统的脉冲响应为h[k],由于存在y[k]=x[k]* h[k],因而可在时域通过解卷积方法求解h[k]。
在实际应用中,进行信号解卷积比较困难。因此,通常从频域来分析系统,这样就可以将时域的卷积转变为频域的乘积,从而通过分析系统输入序列和输出序列的频率特性来确定系统的频率特性 H ( e j? ) ,再由 H ( e ) 得到系统的脉冲响应h[k]。
若该LTI系统输入x[k]的DTFT为 X (e ) ,系统输出y[k]的DTFT为 Y ( e j? ) ,则系统的频率特性可表示为:
j?Y(e)j?H(e)?X(ej?)j?j?有限长序列的DTFT可以利用FFT计算出其在区间
???0,2π?内的
N个等间隔频率点上的样点值:
X = fft(x,N); Y = fft(y,N); 再利用 H = Y./X 和 h = ifft(H,N) 可以得到系统的单位脉冲响应h[k]。 三、 实验内容
1.已知某LTI系统的差分方程为:
y[k]?1.143y[k?1]?0.412y[k?2]?0.0675x[k]?0.1349x[k?1]?0.0675x[k?2]y[?1]?1,y[?2]?2,输入 x[k]?u[k]计(1)初始状态
算系统的完全响应。
??7?x[k]?cos(k)u[k];x[k]?cos(k)u[k];x[k]?cos(k)u[k]23(21)当以下三个信号分别通过系统时,分别计算离散10510系统的 零状态响应:
(3)该系统具有什么特性? (1)
>> N=100; b=[0.0675,0.1349,0.0675]; a=[1,-1.143,0.412]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); >> N=100; b=[0.0675,0.1349,0.0675]; a=[1,-1.143,0.412]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi)
Columns 73 through 78
1.0002 1.0002 1.0002 1.0002 1.0002 1.0002
Columns 79 through 84
1.0002 1.0002 1.0002 1.0002 1.0002 1.0002
Columns 85 through 90
1.0002 1.0002 1.0002 1.0002 1.0002 1.0002
Columns 91 through 96
1.0002 1.0002 1.0002 1.0002 1.0002 1.0002
Columns 97 through 100
1.0002 1.0002 1.0002 1.0002
3.已知LTI系统的输入输出序列分别为
1111(a)x[k]?()ku[k],y[k]?()ku[k]?()ku[k]2424111(b)x[k]?()ku[k],y[k]?()ku[k]?()k?1u[k?1]444(1)利用系统辨识原理从频域分别求解系统的单位脉冲响
h[k]。 应
H(ej?),并求出系统的单位脉冲响(2)利用解析方法确定
h[k]。应 比较解析方法与系统辨识方法得到的系统单位脉冲响应,
分析误差原因。
(1)N=10; (2)k=1:100; K=0:N-1; x1=(1/2).^k; aa=[1,-0.25]; y1=(1/4).*(1/2).^k+(1/4).^k; ab=[5/4,-7/16]; X1=fft(x1);Y1=fft(y1); ba=[1]; bb=[1,-1]; H1=Y1/X1; ay=impz(ab,aa,N); H1=ifft(H1) subplot(2,1,1); x2=(1/4).^k; stem(ay); by=impz(bb,ba,N); y2=(1/4).^k-(1/4).^(k-1); subplot(2,1,2); X2=fft(x2);Y2=fft(y2); stem(by); H2=Y2/X2; h2=ifft(H2)
4.已知某离散系统的输入输出序列。
输入序列:2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,0.001,0.0005,0.0002,0.0001,0.0001,后面的数值均趋于0;
输出序列:0.0056,-0.0259,0.073,-0.1593,0.297, -0.4974,0.7711,-1.1267,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9528,5.7117,-6.3889,6.9034,-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,-29.6298,41.4666,-55.2433,70.7979,-87.7810
(1)绘出输入输出信号的波形。
H(ej?),并绘出其幅频特性。 (2)计算该系统的频率响应
h[k],并绘出其波形。 (3)计算该系统的单位脉冲响应
(1)x=[2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0.0083,0.0041,0.002,0.001,0.0005,0.0002,0.0001,0.0001]; y=[0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.7711,-1.1267,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9528,5.7117,-6.3889,6.9034,-7.1528,7.012,-6.3322,4.9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,-29.6298,41.4666,-55.2433,70.7979,-87.7810]; subplot(2,1,1); stem(x); subplot(2,1,2); stem(y);
(2)x=[2,0.8333,0.3611,0.162,0.0748,0.0354,0.017,0. 0083,0.0041,0.002,0.001,0.0005,0.0002,0.0001,0.000 1]; y=[0.0056,-0.0259,0.073,-0.1593,0.297,-0.4974,0.77 11,-1.1267,1.5702,-2.1037,2.724,-3.4207,4.174,-4.9 528,5.7117,-6.3889,6.9034,-7.1528,7.012,-6.3322,4. 9416,-2.648,-0.7564,5.4872,-11.7557,19.7533,-29.62 98,41.4666,-55.2433,70.7979,-87.7810]; X=fft(x,256);Y=fft(y,256); H=Y./X; H0=abs(H); plot(H0);
(3)h=ifft(H); plot(h0);
[k5.利用load mtlb命令读入一段语音信号得到序列 x ] ,然后k]在该段语音信号中加入500Hz的正弦型干扰信号得到信号y [ ,
利用FFT分析其频谱。
k](1)下列数字滤波器能够滤除信号 y [ 中500Hz正弦型干扰
信号。
0.6877?2.509z?1?3.664z?2?2.509z?3?0.6877z?4H(z)?1?3.594z?1?5.172z?2?3.494z?3?0.945z?4 利用zplane命令做出其零极点分布图,利用命令分析该滤波器的幅频特性和相频特性,比较零极点分布与滤波器频率特性的关系。
(2)利用该数字滤波器滤除信号中的噪声,利用FFT观察其频谱,利用sound函数播放处理前后的信号,比较处理前后的效果。
(1)a=[1,-3.594,5.17,-3.494,0.945];b=[0.6877,-2.509,3.664,-2.509,0.6877]; z=roots(b); p=roots(a); subplot(2,2,1); zplane(b,a); [H,w]=freqz(b,a); subplot(2,2,2); plot(w,abs(H)); xlabel('w(rad)'); ylabel('H0'); title('·ù?μì?D?'); subplot(2,2,3);
(2) load mtlb N=150; x=mtlb(1:N); k=1:N; subplot(2,3,4); plot(k,x); title('′|àí?°'); X=fft(x,150); h=ifft(H); H=fft(H,150); Y=X.*H; y=ifft(Y); subplot(2,3,5); plot((1:N),y); title('′|àíoó');
(2) load mtlb N=150; x=mtlb(1:N); k=1:N; subplot(2,3,4); plot(k,x); title('′|àí?°'); X=fft(x,150); h=ifft(H); H=fft(H,150); Y=X.*H; y=ifft(Y); subplot(2,3,5); plot((1:N),y); title('′|àíoó');
1.系统函数的零极点对系统频率特性有何影响?
答:系统零点越接近1,则幅频响应的波谷越低,即滤波器的滤波衰减最低点越低;
系统极点越接近1,则幅频响应的波峰越高,即滤波器的滤波通带最高点越高。
极点主要影响频率响应的峰值,极点愈靠近单位圆,峰值愈尖锐;零点主要影响频率特性的谷值,零点愈靠近单位圆,谷值愈深,当零点在单位圆上时,频率特性为零,一个传递函数有几个极点幅度响应就有几个峰值,对应出现一些谷值。
2.对于因果稳定实系数的低通、高通、带通、带阻数字滤波器,极点分布有何特点?
答:因为是因果稳定系统,所以极点都在单位圆内。 3.离散系统的系统函数的零极点对系统脉冲响应有何影响? 答:系统函数的极点位置决定序列包络的变化趋势和变化频率,极点的半径决定了序列包络的变化趋势,而极点的幅角决定序列包络的变
化频率,而零点位置只影响冲激响应的幅度大小和相位。 4.若某因果系统不稳定,有哪些主要措施可使之稳定? 答:改变参数,让极点在单位圆内。
5.从频域利用DFT确定离散LTI系统的特性,一般会产生哪些误差,如何改善?
答:频谱混叠,对于带限连续信号,只要提高抽样频率使之满足时域抽样定理;对于非带限信号,更具实际情况对其进行低通滤波,使之成为带限信号;频谱泄露,时域加窗使之成为有限长序列;栅栏现象,在序列后补零,构成新序列后再求频谱
6.若使用DFT对连续LTI系统进行辨识,需要解决哪些问题? 答:频谱混叠的问题,对于带限连续信号,要提高抽样频率使之满足时域抽样定理;对于非带限信号,更具实际情况对其进行低通滤波,使之成为带限信号;
正在阅读:
数字信号处理实验四02-27
高三化学化学反应速率04-13
财务汇总06-02
过程检测及仪表教学大纲修正01-06
英译常用的100个中国成语05-16
数字声波测井仪原理及应用11-17
2010年中考英语试题分类汇编专题九(任务型阅读)12-30
广告文案集萃03-08
2017年注册咨询师继续教育有线通信试题2套94分+94分05-21
反比例函数的定义测试题2210-15
- 计算机试题
- 【2012天津卷高考满分作文】鱼心人不知
- 教育心理学历年真题及答案--浙江教师资格考试
- 20180327-第六届“中金所杯”全国大学生金融知识大赛参考题库
- 洪林兴达煤矿2018年度水情水害预测预报
- 基本要道讲义
- 机电设备安装试运行异常现象分析与对策
- 《有机化学》复习资料-李月明
- 非常可乐非常MC2--非常可乐广告策划提案 - 图文
- 2011中考数学真题解析4 - 科学记数法(含答案)
- 企业人力资源管理师三级07- 09年真题及答案
- 基于单片机的光控自动窗帘控制系统设计说明书1 - 图文
- 20160802神华九江输煤皮带机安装方案001
- (共53套)新人教版一生物必修2(全册)教案汇总 word打印版
- 2014行政管理学总复习
- 中国银监会关于加强地方政府融资平台贷款风险监管的指导意见
- 民宿酒店核心竞争与研究
- 游园活动谜语大全2012
- 河南省天一大联考2016届高三英语5月阶段性测试试题(六)(A卷)
- 小型超市管理系统毕业论文详细设计4
- 信号处理
- 实验
- 数字
- 荷兰和意大利注册会计师行业考察报告
- 视觉传达设计教学大纲
- 房屋建筑与装饰工程预算定额培训资料(doc 76页)
- ACCESS简介
- HPV-DNA标本处理的标准操作程序1
- CorelDRAW12使用详解3
- Adobe Premiere入门操作笔记
- 老顾客维护方案 Word 文档
- 本 科 生 毕 业 实 习 手 册
- 格构柱塔吊基础施工方案 - 图文
- 赵河镇第一中心小学优秀教师、常规管理细则
- 云南省高校建立校园短信平台可行性分析
- 跑步机上的卡路里怎么算
- 用工安全协议书
- 2018届高考英语一轮复习高效学案:必修4Unit2含答案解析
- 课程文本
- 微生物分离纯化技术及其应用
- Ansys作业(有限元)
- 2020版高考生物一轮复习 教学案(打包10套)苏教版必修2
- 新时期企业思想政治工作的实践与思考