tsa
时间同步信号平均。
库::`工程师`
语法
函数调用
-
[参数:ta],[参数:t_out],[参数:p],[参数:rpm]=tsa(___,Name=Value)-使用以下参数为以前的任何语法选项设置附加参数 «名称-值». 例如,可以指定轴每转的转速计脉冲数,或者选择在时域或频域中对信号进行平均。
-
tsa(___;<参数:out>=:情节)-在电流图上绘制对应于每个信号段的时间同步平均信号和时域信号。
争论
输入参数
#
财政司司长 —
抽样调查
费率+
标量,标量
Details
的采样率,设置为正标量。
| 数据类型 |
|
# tp — 脉冲时间
+
标量,标量 | 向量资料
Details
脉冲时间,设置为标量或矢量:
-
标量,标量-旋转发生的恒定时间间隔; -
向量资料-非负的,严格增加的时刻,确定旋转的恒定相位。
使用功能 转速表 以从转速计信号中提取转速计脉冲的采样周期。
| 数据类型 |
|
输入参数 «名称-值»
将可选参数对指定为 名称=值,在哪里 姓名 -参数的名称,以及 价值 -适当的值。 类型参数 «名称-值» 它们应该放在其他参数之后,但对的顺序无关紧要。 您可以指定多个对 «名称-值».
例子: ta=tsa(x,方法="pchip",ResampleFactor=10) 指示信号应在上采样 10 一次并使用分段三次Hermite插值在时域中平均。
# 方法 — 平均算法
+
"线性" (默认情况下)| "样条" | "pchip" | "fft"
Details
由以下值之一定义的插值方案:
-
"线性"-在时域执行线性插值和平均; -
"样条"-执行三次样条插值和时域平均; -
"pchip"-执行分段三次Hermite插值和时域平均; -
"fft"-在频域中执行平均。
例子:
正弦波的同步时间平均
Details
让我们计算噪声正弦波的时间同步平均。
我们将产生一个由沉浸在白高斯噪声中的正弦波组成的信号。 信号以频率采样 500 Hz期间 20 几秒钟。 指定正弦波的频率 10 Hz和噪声色散 0.01. 让我们绘制信号的一个周期。
fs = 500
t = 0:1/fs:20-1/fs
f0 = 10
y = sin.(2π * f0 * t) + randn(length(t)) / 10
plot(t, y, xlims=(0, 1/f0))
让我们计算信号的时间同步平均。 作为同步信号,我们使用一组周期与正弦相同的脉冲。 在功能 tsa 设置参数 [参数:out] 意义 :剧情 来显示结果。
import EngeeDSP.Functions: tsa
tPulse = 0:1/f0:maximum(t)
tsa(y, fs, tPulse, out = :plot)
算法
与指定的输入信号,采样频率和一组转速计脉冲,功能 tsa 执行以下步骤:
-
根据转速计脉冲和为参数指定的值确定周期的开始和结束时间
[参数:脉动]. -
根据为自变量指定的值重新采样输入信号
[参数:重新采样因子]. -
根据为自变量指定的值对过采样信号求平均值
[参数:方法].
*如果为论点 [参数:方法] 值设置 "fft",功能:
-
将信号分成对应于不同周期的段。
-
计算每个段的离散傅立叶变换。
-
截断较长的变换,以便所有变换具有相同的长度。
-
平均光谱。
-
计算平均值的逆离散傅立叶变换以将其转换到时域。
*如果为论点 [参数:方法] 定义了时域方法之一,函数:
-
使用这种方法,它将信号插值到对应于不同周期的均匀分布样本的网格上。
-
根据为自变量指定的值组合通过过采样获得的信号段
[参数:数值]. -
计算所有段的平均值。