まず、scal2frq は、帯域幅を求めるものではなく、疑似周波数に対する、中心周波数を求める機能となります。 残念ながら、疑似周波数に対しての帯域幅を確認する直接的な方法は提供されていません。
お手数をお掛けすることになりますが、下記のようなスクリプト例で帯域幅を確認することができます (原理的には、インパルス信号に対して、疑似周波数に相当するスケールファクタを設定した CWT を適用します。)
下記例は、疑似周波数が 20Hz の場合となります。
fs=2000;
t = 0:1/fs:2-1/fs;
wname = 'cmor1-1';
% すべて0の値を持つ信号の中心にインパルスを付加
sig = zeros(1,length(t));
sig(length(t)/2) = 1;
sf = 100; % scaleファクタ = 100 => Fc = 20Hz
filterd = cwt(sig, sf, wname);
Y = fft(filterd);
F = (0:3999)/4000*fs;
plot(F, abs(Y))
% 中心周波数付近を拡大
xlim([fs-2*scal2frq(sf,'cmor1-1',1/fs) fs])
中心周波数は 1980Hz と現れますが、これは20Hz のミラー成分に相当します. どのあたりまでの振幅の大きさを帯域幅と捉えるかに寄りますが、 例えば、中心周波数に対して、振幅が半分までを帯域幅と捉える場合は、+/-5Hz 程度と解釈できます。
