连续系统离散化分析
1
实验一 离散系统的分析
一 实验目的
1．学习利用采样控制理论；
2．使用MATLAB 理论进行分析；
3. 学习利用z 变换与反变换分析离散控制系统；
二、实验步骤
1．开机执行程序
C ：＼matlab ＼bin ＼matlab.exe (或用鼠标双击图标)进人MATLAB 命令窗口;
2．运用所学自动控制理论z 变换与反变换，使用MATLAB 的基本知识分析离散控制系统的基本性质及进行控制系统的设计。
3. MATLAB 离散系统基本命令
模型转换
1）连续系统离散化
sysd=c2d(sys,T) T 为采样时间
sysd=c2d(sys,T,method)
method 有四种模式：
a. ‘zoh’---采用零阶保持器，
b. ‘foh’---采用一阶保持器，
c. ‘tustin’---采用双线性逼近（tustin ）方法，
d. ‘preqarp’---采用改进的（tustin ）方法，
2）离散系统连续化
sys=d2c(sysd,T,method) T 为采样时间
例 设)
1(1)(+=s s s g ， T=0.1s , 求G(z) 键入命令：sys=tf([1],[1 1 0]);
c2d(sys,0.1) %采样时间0.1s
得到离散传递函数：
当采样时间取T=1s 时：
0.004837 z + 0.004679 G （z ）= ---------------------------- z^2 - 1.905 z + 0.9048 0.3679 z + 0.2642 G （z ）= ---------------------------- z^2 - 1.368 z + 0.3679
2
例 系统脉冲传递函数为6322
.02644.03678.0)()(2+?+=z z z z R z C ，求离散单位阶跃响应。 解. 在MATLAB 窗口键入以下程序
num=[0.3678 0.2644];
den=[1 -1 0.6322];
dstep(num,den)
结果件下图：
三、实验内容
有系统
1）其中K ＝10，T ＝0.25s,求单位阶跃函数r(t)=1 (t)作用下的响应。
从图中求出：上升时间tr,峰值时间tp,超调量σ%，在误差±5%状态下的调整时间ts 各为多少？
2）其中K ＝10，T ＝0.5s,求单位阶跃函数r(t)=1 (t)作用下的响应。
从图中求出：上升时间tr,峰值时间tp,超调量σ%，在误差±5%状态下的调整时间ts 各为多少？
四、实验报告要求
(1) 根据离散控制原理设计并实验，认真作好实验、及记录；
(2) 完成上题给定要求。
实验13 图1 实验系统结构图
离散系统频域分析及matlab实现
《数字信号处理》 课程设计报告 离散系统的频域分析及matlab实现 专业：通信工程 班级：通信11级 组次： 姓名及学号： 姓名及学号：
离散系统的频域分析及matlab 实现 一、设计目的 1.熟悉并掌握matlab 软件的使用； 2.掌握离散系统的频域特性； 3.学会分析离散系统的频域特性的方法； 二、设计任务 1.设计一个系统函数系统的频率响应进行分析； 2.分析系统的频域响应； 3.分析系统的因果稳定性； 4.分析系统的单位脉冲响应； 三、设计原理 1. 系统函数 对于离散系统可以利用差分方程，单位脉冲响应，以及系统函数对系统进行描述。 在本文中利用系统函数H(z)进行描述。若已知一个差分方程为 ∑∑==---=M i N i i i i n y a i n x b n 0 1 )()()(y ，则可以利用双边取Z 变换，最终可以得到系统函数的一 般式H(z)，∑∑=-=-== N i i i M i i i z a z b z X z z H 0 0) () (Y )(。若已知系统的单位脉冲响应，则直接将其进行Z 变换就可以得到系统函数H(z)。系统函数表征系统的复频域特性。 2.系统的频率响应： 利用Z 变化分析系统的频率响应：设系统的初始状态为零，系统对输入为单位脉冲序列 ） （n δ的响应输出称为系统的单位脉冲响应h （n ）。对h(n)进行傅里叶变换，得到： ∑∞ ∞∞-==-)(jw n j |)(|)(e H w j n n j e e H e n h ?ω） （
其中|)(|jwn e H 称为系统的幅频特性函数，)(ω?称为系统的相位特性函数。)(jw e H 表示的是系统对特征序列jwn e 的响应特性。对于一个系统输入信号为n )(ωj e n x =,则系统的输出信号为jwn e )(jw e H 。由上可以知道单频复指数信号jwn e 通过频率响应函数为)(jw e H 后，输出仍为单频复指数信号，其幅度放大了|)(|jw e H ，相移为)(ω?。 对于系统函数H(z)与H(w)之间，若系统函数H(z)的收敛域包含单位圆|z|=1，则有 jw e z jw z H e H ==|)()(，在MATLAB 中可以利用freqz 函数计算系统的频率响应。 （1）[h,w]=freqz(b,a,n) 可得到n 点频率响应，这n 个点均匀地分布在上半单位圆（即 ），并将这n 点频率记录在w 中，相应的频率响应记录在h 中。n 最好能取2的幂次方，如果缺省，则n=512。 （2）[h,w]=freqz(b,a,n,'whole') 在 之间均匀选取n 个点计算频率响应。 （3）[h,w]=freqz(b,a,n,Fs) Fs 为采样频率（以Hz 为单位），在0~Fs/2频率范围内选取n 个频率点，计算相应的频率响应。 （4）[h,w]=freqz(b,a,n,'whole',Fs) 在0~Fs 之间均匀选取n 个点计算频率响应。 （5）freqz(b,a) 可以直接得到系统的幅频和相频特性曲线。其中幅频特性以分贝的形式给出，频率特性曲线的横轴采用的是归一化频率，即Fs/2=1。 3.系统的因果性和稳定性 3.1因果性 因果系统其单位脉冲响应序列h(n)一定是一个因果序列，其z 域的条件是其系统函数H(z)的收敛域一定包含∞，即∞点不是极点，极点 分布在某个圆内，收敛域在某个圆外。 3.2稳定性 系统稳定就要求∞<∑∞ ∞-|h(n)|，由序列的)(jw e H 存在条件和jw e z jw z H e H ==|)()(可以知道 系统稳定的z 域条件就是H(z)的收敛域包含单位圆，即极点全部分布在单位圆内部。 由上3.1和3.2可知，利用系统的零极点分布图可以判断系统的因果性和稳定性。 若在零极点分布图中，若系统的极点都分布在单位圆内，则此系统是因果系统，若有极点分布在单位圆 外，则此系统是非因果系统。在MATLAB 中可以利用zplane 函数画出系统的零极点分布图。系统函数的零极点图的绘制：zplane(b,a)。其中b 为系统函数的分子，a 为系统函数的分母。 4.系统的单位脉冲响应 设系统的初始状态为零，系统对输入为单位脉冲序列）（n δ的响应输出称为系统的单位脉冲响应h （n ）。对于离散系统可以利用差分方程，单位脉冲响应，以及系统函数对系统进行描述。单位脉冲响应是系统的一种描述方法，若已知了系统的系统函数，可以利用系统得出系统的单位脉冲响应。在MATLAB 中利用impz 由函数函数求出单位脉冲响应h(n)。
离散时间系统特性分析
实验五实验报告 实验名称：离散时间系统特性分析
一、实验目的： 1 。深入理解单位样值响应，离散系统的频率响应的概念； 2。 掌握通过计算机进行求得离散系统的单位样值响应，以及离散系统的频率 响应的方法。 二、实验原理： 对于离散系统的单位样值而言，在实际处理过程中，不可能选取无穷多项的取值。往往是选取有限项的取值，当然这里会产生一个截尾误差，但只要这个误差在相对小一个范围里，可以忽略不计。 另外，在一些实际的离散系统中，往往不是事先就能得到描述系统的差分方程的，而是通过得到系统的某些相应值，则此时系统的分析就需借助计算机的数值处理来进行，得到描述系统的某些特征，甚至进而得到描述系统的数学模型。 本实验首先给出描述系统的差分方程，通过迭代的方法求得系统的单位样值响应，进而求得该离散系统的频率响应。限于试验条件，虽然给出了系统方程，但处理的方法依然具有同样的实际意义。 具体的方法是： 1 在给定系统方程的条件下，选取激励信号为δ(n)，系统的起始状态为零 状态，通过迭代法，求得系统的单位样值响应h(n)（n=0,…,N ）。 2 利用公式 其中Ω的取值范围为0～2π 。计算系统的频率响应。 三、实验内容 1 已知系统的差分方程为 利用迭代法求得系统的单位样值响应，取N ＝10。 2 利用公式 其中
#include
实验三___离散时间系统的时域分析
实验三 离散时间系统的时域分析 1.实验目的 （1）理解离散时间信号的系统及其特性。 （2）对简单的离散时间系统进行分析，研究其时域特性。 （3）利用MATLAB对离散时间系统进行仿真，观察结果，理解其时域特性。 2.实验原理 离散时间系统，主要是用于处理离散时间信号的系统，即是将输入信号映射成的输出的某种运算，系统的框图如图所示： （1）线性系统 线性系统就是满足叠加原理的系统。如果对于一个离散系统输入信号为时，输出信号分别为，即：。 而且当该系统的输入信号为时，其中a,b为任意常数，输出为，则该系统就是一个线性离散时间系统。 （2）时不变系统 如果系统的响应与激励加于系统的时刻无关，则该系统是时不变系统。对于一个离散时间系统，若输入，产生输出为，则输入为，产生输出为，即： 若，则。 通常我们研究的是线性时不变离散系统。 3.实验内容及其步骤 （1）复习离散时间系统的主要性质，掌握其原理和意义。 （2）一个简单的非线性离散时间系统的仿真 系统方程为： x = cos(2*pi*0.05*n); x1[n] = x[n+1] x2[n] = x[n] x3[n] = x[n-1] y = x2.*x2-x1.*x3; 或者：y=x*x- x[n+1]* x[n-1] 是非线性。 参考：% Generate a sinusoidal input signal clf; n = 0:200; x = cos(2*pi*0.05*n); % Compute the output signal x1 = [x 0 0]; % x1[n] = x[n+1] x2 = [0 x 0]; % x2[n] = x[n] x3 = [0 0 x]; % x3[n] = x[n-1]
离散控制系统分析方法
实验二 离散控制系统分析方法 一、实验目的 利用MATLAB 对各种离散控制系统进行时域分析。 二、实验指导 1．控制系统的稳定性分析 由前面章节学习的内容可知，对线性系统而言，如果一个连续系统的所有极点都位于s 平面的左半平面，则该系统是一个稳定系统。对离散系统而言，如果一个系统的全部极点都位于z 平面的单位圆内部，则该系统是一个稳定系统。一个连续的稳定系统，如果所有的零点都位于s 平面的左半平面，即所有零点的实部小于零，则该系统是一个最小相位系统。一个离散的稳定系统，如果所有零点都位于z 平面的单位圆内，则称该系统是一个最小相位系统。由于Matlab 提供了函数可以直接求出控制系统的零极点，所以使用Matlab 判断一个系统是否为最小相位系统的工作就变得十分简单。 2．控制系统的时域分析 时域分析是直接在时间域对系统进行分析。它是在一定输入作用下，求得输出量的时域表达式，从而分析系统的稳定性、动态性能和稳态误差。这是一种既直观又准确的方法。 Matlab 提供了大量对控制系统的时域特征进行分析的函数，适用于用传递函数表示的模型。其中常用的函数列入表1，供学生参考。 例1．z z z H 5.05 .1)(2+= 试绘出其单位阶跃响应及单位斜波输入响应。 解：为求其单位阶跃响应及单位斜波输入响应，编制程序如下： num=[1.5]; den=[1 0.5 0];sysd=tf(num,den,0.1) [y,t,x]=step(sysd);
subplot(1,2,1) plot(t,y); xlabel('Time-Sec'); ylabel('y(t)'); gtext('单位阶跃响应') grid; u=0:0.1:1; subplot(1,2,2) [y1,x]=dlsim(num,den,u); plot(u,y1) xlabel('Time-Sec'); ylabel('y(t)'); gtext('单位速度响应') grid 二、实验内容 1、MATLAB在离散系统的分析应用 对于下图所示的计算机控制系统结构图1，已知系统采样周期为T=0.1s，被 控对象的传递函数为 2 () s(0.11)(0.05s1) G s s = ++ ，数字控制器 0.36 () 0.98 z D z z - = + ，试 求该系统的闭环脉冲传递函数和单位阶跃响应。 图1 计算机控制系统结构图 实验步骤： 1）.求解开环脉冲传递函数，运用下面的matlab语句实现：>> T=0.1; >> sys=tf([2],[0.005 0.15 1 0]); %将传函分母展开>> sys1=c2d(sys,T,'zoh'); >> sys2=tf([1 -0.36],[1 0.98],0.1); >> sys3=series(sys2,sys1) 执行语句后，屏幕上显示系统的开环脉冲传递函数为： sys3 = 0.03362 z^3 + 0.05605 z^2 - 0.01699 z - 0.002717 --------------------------------------------------
离散系统稳定性分析
实验一 离散系统稳定性分析 实验学时：2 实验类型：常规 实验要求：必作 一、实验目的： （1）掌握利用MATLAB 绘制系统零极点图的方法； （2）掌握离散时间系统的零极点分析方法； （3）掌握用MATALB 实现离散系统频率特性分析的方法； （4）掌握逆Z 变换概念及MATLAB 实现方法； （5）掌握用MATLAB 分析离散系统稳定性。 二、实验原理： 1、离散系统零极点图及零极点分析； 线性时不变离散系统可用线性常系数差分方程描述，即 ()()N M i j i j a y n i b x n j ==-= -∑∑ （8-1） 其中()y k 为系统的输出序列，()x k 为输入序列。 将式（8-1）两边进行Z 变换的 00 ()()()() () M j j j N i i i b z Y z B z H z X z A z a z -=-== = = ∑∑ （8-2） 将式（8-2）因式分解后有： 11 () ()() M j j N i i z q H z C z p ==-=- ∏∏ （8-3） 其中C 为常数，(1,2,,)j q j M = 为()H z 的M 个零点，(1,2,,)i p i N = 为()H z 的N 个极点。 系统函数()H z 的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。 因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性： ● 系统单位样值响应()h n 的时域特性； ● 离散系统的稳定性；
离散系统的频率特性； 1.1、零极点图的绘制 设离散系统的系统函数为 ()()() B z H z A z = 则系统的零极点可用MA TLAB 的多项式求根函数roots()来实现，调用格式为： p=roots(A) 其中A 为待根求多项式的系数构成的行矩阵，返回向量p 则是包含多项式所有根的列向量。如多项式为231()4 8 B z z z =+ + ，则求该多项式根的MA TLAB 命令为为： A=[1 3/4 1/8]; P=roots(A) 运行结果为： P = -0.5000 -0.2500 需注意的是，在求系统函数零极点时，系统函数可能有两种形式：一种是分子、分母多项式均按z 的降幂次序排列；另一种是分子、分母多项式均按1z -的升幂次序排列。这两种方式在构造多项式系数向量时稍有不同。 （1）()H z 按z 的降幂次序排列：系数向量一定要由多项式最高次幂开始，一直到常数项，缺项要用0补齐；如 3 4 3 2 2()3221 z z H z z z z z += ++++ 其分子、分母多项式系数向量分别为A=[1 0 2 0]、B=[1 3 2 2 1]。 （2）()H z 按1z -的升幂次序排列：分子和分母多项式系数向量的维数一定要相同，不足的要用0补齐，否则0z =的零点或极点就可能被漏掉。如 1 1 2 12()11124 z H z z z ---+= + + 其分子、分母多项式系数向量分别为A=[1 2 0]、B=[1 1/2 1/4]。 用roots()求得()H z 的零极点后，就可以用plot()函数绘制出系统的零极点图。下面是求系统零极点，并绘制其零极点图的MA TLAB 实用函数ljdt()，同时还绘制出了单位圆。 function ljdt(A,B) % The function to draw the pole-zero diagram for discrete system p=roots(A); %求系统极点 q=roots(B); %求系统零点 p=p'; %将极点列向量转置为行向量
离散时间系统的时域分析
第七章离散时间系统的时域分析 §7-1 概述 一、离散时间信号与离散时间系统 离散时间信号：只在某些离散的时间点上有值的 信号。 离散时间系统：处理离散时间信号的系统。 混合时间系统：既处理离散时间信号，又处理连 续时间信号的系统。 二、连续信号与离散信号 连续信号可以转换成离散信号，从而可以用离散时间系统（或数字信号处理系统）进行处理： 三、离散信号的表示方法：
1、 时间函数：f(k)<——f(kT)，其中k 为序号，相当于时间。 例如：)1.0sin()(k k f = 2、 （有序）数列：将离散信号的数值按顺序排列起来。例如： f(k)={1,0.5,0.25,0.125,……,} 时间函数可以表达任意长（可能是无限长）的离散信号，可以表达单边或双边信号，但是在很多情况下难于得到；数列的方法表示比较简单，直观，但是只能表示有始、有限长度的信号。 四、典型的离散时间信号 1、 单位样值函数：? ??==其它001)(k k δ 下图表示了)(n k ?δ的波形。
这个函数与连续时间信号中的冲激函数 )(t δ相似，也有着与其相似的性质。例如： )()0()()(k f k k f δδ=， )()()()(000k k k f k k k f ?=?δδ。 2、 单位阶跃函数：? ??≥=其它001)(k k ε 这个函数与连续时间信号中的阶跃函数)(t ε相似。用它可以产生（或表示）单边信号（这里称为单边序列）。 3、 单边指数序列：)(k a k ε
比较：单边连续指数信号：)()()(t e t e t a at εε=，其 底一定大于零，不会出现负数。 (a) 0.9a = (d) 0.9a =? (b) 1a = (e) 1a =? (c) 1.1a = (f) 1.1a =?
连续传递函数离散化的方法与原理
目录
第一章 模拟化设计基础 数字控制系统的设计有两条道路，一是模拟化设计，一是直接数字设计。如果已经有成熟的模拟控制器，可以节省很多时间和部分试验费用，只要将模拟控制器离散化即可投入应用。如果模拟控制器还不存在，可以利用已有的模拟系统的设计经验，先设计出模拟控制器，再进行离散化。 将模拟控制器离散化，如果用手工进行，计算量比较大。借助数学软件MATLAB 控制工具箱，可以轻松地完成所需要的全部计算步骤。如果需要的话，还可以使用MATLAB 的SIMULINK 工具箱，进行模拟仿真。 第一节 步骤 步骤1 模拟控制器的处理 在数字控制系统中，总是有传输特性为零阶保持器的数模转换器（DAC ），因此，如果模拟控制器尚未设计，则应以下 图的方式设计模拟控制器，即在对象前面加上一个零阶保持器，形成一个新对象Ts 1e G s s ()--，然后针对这个新对象求模拟 控制器D(s)。事实上，模拟控制器一般是已经设计好的，无法或不方便更改了，离散化后的系统只好作为近似设计了。 然而，按照上述思路，可否将已有的控制器除以一个零阶保持器再离散化呢？还没有这方面的实际经验。 以下假设选定的G(s)，D(s)如下图，而且不对G(s)作添加保持器的预处理。 步骤2 离散化模拟控制器 离散化模拟控制器之前，先要确定离散化算法和采样时间。离散化算法有好几种，第二章中有详细的论述，现假定采用双线性变换法。确定采样时间，需要考虑被控对象的特性，计算机的性能，以及干扰信号的影响等，初步可按采样时间T<，Tp 为被控对象时间常数，或T=～τ，为被控对象的纯滞后，初步确定后再综合平衡其它因素，当然这需要一定的经验，现在假定取秒。 假设模拟控制器为s 2 D s 8s 15 +=?+()，在MATLAB 中，用c2d 函数进行离散化，过程为： 转换结果为： 步骤3 检验数字控制器的性能 数字控制器的性能项目比较多，我们仅以直流增益，频率特性，零极点分布说明。 直流增益 dcgain(dz) 返回直流增益 频率特性 bode(ds,'r',dz,'g') 伯德图，见下页左图 零极点分布 pzmap(dz) 零极点分布图，见下页右图 步骤4 离散化控制对象 为了进行模拟仿真，需要对控制对象进行离散化，由于步骤1所说的原因，应把被控对象视为零阶保持器与原对象的串连，即应对 Ts 1e G s s ()--进行离散化，这时可在c2d 函数中使用零阶保持器(zoh)方法，如果认为不需要添加零阶保持器，即直接对G(s)离散化，则应在c2d 函数中使用冲击响应不变法（imp ）。 借用零阶保持器(zoh)方法，将对象20 G s s s 2()() =+带一阶保持器离散化的过程如下： 转换结果为： 步骤5 模拟仿真 求离散系统的闭环传递函数和连续系统的闭环传递函数。 ds=zpk(-2,-15,8) %建立模拟控制器的s 传递函数 dz=c2d(ds,,'tustin') %将模拟控制器按tustin 方法转换为z 传递函数的数字控制器 ...... %模拟控制器D(s)转换为D(z)的过程见前 gs=zpk([ ],[0,-2],20) %建立对象的s 传递函数 g1z=c2d(gs,,'zoh') %借用c2d 函数进行带零阶保持器的对象的离散化
离散时间系统的分析
课程设计报告 课程设计题目：离散时间系统分析学号：201420130206 学生姓名：董晓勇 专业：通信工程 班级：1421301 指导教师：涂其远 2015年12月18日
离散时间系统的分析 一、设计目的和意义 1 . 目的： （1）深刻理解卷积和、相加、相乘运算，掌握求离散序列卷积和、相加相乘的计算方法；（2）加深理解和掌握求离散序列Z变换的方法； （3）加深和掌握离散系统的系统函数零点、函数极点和系统时域特性、系统稳定性的关系。 2 . 意义： 在对《信号与系统》一书的学习中，进行信号与系统的分析是具有十分重要的意义，同时也是必不可少的。利用matlab函数，只需要简单的编程，就可以实现系统的时域、频域分析，对系统特性进行分析，为实际的系统设计奠定了基础。本设计在离散系统Z域分析理论的基础上，利用matlab对离散系统的稳定性和频域响应进行了分析。 二、设计原理
第一部分：对离散时间系统的时域进行分析呈 对离散时间信号的代数运算（相加、相乘、卷积和），是在时域进行分析。相加用“+”来完成，相乘用“·*”来完成，卷积和则用conv 函数来实现，具体形式为y=conv(x1,x2,….)，其中x1,x2,…..为输入的离散序列 ，y 为输出变量。 在零初始状态下，matlab 控制工具箱提供了一个filter 函数，可以计算差分方程描述的系统的响应，其调用形式为： y=filter(b,a,f) 其中，a=[a0,a1,a2,…]、b=[b0,b1,b2,….]分别是系统方程左、右边的系数向量，f 表示输入向量，y 表示输出向量。 第二部分：对离散时间系统的Z 域进行分析 matlab 工具箱提供了计算Z 正变换的函数ztrans,其调用形式为： F=zrtans(f) %求符号函数f 的Z 变换，返回函数的自变量为z 。 Matlab 的zplane 函数用于系统函数的零极点图的绘制，调用方式为： zplane(b,a)其中，b 、a 分别为系统函数分子、分母多项式的系数向量。 matlab 中，利用freqz() 函数可方便地求得系统的频率响应，调用格式为： freqz(b,a,N) 该调用方式将绘制系统在0~PI 范围内N 个频率等分点的幅频特性和相频特性图。 三、 详细设计步骤 1.自己设计两个离散时间序列x1、x2，对其进行相加，相乘，卷积运算，并显示出图形。 2.根据已知的LTI 系统：y[n]-0.7y[n-1]-0.6y[n-2]+y[n-3]=x[n]+0.5[n-1]，得其在Z 域输 入输出的传递函数为： 1 12310.5()10.70.6z H z z z z ----+= --+ 利用matlab 求：（1）系统函数的零点和极点，并在z 平面显示他们的分布；（2）画出幅频响应和相频响应的特性曲线。 四、 设计结果及分析 (1).自行设计产生两个离散序列信号，对其进行相加、乘及卷积运算
离散系统的时域分析实验报告
实验2 离散系统的时域分析 一、实验目的 1、熟悉并掌握离散系统的差分方程表示法； 2、加深对冲激响应和卷积分析方法的理解。 二、实验原理 在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下: 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号, 记系统单位冲激响应，则系统响应为如下的卷积计算式： . .
. . 当 时，h[n]是有限长度的（），称系统为FIR 系统；反之， 称系统为IIR 系统。 三、实验容 1、用MATLAB 求系统响应 1） 卷积的实现 线性移不变系统可由它的单位脉冲响应来表征。若已知了单位脉冲响应和系统激励就可通过卷积运算来求取系统响应，即)(*)()(n h n x n y 程序： x=input(‘Type in the input sequence=’); %输入x h=input(‘Type in the impulse response sequence=’); %输入h y=conv(x,h); % 对x ，h 进行卷积 N=length(y)-1; %求出N 的值 n=0:1:N; %n 从0开始，间隔为1的取值取到N 为止 disp(‘output sequence=’); disp(y); %输出y stem(n,y); %画出n 为横轴，y 为纵轴的离散图 xlabel(‘Time index n ’); ylable(‘Amplitude ’); % 规定x 轴y 轴的标签 输入为： x=[-2 0 1 -1 3] h=[1 2 0 -1] 图形: 2） 单位脉冲响应的求取 线性时不变因果系统可用MA TLAB 的函数filter 来仿真 y=filter(b,a,x);
连续和离散系统分析
实验一连续与离散系统分析 一、实验目得 学习连续系统与离散系统响应得matlab求解方法; 二、实验主要仪器设备与材料 计算机 三、实验方法、步骤及结果测试 实验方法:编程,上机调试,分析实验结果; 步骤: 编程实现上述各实验内容 四、实验结果 1、某系统得传递函数为： 试求系统得冲激响应与阶跃响应。 2、编制程序求解下列两个系统得单位冲激响应与阶跃响应,并绘出其图形。要求
分别用filtｅr、coｎｖ、ｉｍpz三种函数完成。给出理论计算结果与程序计算结果并讨论。 (I) 理论计算结果: 程序计算结果: Ａ：单位冲激响应 （１)用Fiｌter函数(２)用Conｖ函数 (3）用impz函数 单位冲激响应: n 0 1 2 3 4 5 h(n) 1 -1、75 1、19 -0、67 0、355 -0、18 单位阶跃响应: n 0 1 2 3 4 5 y(n) 1 -0、75 0、44 -0、234 0、12 -0、06
Ｂ:单位阶跃响应(1)用Fil ｔer 函数 (2)用Ｃｏnv 函数 (3)用Imp ｚ函数 (II ） 理论计算结果: 程序计算结果： A:单位冲激响应(1)用f ｉlter 函数 单位冲激响应: n 0 1 2 3 4 5 h(n) 0 0、25 0、25 0、25 0、25 单位阶跃响应: N 0 1 2 3 4 5 y(n) 0 0、25 0、5 0、75 1 1
（２)用Conｖ函数 （3)用Impｚ函数 B:单位阶跃响应 （1）用fｉlter函数 （2)用Ｃｏｎv函数 (3)用Impz函数
连续系统离散化处理基本方法
在数字计算机上对连续系统进行仿真时，首先遇到的问题是如何解决数字计算机在数值及时间上的离散性与被仿真系统数值及时间上的连续性这一基本问题。 从根本意义上讲，数字计算机所进行的数值计算仅仅是“数字”计算，它表示数值的精度受限于字长，这将引入舍入误差；另一方面，这种计算是按指令一步一步进行的，因而，还必须将时间离散化，这样就只能得到离散时间点上系统性能。用数字仿真的方法对微分方程的数值积分是通过某种数值计算方法来实现的。任何一种计算方法都只能是原积分的一种近似。因此，连续系统仿真，从本质上是对原连续系统从时间、数值两个方面对原系统进行离散化，并选择合适的数值计算方法来近似积分运算，由此得到的离散模型来近似原连续模型。如何保证离散模型的计算结果从原理上确能代表原系统的行为，这是连续系统数字仿真首先必须解决的问题。 设系统模型为：),,(t u y f y =&，其中u (t )为输入变量，y (t )为系统变量；令仿真时间间隔为h ，离散化后的输入变量为)(?k t u ，系统变量为)(?k t y ，其中k t 表示t=kh 。如果)()(?k k t u t u ≈，)()(?k k t y t y ≈，即0)()(?)(≈-=k k k u t u t u t e ，0)()(?)(≈-=k k k y t y t y t e （对所有k=0,1,2,…），则可认为两模型等价，这称为相似 原理（参见图）。 实际上，要完全保证0)(,0)(==k y k u t e t e 是很困难的。进一步分析离散化引的误差，随着计算机技术的发展，由计算机字长引入的舍入误差可以忽略，关键是数值积分算法，也称为仿真建模方法。相似原理用于仿真时，对仿真建模方法有三个基本要求： （1）稳定性：若原连续系统是稳定的，则离散化后得到的仿真模型也应是稳定的。关于稳定性的详细讨论将在节中进行。 （2）准确性：有不同的准确性评价准则，最基本的准则是： 绝对误差准则：δ≤-=)()(?)(k k k y t y t y t e 相对误差准则：δ≤-= )(?)()(?)(k k k k y t y t y t y t e 其中 规定精度的误差量。 原连续模型 仿真模型 )(≈k y t e 图 相
实验七--离散系统分析的MATLAB实现讲解学习
实验七 离散系统分析的MATLAB 实现 一、实验目的 1、掌握利用MATLAB 绘制系统零极点图的方法； 2、掌握离散时间系统的零极点分析方法； 3、学习离散系统响应的MATLAB 求解方法； 4、掌握用MATALB 实现离散系统频率特性分析的方法； 5、深刻理解离散系统的系统函数零极点对系统频响的影响，可以根据 零极点知识设计简单的滤波器。 二、基本原理 （一）离散系统零极点 线性时不变离散系统可用线性常系数差分方程描述，即 ()()N M i j i j a y n i b x n j ==-=-∑∑ （1） 其中()y k 为系统的输出序列，()x k 为输入序列。 将式（1）两边进行Z 变换， 00 () () ()() () M j j j N i i i b z Y z B z H z X z A z a z -=-== == ∑∑ （2） 将式（2）因式分解后有： 11 () ()() M j j N i i z q H z C z p ==-=-∏∏ （3） 其中C 为常数，(1,2,,)j q j M =L 为()H z 的M 个零点，(1,2,,)i p i N =L 为()H z 的 N 个极点。 系统函数()H z 的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。 （二）离散系统零极点图及零极点分析 1、零极点图的绘制 设离散系统的系统函数为 () ()() B z H z A z =
则系统的零极点可用MATLAB 的多项式求根函数roots()来实现，调用格式为： p=roots(A) 其中A 为待求根多项式的系数构成的行矩阵，返回向量p 则是包含多项式所有根的列向量。如多项式为231 ()48 B z z z =+ +，则求该多项式根的MATLAB 命令为为： A=[1 3/4 1/8]; P=roots(A) 运行结果为： P = -0.5000 -0.2500 需注意的是，在求系统函数零极点时，系统函数可能有两种形式：一种是分子、分母多项式均按z 的降幂次序排列；另一种是分子、分母多项式均按1z -的升幂次序排列。这两种方式在构造多项式系数向量时稍有不同。 （1）()H z 按z 的降幂次序排列：系数向量一定要由多项式最高次幂开始，一直到常数项，缺项要用0补齐。如 34322()3221 z z H z z z z z +=++++ 其分子、分母多项式系数向量分别为A=[1 0 2 0]、B=[1 3 2 2 1]。 （2）()H z 按1z -的升幂次序排列：分子和分母多项式系数向量的维数一定要相同，不足的要用0补齐，否则0z =的零点或极点就可能被漏掉。如 1 1212()11124 z H z z z ---+=++ 其分子、分母多项式系数向量分别为A=[1 2 0]、B=[1 1/2 1/4]。 用roots()求得()H z 的零极点后，就可以用plot()函数绘制出系统的零极点图。下面是求系统零极点，并绘制其零极点图的MATLAB 实用函数ljdt()，同时还绘 制出了单位圆。函数ljdt()的程序如下： function ljdt(A,B) % The function to draw the pole-zero diagram for discrete system p=roots(A) %求系统极点 q=roots(B) %求系统零点 p=p'; %将极点列向量转置为行向量 q=q'; %将零点列向量转置为行向量 x=max(abs([p q 1])); %确定纵坐标范围
离散控制系统的分析与综合
第7章离散控制系统的分析与综合 7.3 离散系统的能控性和能观性 1、离散系统的能控性和能观性判据 ◆能控性和能观性定义： 对有限个采样周期，若能找到控制信号序列，能使任意一个初始状态转移到零状态，则系统是状态完全能控的；若根据有限个采样周期的输出序列，能唯一地确定任意初始状态，则系统是状态完全能观的。 ◆能控性和能观性判据： A B C状态完全能控的充要条件 n阶线性定常离散系统(,,) 是
1 rank rank[,,,]n c Q B AB A B n -== 状态完全能观的充要条件是 1rank rank o n C CA Q n CA -轾犏犏犏==犏犏犏臌 2、连续系统离散化后的能控性与能观性 设具有零阶保持器的n 阶连续系统以采样周期T 离散为离散系统。 定理：若连续系统不能控（不能观），则其离散系统必不能控（不能观）。若连续系统能控（能观），其互异特征值（含 重特征值）为μλλλ，， ， 21，若对一切 μλλ,,2,1,,0][ ==-j i R j i e
的互异特征值满足 ,2,1,2][±±=≠-k T k I j i m πλλ 则其离散系统必保持能控（能观）性。 7.4 离散系统的稳定性 1、离散系统稳定的充要条件 1）赛尔维斯特展开定理 设n 阶系数矩阵A 具有互异特征值n λλλ，，， 21，)(A f 是A 函数，则有 i i n i A f A f )()(1λ∑== 其中 j i i n i j j i I A A λλλ--= ∏≠=,1
2）离散系统稳定的充要条件 线性定常离散系统齐次状态方程 的解为 ()(0)k x k A x = 由系统的特征方程 0zI A -= 可解得系统的特征值。 设A 的特征值n λλλ，，， 21两两互异，则由赛尔维斯特展开定理得 1n k k i i i A λA ==?
实验四-离散时间系统的频域分析(附思考题程序)
实验四 离散时间系统的频域分析 1.实验目的 （1）理解和加深傅里叶变换的概念及其性质。 （2）离散时间傅里叶变换(DTFT)的计算和基本性质。 （3）离散傅里叶变换(DFT)的计算和基本性质。 2.实验原理 对离散时间信号进行频域分析，首先要对其进行傅里叶变换，通过得到的频谱函数进行分析。 离散时间傅里叶变换(DTFT ，Discrete-time Fourier Transform)是傅立叶变换的一种。它将以离散时间nT （其中，T 为采样间隔）作为变量的函数（离散时间信号）f (nT )变换到连续的频域，即产生这个离散时间信号的连续频谱()iw F e ，其频谱是连续周期的。 设连续时间信号f (t )的采样信号为：()()()sp n f t t nT f nT d ￥ =-? = -?，并且其傅里叶变 换为：()()(){}sp n iwt f t f nT t nT dt e d ￥ ￥ -? =-? --= ? òF 。 这就是采样序列f(nT)的DTFT:：()()iwT inwT DTFT n F e f nT e ￥ -=-? = ?,为了方便，通常将采 样间隔T 归一化，则有：()()iw inw DTFT n F e f n e ￥ -=-? = ?，该式即为信号f(n)的离散时间傅 里叶变换。其逆变换为：()1()2iw DTFT inw F e dw f n e p p p -=ò。 离散傅里叶变换（DFT ，Discrete-time Fourier Transform ）是对离散周期信号的一种傅里叶变换，对于长度为有限长信号，则相当于对其周期延拓进行变换。在频域上，DFT 的离散谱是对DTFT 连续谱的等间隔采样。 21 1 20 ()()| ()()DFT k DTFT k w N knT N N i iwT iwnT N n n F w F e f nT e f nT e p p =----==== = 邋 长度为N 的有限长信号x(n)，其N 点离散傅里叶变换为： 1 ()[()]()kn N N n X k DFT x n x n W -=== ?。 X(k)的离散傅里叶逆变换为：10 1()[()]()kn N N k x n IDFT X k X k W N --===?。 DTFT 是对任意序列的傅里叶分析，它的频谱是一个连续函数；而DFT 是把有限长序列作为周期序列的一个周期，对有限长序列的傅里叶分析，DFT 的特点是无论在时域还是频域
离散系统的Z域分析
实验名：离散系统的Z 域分析 一、实验目的 1、掌握离散序列z 变换的计算方法。 2、掌握离散系统系统函数零极点的计算方法和零极点图的绘制方法，并能根据零极点图分析系统的因果性和稳定性。 3、掌握利用MATLAB 进行z 反变换的计算方法。 二、实验原理与计算方法 1、z 变换 离散序列x (n )的z 变换定义为：∑∞ -∞ =-= n n z n x Z X )()(。 在MA TLAB 中可以利用符号表达式计算一个因果序列的z 变换。其命令格式为： syms n; f=(1/2)^n+(1/3)^n; ztrans(f) 2、离散系统的系统函数及因果稳定的系统应满足的条件 一个线性移不变离散系统可以用它的单位抽样响应h (n )来表示其输入与输出关系，即 y (n )= x (n )* h (n ) 对该式两边取z 变换，得： Y (z )= X (z )· H (z ) 则： ) () ()(z X z Y z H = 将H (z )定义为系统函数，它是单位抽样响应h (n )的z 变换，即 ∑∞ -∞ =-= =n n z n h n h Z z H )()]([)( 对于线性移不变系统，若n <0时，h (n )=0，则系统为因果系统；若 ∞<∑∞ -∞ =n n h |)(|，则 系统稳定。由于h (n )为因果序列，所以H (z )的收敛域为收敛圆外部区域，因此H (z )的收敛域为收敛圆外部区域时，系统为因果系统。因为∑∞ -∞ =-= n n z n h z H )()(，若z =1时H (z )收敛，即 ∞<= ∑∞ -∞ ==n z n h z H |)(||)(1，则系统稳定，即H(z)的收敛域包括单位圆时，系统稳定。 因此因果稳定系统应满足的条件为：1,||<∞≤<ααz ，即系统函数H (z )的所有极点全部落在z 平面的单位圆之内。 3、MA TLAB 中系统函数零极点的求法及零极点图的绘制方法 MATLAB 中系统函数的零点和极点可以用多项式求根函数roots ()来实现，调用该函数的命令格式为：p=roots(A)。其中A 为待求根多项式的系数构成的行向量，返回向量p 是包含该多项式所有根位置的列向量。 如：求多项式8 1 43)(2++=z z z A 的根的MA TLAB 命令为： A=[1 3/4 1/8]; p=roots(A) 运行结果为： p= -0.5000 -0.2500 也可以用[z,p,k]=tf2zp(B,A)函数求得。其中z 为由系统的零点构成的向量，p 为由系统的极点构成的向量，k 表示系统的增益；B 、A 分别为系统函数中分子分母多项式的系数向
线性离散系统的数学模型和方法分析
§10-2 线性离散系统的数学模型和分析方法 大多数计算机控制系统可以用线性时不变离散系统的数学模型来描述。对于单输入单输出线性离散系统，人们习惯用线性常系数差分方程或脉冲传递函数来表示。离散系统的线性常系数差分方程和脉冲传递函数，分别和连续系统的线性常系数微分方程和传递函数在结构、性质和运算规则上相类似。对于多变量、时变和非线性系统用状态空间方法处理比较方便。 一、线性离散系统的数学描述 1. 差分方程 对简单的单输入单输出线性离散系统，其输入)(kT u 和输出)(kT y 之间的关系可用下列线性常系数差分方程来表示 )()()()()()(101nT kT u b T kT u b kT u b nT kT y a T kT y a kT y n n -++-+=-++-+ (10.17) (10.17)式也可以写成如下紧缩的形式 ∑∑==-=-+n i n i i i iT kT u b iT kT y a kT y 1 )()()( (10.18) 如果引入后移算子1 -q ，即 )()(1T kT y kT y q -=- (10.19) 则(10.18)式可写成多项式的形式 )()()()(11kT u q B kT y q A --= (10.20) 式中 n n q a q a q A ---+++= 1111)( n n q b q b b q B ---+++= 1101)( 方程(10.17)、(10.18)和(10.20)中假设左右两端阶次相同，这并不失一般性，差分方程中最高和最低指数之差n 被称为差分方程的阶数。如果(10.17)式中右端的系数项i b ，n i ,,1,0 =，不全为零，则此方程被称为非齐次方程。方程右端又被称为驱动项。方程的阶数和系数反映系统的结构特征。用差分方程作为物理系统的数学模型时，方程中各变量代表一定的物理量，其系数有时具有明显的物理意义。如果(10.17)式右端的系数全为零，则被称作齐次方程。齐次差分方程表征了线性离散系统在没有外界作用的情况下，系统的自由运动，它反映了系统本身的物理特性。 2. 差分方程的解 线性常系数差分方程求解方法和线性代数方程的求解相类似，其全解)(kT y 由齐次方程的通解