产生一种类似于汽车发动机发动时产生的振动的信号。信号在1 kHz采样30秒，包含1、2.4和3阶的三个谐波分量，振幅分别为5、4和0.5。将信号嵌入单位方差高斯白噪声中，并存储在MATLAB®时间表中。将瞬时频率乘以60得到RPM配置文件。绘制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：租赁：30001x1持续时间值：Min 0秒中位数15秒最多30秒钟时间表0.001秒变量：RPM：30001x1双值：Min 2.2204E-16中位数4327.2 Max 8879.8

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

抓住上绘图（秒（rpmest.tout），rpmest.rpm，'.-'）绘图（t（ndx），fis（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]）

当风扇叶片关闭后转速下降时，请对其转速进行估算。

关闭20,000 rpm的工业屋顶风扇被关闭。空气阻力（轴承摩擦的贡献可忽略不计）使风扇转子在大约6秒内停止。高速相机测量X- 一个风扇叶片的速度为1 kHz的速率。

FS = 1000;T = 0：1 / FS：6-1 / FS;RPM0 = 20000;

理想化风扇刀片作为点质量，将转子中心以50cm的半径盘旋。刀片体验成比例到速度的拖曳力，从而导致相位角的表达式：

在哪里 ${\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() + randn(size())/10;Y = a*sin(phi) + randn(size(phi))/10;情节(t t, x,, y)

使用rpmtrack函数来确定RPM配置文件。类型

rpmtrack（x，fs）

在命令行打开交互式图。

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

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

% MATLAB代码从rpmtrack GUI%由MATLAB 9.4和信号处理工具箱8.0生成％在18-dec-2017 19:00:35生成％设置采样率FS = 1000.0000;%设定兴趣脊的顺序订单= 1.0000;％在感兴趣的山脊上设定脊点ridgepoints = [......0.4501 179.8246;......0.9944 88.5965;......2.4161 - 11.4035);％估计转/分[RPMOUT，TOUT] = RPMTRACK（X，FS，ORDER，RIDGPOINTS，......'方法'那'stft'那......“FrequencyResolution”, 11.1612,......'PowerGenalty'，inf，......'频率金蛋白'，0.0000，......'开始时间'，0.0000，......'时间结束'，3.0000）;

运行脚本。在半机型绘图中显示RPM配置文件。

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