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

fs = 1000;t = (0:1 / fs: 30) ';适合= @ (a, x)(时距)。^ 6。* exp(-(时距))。*((时距)> = 0)*”;Fis = fit([0.4 1 0.6 1]，[0 6 13 17]);φ= 2 *π* cumtrapz (t, fis);Ol = [1 2.4 3];Amp = [5 4 0.5]';Vib = cos(phi.*ol)*amp + randn(size(t));xt =时间表(秒(t), vib);情节(t, fis * 60)

使用rpmtrack函数从振动信号导出RPM轮廓。以5秒为间隔使用4个点来指定与2.4顺序相对应的脊。显示输出时间表的摘要。

ndx = (5:5:20) * fs;订单= ol (2);P = [t(ndx) order*fis(ndx)];rpm = rpmtrack (xt,秩序,p);总结(转)

RowTimes: tout: 30001x1 duration值:Min 0 sec中位数15 sec Max 30 sec TimeStep 0.001 sec变量:rpm: 30001x1 double值: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,．..“PowerPenalty”正,．..“FrequencyPenalty”, 0.0000,．..“开始时间”, 0.0000,．..“EndTime”, 3.0000);

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

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