ノイズの多い正弦波の時間同期平均を計算します。

ホワイトガウスノイズに含まれ正弦波で構成される信号を生成します。信号は,500 Hzで20秒間サンプリングされます。10Hz の正弦波周波数と 0.01 のノイズ分散を指定します。信号の 1 周期をプロットします。

fs = 500;t = 0:1 / fs: 20:1 / fs;f0 = 10;Y = sin(2*pi*f0*t) + randn(size(t))/10;情节(t, y) xlim ([0 1 / f0])

信号の時間同期平均を計算します。信号を同期させるために,正弦波と同じ周期をもつ一連のパルスを使用します。出力引数なしで运输安全管理局を使用して結果を表示します。

tPulse = 0:1 / f0:马克斯(t);tsa (y, fs, tPulse)

指数的に減衰する2次チャープから構成される信号を生成します。信号は1 kHzで2秒間サンプリングされます。チャープの初期周波数は2赫兹で最初の1秒後に280 Hzまで増大します。減衰の特性時間は1/2秒です。信号をプロットします。

fs = 1 e3;t = 0:1 / fs: 2;x = exp (2 * t”)。*唧唧声(t ', 2, 1, 28日“二次”）;情节(t, x)

時間ベクトルを使用して时间配列を作成します。この时间配列と信号で时间表を構築します。信号のピークの位置を使用してパルス時間を判別します。時間同期平均を表示します。

ts =秒(t) ';tx =时间表(ts, x);(~, lc) = findpeaks (x, t);tsa (tx, lc)

時間同期平均を計算します。出力引数の型を表示します。サンプル時間は duration 配列に格納されます。

[股上扬,xt, xp, xrpm] = tsa (tx, lc);谁x *

名称大小字节类属性x 2001x1 16008 double xp 9x1 1135时间表xrpm 1x1 8 double xt 9x1 74 duration xta 9x1 1133时间表

时间配列をdatetimeベクトルに変換します。datetimeベクトルと信号を使用して时间表を構築します。時間同期平均を計算しますが,今度は,15回転分の平均を求めます。

出力引数の型を表示します。入力的时间表がdatetimeベクトルを使用していても,ここでも,サンプル時間は时间配列に格納されます。

dtb = datetime (datevec (ts));德勤=时间表(防晒霜,x);nr = 15;tsa(德勤,信用证,“NumRotations”nr)

[dta, dt, dp, drpm] = tsa(德勤、lc、“NumRotations”、nr);谁d *

名称大小字节类属性dp 135x1 3151时刻表drpm 1x1 8 double dt 135x1 1082 duration dta 135x1 3149时刻表dtb 2001x1 32016 datetime dtt 2001x1 49001时间表

スイッチオフ後ファンブレードが遅くなっていくときのファンブレードの位置の時間同期平均を計算します。

2400 rpmで回転する卓上扇風機をオフにします。空気抵抗(ベアリングの摩擦からの寄与は無視できるとします)により約5秒でファンのローターが停止します。高速カメラでファンブレードの1つのx座標を1 kHzのレートで測定します。

fs = 1000;t = 0:1 / fs: 5 - 1 / f;rpm0 = 2400;

このファンブレードを,ローターを中心として半径10厘米を回転する質点として理想化します。ブレードには速度に比例する抵抗力がかかり,その結果,位相角は次の式で求められます。

ここで, ${\mathit{f}}_0$は初期周波数, $\mathit{T}＝0．75$秒は減衰時間です。

一个= 0.1;f0 = rpm0/60;T = 0.75;φ= 2 *π* f0 * T * (1-exp (T / T));

x座標とy座標を計算し,プロットします。ホワイトガウスノイズを付加します。

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

同期される信号を求めます。関数tachorpmを使用して,パルス時間を見つけます。検索する時間を2.5秒までに制限します。回転速度をプロットして,その指数関数的減衰を確認します。

(rpm, ~, tp) = tachorpm (x (t < 2.5), fs);tachorpm (x (t < 2.5), fs)

正弦波の周期に対応する時間同期平均信号を計算し,プロットします。周波数領域で平均化を実行します。

clf tsa (x, fs, tp,“方法”，fft的）