光谱

绘图或返回输出功率谱的时间序列模型或线性输入输出模型的干扰频谱

页面上倒塌

句法

谱(系统)

光谱（SYS，{WMIN，WMAX}）

光谱（SYS，W）

光谱（SYS1，...，SYSN，W）

PS =频谱（SYS，W）

[PS，WOUT] =频谱（SYS）

[PS，WOUT，SDP] =频谱（SYS）

描述

绘图结果

例子

光谱(SYS.）绘制识别的时间序列模型的输出功率谱SYS.或已识别的输入-输出模型的干扰谱SYS.．该功能自动选择频率范围和点数。

如果SYS.是一个时间序列模型，然后SYS.代表系统：

$y （ T. ） = H E. （ T. ）$

这里，E.（T.）是高斯白噪音y（T.)为可观察的输出。
光谱绘图|H'H|乘以的方差E.（T.）和采样时间。
如果SYS.是一个投入产出模型，SYS.代表系统：

$y （ T. ） = G 你（ T. ） + H E. （ T. ）$

这里，你（T.)为测量输入，E.（T.为高斯白噪声和y（T.)为可观察的输出。
在这种情况下，光谱绘制扰动组件的光谱他（T.）。

用于采样时间的离散时间模型T._S.那光谱使用转换 $Z. = {E.}^{j ω {T.}_{S.}}$ 将单位圆映射到实际频率轴。该函数仅对小于尼奎斯特频率π/的频率绘制频谱T._S.，当不指定“Ts”时，默认为1个时间单位。

例子

光谱(SYS.，{w那wmax}）为频率产生范围的频谱曲线w到wmax．

光谱(SYS.那W.）使用向量中指定的频率创建频谱图W.．

例子

光谱（SYS1，...，SYSN，W.）在单个绘图上创建几个识别模型的频谱曲线图。这W.参数是可选的。

您可以为每个模型指定颜色，线条样式和标记。例如，谱(sys1‘r’,sys2, y,, sys3, gx)使用红色SYS1.，黄色的划线笔表示SYS2.和绿色X标记SYS3.．

返回结果

PS.=频谱（SYS.那W.）的功率谱振幅SYS.对于指定的频率，W.．没有绘图显示。

PS.=频谱（SYS.那W.）的功率谱振幅SYS.对于指定的频率，W.．

[PS.那wout] =频谱（SYS.）返回频率向量，wout，计算输出功率谱。

[PS.那wout那SDPS.] =频谱（SYS.）返回功率谱的估计标准差。

例子

全部收缩

时间序列模型的绘图输出频谱

打开生活的脚本

加载时间序列估计数据。

加载iddata9Z9.

用最小二乘方法估计四阶AR模型。

sys = ar (z9 4“ls”）;

绘制模型的输出频谱。

谱(系统);

图中包含一个轴对象。标题为“From: e@y1 To: y1”的axes对象包含一个类型为line的对象。该对象表示sys。

绘制SISO线性识别模型的噪声谱

打开生活的脚本

加载估计数据。

加载iddata1Z1;

估计单输入单输出状态空间模型。

sys = n4sid (z1, 2);

绘制模型的噪声频谱。指定频率范围从0.1到50 rad/sec。

光谱（SYS，{0.1,50}）;

图中包含一个轴对象。标题为“From: e@y1 To: y1”的axes对象包含一个类型为line的对象。该对象表示sys。

该函数绘制频谱，但将频率范围限制在约31.4 rad/sec的奈奎斯特频率范围内。

用正方形比较正弦信号的频谱

打开生活的脚本

创建一个由5个正弦波组成的输入，每个正弦波在全频率范围内传播。把这个信号的频谱与其平方的频谱进行比较。

创建一个扩展20个周期的正弦信号之和，每个周期包含100个样本。指定信号组合5个随机相位的正弦波，使用10次试验来找到信号传播最小的集合。有关此步骤的更多信息，请参见id inippul.．

u = id iniNput（[100 1 20]，'正弦'，[]，[]，[5 10 1]）;

创建一个handlder程序iddata目的你它包含输入u，周期为100。

u = iddata ([], u, 1,'每', 100);

平方输入值并将其存储在新的iddata目的u2。

U2 = U.U. ^ 2;U2 = IDDATA（[]，U2,1，'每', 100);

使用etfe来估计经验传递函数模型你和U2.．把这些模型的功率谱画在一起。使用不同的标记颜色和类型来区分光谱来源。

光谱（ETFE（U），的r *，etfe（u2），“+”）

图中包含一个轴对象。具有标题的轴对象：E @ U1至：U1包含2个类型的2个对象。这些对象表示Untitled1，Untitled2。

该曲线显示了一些频率分裂在其中U2.基于频谱没有排列u1,基于频谱，但相反，包含两个侧翼的光谱点U1.的点。这表明了平方系统的非线性。

输入参数

全部收缩

`SYS.`-鉴定的模型
`idpoly`目的|`IDProc.`目的|`IDS.`目的|`idtf.`目的

识别的模型，指定为idpoly对象,一个IDProc.对象,一个IDS.对象，或者idtf.对象。

如果SYS.是一个时间序列模型，然后SYS.代表系统：

$y （ T. ） = H E. （ T. ）$

这里，e (t)是高斯白噪音y (t)是观察到的输出。
如果SYS.是一个输入输出模型，然后SYS.代表系统：

$y （ T. ） = G 你（ T. ） + H E. （ T. ）$

这里，u (t)为测量输入，e (t)是高斯白噪音y (t)是观察到的输出。