离散时间系统的状态空间描述

燕山大学

课程设计说明书题目：离散时间系统的状态空间描述

学院(系）：电气工程学院

年级专业：_11级精仪1班

学号： 110103020058

学生姓名：

指导教师：

教师职称：

电气工程学院《课程设计》任务书

课程名称：数字信号处理课程设计

说明：1、此表一式四份，系、指导教师、学生各一份，报送院教务科一份。

2、学生那份任务书要求装订到课程设计报告前面。

电气工程学院教务科

摘要

摘要：线性时不变离散时间系统是最基本的数字系统，差分方程和系统函数是描述系统的常用数学模型，单位脉冲响应和频率响应是描述系统特性的主要特征参数，零状态响应和因果稳定性是系统分析的重要内容。文章从系统的分析流程、系统模型的创建、时域分析、频域分析和因果稳定性分析等方面，介绍了线性时不变离散时间系统的基本分析方法，并以实例形式列举了MATLAB实现程序。

关键词：MATLAB；离散时间系统；系统分析；传输函数

目录

第一章离散时间系统与状态空间描述 (1)

1.1 离散时间系统 (1)

1.2 状态空间描述 (3)

1.3 LSI系统的求解方法 (5)

第二章软件仿真设计 (5)

2.1状态方程 (5)

2.2输出方程 (6)

2.3 LSI系统的单位冲击响应 (7)

第三章仿真结果分析 (10)

3.1状态方程 (10)

3.2 输出方程 (10)

3.3 LSI系统的单位冲击响应 (11)

第四章学习心得 (11)

第五章设计与实验过程中遇到的问题和分析 (12)

第一章相关离散时间系统的知识

1.1离散时间系统

离散时间系统离散时间系统是将一个序列变换成另一序列的系统，它有多种类型，其中线性时不变离散时间系统是最基本、最重要的系统。差分方程反映了系统输入与输出的运动状态，是在时域描述系统的通用数学模型；系统函数是零状态下系统输出与输入的Z变换之比，在时域与频域之间起桥梁作用。分析系统就是在已知系统结构或系统模型条件下，从时域和频域两方面分析系统输入与输出的关系，前者重点研究系统的时间特性，后者主要研究系统的频率特性。下面从系统分析流程、系统模型创建、系统时域分析、系统频域分析和因果稳定性分析等方面，介绍线性时不变离散时间系统的基本分析方法，并以实例形式列举MATLAB在系统分析过程中的具体应用。

二、单位脉冲响应的计算根据差分方程求解单位脉冲激励下系统的零状态响应，或将系统函数进行Z反变换都可算出系统的单位脉冲响应，具体算法可参见参考文献[3]。在MATLAB中描述系统的差分方程或系统函数都是用系数向量表示，调用impz函数就可直接算出系统的单位脉冲响应。如实例1描述的系统，其单位脉冲响应的计算及显示程序如下：b=[0.3,0.06,0,0]; %系数向量不齐后面补0

a=[1,-1.1,0.55,-0.125]; %系数向量不齐后面补0

[hn,n]=impz(b,a,16), %列向求出16点单位脉冲响应

stem(n,hn,'.'); grid; %绘制点状图并加网格

xlabel('n');ylabel('hn');title('单位脉冲响应');

若要写出闭环形式，可调用residuez函数将系统函数展开成部分分式形式，再通过查表求Z反变换即可。

三、系统输出的时域计算

在时域上计算离散时间系统的输出，实际上就是直接求解差分方程或作卷积运算。参考文献[3]列举了迭代法、时域经典法、卷积法等常用方法及应用实例。考虑到分析系统的目的在于综合，系统设计时不存在初始问题，因此，分析系统响应重点分析零状态响应。只要掌握了分析系统的概念、原理和方法，繁杂的计算可由MATLAB完成。

实例2：试计算实例1中，当输入序列分别为单位脉冲、单位阶跃和一般序列时，系统的输出响应。

方法1：调用filter函数实现

b=[0.3,0.06,0,0]; a=[1,-1.1,0.55,-0.125];

x1=[1,zeros(1,15)]; %产生16点单位脉冲序列

x2=ones(1,16); %产生16点单位阶跃序列

x3=exp(-0.2*[0:15]); %用指数序列代表一般序列

y1=filter(b,a,x1), %计算单位脉冲响应

y2=filter(b,a,x2), %计算单位阶跃响应

y3=filter(b,a,x3), %计算一般序列响应

方法2：调用conv函数实现

b=[0.3,0.06,0,0]; a=[1,-1.1,0.55,-0.125];

[hn,n]=impz(b,a,16); %求出16点单位脉冲响应

x=exp(-0.2*(0:15)); %输入或产生一般序列

y1=conv(hn,x), %用线性卷积求系统响应

y2=filter(b,a,x), %用系统函数求系统响应

k=1:16;dy=y1(k)-y2(k), %两种计算的误差对比

结果表明，用有限长单位脉冲响应序列代替无限长单位脉冲响应系统会有一定的误差，但可通过增加单位脉冲响应的长度逼近。

四、频率响应的计算

稳定系统的频率响应就是系统函数在单位圆上的取值，计算系统的频率响应，可将系统函数中的Z变量用 j e代入即可得到。频率响应是一个复函数，其模叫幅度响应，其相角叫相位响应，它反映了输入序列的频谱经系统后所发生的变化规律。从幅频曲线上可直观看到各频率分量的幅度变化情况，从相频曲线上可直观看到各频率分量的相移情况。根据频响曲线分析系统对信号频谱的影响，概念清楚、简单直观，对信号综合也意义重大，但要将一个较复杂的频率响应复函数转化成幅度响应和相位响应并图示，计算量大且容易出错，图示结果也不一定精确。利用MATLAB函数这些问题都迎刃而解。

实例3：利用MATLAB函数计算实例1中离散系统的频率响应并图示。

由系统函数绘制频响曲线的程序如下：

N=100; w=[0:(N-1)]*2*pi/N; %确定频点

z=exp(j*w); %求频点对应的z点

b=[0.3,0.06,0,0]; a=[1,-1.1,0.55,-0.125];