600hzで5秒間サンプリングされたシミュレ。テスト中のシステムは,観察期間中に回転速度を1秒あたり10回転から40回転に上昇させます。

タコメタの読み取り値を生成します。

Fs = 600;T1 = 5;T = 0:1/fs:t1;F0 = 10;F1 = 40;RPM = 60*linspace(f0,f1,length(t));

信号は、1、0.5、4および6の次数をもつ,調和的に関連した4つのチャープで構成されています。次数4のチャプの振幅は他のチャプの振幅の2倍です。チャプを生成するために，台形則を使用して位相を回転速度の積分で表します。

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

信号の次数-rpmマップを可視化します。

rpmordermap (x, fs, rpm)

信号の平均次数スペクトルを可視化します。スペクトルのピクは次数-rpmマップのリッジに対応します。

orderspectrum (x, fs, rpm)

ヘリコプタのコクピットに設置された加速度計から取得したシミュレションデタを解析します。

ヘリコプタのデタを読み込みます。振動測定値vibは，500 Hzのレトで10秒間サンプリングされています。デタには線形トレンドがあります。トレンドを除去して，トレンドによる次数推定の質の低下を防ぎます。

负载(“helidata.mat”) vib =趋势(vib);

非線形RPMプロファereplicationルをプロットします。ローターは,最大回転速度である1分あたり27600回転に達するまで回転数を上昇させ,その後降下させます。

情节(t, rpm)包含(“时间(s)”) ylabel (“转”）

信号の平均次数スペクトルを計算します。既定の次数分解能を使用します。

orderspectrum (vib fs, rpm)

rpmordermapを使用して，より細かい次数分解能で計算を繰り返します。次数が低いほど，より明確に分解されます。

[map,order] = rpmordermap(vib,fs,rpm,0.005);orderspectrum(地图,顺序)

各推定次数のパワレベルを計算します。結果をデシベルで表示します。

[map,order] = rpmordermap(vib,fs,rpm,0.005，“振幅”，“权力”）;Spec = orderspectrum(地图，顺序);情节(顺序,pow2db(规范))包含(“订单号”) ylabel (“阶功率振幅(dB)”网格)在