产生一种类似于汽车引擎旋转时产生的振动的信号。信号以1 kHz采样30秒，包含三个阶数为1、2.4和3的谐波分量，振幅分别为5、4和0.5。将信号嵌入到单位方差高斯白噪声中，并存储在MATLAB®时间表中。将瞬时频率乘以60得到转速剖面。绘制RPM配置文件。

fs = 1000;t =（0：1 / fs：30）';fit = @（a，x）（t-x）。^ 6. * exp（ - （t-x））。*（（t-x）> = 0）* a';FIS = FIT（[0.4 1 0.6 1]，[0 6 13 17]）;phi = 2 * pi * cumtrapz（t，fis）;ol = [1 2.4 3];amp = [5 4 0.5]';VIB = COS（PHI。* OL）* AMP + RANDN（尺寸（t））;xt =时间表（秒（t），vib）;情节（T，FIS * 60）

使用rpmtrack功能从振动信号推导RPM配置文件。使用5秒间隔的四个点来指定对应于订单2.4的脊。显示输出时间表的摘要。

ndx =（5：5：20）* fs;订单= ol（2）;p = [t（ndx）顺序* fis（ndx）];rpmest = rpmtrack（xt，order，p）;摘要（rpmest）

rowtimes：Tout：30001x1持续时间值：Min 0秒中位数15秒MAX 30秒SIMESTEP 0.001秒变量：RPM：30001x1双值：MIN 2.2204E-16中位数4327.2 MAX 8879.8

绘制重建的RPM配置文件和重建中使用的点。

抓住在rpmest.rpm情节(秒(rpmest.tout),'.-')情节(t (ndx),金融中间人(ndx) * 60岁'好的')举行从传奇(“原始”那“重建”那“脊点”那“位置”那“西北”的）

使用提取的RPM配置文件生成信号的订单-RPM映射。

rpmordermap (vib fs, rpmest.rpm)

重构并绘制构成信号的时域波形。放大瞬态衰减后发生的时间间隔。

xrc = orderwaveform (vib、fs、rpmest.rpm ol);图绘制(t, xrc)传说([repmat (订单= '，[3 1]）NUM2STR（OL'）]）XLIM（[5 20]）

估计风扇叶片的RPM轮廓，因为它在开关后减速。

每分钟2万转的工业用屋顶风扇被关闭。空气阻力(来自轴承摩擦的微不足道的贡献)导致风扇转子在大约6秒内停止。一台高速摄像机测量X-坐标的一个风扇叶片的速率为1khz。

fs = 1000;t = 0:1 / fs: 6 - 1 / f;rpm0 = 20000;

将风扇叶片理想化为以50厘米半径绕转子中心旋转的质点。叶片受到的阻力与速度成正比，相位角的表达式如下:

在哪里 ${\mathit{F}}_{0.}$初始频率是和吗 $\mathit{T.}=0.．75$二是衰变时间。

一个= 0.5;f0 = rpm0/60;T = 0.75;φ= 2 *π* f0 * T * (1-exp (T / T));

计算并绘制X- - -y- 刀片的控制。添加白色高斯噪声的方差 $0.．1^2$．

X = a*cos(phi) + randn(size(phi))/10;Y = a*sin(phi) + randn(size(phi))/10;情节(t t, x,, y)

使用rpmtrack用于确定RPM配置文件的功能。类型

rpmtrack（x，fs）

在命令行打开交互式图。

使用滑块调整时频映射的频率分辨率到大约11hz。假设信号分量对应于order 1，并将脊提取的结束时间设置为3.0秒。使用十字准线光标在时频图和添加按钮添加三分位于山脊上。（或者，双击光标以在您选择的位置添加点。）单击估计跟踪和提取RPM配置文件。

确认RPM配置文件呈指数衰减。在出口选项卡，单击出口并选择生成matlab脚本．脚本将出现在编辑器中。

% MATLAB代码从rpmtrack GUI%由matlab9.4和信号处理工具箱8.0生成%生成于2017年12月18日19:00:35％设置采样率fs = 1000.0000;%设置兴趣岭的顺序订单= 1.0000;在感兴趣的山脊上设置山脊点ridgepoints = [．．.0.4501 - 179.8246;．．.0.9944 - 88.5965;．．.2.4161 - 11.4035);％估计转/分[rpmOut,兜售]= rpmtrack (x, fs,秩序,ridgePoints,．．.“方法”那“stft”那．．.“FrequencyResolution”, 11.1612,．．.'PowerGenalty'，inf，．．.“FrequencyPenalty”0.0000，．．.'开始时间'0.0000，．．.'时间结束'，3.0000）;

运行脚本。在半对数图中显示RPM配置文件。

semilogy(吹捧,rpmOut) ylim (20000 [500])