シーケンス / 関数top
オクターブに関するFAQ
Qオクターブ解析結果が異なります
FAQ ID:c053
ABCRating関数のA特性フィルタを使用した場合と、OCTI関数のA特性を適用した場合で解析した結果が異なります。
具体的には以下のシーケンスを参照してください。
;;-------------------------------------------------------
ABCRating_data = ABCRating(data,1,0,1,0) ; A-Rating
OctI(0,-1,0,-1,0,0,0)
Oct_1 = dB(OctA(ABCRating_data,20,20000)/(20*10^-6)); 1/1 Octave
OctI(0,-1,0,-1,1,0,0)
Oct_2 = dB(OctA(data,20,20000)/(20*10^-6)) ; A-Rating & 1/1 Octave
;-------------------------------------------------------
計算アルゴリズムの違いにより差異が生じます。
1) ABCRating()関数では、Aウェイティングに対してスムーズな伝達関数を持ちます。
e.g.
-30.2dB @ 50Hz
-29.9dB @ 51Hz
-29.6dB @ 52Hz
...
-26.2dB @63Hz
...
-22.5dB @ 80Hz
(これらの値はサンプルで示していますので、正確な値とは少し異なります。)
これに対して OctI(0,-1,0,-1,1,0,0)では、オクターブバンドと等しい減衰を持つバンドパスフィルタを適用します。例えば63Hzでは、バンド全体に対するRMS結果に対して-26.2dBの重みをかけます。
50Hz、0dB RMSの正弦波を入力した場合、使用する関数により -30.2dB or-26.2dBとなります。
2) DIN EN 61260や他の規格でバンド誤差を持つオクターブフィルタを定義しています。これらは実際に使用されるバンドパスフィルタの組み込みに対してさまざまな方法を許容しています。
e.g. 4, 6, 8次のバターワース or チェビシェフフィルタを使用することができます。それぞれは異なる伝達関数を持ちます。
差異は中心バンドでは明らかではありませんが、ストップバンドで大きくなります。
Qオクターブ分析結果の3D表示
FAQ ID:c030
測定データに対してオクターブ分析を行い、これを3D表示させるため、以下のシーケンス作成しました。
data = sin ( ramp ( 0, 1/40000, 1000 ) * (6.28 * 300) )
OctI(0, -1, 1, -1, 0, 0, 0) ; 1/1オクターブ、Fast時間重み
res = OctA( data, 16, 8000) ; 解析周波数16~8000Hz
1/1オクターブ解析のため16, 31.5, 63, 125, 250, 500, 1000, 2000, 4000, 8000Hz、すなわち変数resには10点ごとにFast時間重みでの1/1オクターブ解析が格納されています。
変数resを周期ごとに3D表示するには下記の2通りがあります。((1)がおすすめです。)
(1)
SetSegLen(res, 10)
のようにSetSegLen()関数で変数resをセグメント化してカーブウィンドウでウォーターフォールまたはカラーマップ表示を行ってください。
(2)
カーブウィンドウでウォーターフォールまたはカラーマップ表示を行った後、[設定]>[イベント、セグメント、間隔]をクリックして、下記のように[周期]タブで[周期の選択]に[全て]、時間の配置に[最初の周期のx0]、[周期長さ]に[10]を指定して[OK]をクリックしてください。
結果は下記のような表示になります。
Qオクターブ分析データの表示
FAQ ID:c021
オクターブ分析されたデータのX軸は周波数スケールで等間隔になるように定義されています。但し、これらの値を以下のようにそのまま表示すると端数が出てしまいます。
このように単なる対数表示ではX/軸ラベルが一般的なオクターブ表示と多少異なりますし、軸の設定も煩雑です。通常、これらの端数を丸めた値を使用しています。
下表は1/3オクターブの中心周波数です。
Center [Hz] |
---|
1000 |
1250 |
1600 |
2000 |
2500 |
3150 |
4000 |
5000 |
6300 |
8000 |
10000 |
そこで、オクターブ分析されたデータのX軸をスケーリングします。スケーリングは以下の式に従います。
中心周波数=10^(スケールパラメータ/10)
X-Scaling | 中心周波数[Hz] |
---|---|
.. | .. |
-3 | 0.5 |
-2 | 0.63 |
-1 | 0.8 |
0 | 1 |
1 | 1.25 |
2 | 1.6 |
3 | 2 |
4 | 2.5 |
5 | 3.15 |
6 | 4 |
7 | 5 |
8 | 6.3 |
9 | 8 |
10 | 10 |
11 | 12.5 |
.. | .. |
20 | 100 |
30 | 1000 |
40 | 10000 |
41 | 12500 |
42 | 16000 |
43 | 20000 |
.. | .. |
このようにスケーリングすることで、X軸データは周波数スケールからリニアスケールに変換され、1/3オクターブではXスケールは0, 1, 2...となります。1/1オクターブでは0, 3, 6, 9, 12...、1/12-オクターブでは0, 0.25, 0.5, 0.75, 1, 1.25...です。
各周波数バンドのX間隔は次のようになります。
Width of Filter Band | Delta-x |
---|---|
1/1オクターブ | 3 |
1/3オクターブ | 1 |
1/12オクターブ | 0.25 |
1/24オクターブ | 0.125 |
FAMOSはこのスケーリングに対応した表示モードを備えていますので、カーブウィンドウに表示させる際にはモードを切り替えるようにしてください。
標準表示モード
オクターブ表示モード