rpmfreqmap的
用于分析订单的频率反转图。
库::`工程师`
语法
函数调用
-
[参数:map],[参数:freq],[参数:rpm_out],[参数:time],[参数:res_out]=rpmfreqmap(<参数:x>>,<参数:fs>>,<参数:rpm_in>>)-返回频率周转卡的矩阵[参数:地图]输入矢量的频率分析的结果而获得的[参数:x]，以及具有频率的矢量[参数：freq]、转速[参数：rpm_out]和时间的瞬间[参数：时间]其中计算频率图[参数:地图]. 该函数还返回使用的分辨率带。[参数：res_out].输入参数
[参数:x]定的速度值进行测量。[参数：rpm_in]以rpm表示。 论点[参数:fs]-以Hz为单位的采样率。 矩阵的每一列[参数:地图]包含每个值的频谱分量幅度的RMS估计值[参数：rpm_in]. 功能rpmfreqmap使用短时傅立叶变换来分析信号的频谱成分[参数:x].
-
map,freq,rpm,time,res=rpmfreqmap(x,fs,rpm,[Argument:res_in>)-设置输入分辨率波段<<参数：res_in]以Hz为单位的地图。
-
rpmfreqmap（___,[参数:out]=:情节)-建立一个频率图作为旋转速度和时间的函数在交互式图形。 该图也称为坎贝尔图。
争论
输入参数
#
财政司司长 —
抽样调查
费率+
标量,标量
Details
采样率，设置为正标量，以Hz表示。
# res — 分辨率波段
+
财政司司长/128 （默认情况下）| 标量,标量
Details
频率反转卡的分辨率带，设置为正标量。 如果论点是 资源,资源 如果未指定，则函数 rpmfreqmap 设置其值等于采样频率除以 128. 如果信号持续时间不足，则该函数使用信号的整个长度来计算单个频率估计。
| 数据类型 |
|
输入参数 «名称-值»
将可选参数对指定为 名称，值，在哪里 姓名 -参数的名称，以及 价值 -适当的值。 类型参数 «名称-值» 它们应该放在其他参数之后，但对的顺序并不重要。 您可以指定多个对 «名称-值».
使用逗号分隔名称和值，以及 姓名 把它放在引号里。
例子: rpmfreqmap(x,fs,rpm,"Scale","dB","Window","hann") 以dB为单位设置频率图评级比例，并使用Hanna窗口。
# 振幅 — 振幅类型
+
"rms" （默认情况下）| "山顶" | "权力"
Details
频率反转卡的振幅类型，由以下值之一设置:
-
"rms"-返回每个估计频率的均方根振幅; -
"山顶"-返回每个估计频率的峰值幅度; -
"权力"-返回每个估计频率的功率电平。
# 重叠 — 相邻段之间的重叠百分比
+
50 （默认情况下）| 标量,标量
Details
相邻段之间的重叠百分比，以 0 以前 100. 意义 0 意味着相邻段不重叠。 意义 100 表示相邻段偏移一个计数。 更高的重叠百分比提供了更平滑的地图，但增加了计算时间。 有关详细信息，请参阅功能说明 rpmordermap.
| 数据类型 |
|
# 比例尺 — 频率-周转卡的规模
+
"线性" （默认情况下）| "dB"
Details
频率反转卡的刻度，由以下值之一设置:
-
"线性"-返回具有线性比例的地图; -
"dB"-返回具有对数刻度的地图，其值以dB表示。
# 窗口 — 分析窗口
+
"汉恩" （默认情况下）| "切布温" | "弗拉托普温" | "汉明" | "凯撒" | "直肠"
Details
分析窗口，设置为以下值之一:
-
"汉恩"-汉娜设置窗口。 有关详细信息，请参阅功能说明汉恩. -
"切布温"-设置切比雪夫窗口。 使用元组以dB表示旁瓣的衰减。 旁瓣的减弱应该更大45分贝。 如果未设置衰减值，则默认值为100分贝。 有关详细信息，请参阅功能说明切布温.例子:
rpmfreqmap（x，fs，rpm，"幅度"，"峰值"，"窗口"，（"chebwin"，80））设置峰值幅度并使用切比雪夫窗口与旁瓣的衰减80分贝。 -
"弗拉托普温"-设置一个平顶的窗口。 有关详细信息，请参阅功能说明弗拉托普温. -
"汉明"-设置汉明窗。 有关详细信息，请参阅功能说明汉明. -
"凯撒"-设置凯撒窗口。 使用元组指定形状因子。 . 形状因子必须是正标量。 如果未指定form参数的值，则默认使用该值。0.5. 有关详细信息，请参阅功能说明凯撒. -
"直肠"-设置一个矩形窗口。 有关详细信息，请参阅功能说明[医]直肠.
| 数据类型 |
|
例子：
频率-周转卡的表面图
Details
我们将生成一个由两个线性和一个二次线性调频组成的信号，所有这些都具有采样频率。 600 Hz期间 15 几秒钟。 产生信号的系统将其旋转速度从 10 以前 40 测试期间的rpm。
让我们生成转速计读数。
fs = 600
t1 = 15
t = 0:1/fs:t1
f0 = 10
f1 = 40
rpm = collect(60 * range(f0, f1, length=length(t)))线性啁啾有数量级 1 和 2.5. 顺序的组件 1 其具有等于另一分量振幅的一半的振幅。 二次线性调频从阶数开始 6 并在测量结束时返回到该顺序。 其振幅为 0.8. 让我们使用这些信息创建一个信号。
o1 = 1
o2 = 2.5
o6 = 6
import EngeeDSP.Functions: chirp
x = 0.5 * chirp(t, o1 * f0, t1, o1 * f1) + chirp(t, o2 * f0, t1, o2 * f1) +
0.8 * chirp(t, o6 * f0, t1, o6 * f1, "quadratic")让我们计算信号频率对旋转速度的依赖性的映射。 我们在每个测量单元中使用峰值幅度。 指定分辨率 6 赫兹。 让我们使用平顶窗口处理数据。
import EngeeDSP.Functions: rpmfreqmap
map, fr, rp = rpmfreqmap(x, fs, rpm, 6, "Amplitude", "peak", "Window", "flattopwin")让我们以曲面图的形式构建一个频率-周转图。
FR = repeat(fr', length(rp), 1)
RP = repeat(rp, 1, length(fr))
surface(FR, RP, map';
camera = (-6, 60),
xlabel = "Frequency (Hz)",
ylabel = "RPM",
zlabel = "Amplitude")