周波数対ウェーブレット変換の大きさをプロットするにはどうすればよいですか？

Morlet連続ウェーブレット変換を実行しています。私はwscalogram信号を受け取ったので、次の図のように周波数の大きさをプロットしたいのですが、どうすればいいのかわかりません。

scal2freqMATLAB関数を使用して、スケールを擬似周波数に変換しました。また、信号にいくつかの周波数があり、それらは大きな減衰比（4％）を持っているため、プロットではよく見えません。これらの非常に減衰したモードを誇張するにはどうすればよいですか？

私はMATLABを使用しています、ここに私のコードがあります：

% Import the text4.txt to matlab workspace. and save it under name "data"
t=linspace(0,30,301);
Fs=ceil(inv(t(2)-t(1)));
x=data(:,4); % use x=data(:,3),x=data(:,5) too. first column is time,second is refrence
wname = 'morl';
scales = 1:1:256;
coefs = cwt(x,scales,wname,'lvlabs');
freq = scal2frq(scales,wname,1/Fs);
surf(t,freq,abs(coefs));shading('interp');
axis tight; xlabel('Seconds'); ylabel('Pseudo-Frequency (Hz)');
axis([0 30 0 1 0 60])
xlabel('Time'); ylabel('Frequency');
figure;
sc=wscalogram('image',coefs,'scales',freq,'ydata',x);
hold on 

私の現在のプロットは次のようになります。

コードを少し変更しますが、大きな変更はありません。どちらにしても正しい結果が得られます。ここでこのテンプレートコードを使用すると、問題は発生しません。正しい結果が得られます。

clear all; 
t=linspace(0,30,301); 
Fs = (inv(t(2)-t(1))); 
x=randn(100,1);  
wname = 'morl'; 
scales = 1:1:256; 
chefs = cwt(x,scales,wname,'lvlabs');

freq = scal2frq(scales,wname,1/Fs);

figure; 
coefsSquared = abs(coefs).^2; 
imagesc(coefsSquared); 
grid off;

%Pick one of the columsn to plot: 
figure; 
plot(coefsSquared(:,47))

最初に表示した画像の右側にある図を取得するには、それは単にタイムスライススペクトルです。つまり、周波数と時間のマトリックスの1列です。

たとえば、周波数と時間の関係で次の2乗係数行列を取得できます。

次に、47列目のタイムスライスのスペクトルを確認するとします。次に、これを示します。

ここで、列47に対応するタイムスライスの係数のパワー対周波数を確認できます。

連続ウェーブレット変換（CWT）は、時間スケール分析の方法です。はい、周波数ではなく、スケールを正しく読みます。ただし、スケールを周波数にマップすることは可能であり、非常に簡単です。MATLABユーザーであるため、次のことを行うこの関数を使用することをお勧めします。

F = scal2frq（A、 'wname'、DELTA）は、Aで指定されたスケール、ウェーブレット関数 'wname'（詳細についてはwavefunを参照）、およびサンプリング周期DELTAに対応する擬似周波数を返します。

ご覧のとおり、周波数は信号の実際の周波数ではなく、近似値です。実周波数が必要な場合は、短時間フーリエ変換（STFT）を使用する必要があります。

それでは、いつCWTを使用すべきでしょうか？一般的な原則として、ウェーブレットは、信号が過渡的な場合（つまり、変化が速い場合）にうまく機能します。

編集：笑、私はここで何をしているのかわかりません。すみません:-Dはあなたの質問のポイントを完全に逃しました。タイトルを読んでください。

まだ助けが必要かどうかわからない...

周波数対振幅を求める場合は、「mesh / surf」コマンドを使用して係数をプロットします。

次に、求める結果が得られるようにプロットの方向を変更します。これについては、「view」コマンドを参照してください。