主要内容gydF4y2Ba

水疗中心gydF4y2Ba

估计使用光谱分析固定频率分辨率的频率响应gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

GgydF4y2Ba= spa (gydF4y2Ba数据gydF4y2Ba)gydF4y2Ba估计频率响应,以及不确定性,和时间——或者频域数据的噪声谱gydF4y2Ba数据gydF4y2Ba。gydF4y2Ba

数据gydF4y2Ba包含时域或频域的输入/输出数据或时间序列数据。gydF4y2Ba数据gydF4y2Ba能的形式吗gydF4y2Ba时间表gydF4y2Ba逗号分隔的数字矩阵,一个数字矩阵,一个gydF4y2BaiddatagydF4y2Ba对象,或一个gydF4y2BaidfrdgydF4y2Ba对象可以包含复杂的值。gydF4y2Ba

如果gydF4y2Ba数据gydF4y2Ba是一个时间序列,即不包含输入,然后呢gydF4y2Ba水疗中心gydF4y2Ba(gydF4y2Ba数据)gydF4y2Ba返回输出功率谱以及不确定性。gydF4y2Ba水疗中心gydF4y2Ba计算光谱在128 0(排除)之间的等距的频率值和π,使用损害窗口。gydF4y2Ba

如果gydF4y2Ba数据gydF4y2Ba的形式是一个时间表,软件解释最后一个变量是单一的输出变量。改变这个解释,使用gydF4y2BaInputNamegydF4y2Ba和gydF4y2BaOutputNamegydF4y2Ba名称变量参数。gydF4y2Ba

例子gydF4y2Ba

GgydF4y2Ba= spa (gydF4y2Ba数据gydF4y2Ba,gydF4y2BawinSizegydF4y2Ba,gydF4y2Ba频率gydF4y2Ba)gydF4y2Ba估计频率响应中包含的频率gydF4y2Ba频率gydF4y2Ba用损害窗口大小gydF4y2BawinSizegydF4y2Ba。gydF4y2Ba

GgydF4y2Ba= spa (gydF4y2Ba数据gydF4y2Ba,gydF4y2BawinSizegydF4y2Ba,gydF4y2Ba频率gydF4y2Ba,gydF4y2Ba最大尺寸gydF4y2Ba)gydF4y2Ba将输入/输出数据分为段,每一段包含少于gydF4y2Ba最大尺寸gydF4y2Ba元素。使用这个语法来提高计算性能。gydF4y2Ba

sysgydF4y2Ba= spa (gydF4y2Ba___gydF4y2Ba,gydF4y2Ba名称,值gydF4y2Ba)gydF4y2Ba使用额外的模型选项指定一个或多个名称参数。gydF4y2Ba

例如,指定变量名称使用输入和输出信号gydF4y2Basys = spa(数据、“InputName”(“u1”、“u3”),“OutputName”,(“日元”、“y4))gydF4y2Ba。gydF4y2Ba

您可以使用该语法与任何以前的输入参数组合。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

估计频率响应的输入/输出数据gydF4y2BaiddatagydF4y2Ba对象gydF4y2Baz3gydF4y2Ba。使用默认的固定分辨率为128等距的对数频率值在0(排除)和之间gydF4y2Ba πgydF4y2Ba 。gydF4y2Ba

负载gydF4y2Baiddata3gydF4y2Baz3gydF4y2Ba;g = spa (z3);波德(g)gydF4y2Ba

图包含2轴对象。坐标轴对象1的标题:u1:日元ylabel级(dB)包含一个类型的对象。这个对象表示g。坐标轴对象2 ylabel阶段(度)包含一个类型的对象。这个对象表示g。gydF4y2Ba

生成对数间隔向量gydF4y2BafgydF4y2Ba。gydF4y2Ba

f = logspace(2π,128);gydF4y2Ba

估计的频率响应输入/输出数据gydF4y2Baz3gydF4y2Ba。指定窗口大小gydF4y2Ba[]gydF4y2Ba获取默认滞后窗口大小。gydF4y2Ba

负载gydF4y2Baiddata3gydF4y2Baz3gydF4y2Ba;g = spa (z3, [], f);gydF4y2Ba

情节预示反应和光谱干扰的置信区间3个标准差。gydF4y2Ba

h = bodeplot (g);showConfidence (h, 3)gydF4y2Ba

图包含2轴对象。坐标轴对象1的标题:u1:日元ylabel级(dB)包含一个类型的对象。这个对象表示g。坐标轴对象2 ylabel阶段(度)包含一个类型的对象。这个对象表示g。gydF4y2Ba

图h = spectrumplot (g);showConfidence (h, 3)gydF4y2Ba

图包含一个坐标轴对象。坐标轴对象的标题:e@y1:日元ylabel权力(dB)包含一个类型的对象。这个对象表示g。gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

均匀采样估计数据,指定为一个时间表,逗号分隔矩阵对,或数据对象,如以下部分描述。gydF4y2Ba

时间表gydF4y2Ba

指定gydF4y2Ba数据gydF4y2Ba作为一个gydF4y2Ba时间表gydF4y2Ba定期使用间隔时间向量。gydF4y2Ba数据gydF4y2Ba包含变量代表输入和输出通道。gydF4y2Ba

当你使用时间表估计,您可以使用所有的变量或指定要使用通道的一个子集。gydF4y2Ba

对于multiexperiment数据,指定数据作为gydF4y2Ba不gydF4y2Ba1单元阵列的时间表,gydF4y2Ba不gydF4y2Ba是实验的数量。所有的样品时间实验必须匹配。gydF4y2Ba

选择单个输入和输出通道用于估计,使用gydF4y2BaInputNamegydF4y2Ba和gydF4y2BaOutputNamegydF4y2Ba名称-值参数。gydF4y2Ba

对于时间序列数据,要么指定gydF4y2Ba数据gydF4y2Ba作为一个变量时间表只包含一个输出变量,或使用gydF4y2BaInputNamegydF4y2Ba和gydF4y2BaOutputNamegydF4y2Ba名称-值参数指定只输出变量估计。gydF4y2Ba

例如,gydF4y2Basys = spa(数据、“OutputName”、“日元”,“InputName”, [])gydF4y2Ba估计时间序列模型gydF4y2BasysgydF4y2Ba从时间表gydF4y2Ba数据gydF4y2Ba使用输出变量gydF4y2Bay2gydF4y2Ba没有输入变量。gydF4y2Ba

逗号分隔矩阵对gydF4y2Ba

指定gydF4y2Ba数据gydF4y2Ba作为一个逗号分隔的矩阵gydF4y2BaugydF4y2Ba,gydF4y2BaygydF4y2Ba含有均匀采样输入和输出时域信号值。依赖于数据没有提供采样时间信息。软件假定样本时间是1秒。不推荐使用依赖于数据为连续时间系统。gydF4y2Ba

对于输出系统,指定gydF4y2BaugydF4y2Ba,gydF4y2BaygydF4y2Ba列向量的长度gydF4y2BaNsgydF4y2Ba,在那里gydF4y2BaNsgydF4y2Ba是样品的数量。gydF4y2Ba

对于MIMO系统,指定gydF4y2BaugydF4y2Ba,gydF4y2BaygydF4y2Ba作为一个矩阵对以下维度:gydF4y2Ba

  • ugydF4y2Ba- - - - - -gydF4y2BaNgydF4y2Ba年代gydF4y2Ba——- - - - - -gydF4y2BaNgydF4y2BaugydF4y2Ba,在那里gydF4y2BaNgydF4y2BaugydF4y2Ba输入的数量吗gydF4y2Ba

  • ygydF4y2Ba- - - - - -gydF4y2BaNgydF4y2Ba年代gydF4y2Ba——- - - - - -gydF4y2BaNgydF4y2BaygydF4y2Ba,在那里gydF4y2BaNgydF4y2BaygydF4y2Ba是输出的数量吗gydF4y2Ba

对于multiexperiment数据,指定gydF4y2BaugydF4y2Ba,gydF4y2BaygydF4y2Ba作为一对1×-gydF4y2BaNgydF4y2BaegydF4y2Ba细胞阵列,gydF4y2BaNgydF4y2BaegydF4y2Ba是实验的数量。gydF4y2Ba

对于时间序列系统,指定一个空gydF4y2BaugydF4y2Ba,也就是说,gydF4y2Ba[]gydF4y2Ba,gydF4y2BaygydF4y2Ba。gydF4y2Ba

数据对象gydF4y2Ba

一个估计的数据对象,指定为时域或频域gydF4y2BaiddatagydF4y2Ba对象或一个gydF4y2BaidfrdgydF4y2Ba对象包含均匀采样输入和输出值。数据对象可以有一个或多个输出渠道和零个或多个输入通道。默认情况下,软件集的样本时间模型的样本时间估计数据。gydF4y2Ba

关于处理估计数据类型的更多信息,见gydF4y2Ba数据域和数据类型在系统辨识工具箱gydF4y2Ba。gydF4y2Ba

损害窗口大小,也称为延迟大小,指定为一个标量整数。默认情况下,窗口大小的函数集gydF4y2Ba分钟(长度(数据)/ 10、30)gydF4y2Ba。gydF4y2Ba

频率估计光谱响应,指定为一个行向量在rad /单位gydF4y2BaTimeUnitgydF4y2Ba,在那里gydF4y2BaTimeUnitgydF4y2Ba指的是gydF4y2BaTimeUnitgydF4y2Ba的属性gydF4y2Ba数据gydF4y2Ba。默认情况下,函数集gydF4y2Ba频率gydF4y2Ba128矢量值的范围(0,π),均匀间隔的对数。对于离散时间模型,设置gydF4y2Ba频率gydF4y2Ba在奈奎斯特频率。gydF4y2Ba

段内的最大大小gydF4y2Ba数据gydF4y2Ba,指定为一个正整数。如果你忽略这个参数,函数执行评估中使用完整的数据集gydF4y2Ba数据gydF4y2Ba而不是细分数据。gydF4y2Ba

名称-值参数gydF4y2Ba

指定可选的双参数作为gydF4y2BaName1 = Value1,…,以=家gydF4y2Ba,在那里gydF4y2Ba的名字gydF4y2Ba参数名称和吗gydF4y2Ba价值gydF4y2Ba相应的价值。名称-值参数必须出现在其他参数,但对的顺序无关紧要。gydF4y2Ba

R2021a之前,用逗号来分隔每一个名称和值,并附上gydF4y2Ba的名字gydF4y2Ba在报价。gydF4y2Ba

例子:gydF4y2Basys = spa(数据、“InputName”、“u2”)gydF4y2Ba

输入通道名称,指定为一个字符串,特征向量,特征向量的字符串数组,数组或单元。gydF4y2Ba

如果您正在使用数据源,时间表中的名称gydF4y2BaInputNamegydF4y2Ba必须的一个子集的时间表变量。gydF4y2Ba

例子:gydF4y2Basys = spa (tt, __, InputName, [u1”“u2”])gydF4y2Ba选择的变量gydF4y2Bau1gydF4y2Ba和gydF4y2Bau2gydF4y2Ba作为输入通道的时间表gydF4y2BattgydF4y2Ba用于估计。gydF4y2Ba

输出通道名称,指定为一个字符串,特征向量,特征向量的字符串数组,数组或单元。gydF4y2Ba

如果您正在使用数据源,时间表中的名称gydF4y2BaOutputNamegydF4y2Ba必须的一个子集的时间表变量。gydF4y2Ba

例子:gydF4y2Basys = spa (tt, __“OutputName”,(“日元”“y3))gydF4y2Ba选择的变量gydF4y2Ba日元gydF4y2Ba和gydF4y2Bay3gydF4y2Ba作为输出通道的时间表gydF4y2BattgydF4y2Ba用于估计。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

频率响应的不确定性和噪声频谱,指定为一个gydF4y2BaidfrdgydF4y2Ba对象。对于时间序列数据,gydF4y2BaGgydF4y2Ba估计谱和标准偏差。gydF4y2Ba

评估结果和信息存储在选择使用gydF4y2Ba报告gydF4y2Ba模型的属性。gydF4y2Ba报告gydF4y2Ba有以下字段。gydF4y2Ba

报告字段gydF4y2Ba 描述gydF4y2Ba
状态gydF4y2Ba

总结模型的状态,这表明模型是由建设或是否得到评估。gydF4y2Ba

方法gydF4y2Ba

评估使用的命令。gydF4y2Ba

windowSizegydF4y2Ba

损害的大小窗口。gydF4y2Ba

DataUsedgydF4y2Ba

属性的数据用于估计。与以下领域的结构。gydF4y2Ba

场gydF4y2Ba 描述gydF4y2Ba
的名字gydF4y2Ba

数据集的名称。gydF4y2Ba

类型gydF4y2Ba

数据类型。gydF4y2Ba

长度gydF4y2Ba

数据样本的数量。gydF4y2Ba

TsgydF4y2Ba

样品时间。这相当于gydF4y2Badata.TsgydF4y2Ba。gydF4y2Ba

InterSamplegydF4y2Ba

输入intersample行为。下列值之一:gydF4y2Ba

  • “zoh”gydF4y2Ba——零维护样本之间的分段常数输入信号。gydF4y2Ba

  • “呸”gydF4y2Ba——一阶保持维护一个分段线性输入信号之间的样本。gydF4y2Ba

  • “提单”gydF4y2Ba——带限行为指定连续时间输入信号零功率高于奈奎斯特频率。gydF4y2Ba

的价值gydF4y2BaInterSamplegydF4y2Ba离散时间模型的估计结果没有影响。gydF4y2Ba

InputOffsetgydF4y2Ba

在估计抵消从时域输入数据。gydF4y2Ba

OutputOffsetgydF4y2Ba

在估计抵消从时域输出数据。gydF4y2Ba

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

频率响应函数gydF4y2Ba

一个gydF4y2Ba频率响应函数gydF4y2Ba描述了系统的稳态响应正弦输入。对于一个线性系统,特定频率的正弦输入结果的输出也是一个正弦信号具有相同的频率,但用不同的振幅和相位。频率响应函数描述了振幅变化和相移作为频率的函数。gydF4y2Ba

为了更好地理解频率响应函数,线性动态系统的考虑以下描述:gydF4y2Ba

ygydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba =gydF4y2Ba GgydF4y2Ba (gydF4y2Ba 问gydF4y2Ba )gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba +gydF4y2Ba vgydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba

在这里,gydF4y2BaugydF4y2Ba(gydF4y2BatgydF4y2Ba),gydF4y2BaygydF4y2Ba(gydF4y2BatgydF4y2Ba)输入和输出信号,分别。gydF4y2BaGgydF4y2Ba(gydF4y2Ba问gydF4y2Ba)称为传递函数的系统相分离——捕获的系统动力学的输入输出。的符号gydF4y2BaGgydF4y2Ba(gydF4y2Ba问gydF4y2Ba)gydF4y2BaugydF4y2Ba(gydF4y2BatgydF4y2Ba)代表以下操作:gydF4y2Ba

GgydF4y2Ba (gydF4y2Ba 问gydF4y2Ba )gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba ∞gydF4y2Ba ggydF4y2Ba (gydF4y2Ba kgydF4y2Ba )gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba −gydF4y2Ba kgydF4y2Ba )gydF4y2Ba

问gydF4y2Ba是gydF4y2Ba运营商转变gydF4y2Ba,由以下方程:gydF4y2Ba

GgydF4y2Ba (gydF4y2Ba 问gydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba ∞gydF4y2Ba ggydF4y2Ba (gydF4y2Ba kgydF4y2Ba )gydF4y2Ba 问gydF4y2Ba −gydF4y2Ba kgydF4y2Ba 问gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba )gydF4y2Ba

GgydF4y2Ba(gydF4y2Ba问gydF4y2Ba)是gydF4y2Ba频率特性函数gydF4y2Ba在单位圆上评估时,gydF4y2BaGgydF4y2Ba(gydF4y2Ba问gydF4y2Ba=gydF4y2BaegydF4y2BaiωgydF4y2Ba)。gydF4y2Ba

在一起,gydF4y2BaGgydF4y2Ba(gydF4y2Ba问gydF4y2Ba=gydF4y2BaegydF4y2BaiωgydF4y2Ba)和输出噪声谱gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba 组成系统的频域描述。gydF4y2Ba

使用Blackman-Tukey频率特性函数估计方法是由以下方程:gydF4y2Ba

GgydF4y2Ba ^gydF4y2Ba NgydF4y2Ba (gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba

在这种情况下,gydF4y2Ba^gydF4y2Ba代表近似数量。这个方程的推导,看到这一章非参数时间和频域方法gydF4y2Ba[1]gydF4y2Ba。gydF4y2Ba

输出噪声频谱gydF4y2Ba

输出噪声频谱(频谱gydF4y2BavgydF4y2Ba(gydF4y2BatgydF4y2Ba)是由以下方程:gydF4y2Ba

ΦgydF4y2Ba ^gydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba −gydF4y2Ba |gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba |gydF4y2Ba 2gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba

这个方程的噪声频谱派生假设线性关系gydF4y2Ba ygydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba =gydF4y2Ba GgydF4y2Ba (gydF4y2Ba 问gydF4y2Ba )gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba +gydF4y2Ba vgydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba 持有,gydF4y2BaugydF4y2Ba(gydF4y2BatgydF4y2Ba)是独立的gydF4y2BavgydF4y2Ba(gydF4y2BatgydF4y2Ba),以下光谱之间的关系持有:gydF4y2Ba

ΦgydF4y2Ba ygydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba |gydF4y2Ba GgydF4y2Ba (gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba |gydF4y2Ba 2gydF4y2Ba ΦgydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba +gydF4y2Ba ΦgydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba

ΦgydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba GgydF4y2Ba (gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba ΦgydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba

这里的噪声频谱是由以下方程:gydF4y2Ba

ΦgydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba ≡gydF4y2Ba ∑gydF4y2Ba τgydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba RgydF4y2Ba vgydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba egydF4y2Ba −gydF4y2Ba 我gydF4y2Ba wgydF4y2Ba τgydF4y2Ba

ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba 出力互谱和吗gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba 是输入频谱。gydF4y2Ba

另外,干扰gydF4y2BavgydF4y2Ba(gydF4y2BatgydF4y2Ba)可以被描述为过滤白噪声:gydF4y2Ba

vgydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba =gydF4y2Ba HgydF4y2Ba (gydF4y2Ba 问gydF4y2Ba )gydF4y2Ba egydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba

在这里,gydF4y2BaegydF4y2Ba(gydF4y2BatgydF4y2Ba)是白噪声的方差gydF4y2Ba λgydF4y2Ba 噪声功率谱是由以下方程:gydF4y2Ba

ΦgydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba λgydF4y2Ba |gydF4y2Ba HgydF4y2Ba (gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba |gydF4y2Ba 2gydF4y2Ba

算法gydF4y2Ba

水疗中心gydF4y2Ba应用Blackman-Tukey光谱分析方法通过以下步骤:gydF4y2Ba

  1. 计算协方差和cross-covariancegydF4y2BaugydF4y2Ba(gydF4y2BatgydF4y2Ba),gydF4y2BaygydF4y2Ba(gydF4y2BatgydF4y2Ba):gydF4y2Ba

    RgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba NgydF4y2Ba ∑gydF4y2Ba tgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba NgydF4y2Ba ygydF4y2Ba (gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τgydF4y2Ba )gydF4y2Ba ygydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba RgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba NgydF4y2Ba ∑gydF4y2Ba tgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba NgydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τgydF4y2Ba )gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba RgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba NgydF4y2Ba ∑gydF4y2Ba tgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba NgydF4y2Ba ygydF4y2Ba (gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τgydF4y2Ba )gydF4y2Ba ugydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba

  2. 计算协方差和cross-covariance的傅里叶变换:gydF4y2Ba

    ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba τgydF4y2Ba =gydF4y2Ba −gydF4y2Ba 米gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba WgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba egydF4y2Ba −gydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba τgydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba τgydF4y2Ba =gydF4y2Ba −gydF4y2Ba 米gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba WgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba egydF4y2Ba −gydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba τgydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba τgydF4y2Ba =gydF4y2Ba −gydF4y2Ba 米gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba WgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba egydF4y2Ba −gydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba τgydF4y2Ba

    在哪里gydF4y2Ba WgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba 是损害窗口的宽度(滞后大小)gydF4y2Ba米gydF4y2Ba。您可以指定gydF4y2Ba米gydF4y2Ba控制的频率分辨率估计,约等于2π/gydF4y2Ba米gydF4y2Barad /样品时间。gydF4y2Ba

    默认情况下,这个操作使用128 0(排除)之间的等距的频率值和πgydF4y2BawgydF4y2Ba=gydF4y2Ba(1:128)/ 128 *π/ TsgydF4y2Ba和gydF4y2BaTsgydF4y2Ba是数据集的样本时间。默认延迟损害窗口的大小是什么gydF4y2BaM = min(长度(数据)/ 10、30)gydF4y2Ba。违约频率,操作使用快速傅里叶变换(FFT),这是更有效的比定义的频率。gydF4y2Ba

    请注意gydF4y2Ba

    M =gydF4y2Baγ是在表6.1gydF4y2Ba[1]gydF4y2Ba。标准差是在184和188页gydF4y2Ba[1]gydF4y2Ba。gydF4y2Ba

  3. 计算频率特性函数gydF4y2Ba GgydF4y2Ba ^gydF4y2Ba NgydF4y2Ba (gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba 和输出噪声谱gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba 。gydF4y2Ba

    GgydF4y2Ba ^gydF4y2Ba NgydF4y2Ba (gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ygydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba ΦgydF4y2Ba ^gydF4y2Ba ugydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba

    ΦgydF4y2Ba vgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba ≡gydF4y2Ba ∑gydF4y2Ba τgydF4y2Ba =gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba RgydF4y2Ba vgydF4y2Ba (gydF4y2Ba τgydF4y2Ba )gydF4y2Ba egydF4y2Ba −gydF4y2Ba 我gydF4y2Ba wgydF4y2Ba τgydF4y2Ba

光谱gydF4y2Ba的光谱矩阵输出和输入通道。也就是说,如果gydF4y2Baz = [data.OutputDatagydF4y2Ba,gydF4y2Badata.InputData]gydF4y2Ba,gydF4y2Ba光谱gydF4y2Ba包含了光谱数据的矩阵值功率谱gydF4y2BazgydF4y2Ba。gydF4y2Ba

年代gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba =gydF4y2Ba −gydF4y2Ba 米gydF4y2Ba 米gydF4y2Ba EgydF4y2Ba zgydF4y2Ba (gydF4y2Ba tgydF4y2Ba +gydF4y2Ba 米gydF4y2Ba )gydF4y2Ba zgydF4y2Ba (gydF4y2Ba tgydF4y2Ba )gydF4y2Ba ′gydF4y2Ba WgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba TgydF4y2Ba 年代gydF4y2Ba )gydF4y2Ba 经验值gydF4y2Ba (gydF4y2Ba −gydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba 米gydF4y2Ba )gydF4y2Ba

在这里,gydF4y2Ba”gydF4y2Ba是一个复共轭转置。gydF4y2Ba

引用gydF4y2Ba

[1]Ljung Lennart。gydF4y2Ba为用户系统标识:理论gydF4y2Ba。第二版。普伦蒂斯霍尔信息和系统科学系列。上台北:Prentice Hall PTR, 1999。gydF4y2Ba

版本历史gydF4y2Ba

之前介绍过的R2006agydF4y2Ba

全部展开gydF4y2Ba