主要内容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到奈奎斯特频率。对于频域数据，频率范围从数据集中的最小频率到最大频率。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

Baidu
map