产生一个类似于汽车发动机旋转引起的振动的信号。信号以1 kHz采样30秒，包含阶为1、2.4和3的三个谐波分量，振幅分别为5、4和0.5。将信号嵌入单位方差高斯白噪声中，并将其存储在MATLAB®时间表中。将瞬时频率乘以60得到RPM剖面。绘制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*pi*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)

从振动信号中推导出转速曲线。以5秒间隔使用4个点来指定对应于2.4阶的脊。显示输出时间表的摘要。

NDX = (5:5:20)*fs;Order = ol(2);P = [t(ndx) order*fis(ndx)];Rpmest = rpmtrack(xt,order,p);总结(转)

RowTimes: tout: 30001x1 duration值:Min 0 sec Median 15 sec Max 30 sec TimeStep 0.001 sec变量:rpm: 30001x1 double值:Min 2.2204e-16 Median 4327.2 Max 8879.8

绘制重建的RPM剖面和重建中使用的点。

持有在rpmest.rpm情节(秒(rpmest.tout),“。”)情节(t (ndx),金融中间人(ndx) * 60岁“好吧”)举行从传奇(“原始”，“重建”，“脊点”，“位置”，“西北”）

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

rpmordermap (vib fs, rpmest.rpm)

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

XRC =顺序波形(vib,fs,rpme .rpm,ol);图(t,xrc)图例([repmat('Order = '，[3 1]) num2str(ol')]) xlim([5 20])

估计风机叶片在关闭后减速时的RPM分布。

一台转速为每分钟20000转的工业屋顶风机被关闭。空气阻力(轴承摩擦的贡献可以忽略不计)导致风扇转子在大约6秒内停止。一个高速摄像机测量x-一个风扇叶片的坐标，速率为1khz。

Fs = 1000;T = 0:1/fs:6-1/fs;Rpm0 = 20000;

理想的风扇叶片作为一个点质量围绕转子中心半径为50厘米。叶片受到与速度成比例的阻力，导致相位角的表达式如下

在哪里 ${\mathit{f}}_0$初始频率和 $\mathit{T}＝0．75$第二是衰减时间。

A = 0.5;F0 = rpm0/60;T = 0.75;= 2*pi*f0*T*(1-exp(-t/T));

计算并绘制x- - -y-刀片的坐标。加上方差的高斯白噪声 $0．1^2$．

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

使用rpmtrack函数来确定RPM剖面。类型

rpmtrack (x, fs)

在命令行中打开交互式图形。

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

验证RPM配置文件呈指数衰减。在出口选项卡上,单击出口并选择生成MATLAB脚本．脚本出现在编辑器中。

% MATLAB代码从rpmtrack GUI由MATLAB 9.12和信号处理工具箱8.7生成%生成于2021年10月12日09:36:49设定抽样率Fs = 1000.0000;设置感兴趣的山脊的顺序订单= 1.0000;在感兴趣的山脊上设置山脊点ridgePoints = [.．.0.4077 - 194.6023;.．.0.9781 - 89.4886;.．.2.3678 - 15.6250);估计转速[rpmOut,tOut] = rpmtrack(x,fs,order,ridgePoints，.．.“方法”，“stft”，.．.“FrequencyResolution”, 11.1612,.．.“PowerPenalty”正,.．.“FrequencyPenalty”, 0.0000,.．.“开始时间”, 0.0000,.．.“EndTime”, 3.0000);

运行脚本。在半对数图中显示RPM剖面。

符号(tOut,rpmOut) ylim([500 20000])