rpmtrack

产生一个信号,像踩油门发动机引起的振动。信号采样1 kHz 30秒和包含三个谐波组件订单1,2.4,3,5,振幅与4和0.5,分别。将信号嵌入unit-variance高斯白噪声,并将它存储在一个MATLAB®的时间表。瞬时频率乘以60获得RPM概要文件。情节RPM概要文件。

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

RPM文件来自振动信号。使用4分5秒间隔指定脊对应订单2.4。显示一个总结输出的时间表。

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

RowTimes:吹捧:30001 x1持续时间值:分0秒值15秒马克斯30秒步伐0.001秒变量:转:30001 x1双重价值:8879.8最小2.2204 e-16值4327.2最大

绘制重建RPM概要和点用于重建。

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

使用RPM的概要文件生成order-RPM地图中提取信号。

rpmordermap (vib fs, rpmest.rpm)

重建和情节构成信号的时域波形。放大后的时间间隔发生瞬变腐朽。

xrc = orderwaveform (vib、fs、rpmest.rpm ol);图绘制(t, xrc)传说([repmat (订单= ',[3])num2str (ol)]) xlim (20 [5])

估计一个风扇叶片的RPM概要switchoff后减慢。

一个工业屋顶风机旋转在20000 rpm是关闭的。空气阻力(与轴承摩擦的贡献可以忽略不计)导致风扇转子停止在大约6秒。高速相机的措施x协调的风机叶片的速度1 kHz。

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

理想化的风扇叶片作为一个质点绕转子中心半径为50厘米。叶片阻力与速度成正比,导致以下相位角的表达式

在哪里 ${\mathit{f}}_0$初始频率和吗 $\mathit{T}=0。\mathrm{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(大小(φ))/ 10;y = a * sin(φ)+ randn(大小(φ))/ 10;情节(t t, x,, y)

使用rpmtrack函数来确定RPM概要文件。类型

rpmtrack (x, fs)

在命令行打开交互图。

使用滑块调整的频率分辨率时频映射到大约11赫兹。假定信号分量对应订单1和脊提取设置结束时间3.0秒。使用十字光标在时频图添加按钮以添加三个点躺在山脊上。另外,双击鼠标添加点在您选择的位置。点击估计跟踪和提取RPM概要文件。

验证RPM剖面指数衰减。在出口选项卡上,单击出口并选择生成MATLAB脚本。脚本出现在编辑器中。

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

运行脚本。在半对数图显示RPM概要文件。

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