主要内容

rpmfreqmap.

用于订单分析的频率-转速图

描述

例子

地图= rpmfreqmap(XFS.rpm.返回频率- rpm映射矩阵,地图,这是对输入矢量进行频率分析的结果,XX是在一个设定范围内测量的吗rpm.以每分钟转数表示的转速。FS.是Hz中的采样率。每列地图包含在每个值的频谱内容的根均方(RMS)幅度估计rpm.rpmfreqmap.利用短时傅里叶变换分析X

例子

地图= rpmfreqmap(XFS.rpm.res.指定映射的分辨率带宽,以Hz为单位。

例子

地图= rpmfreqmap(___名称,值指定的选择使用名称,值对,以及前面语法中的输入参数。

[地图弗里克rpm.时间res.] = rpmfreqmap(___返回具有频率映射的频率,转速和时间速度的向量。它还返回使用的分辨率带宽。

例子

rpmfreqmap(___没有输出参数将频率贴图绘制在交互式上的转速和时间的函数。剧情也被称为一个坎贝尔图

例子

全部折叠

创建一个以600hz采样5秒的模拟信号。正在测试的系统在观测期间将其转速从每秒10转提高到每秒40转。

生成转速表读数。

fs = 600;t1 = 5;t = 0:1 / fs: t1;f0 = 10;f1 = 40;rpm = 60 * linspace (f0 f1长度(t));

该信号由四个阶为1、0.5、4和6的谐波相关啁啾组成。阶4啁啾的振幅是其他啁啾的两倍。为了产生啁啾,使用梯形规则来表示相位为转速的积分。

O1 = 1;O2 = 0.5;O3 = 4;O4 = 6;pH = 2 * pi * cumtrapz(rpm / 60)/ fs;x = [111 2 1] * cos([o1 o2 o3 o4]'* pH);

可视化信号的频率-rpm映射。

RPMFreqMap(x,fs,rpm)

图频率映射包含UIMEnu,UIToolbar,UiflowContainer类型的对象。

分析放置在直升机驾驶舱内的加速度计的模拟数据。

加载直升机数据。振动测量,VIB.,以500Hz的速率进行取样10秒。检查数据显示它具有线性趋势。删除趋势以防止其降低频率估计的质量。

加载('helidata.mat') vib = detrend(vib);

绘制非线性RPM配置文件。转子延伸,直到它达到每分钟约27,600转的最大转速,然后向下沿海。

绘制(t,rpm)xlabel(“时间(s)”)ylabel('rpm'

图中包含一个轴。轴包含一个类型为line的对象。

计算频率-RPM映射。指定2.5 Hz的分辨率带宽。

[地图,频率,rpmout,时间] = rpmfreqmap(vib,fs,rpm,2.5);

可视化地图。

ImageC(时间,频率,地图)AX = GCA;ax.ydir ='普通的';包含(“时间(s)”)ylabel(的频率(赫兹)

图中包含一个轴。轴包含一个图像类型的对象。

使用更精细的分辨率带宽重复计算。使用内置功能绘制地图rpmfreqmap..频率分辨率的增加是以牺牲时间分辨率为代价的。

RPMFreqMap(VIB,FS,RPM,1.5);

图频率映射包含UIMEnu,UIToolbar,UiflowContainer类型的对象。

产生一个信号,包括两个线性啁啾和一个二次啁啾,所有在600hz采样15秒。在测试期间,产生信号的系统将其转速从每秒10转提高到每秒40转。

生成转速表读数。

fs = 600;t1 = 15;t = 0:1 / fs: t1;f0 = 10;f1 = 40;rpm = 60 * linspace (f0 f1长度(t));

线性啁啾具有1和2.5的订单。具有订单1的组件具有另一个的幅度。二次啁啾以订单6开始,并在测量结束时返回此顺序。它的幅度为0.8。使用此信息创建信号。

O1 = 1;O2 = 2.5;O6 = 6;x = 0.5 * chirp(t,O1 * f0,t1,O1 * f1)+ chirp(t,O2 * f0,t1,O2 * f1)+......0.8 * chirp(t,o6 * f0,t1,o6 * f1,'二次');

计算信号的频率-rpm映射。在每个测量单元中使用峰值幅度。指定6 Hz的分辨率。窗口与平顶窗口的数据。

[Map,FR,RP] = RPMFreqmap(x,fs,rpm,6,......“振幅”“高峰”“窗口”'flattopwin');

将频率- rpm图绘制成瀑布图。

[FR, RP] = meshgrid (FR, RP);瀑布(FR, RP,地图”)视图(6、60)包含(的频率(赫兹))ylabel('rpm') zlabel (“振幅”

图中包含一个轴。轴包含类型贴片的对象。

通过呼叫绘制交互式频率-RPM地图rpmfreqmap.没有输出参数。

加载包含放置在直升机的驾驶舱内的加速度计的模拟振动数据的文件。数据以500 Hz的速率进行采样10秒。删除数据中的线性趋势。称呼rpmfreqmap.生成频率-RPM映射的交互式图。指定2 Hz的频率分辨率。

加载Helidata.matrpmfreqmap(去趋势(vib)、fs、rpm, 2)

在图5秒后移动图中的十字线游标以确定频率为25Hz的RPM和RMS幅度。

单击zoom x.按钮在工具栏中放大2到4秒之间的时间区域。窗格出现在底部绘图中。

单击瀑布图按钮在工具栏中以瀑布图的形式显示频率- rpm映射。为提高可见性,请使用向左旋转按钮三次。移动平移器到4到6秒之间的间隔。

输入参数

全部折叠

输入信号,指定为行或列向量。

例子:COS(PI / 4 *(0:159))+ RANDN(1,160)指定嵌入在白色高斯噪声中的正弦曲线。

采样率,指定为Hz中表达的正标量。

转速,指定为以每分钟转数表示的正值的向量。rpm.必须具有与X

  • 如果您有转速表脉冲信号,请使用tachorpm提取rpm.直接。

  • 如果你没有转速计脉冲信号,使用rpmtrack提取rpm.从振动信号。

例子:100:10:3​​000指定系统最初以每分钟100转的速度旋转,最高以每分钟3000转的速度运行,以10为增量。

频率-RPM地图的分辨率带宽指定为正标量。如果res.没有指定,然后rpmfreqmap.将其设置为除以128的采样率。如果信号不够长,则该功能使用整个信号长度来计算单个频率估计。

数据类型:单身的|

名称-值对的观点

指定可选的逗号分离对名称,值参数。名称是参数名称和价值是相应的价值。名称必须出现在引号内。您可以按如下顺序指定几个名称和值对参数name1,value1,...,namen,valuen

例子:“规模”、“数据库”、“窗口”、“损害”指定频率映射估计将以分贝缩放并使用HANN窗口确定。

频率-RPM地图幅度,指定为逗号分隔对“振幅”一之一'rms'“高峰”, 或者'力量'

  • 'rms'- 返回每个估计频率的根均方幅度。

  • “高峰”-返回每个估计频率的峰值振幅。

  • '力量'-返回功率级别为每个估计频率。

相邻段之间的重叠百分比,指定为逗号分隔的对,由'重叠的'并且标量从0到100.值为0表示相邻的段不重叠。值100表示​​相邻的段由一个样品移位。更大的重叠百分比产生更平滑的地图,但增加了计算时间。看rpmordermap.为更多的信息。

数据类型:|单身的

频率-RPM映射缩放,指定为逗号分隔对组成“规模”,要么'线性'或者'D b'

  • '线性'-返回线性比例地图。

  • 'D b'- 返回具有在分贝中表示的值的对数映射。

分析窗口,指定为逗号分隔的对,由“窗口”其中一个值是:

  • “损害”指定HANN窗口。看为更多的细节。

  • 'chebwin'指定Chebyshev窗口。使用单元格数组在分贝中指定Sidelobe衰减。侧链衰减必须大于45 dB。如果未指定,则它默认为100 dB。看chebwin为更多的细节。

  • 'flattopwin'指定一个平顶窗口。看扁平翅膀为更多的细节。

  • '汉明'指定汉明窗口。看汉字为更多的细节。

  • “皇帝”指定kaiser窗口。使用单元格数组指定形状参数,β.形状参数必须是一个正标量。如果没有指定,则默认为0.5。看凯泽为更多的细节。

  • 'etcwin'指定矩形窗口。看retwin.为更多的细节。

例子:'窗口','chebwin'指定具有100 dB的Sidelobe衰减的Chebyshev窗口。

例子:'窗口',{'chebwin',60}指定具有60 dB的Sidelobe衰减的Chebyshev窗口。

例子:'窗口','kaiser'指定具有0.5的形状参数的Kaiser窗口。

例子:'窗口',{'kaiser',1}指定具有1的形状参数的kaiser窗口。

数据类型:字符|细绳|细胞

输出参数

全部折叠

频率- rpm映射,以矩阵形式返回。

频率,作为向量返回。

转速,作为向量返回。

时间瞬间,作为一个矢量返回。

分辨率带宽,作为标量返回。

参考

[1]布兰德,安德斯。噪声和振动分析:信号分析和实验程序.奇切斯特,英国:John Wiley&Sons,2011年。

扩展能力

C / c++代码生成
使用MATLAB®Coder™生成C和c++代码。

在R2015B中介绍