如何绘制频率与小波变换幅度的关系？

我正在运行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))

为了使该图位于您放置的第一张图像的右侧，这就是时间切片频谱。换句话说，频率时间矩阵的一列。

例如，我们可以得到以下频率与时间的平方系数矩阵：

然后让我们假设我们希望在第47列中查看时间片的频谱。然后，在此处给出：

在这里，您将看到对应于第47列的时间片的系数幂与频率的关系。

连续小波变换（CWT）是一种时标分析方法。是的，您正确地阅读了它，按比例而不是频率。但是，可以将比例尺映射到频率，甚至非常容易。由于您是MATLAB用户，因此您可能需要使用此函数，该函数将执行以下操作：

F = scal2frq（A，'wname'，DELTA）返回与A给定的小数位对应的伪频率，小波函数'wname'（有关更多信息，请参见wavefun）和采样周期DELTA。

如您所见，频率不是信号中的实际频率，而是近似值。如果需要实际频率，则应使用短时傅立叶变换（STFT）。

那么，什么时候应该使用CWT？作为一般原则，当信号是瞬态的（即具有快速变化）时，小波效果很好。

编辑：哈哈，我不知道我在这里做什么。抱歉：-D完全没有回答您的问题。只需阅读标题即可。

不知道您是否仍然需要帮助...

如果您想要频率与幅度的关系，请使用“网格/冲浪”命令绘制系数。

然后更改绘图的方向，以便获得所需的结果。为此，请参见“查看”命令。