冲动
动态系统的脉冲响应图;脉冲响应数据
语法
描述
脉冲响应图
冲动(绘制动态系统模型对脉冲输入的响应。该模型sys）sys可以是连续的，也可以是离散的。
对连续时间
sys，脉冲输入为狄拉克脉冲δ(t)．对连续时间sys通过直接馈通，冲动忽略无限脉冲在T = 0．对于离散时间
sys用采样时间Ts，该函数计算对长度单位面积脉冲的响应Ts和高度1 / Ts．这个脉冲接近连续时间的狄拉克脉冲δ(t)作为Ts趋于0。
sys可以是SISO或MIMO。对于MIMO系统，该图显示每个I/O通道的脉冲响应。冲动根据系统动力学自动确定仿真的时间步长和持续时间。
例子
动力系统的脉冲响应
用下列传递函数表示连续时间系统的脉冲响应。
对于本例，创建一个特遣部队表示传递函数的模型。您可以类似地绘制其他动态系统模型类型的脉冲响应，例如零极增益(zpk)或状态空间(党卫军)模型。
Sys = tf(4，[1 2 10]);
画出脉冲响应。
冲动(系统)
的冲动图中自动包含一条虚线，表示稳态响应。在MATLAB®图形窗口中，您可以右键单击图形以查看其他脉冲响应特征，如峰值响应和瞬态时间。
离散时间系统的脉冲响应
画出离散时间系统的脉冲响应。该系统的采样时间为0.2 s，由以下状态空间矩阵表示。
A = [1.6 -0.7;1 0];B = [0.5;0);C = [0.1 0.1];D = 0;
创建状态空间模型并绘制其脉冲响应。
sys = ss(A,B,C,D,0.2);冲动(系统)
脉冲响应反映了模型的离散化，因为它显示了每0.2秒计算一次的响应。
特定时间的脉冲响应
检查以下零极增益模型的脉冲响应。
sys = zpk (1, (-0.2 + 3 j, -0.2 3 j], 1) *特遣部队([1],0.05 [1])
sys = (s + 1) ^ 2 ---------------------------- ( s + 0.05) (s ^ 2 + 0.4 + 9.04)连续时间零/钢管/增益模型。
冲动(系统)
默认情况下,冲动选择一个显示响应趋向于的稳态的结束时间。为了更仔细地观察瞬态响应，将脉冲图限制为t= 20秒。
冲动(sys, 20)
或者，您可以指定您想要检查脉冲响应的确切时间，前提是它们被一个恒定的间隔隔开。例如，检查从暂态结束到系统达到稳态的响应。
T = 20:0.2:120;冲动(sys, t)
尽管这个情节开始于t= 20,冲动总是应用脉冲输入在t= 0。
MIMO系统的脉冲响应图
考虑以下二阶状态空间模型:
A = [-0.5572，-0.7814;0.7814,0];B = [1，-1;0,2];C = [1.9691,6.4493];sys = ss(A,B,C,0);
这个模型有两个输入和一个输出，所以它有两个通道:从第一个输入到输出，从第二个输入到输出。每个通道都有自己的脉冲响应。
当你使用冲动，计算各通道的响应。
冲动(系统)
左图为第一个输入通道的脉冲响应，右图为第二个输入通道的脉冲响应。当你使用冲动为了绘制MIMO模型的响应，它生成了一个表示模型的所有I/O通道的图数组。例如，创建一个随机状态空间模型，它有五个状态、三个输入和两个输出，并绘制它的脉冲响应。
Sys = rss(5,2,3);冲动(系统)
在MATLAB图形窗口中，可以通过右键单击绘图并选择将绘图限制为通道的子集I / O选择器．
比较多个系统的脉冲响应
冲动允许您在同一轴上绘制多个动态系统的响应。例如，比较PI控制器和PID控制器系统的闭环响应。创建系统的传递函数并调优控制器。
H = tf(4，[1 2 10]);C1 = pidtune(H，“π”）;C2 = pidtune(H，“PID”）;
建立闭环系统并绘制它们的脉冲响应。
sys1 = feedback(H*C1,1);sys2 = feedback(H*C2,1);冲动(sys1 sys2)传说(“π”，“PID”，“位置”，“东南”）
默认情况下,冲动为您绘制的每个系统选择不同的颜色。属性指定颜色和线条样式LineSpec输入参数。
冲动(sys1“r——”sys2,“b”)传说(“π”，“PID”，“位置”，“东南”）
第一个LineSpec“r——”为PI控制器的响应指定一条红色虚线。第二个LineSpec“b”为PID控制器的响应指定一条蓝色实线。图例反映了指定的颜色和线条样式。有关更多情节定制选项，请使用impulseplot．
模型阵列中系统的脉冲响应
这个例子比较多个系统的脉冲响应展示如何在单个轴上绘制多个单独系统的响应。当你有多个动态系统排列在一个模型数组中，冲动同时画出他们所有的反应。
创建一个模型数组。对于这个例子，使用具有不同固有频率的二阶传递函数的一维数组。首先，为模型数组预分配内存。下面的命令创建一个1 × 5的零增益SISO传递函数行。前两个维度表示模型输出和输入。剩下的维度是数组维度。
Sys = tf(0 (1,1,1,5));
填充数组。
W0 = 1.5:1:5.5;%固有频率Zeta = 0.5;%阻尼常数为i = 1:长度(w0)系统(:,:,1,我)=特遣部队(w0 (i) ^ 2,[1 2 *ζ* w0 (i) w0 (i) ^ 2]);结束
(有关模型数组和如何创建它们的更多信息，请参见模型阵列)。绘制数组中所有模型的脉冲响应。
冲动(系统)
冲动对数组中所有条目的响应使用相同的线型。项之间的区分方法之一是使用SamplingGrid属性将数组中的每个项与相应的项关联w0价值。
sys。SamplingGrid = struct(“频率”, w0);
现在，当您在MATLAB图形窗口中绘制响应时，您可以单击一条跟踪来查看它对应的频率值。
脉冲响应数据
当你给它一个输出参数时，冲动返回一个响应数据数组。对于SISO系统，响应数据作为长度等于响应采样的时间点数量的列向量返回。你可以提供时间点的向量t，或者允许冲动根据系统动力学为您选择时间点。例如，提取SISO系统在101个时间点之间的脉冲响应t= 0和t= 5秒。
Sys = tf(4，[1 2 10]);T = 0:05:5;Y =冲量(sys,t);大小(y)
ans =1×2101年1
对于MIMO系统，响应数据以维数组的形式返回N——- - - - - -纽约——- - - - - -ν,在那里纽约而且ν是动态系统的输出和输入的数量。例如，考虑下面的状态空间模型，它表示一个双输入一输出系统。
A = [-0.5572，-0.7814;0.7814,0];B = [1，-1;0,2];C = [1.9691,6.4493];sys = ss(A,B,C,0);
提取系统之间200个时间点的脉冲响应t= 0和t= 20秒。
T = linspace(0,20,200);Y =冲量(sys,t);大小(y)
ans =1×3200 1 2
y (:, i, j)列向量包含脉冲响应吗j的Th输入我泰晤士报的输出t．例如，提取第二个输入到输出的脉冲响应。
Y12 = y(:，1,2);情节(t,日元)
具有置信区域的已识别模型的脉冲响应
比较一个参数识别模型的脉冲响应与一个非参数(经验)模型。同时查看他们的3个 地区的信心。
加载数据。
负载iddata1z1
估计一个参数模型。
Sys1 = sest(z1,4);
估计一个非参数模型。
Sys2 =冲量(z1);
绘制脉冲响应图以便比较。
T = (0:0.1:10)';[y1， ~， ~， ysd1] =脉冲(sys1,t);[y2， ~， ~， ysd2] =脉冲(sys2,t);情节(t, y1,“b”， t, y +3*ysd1，”乙:“， t, y -3*ysd1，”乙:“)举行在情节(t, y2,‘g’t, y2+3*ysd2，“旅客:”， t, y2-3*ysd2，“旅客:”）
已识别时间序列模型的脉冲响应
计算一个确定的时间序列模型的脉冲响应。
时间序列模型，也称为信号模型，是没有测量输入信号的模型。该模型的脉冲图使用其(未测量的)噪声通道作为脉冲信号应用的输入通道。
加载数据。
负载iddata9；
估计一个时间序列模型。
Sys = ar(z9,4);
sys是模型的形式吗A y(t) = e(t),在那里e (t)表示噪声通道。为了计算脉冲响应，e (t)作为输入通道，并命名为e@y1．
画出脉冲响应。
冲动(系统)
输入参数
输出参数
限制
非零连续系统的脉冲响应D矩阵是无穷大的t＝0．
冲动忽略这个不连续并返回较低的连续性值Cb在t＝0．的
冲动命令在具有内部延迟的连续时间模型上不起作用。对于这样的模型，请使用pade为了在计算脉冲响应之前近似地计算时间延迟。
提示
当您需要额外的情节定制选项时，请使用
impulseplot代替。要模拟系统对任意输入信号的响应，使用
lsim．
算法
连续时间模型首先转换为状态空间形式。单输入状态空间模型的脉冲响应
是否等价于以下具有初始状态的非强制响应b．
为了模拟这种响应，在输入端使用零阶保持器对系统进行离散化。采样时间是根据系统动力学自动选择的，除非是一个时间向量t = 0:dt:Tf供应。因此,dt作为采样时间。
版本历史
R2006a之前介绍
您也可以从以下列表中选择一个网站: