主要内容gydF4y2Ba

估计时间序列功率谱gydF4y2Ba

一个gydF4y2Ba频率特性gydF4y2Ba模型封装了一个线性系统在一定频率范围内的频率响应。当数据同时包含输入和输出通道时，频率响应模型描述了系统对正弦输入的稳态响应。时间序列数据不包含输入通道。时间序列模型的频率响应减少到输出数据的频谱表示。输出数据隐式地包含了未测量输入噪声的影响。gydF4y2Ba

对于一个离散时间系统，输入和输出都按时间间隔采样gydF4y2BaTgydF4y2Ba，传递函数gydF4y2BaG (z)gydF4y2Ba表示输入的z变换gydF4y2BaU (z)gydF4y2Ba和输出gydF4y2BaY (z)gydF4y2Ba：gydF4y2Ba

Y (z)gydF4y2Ba＝gydF4y2BaGgydF4y2Ba（gydF4y2BazgydF4y2Ba）gydF4y2BaUgydF4y2Ba（gydF4y2BazgydF4y2Ba) +gydF4y2BaHgydF4y2Ba（gydF4y2BazgydF4y2Ba）gydF4y2BaEgydF4y2Ba（gydF4y2BazgydF4y2Ba）gydF4y2Ba

HgydF4y2Ba（gydF4y2BazgydF4y2Ba)表示每个输出的噪声传递函数和gydF4y2BaEgydF4y2Ba（gydF4y2BazgydF4y2Ba)为加性扰动的z变换gydF4y2BaegydF4y2Ba（gydF4y2BatgydF4y2Ba)与方差gydF4y2BaΛgydF4y2Ba．gydF4y2Ba

对于时间序列模型，此方程可简化为:gydF4y2Ba

Y (z)gydF4y2Ba＝gydF4y2BaHgydF4y2Ba（gydF4y2BazgydF4y2Ba）gydF4y2BaEgydF4y2Ba（gydF4y2BazgydF4y2Ba）gydF4y2Ba

在这种情况下，gydF4y2BaEgydF4y2Ba（gydF4y2BazgydF4y2Ba)表示假设但未测量的白噪声输入扰动。单输出噪声谱gydF4y2BaΦgydF4y2Ba在存在标量方差干扰噪声的情况下gydF4y2BaλgydF4y2Ba定义为:gydF4y2Ba

ΦgydF4y2Ba vgydF4y2Ba （gydF4y2Ba ωgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba λgydF4y2Ba TgydF4y2Ba |gydF4y2Ba HgydF4y2Ba （gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba TgydF4y2Ba ）gydF4y2Ba |gydF4y2Ba 2gydF4y2Ba

等效多输出噪声功率谱为:gydF4y2Ba

ΦgydF4y2Ba vgydF4y2Ba （gydF4y2Ba ωgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba TgydF4y2Ba HgydF4y2Ba （gydF4y2Ba egydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba TgydF4y2Ba ）gydF4y2Ba ΛgydF4y2Ba HgydF4y2Ba ”gydF4y2Ba （gydF4y2Ba egydF4y2Ba −gydF4y2Ba 我gydF4y2Ba ωgydF4y2Ba TgydF4y2Ba ）gydF4y2Ba

在这里,gydF4y2BaΛgydF4y2Ba长度等于输出数的方差向量。gydF4y2Ba

在命令行中估计时间序列功率谱gydF4y2Ba

您可以使用gydF4y2BaetfegydF4y2Ba，gydF4y2Ba水疗中心gydF4y2Ba,gydF4y2BaspafdrgydF4y2Ba用于估计时域和频域数据的时间序列功率谱的命令。这些函数返回由表示的估计模型gydF4y2BaidfrdgydF4y2Ba模型对象，在属性中包含光谱数据gydF4y2BaSpectrumDatagydF4y2Ba属性中的谱方差gydF4y2BaNoiseCovariancegydF4y2Ba．对于多输出数据，gydF4y2BaSpectrumDatagydF4y2Ba包含每个输出的功率谱和每个输出对之间的交叉谱。gydF4y2Ba

估计和比较时间序列频率响应的命令gydF4y2Ba

命令gydF4y2Ba 描述gydF4y2Ba
etfegydF4y2Ba

利用傅里叶分析估计周期图。gydF4y2Ba
水疗中心gydF4y2Ba

使用谱分析估计功率谱的标准偏差。gydF4y2Ba
spafdrgydF4y2Ba

使用可变频率分辨率估计功率谱的标准偏差。gydF4y2Ba
光谱gydF4y2Ba 估计和绘制时间序列模型的输出功率谱。gydF4y2Ba

例如，假设gydF4y2BaygydF4y2Ba是时间序列数据。估计功率谱gydF4y2BaggydF4y2Ba还有周期图gydF4y2BapgydF4y2Ba使用gydF4y2Ba水疗中心gydF4y2Ba而且gydF4y2BaetfegydF4y2Ba．将模型与三个标准差置信区间一起绘制gydF4y2Ba光谱gydF4y2Ba．gydF4y2Ba

G = spa(y);P = etfe(y);谱(g、p);gydF4y2Ba

有关光谱估计的更详细示例，请参见gydF4y2Ba在命令行识别时间序列模型gydF4y2Ba．有关各个命令的更多信息，请参阅相应的参考页面。gydF4y2Ba

使用App估计时间序列功率谱gydF4y2Ba

您必须已经将数据导入应用程序。gydF4y2Ba

在System Identification app中估计时间序列谱模型:gydF4y2Ba

  1. 在“系统识别”应用程序中，选择gydF4y2Ba估计gydF4y2Ba>gydF4y2Ba谱模型gydF4y2Ba打开光谱模型对话框。gydF4y2Ba

  2. 在gydF4y2Ba方法gydF4y2Ba列表中，选择要使用的光谱分析方法。有关每种方法的信息，请参见gydF4y2Ba谱模型计算方法的选择gydF4y2Ba．gydF4y2Ba

  3. 指定计算频谱模型的频率，方法如下:gydF4y2Ba

    • 在gydF4y2Ba频率gydF4y2Ba字段，输入一个矢量的值，一个MATLABgydF4y2Ba®gydF4y2Ba求值为向量的表达式，或MATLAB工作区中向量的变量名。例如,gydF4y2Balogspace (2500)gydF4y2Ba．gydF4y2Ba

    • 使用以下组合gydF4y2Ba频率间隔gydF4y2Ba而且gydF4y2Ba频率gydF4y2Ba构造值的频率向量:gydF4y2Ba

      • 在gydF4y2Ba频率间隔gydF4y2Ba列表中,选择gydF4y2Ba线性gydF4y2Ba或gydF4y2Ba对数gydF4y2Ba频率间隔。gydF4y2Ba

        请注意gydF4y2Ba

        为gydF4y2BaetfegydF4y2Ba，只有gydF4y2Ba线性gydF4y2Ba选项可用。gydF4y2Ba

      • 在gydF4y2Ba频率gydF4y2Ba字段，输入频率点的个数。gydF4y2Ba

      对于时域数据，频率范围为0 ~ Nyquist频率。对于频域数据，频率范围从数据集中的最小频率到最大频率。gydF4y2Ba

  4. 在gydF4y2Ba频率分辨率gydF4y2Ba字段，输入频率分辨率，如中所述gydF4y2Ba控制频谱模型的频率分辨率gydF4y2Ba．如需使用默认值，请输入gydF4y2Ba默认的gydF4y2Ba或者让这个区域空着。gydF4y2Ba

  5. 在gydF4y2Ba模型名称gydF4y2Ba字段中，输入相关分析模型的名称。模型名称在模型板中应该是唯一的。gydF4y2Ba

  6. 点击gydF4y2Ba估计gydF4y2Ba将该模型添加到系统识别应用程序中的模型板。gydF4y2Ba

  7. 在“光谱模型”对话框中，单击gydF4y2Ba关闭gydF4y2Ba．gydF4y2Ba

  8. 要查看估计的扰动谱，请选择gydF4y2Ba噪声频谱gydF4y2Ba复选框。有关使用此绘图的更多信息，请参见gydF4y2Ba噪声谱图gydF4y2Ba．gydF4y2Ba

若要将模型导出到MATLAB工作区，请将其拖到gydF4y2Ba到工作空间gydF4y2Ba在“系统识别”应用程序中，可以查看输出的功率谱和置信区间gydF4y2BaidfrdgydF4y2Ba模型对象使用gydF4y2Ba光谱gydF4y2Ba命令。gydF4y2Ba

另请参阅gydF4y2Ba

|gydF4y2Ba|gydF4y2Ba|gydF4y2Ba

相关的话题gydF4y2Ba