synthesize_bearing_signal — ACOUSTICS synthesis op

資料種類:signal(僅由參數決定的運算子 —— 不接受影像或資料輸入)

呼叫:import acoustics; acoustics.synthesize_bearing_signal(rate=25600.0, duration=1.0, carrier_hz=3000.0, defect_hz=107.0, modulation=0.5, mode='am', damping=0.05, noise_sigma=0.0, seed=None)(或 opsacoustics.get("synthesize_bearing_signal"))

用法

以已知缺陷頻率調幅的共振 —— 答案已知的真值。

> 以下的詳細說明為原文 —— 摘要與標題已翻譯。

This is the whole reason envelope analysis exists, built forwards so the

answer is known before the measurement. A spall on a bearing race does not

radiate at the defect rate; it strikes a structure that rings at a much

higher resonance, once per defect passage. What reaches the microphone is a

carrier at the resonance, modulated at the defect rate, and the defect

rate itself is not present in the signal as a frequency component at all.

`mode="am"` gives the exactly analysable case,

`x(t) = (1 + m cos(2 pi f_d t)) sin(2 pi f_c t)`. Its analytic envelope is

exactly `1 + m cos(2 pi f_d t) for m < 1`, so the single-sided envelope

spectrum has a line of amplitude exactly m at `f_d` and nothing else.

Measured with `m = 0.5: :func:envelope_spectrum` returns a peak at

107.000000 Hz of amplitude 0.499677 (the 0.06 % shortfall is the band-pass

filter rolling off across the two sidebands, not the demodulation).

`mode="impulse"` gives the physically shaped case: an impulse train at

`f_d, each impulse ringing down as `exp(-2 pi zeta f_c t) sin(2 pi f_c

t)`. The envelope spectrum then shows f_d` and its harmonics, which

is what a real record looks like. Measured with `f_d = 107` Hz: the

envelope-spectrum peak is at 107.000000 Hz and the harmonics at 214 and

321 Hz carry 0.6542 and 0.4748 of the fundamental's amplitude.

In am mode the raw spectrum has nothing at `f_d`: measured, the raw

single-sided amplitude at 107 Hz is 4.3e-16, while the carrier reads

1.000000 and each sideband at 2893 and 3107 Hz reads 0.250000 — exactly

`m/2`, as amplitude modulation requires. In impulse mode the raw amplitude

at 107 Hz is 0.01165, not zero (an impulse train is not a pure product), but

still 18x below what the envelope spectrum recovers from the same record.

Raises `ValueError`: any non-real / non-finite / string / bool scalar,

`rate <= 0, duration <= 0, modulation outside [0, 1)` in am

mode (at `m >= 1 the envelope is |1 + m cos|`, which folds and puts

energy at `2 f_d — a rectified envelope, not the modulation), damping`

outside `(0, 1), a total length over :data:MAX_SAMPLES`, and — the one

that matters — **any requested frequency at or above Nyquist, including the

upper modulation sideband** `f_c + f_d`. An aliased carrier would come

back as a plausible signal at the wrong frequency with no error.

詳細使用指南

acoustic_condition_monitoring 族使用指南

參考(範例資料・文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者・年份)與用途見上面的族使用指南

可執行的範例(實際呼叫該運算子並已驗證的樣例)

acoustic_condition_monitoringpy -3.11 examples/acoustic_condition_monitoring.py

poc_bearing_diagnosispy -3.11 examples/poc_bearing_diagnosis.py

型別可銜接的下一個運算子(可接受 signal 作為輸入)

stft · envelope_spectrum · spectral_kurtosis · cepstrum · angular_resample · order_spectrum · octave_spectrum · weighting_response

同類別(synthesis)

synthesize_speed_ramp


*Provenance: acoustics.py — ACOUSTICS 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.